Description

有一个n行m列的黑白棋盘，你每次可以交换两个相邻格子（相邻是指有公共边或公共顶点）中的棋子，最终达到目标状态。要求第i行第j列的格子只能参与mi,j次交换。

Input
第一行包含两个整数n，m（1<=n, m<=20）。以下n行为初始状态，每行为一个包含m个字符的01串，其中0表示黑色棋子，1表示白色棋子。以下n行为目标状态，格式同初始状态。以下n行每行为一个包含m个0~9数字的字符串，表示每个格子参与交换的次数上限。

Output
输出仅一行，为最小交换总次数。如果无解，输出-1。

Sample Input
3 3
110
000
001
000
110
100
222
222
222

Sample Output
4

分析：
题面让人觉得很奇怪
实际上我们可以看做只有黑棋子，这样就好处理多了

说实话，这道题的建图是我见过的最zz的
显然初始状态和最终状态该位置都是黑棋时，
这个黑棋是没有必要移动的，对答案没有贡献
所以我们就可以把ta变成白棋

题目的限制不是在棋子上，而是在格子上
我们把每个格子拆成三个点a,b,c

(1)格子的初始状态是1，连接(s,a,1,0),(b,a,m[i][j]/2,0),(a,c,(m[i][j]+1)/2,0)
(2)格子末状态是1，连接(a,t,1,0),(b,a,(m[i][j]+1)/2,0),(a,c,m[i][j]/2,0)
(3)始末都是0，连接(b,a,m[i][j]/2,0),(a,c,m[i][j]/2,0)
(4)相邻格子x,y,连接(xc,yb,INF,1)

解释
对于一次交换，会使用相邻两个格子各一次，一共会使用两次，
所以我们将边权除以二
对于初始黑子在的点，与末尾黑子在的点，
其实对于这个黑子只用交换一次，所以只会使用1的流量，
所以我们把ta的边权令为(m[i][j]+1)/2，这个细节比较重要

附测试样例及图
input
2 2
10
01
01
10
11
11

output
2

tip

发现自己的spfa都写不好了

注意所谓的相邻格子包括对角线

在获取点的编号的时候，get函数中n和m不要弄混了
手残党。。。STO

这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>using namespace std;const int INF=0x33333333;
const int N=50010;
int st[N],tot=-1,s,t;
struct node{int x,y,v,c,nxt;
};
node way[N<<2];
int pre[N],dis[N],q[N],tou,wei;
bool p[N];
int mp[100][100],num[100][100],n,m;int get(int x,int y){return (x-1)*m+y;}void add(int u,int w,int v,int cc)
{tot++;way[tot].x=u;way[tot].y=w;way[tot].v=v;way[tot].c=cc;way[tot].nxt=st[u];st[u]=tot;tot++;way[tot].x=w;way[tot].y=u;way[tot].v=0;way[tot].c=-cc;way[tot].nxt=st[w];st[w]=tot;
}int spfa(int s,int t)
{memset(dis,0x33,sizeof(dis));memset(pre,-1,sizeof(pre));memset(p,1,sizeof(p));wei=tou=0;q[++wei]=s;dis[s]=0;p[s]=0;do{int r=q[++tou];for (int i=st[r];i!=-1;i=way[i].nxt)if (way[i].v&&way[i].c+dis[r]<dis[way[i].y]){dis[way[i].y]=dis[r]+way[i].c;pre[way[i].y]=i;if (p[way[i].y]){q[++wei]=way[i].y;p[way[i].y]=0;} }p[r]=1;}while (tou<wei);return dis[t]!=INF;
}void doit()
{int ans=0;while (spfa(s,t)){int sum=INF;for (int i=t;i!=s;i=way[pre[i]].x)sum=min(sum,way[pre[i]].v);   //pre[i]ans+=sum*dis[t];for (int i=t;i!=s;i=way[pre[i]].x)way[pre[i]].v-=sum,way[pre[i]^1].v+=sum;}printf("%d",ans);
}void lianbian()
{int i,j;s=0; t=n*m*3+1;for (i=1;i<=n;i++)for (j=1;j<=m;j++){if (mp[i][j]==1){add(s,get(i,j),1,0);add(get(i,j)+n*m,get(i,j),num[i][j]/2,0);add(get(i,j),get(i,j)+n*m*2,(num[i][j]+1)/2,0);}else if (mp[i][j]==2){add(get(i,j),t,1,0);add(get(i,j)+n*m,get(i,j),(num[i][j]+1)/2,0);add(get(i,j),get(i,j)+n*m*2,num[i][j]/2,0);}else{add(get(i,j)+n*m,get(i,j),num[i][j]/2,0);add(get(i,j),get(i,j)+n*m*2,num[i][j]/2,0);}if (i+1<=n) add(get(i,j)+n*m*2,get(i+1,j)+n*m,INF,1);if (j+1<=m) add(get(i,j)+n*m*2,get(i,j+1)+n*m,INF,1);if (i-1>0) add(get(i,j)+n*m*2,get(i-1,j)+n*m,INF,1);if (j-1>0) add(get(i,j)+n*m*2,get(i,j-1)+n*m,INF,1);if (i+1<=n&&j+1<=m) add(get(i,j)+n*m*2,get(i+1,j+1)+n*m,INF,1);if (i-1>0&&j-1>0) add(get(i,j)+n*m*2,get(i-1,j-1)+n*m,INF,1);if (i+1<=n&&j-1>0) add(get(i,j)+n*m*2,get(i+1,j-1)+n*m,INF,1);if (i-1>0&&j+1<=m) add(get(i,j)+n*m*2,get(i-1,j+1)+n*m,INF,1);}
}int main()
{scanf("%d%d",&n,&m);memset(st,-1,sizeof(st));char ch[30];int cnt=0;for (int i=1;i<=n;i++){scanf("%s",&ch);for (int j=0;j<m;j++){mp[i][j+1]=ch[j]-'0';if (ch[j]-'0'==1) cnt++;}}for (int i=1;i<=n;i++){scanf("%s",&ch);for (int j=0;j<m;j++){int x=ch[j]-'0';if (x) cnt--;if (mp[i][j+1]==1&&x==0) mp[i][j+1]=1;else if (mp[i][j+1]==0&&x==1) mp[i][j+1]=2;else mp[i][j+1]=0;}}for (int i=1;i<=n;i++){scanf("%s",&ch);for (int j=0;j<m;j++)num[i][j+1]=ch[j]-'0';}if (cnt!=0){printf("-1");return 0;}lianbian();doit();return 0;
}