标签：
题解：
二元组建图：http://www.cnblogs.com/chenyushuo/p/5146626.html
这看上去是一个多元关系，但我们可以转化为二元关系
我们对每个节点x0新建两个点x1，x2分别表示x0周围是否都是选文和是否都选理，然后有以下四种二元关系
x0选文且x1为是，那么有v1的收入
x0选理且x2为是，那么有v2的收入
x1为是，与x0相邻的点必须选文，如果选了理，则造成inf的代价
x2为是，与x0相邻的点必须选理，如果选了文，则造成inf的代价
然后二元组建图即可
code：
1 #include<cstdio> 2 #include<iostream> 3 #include<cmath> 4 #include<cstring> 5 #include<algorithm> 6 #define maxn 30005 7 #define maxm 1200000 8 #define inf 1061109567 9 using namespace std; 10 char ch; 11 bool ok; 12 void read(int &x){ 13 for (ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if (ch==‘-‘) ok=1; 14 for (x=0;isdigit(ch);x=x*10+ch-‘0‘,ch=getchar()); 15 if (ok) x=-x; 16 } 17 int n,m,sum,pos[105][105][3],x; 18 struct flow{ 19 int s,t,tot,now[maxn],son[maxm],pre[maxm],val[maxm]; 20 int dis[maxn],head,tail,list[maxn]; 21 bool bo[maxn]; 22 void init(){s=0,t=n*m*3+1,tot=1,memset(now,0,sizeof(now));} 23 void put(int a,int b,int c){pre[++tot]=now[a],now[a]=tot,son[tot]=b,val[tot]=c;} 24 void add(int a,int b,int c){put(a,b,c),put(b,a,0);} 25 bool bfs(){ 26 memset(bo,0,sizeof(bo)); 27 head=0,tail=1,list[1]=s,dis[s]=0,bo[s]=1; 28 while (head<tail){ 29 int u=list[++head]; 30 for (int p=now[u],v=son[p];p;p=pre[p],v=son[p]) 31 if (val[p]&&!bo[v]) bo[v]=1,dis[v]=dis[u]+1,list[++tail]=v; 32 } 33 return bo[t]; 34 } 35 int dfs(int u,int rest){ 36 if (u==t) return rest; 37 int ans=0; 38 for (int p=now[u],v=son[p];p&&rest;p=pre[p],v=son[p]) 39 if (val[p]&&dis[v]==dis[u]+1){ 40 int d=dfs(v,min(rest,val[p])); 41 val[p]-=d,val[p^1]+=d,ans+=d,rest-=d; 42 } 43 if (!ans) dis[u]=-1; 44 return ans; 45 } 46 int dinic(){ 47 int ans=0; 48 while (bfs()) ans+=dfs(s,inf); 49 return ans; 50 } 51 }f; 52 const int dx[4]={1,0,-1,0}; 53 const int dy[4]={0,1,0,-1}; 54 int main(){ 55 read(n),read(m),f.init(); 56 for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) 57 pos[i][j][0]=(i-1)*m+j,pos[i][j][1]=pos[i][j][0]+n*m,pos[i][j][2]=pos[i][j][1]+n*m; 58 for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) read(x),sum+=x,f.add(pos[i][j][0],f.t,x<<1); 59 for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) read(x),sum+=x,f.add(f.s,pos[i][j][0],x<<1); 60 for (int i=1;i<=n;i++) for (int j=1;j<=m;j++){ 61 read(x),sum+=(x<<1); 62 f.add(pos[i][j][0],pos[i][j][1],x),f.add(pos[i][j][1],pos[i][j][0],x); 63 f.add(f.s,pos[i][j][0],x),f.add(f.s,pos[i][j][1],x); 64 f.add(pos[i][j][0],f.t,x<<1),f.add(pos[i][j][1],f.t,x<<1); 65 } 66 for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) for (int k=0;k<4;k++){ 67 int x=i+dx[k],y=j+dy[k]; 68 if (x<=0||x>n||y<=0||y>m) continue; 69 f.add(pos[x][y][0],pos[i][j][1],inf); 70 } 71 for (int i=1;i<=n;i++) for (int j=1;j<=m;j++){ 72 read(x),sum+=(x<<1); 73 f.add(pos[i][j][0],pos[i][j][2],x),f.add(pos[i][j][2],pos[i][j][0],x); 74 f.add(pos[i][j][0],f.t,x),f.add(pos[i][j][2],f.t,x); 75 f.add(f.s,pos[i][j][0],x<<1),f.add(f.s,pos[i][j][2],x<<1); 76 } 77 for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) for (int k=0;k<4;k++){ 78 int x=i+dx[k],y=j+dy[k]; 79 if (x<=0||x>n||y<=0||y>m) continue; 80 f.add(pos[i][j][2],pos[x][y][0],inf); 81 } 82 printf("%d\n",sum-(f.dinic()>>1)); 83 return 0; 84 }
标签：
原文地址：http://www.cnblogs.com/chenyushuo/p/5147375.html
