BZOJ_4439_[Swerc2015]Landscaping_最小割
Description
FJ有一块N*M的矩形田地,有两种地形高地(用‘#’表示)和低地(用‘.’表示)
FJ需要对每一行田地从左到右完整开收割机走到头,再对每一列从上到下完整走到头,如下图所示
对于一个4*4的田地,FJ需要走8次。
收割机是要油的,每次从高地到低地或从低地到高地需要支付A的费用。
但是FJ有黑科技,可以高地与低地的互变,都只需要一个支付B的费用。
询问FJ需要支付最小费用。
Input
第一行包含四个整数N,M,A,B,意义如上文所述。
接下来是一个N*M的字符串矩阵,表示农田的地形,’#’表示高地,’.’表示低地。
Output
只包含一个正整数,表示最小费用。
1<=N,M<=50
1<=A,B<=100000
1<=A,B<=100000
Sample Input
5 4 1000 2000
...#
#..#
...#
##..
###.
...#
#..#
...#
##..
###.
Sample Output
11000
样例解释:
把(2,1)的高地变成低地花费2000,燃料花费9000
用最小割的思想。
S->高地(B) 低地->T(B),割这些边表示一开始将高低互换。
然后对于每个点向四周连边(A),割这个表示支付A从而改变与四周的联系,可以理解为边是双向的。
求最小割即可。
代码:
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define N 3050
#define M 300050
#define S (n*m+1)
#define T (n*m+2)
#define inf 100000000
#define p(i,j) ((i-1)*m+j)
int head[N],to[M],nxt[M],flow[M],cnt=1,n,m,dep[N],Q[N],l,r;
inline void add(int u,int v,int f) {
to[++cnt]=v; nxt[cnt]=head[u]; head[u]=cnt; flow[cnt]=f;
to[++cnt]=u; nxt[cnt]=head[v]; head[v]=cnt; flow[cnt]=0;
}
bool bfs() {
memset(dep,0,sizeof(dep));
dep[S]=1; l=r=0; Q[r++]=S;
while(l<r) {
int x=Q[l++],i;
for(i=head[x];i;i=nxt[i]) {
if(!dep[to[i]]&&flow[i]) {
dep[to[i]]=dep[x]+1;
if(to[i]==T) return 1;
Q[r++]=to[i];
}
}
}
return 0;
}
int dfs(int x,int mf) {
if(x==T) return mf;
int nf=0,i;
for(i=head[x];i;i=nxt[i]) {
if(dep[to[i]]==dep[x]+1&&flow[i]) {
int tmp=dfs(to[i],min(mf-nf,flow[i]));
if(!tmp) dep[to[i]]=0;
nf+=tmp;
flow[i]-=tmp;
flow[i^1]+=tmp;
if(nf==mf) break;
}
}
return nf;
}
void dinic() {
int ans=0,f;
while(bfs()) while(f=dfs(S,inf)) ans+=f;
printf("%d\n",ans);
}
char s[60];
int main() {
int A,B,i,j;
scanf("%d%d%d%d",&n,&m,&A,&B);
for(i=1;i<=n;i++) {
scanf("%s",s+1);
for(j=1;j<=m;j++) {
if(s[j]=='#') {
add(S,p(i,j),B);
}else {
add(p(i,j),T,B);
}
}
}
for(i=1;i<=n;i++) {
for(j=1;j<=m;j++) {
if(i>1) add(p(i,j),p(i-1,j),A);
if(i<n) add(p(i,j),p(i+1,j),A);
if(j>1) add(p(i,j),p(i,j-1),A);
if(j<m) add(p(i,j),p(i,j+1),A);
}
}
dinic();
}