[CQOI2012]交换棋子

传送门
此博客无讲解，讲解可以看hzwer
这个题还是比较巧妙的，把点权转换成边权，这个没什么，常用技巧。
主要觉得比较妙的是那个+1然后/2，这一步感觉很妙，而且这次拆成3个点，很妙。
（其实我这个博客就是来吐槽一波。。）
代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstdlib>
#define ll long long
using namespace std;
inline int read(){
    int x=0;char ch=' ';int f=1;
    while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
    if(ch=='-')f=-1,ch=getchar();
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
const int N=1251,M=N*50;
const int B[8][2]={{0,1},{1,0},{0,-1},{-1,0},{1,1},{-1,-1},{1,-1},{-1,1}};
int n,m,s,t,tot=-1;
char ch[21];
int a[21][21],b[21][21],tm[21][21];
int head[N],cur[N],mark[N],inq[N],d[N],q[M],to[M],Next[M],flow[M],cost[M];
inline void addedge(int x,int y,int l,int c){
    to[++tot]=y;Next[tot]=head[x];flow[tot]=l;cost[tot]=c;head[x]=tot;
    to[++tot]=x;Next[tot]=head[y];flow[tot]=0;cost[tot]=-c;head[y]=tot;
}
inline bool spfa(){
    for(int i=s;i<=t;i++)inq[i]=0,d[i]=0x3f3f3f3f;
    int l=M/2,r=M/2;q[l]=s;d[s]=0;inq[s]=1;
    while(l<=r){
        int x=q[l++];
        for(int i=head[x];i!=-1;i=Next[i]){
            int u=to[i];
            if(flow[i]&&d[u]>d[x]+cost[i]){
                d[u]=d[x]+cost[i];
                if(!inq[u]){
                    inq[u]=1;
                    if(d[u]<d[q[l]])q[--l]=u;
                    else q[++r]=u;
                }
            }
        }
        inq[x]=0;
    }
    return d[t]!=0x3f3f3f3f;
}
int C;
inline int dfs(int x,int a){
    mark[x]=1;
    if(x==t||!a)return a;
    int F=0,f;
    for(int &i=cur[x];i!=-1;i=Next[i]){
        int u=to[i];
        if(!mark[u]&&flow[i]&&d[u]==d[x]+cost[i]&&(f=dfs(u,min(a,flow[i])))>0){
            flow[i]-=f;
            flow[i^1]+=f;
            F+=f;
            a-=f;
            C+=f*cost[i];
            if(!a)return F;
        }
    }
    return F;
}
inline int MCMF(){
    int F=0;
    while(spfa()){
        mark[t]=1;
        while(mark[t]){
            for(int i=s;i<=t;i++)cur[i]=head[i],mark[i]=0;
            F+=dfs(s,0x3f3f3f3f);
        }
    }
    return F;
}
int main(){
    memset(head,-1,sizeof(head));
    n=read();m=read();s=0;t=n*m*3+1;
    int all=n*m;int sum1=0,sum2=0,sum=0;
    for(int i=1;i<=n;i++){
        scanf("%s",ch+1);
        for(int j=1;j<=m;j++){
            a[i][j]=ch[j]-'0';
            if(a[i][j])sum1++;
        }
    }
    for(int i=1;i<=n;i++){
        scanf("%s",ch+1);
        for(int j=1;j<=m;j++){
            b[i][j]=ch[j]-'0';
            if(b[i][j])sum2++;
            if(a[i][j]&&b[i][j])sum++;
        }
    }
    for(int i=1;i<=n;i++){
        scanf("%s",ch+1);
        for(int j=1;j<=m;j++){
            tm[i][j]=ch[j]-'0';
        }
    }
    if(sum1!=sum2){printf("-1");return 0;}
    int add1,add2,tx,ty,id;
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            id=(i-1)*m+j;
            if(a[i][j])addedge(s,id+all,1,0);
            if(b[i][j])addedge(id+all,t,1,0);
            add1=add2=0;
            if(!a[i][j]&&b[i][j])add1=1;
            if(a[i][j]&&!b[i][j])add2=1;
            addedge(id,id+all,(tm[i][j]+add1)/2,1);
            addedge(id+all,id+all+all,(tm[i][j]+add2)/2,0);

            for(int pos=0;pos<8;pos++){
                tx=i+B[pos][0];ty=j+B[pos][1];
                if(tx<1||tx>n||ty<1||ty>m)continue;
                addedge(id+all+all,(tx-1)*m+ty,0x3f3f3f3f,0);
            }
        }
    }
    MCMF();
    printf("%d",C);
    return 0;
}