[bzoj2668] [cqoi2012]交换棋子

　　到现在还不是很懂。网络流这块基础太弱了QAQ

　　一开始先把初始状态黑白棋子数目不同的情况去掉。

　　然后只要关注黑色棋子就好了。

　　每个棋子拆点，入点往出点流量限制是m(i,j)/2，费用为2。但如果这个棋子入点是黑或者出点是黑的话就加1/2

　　棋子初始状态为黑的话，就从S往它入点连一条流量1，费用-1的边，结束状态为黑的话，就从它出点往T连一条同样的边。

　　每个棋子从出点往相邻棋子的入点连一条流量正无穷，费用为0的边。

　　最小交换次数就是最小费用的1/2

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=23*23*2,inf=1002333333;
const int xx[8]={-1,-1,-1,0,0,1,1,1},yy[8]={-1,0,1,-1,1,-1,0,1};
struct zs{int too,pre,flow,dis;}e[233333];int tot,last[maxn];
int dis[maxn],dl[1002333];
int id[23][23],mx[maxn];
bool u[maxn],ins[maxn],st[maxn],ed[maxn];
int i,j,k,n,m,s,t,ans;
char ch[23];
 
int ra;char rx;
inline int read(){
    rx=getchar(),ra=0;
    while(rx<'0'||rx>'9')rx=getchar();
    while(rx>='0'&&rx<='9')ra*=10,ra+=rx-48,rx=getchar();return ra;
}
inline bool spfa(){
    memset(dis,60,(t+1)<<2);
    int l=0,r=1,i,now;dl[1]=s,dis[s]=0,u[s]=1;
    while(l<r)
        for(i=last[now=dl[++l]],u[now]=0;i;i=e[i].pre)if(e[i].flow&&dis[e[i].too]>dis[now]+e[i].dis){
            dis[e[i].too]=dis[now]+e[i].dis;
            if(!u[e[i].too])dl[++r]=e[i].too,u[e[i].too]=1;
        }
//  for(i=1;i<=t;i++)printf("0-->%d  %d
",i,dis[i]);
    return dis[t]<inf;
}
inline int min(int a,int b){return a<b?a:b;}
int dfs(int x,int mx){
    if(x==t)return mx;
    int i,used=0,w;ins[x]=1;
    for(i=last[x];i;i=e[i].pre)if(e[i].flow&&dis[e[i].too]==dis[x]+e[i].dis&&!ins[e[i].too]){
        w=dfs(e[i].too,min(mx-used,e[i].flow));if(w){
            e[i].flow-=w,e[i^1].flow+=w,used+=w,ans+=e[i].dis*w;
            if(used==mx){ins[x]=0;return mx;}
        }
    }
    dis[x]=inf,ins[x]=0;return used;
}
inline void insert(int a,int b,int c,int d){
//  printf("   %d-->%d %d %d
",a,b,c,d);
    e[++tot].too=b,e[tot].flow=c,e[tot].dis= d,e[tot].pre=last[a],last[a]=tot,
    e[++tot].too=a,e[tot].flow=0,e[tot].dis=-d,e[tot].pre=last[b],last[b]=tot;
}
 
int main(){
    n=read(),m=read();int num1=0,num2=0,tmp=0;
    for(i=1;i<=n;i++)for(j=1;j<=m;j++)id[i][j]=++tmp;
    for(i=1;i<=n;i++){
        scanf("%s",ch+1);
        for(j=1;j<=m;j++)
            if(ch[j]=='1')st[id[i][j]]=1,num1++;
    }
    for(i=1;i<=n;i++){
        scanf("%s",ch+1);
        for(j=1;j<=m;j++)
            if(ch[j]=='1')ed[id[i][j]]=1,num2++;
    }
    if(num1!=num2){puts("-1");return 0;}
    for(i=1;i<=n;i++){
        scanf("%s",ch+1);
        for(j=1;j<=m;j++)
            mx[id[i][j]]=ch[j]-'0';
    }
    s=0,t=tmp<<1|1,tot=1;
    for(i=1;i<=n;i++)for(j=1;j<=m;j++){
        int now=id[i][j],x,y;
        if(st[now])insert(s,now,1,-1);
        if(ed[now])insert(now+tmp,t,1,-1);
        insert(now,now+tmp,(mx[now]+st[now]+ed[now])>>1,2);
        for(k=0;k<8;k++){
            x=i+xx[k],y=j+yy[k];
            if(x<1||y<1||x>n||y>m)continue;
            insert(now+tmp,id[x][y],inf,0);
        }
    }int flow=0;
    while(spfa())flow+=dfs(s,inf);
    if(flow<num1)puts("-1");else printf("%d
",ans>>1);
}