GYM 101128 F.Landscaping(网络流)

链接：http://codeforces.com/gym/101128

题意：给定一个n*m的地，每块有两种高度，可以花费b将某一块高变低或低变高，有n+m辆车，分别从左往右、从上往下走，经过不同高度的地交界处花费为a，求最小花费。

分析：想到网络流，没想到怎么构造。。。实际上就是要把高和低分开，同时存在高变低和低变高，求最小费用，因此相邻连流量为a的边，源点向低地连流量为b的边，高地向汇点连流量为b的边，最后求一下最小割即可。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
const int delta=50;
const int maxn=2505,inf=1e9+5;
char field[55][55];
struct Edge{
    int from,to,cap,flow;
    Edge(int u,int v,int c,int f):from(u),to(v),cap(c),flow(f){}
};
struct EdmodsKarp{
    int n;
    vector<Edge> edges;
    vector<int> G[maxn];
    int a[maxn];
    int p[maxn];
    void init(int n){
        for(int i=0;i<n;i++)G[i].clear();
        edges.clear();
        this->n=n;
    }
    void AddEdge(int from,int to,int cap){
        edges.push_back(Edge(from,to,cap,0));
        edges.push_back(Edge(to,from,0,0));
        G[from].push_back(edges.size()-2);
        G[to].push_back(edges.size()-1);
    }
    int Maxflow(int s,int t){
        int flow=0;
        while(1){
            memset(a,0,sizeof(a));
            queue<int> Q;
            Q.push(s);
            a[s]=inf;
            while(!Q.empty()){
                int x=Q.front();Q.pop();
                for(int i=0;i<G[x].size();i++){
                    Edge &e=edges[G[x][i]];
                    if(!a[e.to]&&e.cap>e.flow){
                        p[e.to]=G[x][i];
                        a[e.to]=min(a[x],e.cap-e.flow);
                        Q.push(e.to);
                    }
                }
                if(a[t])break;
            }
            if(!a[t])break;
            for(int u=t;u!=s;u=edges[p[u]].from){
                edges[p[u]].flow+=a[t];
                edges[p[u]^1].flow-=a[t];
            }
            flow+=a[t];
        }
        return flow;
    }
}ek;
void add(int s,int t,int cap){
    ek.AddEdge(s,t,cap);
    ek.AddEdge(t,s,cap);
}
int main(){
//    freopen("e:\in.txt","r",stdin);
    int n,m,a,b;
    scanf("%d%d%d%d",&n,&m,&a,&b);
    ek.init(m*n+2);
    for(int i=0;i<n;i++){
        for(int j=0;j<m;j++){
            cin>>field[i][j];
        }
    }
    for(int i=0;i<n;i++){
        for(int j=0;j<m;j++){
            if(i<n-1){
                add(i*delta+j,(i+1)*delta+j,a);
            }
            if(j<m-1){
                add(i*delta+j,i*delta+j+1,a);
            }
            if(field[i][j]=='#')
                ek.AddEdge(i*delta+j,delta*delta+1,b);
            else
                ek.AddEdge(delta*delta,i*delta+j,b);
        }
    }
    printf("%d
",ek.Maxflow(delta*delta,delta*delta+1));
    return 0;
}