bzoj2132【圈地计划】

题面

思路：

一开始以为和为了博多一样，两边连一样的，后来发现中间连负边的话根本不会割，即割断两块收益为负，所以WA的起飞……

正解是先黑白染色，每个点和它周围的点连边方式不同。对于黑点A，S-->A表示商业区的价值，A-->T表示工业区的价值，白点相反，对于相邻的情况，边权表示相邻的价值和，这样割断相当于选择同一用途，代价为正。

 #include <map>
 #include <set>
 #include <cmath>
 #include <queue>
 #include <stack>
 #include <cstdio>
 #include <string>
 #include <vector>
 #include <cstring>
 #include <complex>
 #include <cstdlib>
 #include <iostream>
 #include <algorithm>
 #define rg register
 #define ll long long
 using namespace std;

 inline int gi()
 {
     rg int r = ; rg bool b = ; rg char c = getchar();
     while ((c < '' || c > '') && c != '-') c = getchar();
     if (c == '-') b = ;
     while (c >= '' && c <= '') { r = (r << ) + (r << ) + c - '', c = getchar(); }
     if (b) return r; return -r;
 }

 const int inf = 2e8+, N = , M = 5e5+;
 int n,m,num,S,T;
 int dx[]={,,-,},dy[]={-,,,},dep[N*N],f[N*N],c[N][N];
 bool col[N][N];
 queue <int> q;
 struct Edge
 {
     int to,nx,cap;
 }eg[M];

 inline void add(int x,int y,int z)
 {
     eg[++num]=(Edge){y,f[x],z}, f[x]=num;
 }

 inline void link(int x,int y,int z)
 {
     add(x,y,z), add(y,x,);
 }

 inline int dfs(int o,int flow)
 {
     if (o == T || !flow) return flow;
     rg int to,cap,fl,tmp,i;
     fl=;
     for (i=f[o]; i; i=eg[i].nx)
         {
             to=eg[i].to, cap=eg[i].cap;
             if (dep[o] != dep[to]- || cap <= )
                 continue;
             tmp=dfs(to,min(flow,cap));
             if (!tmp) continue;
             fl+=tmp, flow-=tmp;
             eg[i].cap-=tmp, eg[i^].cap+=tmp;
             if (!flow) break;
         }
     if (!fl) dep[o]=;
     return fl;
 }

 inline int bfs()
 {
     memset(dep,,sizeof(dep));
     rg int o,to,cap,i;
     q.push(S); dep[S]=;
     while (!q.empty())
         {
             o=q.front(), q.pop();
             for (i=f[o]; i; i=eg[i].nx)
                 {
                     to=eg[i].to, cap=eg[i].cap;
                     if (dep[to] || cap <= )
                         continue;
                     dep[to]=dep[o]+, q.push(to);
                 }
         }
     return dep[T];
 }

 inline int dinic()
 {
     rg int flow;
     flow=;
     while (bfs()) flow+=dfs(S,inf);
     return flow;
 }

 int main()
 {
     rg int i,j,k,x,y,dis,ans;
     n=gi(), m=gi();
     ans=S=, T=n*m+, num=;
     for (i=; i<=n; ++i)
         for (j=; j<=m; ++j)
             if ((i+j)&)
                 col[i][j]=;
     for (i=; i<=n; ++i)
         for (j=; j<=m; ++j)
             {
                 dis=gi(), ans+=dis;
                 if (col[i][j])
                     link(S,(i-)*m+j,dis);
                 else
                     link((i-)*m+j,T,dis);
             }

     for (i=; i<=n; ++i)
         for (j=; j<=m; ++j)
             {
                 dis=gi(), ans+=dis;
                 if (col[i][j])
                     link((i-)*m+j,T,dis);
                 else
                     link(S,(i-)*m+j,dis);
             }
     for (i=; i<=n; ++i) for (j=; j<=m; ++j) c[i][j]=gi();
     for (i=; i<=n;++i)
         for (j=; j<=m; ++j)
             {
                 if (!col[i][j]) continue;
                 for (k=; k<; ++k)
                     {
                         x=i+dx[k], y=j+dy[k];
                         if (x > n || x <  || y > m || y < )
                             continue;
                         add((i-)*m+j,(x-)*m+y,c[i][j]+c[x][y]);
                         add((x-)*m+y,(i-)*m+j,c[i][j]+c[x][y]);
                         ans+=c[i][j]+c[x][y];
                     }
             }
     printf("%d\n",ans-dinic());
     return ;
 }