可以只看做把黑点移动到目标位置。对每个点的交换次数有限制。我们容易想到拆点。但是这题比较妙的地方在于拆成了三个点。把一个黑点换到目标位置上的一条路径上,除了起终点交换了一次之外,其他点都交换了两次。我们如何体现这个差异呢?拆成三个点x1,x0,x2。
如果起始和目标状态相同,则建边x1->x0,c/2,0 x0->x2,c/2,0
如果只有起始是黑的 ,则建边x1->x0,c/2,0 x0->x2,(c+1)/2,0 S->x0,1,0
如果只有目标是黑的 ,则建边x1->x0,(c+1)/2,0 x0->x2,c/2,0 x0->T,1,0
对于所有联通的格子x,y 建边 x2->y1,inf,1.
我们发现从S->T的一个流就对应了这样的一条路径。最小费用就是答案。
#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 1220
inline int read(){
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
return x*f;
}
int n,m,id[21][21],h[N],num=1,T=1210,ans=0,mxflow=0,path[N],dis[N];
int dx[8]={0,0,1,1,1,-1,-1,-1},dy[8]={1,-1,-1,0,1,-1,0,1};bool inq[N];
struct edge{
int to,next,w,c;
}data[10000];
char mp1[30][30],mp2[30][30],mp3[30][30];
inline void add(int x,int y,int w,int c){
data[++num].to=y;data[num].next=h[x];h[x]=num;data[num].w=w;data[num].c=c;
data[++num].to=x;data[num].next=h[y];h[y]=num;data[num].w=0;data[num].c=-c;
}
inline bool spfa(){
deque<int>q;memset(path,0,sizeof(path));memset(dis,inf,sizeof(dis));
q.push_back(0);dis[0]=0;inq[0]=1;
while(!q.empty()){
int x=q.front();q.pop_front();inq[x]=0;
for(int i=h[x];i;i=data[i].next){
int y=data[i].to;if(!data[i].w) continue;
if(data[i].c+dis[x]<dis[y]){
dis[y]=dis[x]+data[i].c;path[y]=i;
if(!inq[y]){
if(!q.empty()&&dis[q.front()]>dis[y]) q.push_front(y);
else q.push_back(y);inq[y]=1;
}
}
}
}return path[T];
}
int main(){
// freopen("a.in","r",stdin);
n=read();m=read();int tot=0,tot1=0,tot2=0;
for(int i=1;i<=n;++i)
for(int j=1;j<=m;++j) id[i][j]=++tot;
for(int i=1;i<=n;++i) scanf("%s",mp1[i]+1);
for(int i=1;i<=n;++i) scanf("%s",mp2[i]+1);
for(int i=1;i<=n;++i) scanf("%s",mp3[i]+1);
for(int i=1;i<=n;++i)
for(int j=1;j<=m;++j){
if(mp1[i][j]==mp2[i][j]){
add(id[i][j]+tot,id[i][j],mp3[i][j]-'0'>>1,0);
add(id[i][j],id[i][j]+tot*2,mp3[i][j]-'0'>>1,0);continue;
}if(mp1[i][j]=='1'){
add(0,id[i][j],1,0),++tot1;
add(id[i][j]+tot,id[i][j],mp3[i][j]-'0'>>1,0);
add(id[i][j],id[i][j]+tot*2,mp3[i][j]-'0'+1>>1,0);
}if(mp2[i][j]=='1'){
add(id[i][j],T,1,0),++tot2;
add(id[i][j]+tot,id[i][j],mp3[i][j]-'0'+1>>1,0);
add(id[i][j],id[i][j]+tot*2,mp3[i][j]-'0'>>1,0);
}
}if(tot2!=tot1){puts("-1");return 0;}
for(int x=1;x<=n;++x)
for(int y=1;y<=m;++y)
for(int k=0;k<8;++k){
int xx=x+dx[k],yy=y+dy[k];
if(xx<1||xx>n||yy<1||yy>m) continue;
add(id[x][y]+tot*2,id[xx][yy]+tot,inf,1);
}
while(spfa()){
int low=inf,now=T;
while(path[now]) low=min(low,data[path[now]].w),now=data[path[now]^1].to;
ans+=low*dis[T];mxflow+=low;now=T;
while(path[now]) data[path[now]].w-=low,data[path[now]^1].w+=low,now=data[path[now]^1].to;
}if(mxflow!=tot1) puts("-1");
else printf("%d\n",ans);
return 0;
}