记录编号 484226 评测结果 AAAAAAAAAAA
题目名称 [网络流24题] 方格取数问题 最终得分 100
用户昵称 Gravatarpanda_2134 是否通过 通过
代码语言 C++ 运行时间 0.023 s
提交时间 2018-01-22 17:08:31 内存使用 1.08 MiB
显示代码纯文本
  1. //发现要用到复杂(超过6个符号)的位运算时马上住手,改用循环预处理!虽然慢一点,但是不会错!
  2. #include <bits/stdc++.h>
  3.  
  4. using namespace std;
  5.  
  6. struct Edge {
  7. int v, flow, cap, next;
  8. };
  9.  
  10. const int MAXN = 10002, MAXM = 4e4, INF = 0x3f3f3f3f;
  11. const int di[] = { 1, -1, 0, 0 }, dj[] = { 0, 0, 1, -1 };
  12. int R, C, N, M, S, T, Tot, e_ptr=1, head[MAXN+10]; Edge E[(MAXM<<1)+10];
  13. int ID[110][110];
  14. void AddEdge(int u, int v, int cap) {
  15. E[++e_ptr] = (Edge) { v, 0, cap, head[u] }; head[u] = e_ptr;
  16. E[++e_ptr] = (Edge) { u, 0, 0, head[v] }; head[v] = e_ptr;
  17. }
  18.  
  19. int d[MAXN+10], cur[MAXN+10];
  20. bool BFS() {
  21. queue<int> Q;
  22. memset(d, 0xff, sizeof(d));
  23. Q.push(S); d[S] = 0;
  24. while(!Q.empty()) {
  25. int u = Q.front(); Q.pop();
  26. for(int j=head[u]; j; j=E[j].next) {
  27. int v = E[j].v, f = E[j].flow, c = E[j].cap;
  28. if(f<c && d[v] == -1) {
  29. d[v] = d[u] + 1;
  30. if(v == T) return true;
  31. else Q.push(v);
  32. }
  33. }
  34. }
  35. return false;
  36. }
  37.  
  38. int DFS(int u, int flow) {
  39. if(u == T || flow == 0) return flow;
  40. int res = flow;
  41. for(int &j=cur[u]; j; j=E[j].next) {
  42. int v = E[j].v, f = E[j].flow, c = E[j].cap;
  43. if(d[v] == d[u] + 1) {
  44. int tmp = DFS(v, min(res, c-f));
  45. res -= tmp;
  46. E[j].flow += tmp;
  47. E[j^1].flow -= tmp;
  48. }
  49. }
  50. return flow - res;
  51. }
  52.  
  53. int Dinic() {
  54. int MaxFlow = 0, CurFlow = 0;
  55. while(BFS()) {
  56. for(int i=1; i<=N; i++)
  57. cur[i] = head[i];
  58. while( (CurFlow = DFS(S, INF)) != 0 )
  59. MaxFlow += CurFlow;
  60. }
  61. return MaxFlow;
  62. }
  63.  
  64. inline bool valid(int i, int j) {
  65. return i>=1 && i<=R && j>=1 && j<=C;
  66. }
  67.  
  68. void Label() {
  69. int idx = 0;
  70. for(int i=1; i<=R; i++)
  71. for(int j = 2 - (i&1); j<=C; j+=2)
  72. ID[i][j] = ++idx;
  73. for(int i=1; i<=R; i++)
  74. for(int j = 1 + (i&1); j<=C; j+=2)
  75. ID[i][j] = ++idx;
  76. }
  77.  
  78. void Init() {
  79. int Cur;
  80. scanf("%d%d", &R, &C);
  81. N = R*C + 2; S = R*C + 1; T = R*C + 2;
  82. Label();
  83. for(int i=1; i<=R; i++)
  84. for(int j=1; j<=C; j++) {
  85. scanf("%d", &Cur); Tot += Cur;
  86. if(((i^j)&1) == 0) {
  87. AddEdge(S, ID[i][j], Cur);
  88. for(int t=0; t<4; t++) {
  89. int ni = i + di[t], nj = j + dj[t];
  90. if(!valid(ni, nj)) continue;
  91. AddEdge(ID[i][j], ID[ni][nj], INF);
  92. }
  93. } else
  94. AddEdge(ID[i][j], T, Cur);
  95. }
  96. }
  97.  
  98. void Work() {
  99. printf("%d", Tot - Dinic());
  100. }
  101.  
  102. int main() {
  103. freopen("grid.in", "r", stdin);
  104. freopen("grid.out", "w", stdout);
  105. Init(); Work();
  106. return 0;
  107. }