【最小割】[JZOJ4020] Revolution

Description

地图是个n*m的网格。
选择投资格子(i,j)会消耗a[i][j]的金钱
如果一个格子被投资，或者四连通的格子都被投资（边界可能不需要四个），就可以获得该格子的收益。
求收益-花费的最大值，当然也可以不投资。
第一行两个数 n,m(n,m ≤ 20)，表示矩形的长和宽。
接下来 n 行，每行是 m 个字符组成的字符串，描述投资的花费。
接下来 n 行，每行是 m 个字符组成的字符串，表示该格子的收益。
花费和收益按照一种奇葩的方式给出：
字符 数
‘0’ -’ 9’ 0-9
‘a’ -’ z’ 10-35
‘A’ -’ Z’ 36-61

Solution

典型的最小割问题
考虑建图

连边的时候为了避免负权边，我们用总代价减去最小割，割掉花费说明投资了这个格子，割掉代价说明放弃了这个格子

考虑第一种关系，投资这个格子即有收益
那么应该是这样的

要么割cost，要么割profit（cost和profit反过来也可以）

第二种与第一种是或的关系，不会既选了这个格子，又选了它四相邻的所有格子（明显不优）

如图，上图表示了选了上面那个点代表的格子，下面那个点的格子就会有收益
那么两种关系综合起来，就是

我们将一个格子拆成两个点in和out，令源点为st，汇点为ed

选了i，那么割掉cost[i]，这时发现profit[j]也保留了下来
选了j，那么割掉cost[j]，profit[i]也不用被割掉

回到原题目，题目要求是四连通的点都要选，那么我们将整个矩阵黑白染色，规定INF边只从白店连向黑点
白点是上图i的方式连边，黑点是上图j的方式连边，一个白点需要对周围四个黑点如上图连INF的边
这样就把图够好了，直接套用最大流算法即可

Code

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fod(i,a,b) for(int i=a;i>=b;i--)
#define N 23
#define M 805
#define INF 1802201963
using namespace std;
int m1,n,m,f[M][M],nt[10*M],dt[10*M],fs[M],in[N][N],out[N][N],cost[N][N],pr[N][N],n1,st,ed,fx[4][2]={{0,1},{1,0},{0,-1},{-1,0}},h[M],d[M];
void link(int x,int y,int z)
{
    f[x][y]+=z;
    nt[++m1]=fs[x];
    dt[fs[x]=m1]=y;
}
bool bfs()
{
    memset(h,0,sizeof(h));
    int l=0,r=1;
    h[d[1]=st]=1;
    while(l<r)
    {
        int k=d[++l];
        for(int i=fs[k];i;i=nt[i])
        {
            int p=dt[i];
            if(!h[p]&&f[k][p]) h[p]=h[k]+1,d[++r]=p;
        }
    }
    return (h[ed]>0);
}
int dinic(int k,int s)
{
    if(k==ed) return s;
    int sl=0;
    for(int i=fs[k];i;i=nt[i])
    {
        int p=dt[i];
        if(h[p]==h[k]+1&&f[k][p])
        {
            int v=dinic(p,min(s,f[k][p]));
            if(v)
            {
                sl+=v,s-=v;
                f[k][p]-=v,f[p][k]+=v;
                if(!s) break;
            }
        }   
    }
    if(sl==0) h[k]=-1;
    return sl;
}
int main()
{
    cin>>n>>m;
    st=1,ed=2,n1=2;
    fo(i,1,n)
    {
        scanf("\n");
        fo(j,1,m)
        {
            char ch=getchar();
            if(ch>='0'&&ch<='9') cost[i][j]=ch-'0';
            else if(ch>='a'&&ch<='z') cost[i][j]=ch-'a'+10;
            else cost[i][j]=ch-'A'+36;
        }
    }
    int s=0;
    fo(i,1,n)
    {
        scanf("\n");
        fo(j,1,m)
        {
            char ch=getchar();
            if(ch>='0'&&ch<='9') pr[i][j]=ch-'0';
            else if(ch>='a'&&ch<='z') pr[i][j]=ch-'a'+10;
            else pr[i][j]=ch-'A'+36;
            s+=pr[i][j];
            in[i][j]=++n1,out[i][j]=++n1;
            link(in[i][j],out[i][j],pr[i][j]),link(out[i][j],in[i][j],0);
            if((i+j)%2==0) link(st,in[i][j],cost[i][j]),link(in[i][j],st,0);
            else link(out[i][j],ed,cost[i][j]),link(ed,out[i][j],0);
        }
    }
    fo(i,1,n) 
    {
        fo(j,1,m)
        {
            if((i+j)%2==0)
            {
                fo(k,0,3)
                {
                    int x=i+fx[k][0],y=j+fx[k][1];
                    if(x&&y&&x<=n&&y<=m) link(in[i][j],in[x][y],INF),link(in[x][y],in[i][j],0),link(out[i][j],out[x][y],INF),link(out[x][y],out[i][j],0);
                }
            }
        }
    }
    n++;
    n--;
    while(bfs()) 
        s-=dinic(st,INF);
    printf("%d\n",s);
}