bzoj4439[Swerc2015]:Landscaping （最大流、最小割）

Description
秋天到了，有个叫John的农民准备丰收了。由于拖拉机频繁地上山下山，他去年丰收的开支超出了他的预期。

丰收的时候，他的$n+m$台拖拉机需要横向把土地扫一遍，再纵向把土地扫一遍。在下面这幅图中，你可以看到嫩绿色的部分是海拔低的地方，而深绿色的部分是海拔高的地方。收获的时候，他的拖拉机需要上山和下山，上山和下山的位置用红色线标出，这些拖拉机一共要上山或下山8次。

这一年，他准备把一些海拔高的土地铲平变成海拔低的，把一些海拔低的土地填充变成海拔高的来减少拖拉机上山和下山的次数。为了最小化花费，你能帮他决定一下如何铲平/填充使得花费最小吗？

John知道拖拉机每次上山或者下山需要花费 $a$ 金币，而把一块海拔高的土地变成海拔低的土地或者把一块海拔低的土地变成海拔高的土地需要花费 $b$ 金币。

Input
第一行包含四个整数 $n,m,a,b$ ，$n,m$ 表示农场是一个 $n×m$的网格 。
接下来N行，每行M个字符表示John的农场。$′.′$ 表示海拔低的土地，$′\#′$ 表示海拔高的土地。
Output
输出一个整数表示最少花费。

Sample Input
5 4 1000 2000
…#
#..#
…#
##..
###.

Sample Output
11000

Hint
【样例解释】
John有一个5*4的土地，上山下山费用为1000，改变高度费用为2000。
最少花费是11000：将第二行第一个$′\#′$ 改为 $′.′$ ，花费2000，另外上山下山费用为9000。
而不改变任何高度需要花费12000，把所有高土地变低土地需花费18000，把所有低土地变高土地需花费22000。
【数据范围与约定】
对于25%的数据， $n,m≤5$
有10%的数据满足 $n=1$ 或者 $m=1$
对于所有数据， $1≤n,m≤50,1≤a,b≤10^{5}$

题解：最小割

将图放到一个最小割图里面
新建源点 S，汇点T，设更改土地费用为cutcost，登上高地的费用为climbcost。
将S向所有低地连一条容量为cutcost的边，所有高地向T连一条容量为cutcost的边。
接下来所有点向相邻点连一条容量为climbcost的边。
对该图求一遍最大流即可（最大流等于最小割）。

这样做的原因是将S集合与T集合分开求最小割。因为S连向了所有低地，T连向了所有高地，而将一块地的高度变化等价于将该地与相邻的属于不同集合的地的边舍弃。

在上图中，若$\,\,\,cutc(4->t)<climbc(3->4)+climbc(2->4)\,\,\,$则最大流为$\,\,\,cutc(4->t)$
$\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,$若$cutc(4->t)>climbc(3->4)+climbc(2->4)\,\,\,$则最大流为$\,\,\,\,\,climbc(3->4)+climbc(2->4)$
相同集合的1、2、3之间不会产生流量，因为他们同属S集合，若断开的是$cutc(s->2)$（把2变为高地），那么就必须断开$climbc(1->2)$（若断开$climbc(2->4)$则不符合最小割性质（同时断开不是最小割））因此求出来的最大流是原图的最小割。

代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
const int INF=0x3f3f3f3f;
const int Maxn=55;
const int Maxe=Maxn*Maxn*55;
int map[Maxn][Maxn];
char ch[Maxn];
int n,m,cutc,clic,src,des,ecnt=1,ans;
int cap[Maxe],lev[Maxe],cur[Maxe],last[Maxe],before[Maxe],to[Maxe];
const int fx[5]={0,0,0,1,-1};
const int fy[5]={0,1,-1,0,0};
inline void Add(int x,int y,int c)
{
    before[++ecnt]=last[x];
    last[x]=ecnt;
    to[ecnt]=y;
    cap[ecnt]=c;
}
inline void Insert(int x,int y,int z)
{
    Add(x,y,z);
    Add(y,x,0);
}
inline bool Bfs()
{
    static int que[Maxe],qn;
    for(int i=src;i<=des;i++)lev[i]=-1,cur[i]=last[i];
    que[qn=1]=src,lev[src]=0;
    for(int ql=1;ql<=qn;ql++)
    {
        int u=que[ql],v,e;
        for(e=last[u];e,v=to[e];e=before[e])
        {
            if(cap[e]&&lev[v]==-1)
            {
                que[++qn]=v;
                lev[v]=lev[u]+1;
                if(v==des)return true;
            }
        }
    }
    return false;
}
inline int Dinic(int now,int Maxflow)
{
    if(now==des)return Maxflow;
    int v,delta,res=0;
    for(int &e=cur[now];e,v=to[e];e=before[e])
    {
        if(cap[e]&&lev[v]==lev[now]+1)
        {
            delta=Dinic(v,min(Maxflow-res,cap[e]));
            if(delta)
            {
                cap[e]-=delta;
                cap[e^1]+=delta;
                res+=delta;
                if(res==Maxflow)return res;
            }
        }
    }
    if(res!=Maxflow)lev[now]=-1;
    return res;
}
inline void Maxflow()
{
    while(Bfs())
    ans+=Dinic(src,INF);
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&clic,&cutc);
    src=0,des=n*m+1;
    for(int i=1;i<=n;i++)
    {
        scanf("%s",ch+1);
        for(int j=1;j<=m;j++)
        {
            if(ch[j]=='.')map[i][j]=1;
            else map[i][j]=2;
        }
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            if(map[i][j]==1)Insert(src,(i-1)*m+j,cutc);
            else Insert((i-1)*m+j,des,cutc);
            for(int k=1;k<=4;k++)
            {
                int x=i+fx[k],y=j+fy[k];
                if(x>=1&&x<=n&&y>=1&&y<=m)
                {
                    Insert((i-1)*m+j,(x-1)*m+y,clic);
                }
            }
        }
    }
    Maxflow();
    printf("%d\n",ans);
    return 0;
}