bzoj2668 [cqoi2012]交换棋子（费用流）

可以只看做把黑点移动到目标位置。对每个点的交换次数有限制。我们容易想到拆点。但是这题比较妙的地方在于拆成了三个点。把一个黑点换到目标位置上的一条路径上，除了起终点交换了一次之外，其他点都交换了两次。我们如何体现这个差异呢？拆成三个点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;
}