[bzoj2668][网络流-费用流]交换棋子

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

题解

这题超强
被波老师带着做题
如果没有这个交换次数限制的话，那就超级好做对不对。st向起始位置为黑色的点，结束位置为黑色的点向ed连边，跑费用流验证
我们考虑有了怎么办
一个棋子交换到另一个位置的路径上，只有起点和终点只交换了一次，其他均交换了两次，我们要体现这种差异
那么每个位置拆点x1,x0,x2，设限制次数为c
如果这个位置的起始与结束状态相同那么
x1->x0，流量c/2，x0->x2，流量c/2，费用均为1
如果只有起始位置是黑色
x1->x0，流量c/2，费用1，x0->x2，流量c/2+1，费用1，st->x0 流量1 费用0
如果只有结束位置是黑色
x1->x0，流量c/2+1，费用1，x0->x2，流量c/2，费用1，x0->ed 流量1 费用0
再跑费用流验证即可

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
const int dx[8]={-1,-1,0,1,1,1,0,-1};
const int dy[8]={0,-1,-1,-1,0,1,1,1};
struct node
{
    int x,y,c,d,next,other;
}a[210000];int len,last[1500];
void ins(int x,int y,int c,int d)
{
    int k1,k2;
    k1=++len;
    a[len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;
    a[len].next=last[x];last[x]=len;
    k2=++len;
    a[len].x=y;a[len].y=x;a[len].c=0;a[len].d=-d;
    a[len].next=last[y];last[y]=len;
    a[k1].other=k2;a[k2].other=k1;
}
bool v[1500];
int d[1500],pre[1500],tmp[1500];
int list[110000],head,tail,n,m,st,ed;
bool spfa()
{
    memset(d,63,sizeof(d));d[st]=0;
    memset(v,false,sizeof(v));v[st]=true;
    list[1]=st;head=1;tail=2;
    while(head!=tail)
    {
        int x=list[head];
        for(int k=last[x];k;k=a[k].next)
        {
            int y=a[k].y;
            if(d[y]>d[x]+a[k].d && a[k].c>0)
            {
                d[y]=d[x]+a[k].d;
                pre[y]=x;tmp[y]=k;
                if(v[y]==false)
                {
                    v[y]=true;
                    list[tail++]=y;
                    if(tail==ed+1)tail=1;   
                }
            }
        }
        head++;if(head==ed+1)head=1;
        v[x]=false;
    }
    if(d[ed]<1061109567)return true;
    return false;
}
int sum;
int mcf()
{
    int ret=0,cnt=0,minn;bool bk=false;
    while(spfa())
    {
        int x=ed;minn=999999999;
        while(x!=st)
        {
            minn=min(minn,a[tmp[x]].c);
            x=pre[x];
        }
        ret+=d[ed]*minn;
        cnt+=minn;if(cnt==sum){bk=true;}
        x=ed;
        while(x!=st)
        {
            a[tmp[x]].c-=minn;a[a[tmp[x]].other].c+=minn;
            x=pre[x];
        }
    }
    if(bk==false)return -1;
    return ret;
}
char p[35][35],q[35][35],ch[35];
int pt(int x,int y){return (x-1)*m+y;}
bool chk(int x,int y)
{
    if(x<1 || x>n || y<1 || y>m)return false;
    return true;
}
int main()
{
    scanf("%d%d",&n,&m);
    st=n*m*3+1;ed=n*m*3+2;
    for(int i=1;i<=n;i++)scanf("%s",p[i]+1);
    for(int i=1;i<=n;i++)scanf("%s",q[i]+1);
    sum=0;
    for(int i=1;i<=n;i++)
    {
        scanf("%s",ch+1);
        for(int j=1;j<=m;j++)
        {
            if(q[i][j]==p[i][j])ins(pt(i,j),pt(i,j)+n*m,(ch[j]-'0')/2,1),ins(pt(i,j)+n*m,pt(i,j)+n*m*2,(ch[j]-'0')/2,1);
            else if(p[i][j]=='1')ins(pt(i,j),pt(i,j)+n*m,(ch[j]-'0')/2,1),ins(pt(i,j)+n*m,pt(i,j)+n*m*2,(ch[j]-'0')/2+1,1);
            else ins(pt(i,j),pt(i,j)+n*m,(ch[j]-'0')/2+1,1),ins(pt(i,j)+n*m,pt(i,j)+n*m*2,(ch[j]-'0')/2,1);
        }
    }
    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)for(int k=0;k<=7;k++)
        if(chk(i+dx[k],j+dy[k]))ins(pt(i,j)+n*m*2,pt(i+dx[k],j+dy[k]),999999999,0);

    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)if(p[i][j]=='1')ins(st,pt(i,j)+n*m,1,0),sum++;
    for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)if(q[i][j]=='1')ins(pt(i,j)+n*m,ed,1,0);
    int tmp=mcf();
    if(tmp!=-1)printf("%d\n",tmp/2);
    else printf("-1\n");
    return 0;
}