UVa1515——最小割最大流+灵活的建图方式

题目来源：https://vjudge.net/problem/UVA-1515

题目思路：按照刘汝佳的思路，最后目的是将草和洞相邻的边都封闭起来，那么就可以用割这个概念，而草变洞，洞变草，草和洞相邻的边构建栅栏都要费用，所以可以用最小割最大流的模型来解决。先将所有边界的洞变为草，这是必须的费用，再建立超级源点和汇点（这个操作在网络流中是常见模型），与源点连通的的是草，与汇点连通的都是洞，则S与草格子节点的容量就是d，表示将这条弧满载时，草才会变为洞，花费为d，但S与边界的草格子连接的容量为INF，表示这个费用为无穷大，也就是不能切这条边；同理，将T汇点与所有洞节点连一条容量为f的边，表示将洞变为草需花费f。最后，相邻的节点之间建一条容量为b的边，表示如果是草和洞相连，这条边需要切断。这样，根据最小割最大流定理，求个S到T的最大流，加上之前必须的花费，就是答案。

总结：还是建模的方法找不到，为什么这样建模有时候还无法理解，看了许久，参考了大佬博客，才能弄个半懂，果然图论题还是要多刷才能有经验，加油！

思路代码均参考大佬博客：https://blog.csdn.net/jijijix/article/details/51941014

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#define mem(a,x) memset(a,x,sizeof(a))
#define Rep(i,x,n)  for(int i= x;i<=n;i++)
using namespace std;
const int maxn = 3000;
const int inf = 100000000;
char mp[55][55];
struct Edge
{
    int from,to,cap,flow;
};
int drx[]={-1,1,0,0};
int dry[]={0,0,-1,1};
struct Dicnic
{
    int n,m,s,t;
    int d[maxn],cur[maxn],vis[maxn];
    vector<Edge>edges;
    vector<int>G[maxn];


    void init(int n)
    {
        this -> n = n;
        for(int i= 0;i<=n;i++)
            G[i].clear();
        edges.clear();
    }


    void AddEdge(int from,int to,int cap)
    {
        Edge e1 = {from,to,cap,0};
        Edge e2 = {to,from,0,0};
        edges.push_back(e1);
        edges.push_back(e2);
        int m = edges.size();
        G[from].push_back(m-2);
        G[to].push_back(m-1);
    }


    bool BFS()
    {
        mem(vis,0);
        queue<int >Q;
        vis[s] = 1;
        d[s] = 0;
        Q.push(s);
        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(!vis[e.to]&&e.cap>e.flow){
                    d[e.to] = d[x] + 1;
                    vis[e.to] = 1;
                    Q.push(e.to);
                }
            }
        }
        return vis[t];
    }


    int DFS(int x,int a)
    {
        if(x==t||a==0)  return a;
        int f,flow = 0;
        for(int &i= cur[x];i<G[x].size();i++){
            Edge &e = edges[G[x][i]];
            if(d[e.to] == d[x]+1&&(f=(DFS(e.to,min(e.cap-e.flow,a))))>0){
                flow += f;
                a -= f;
                e.flow += f;
                edges[G[x][i]^1].flow -= f;
                if(a==0)    break;
            }
        }
        return flow;
    }


    int MaxFlow(int s,int t)
    {
        this -> s = s;  this -> t = t;
        int flow = 0;
        while(BFS()){
            mem(cur,0);
            flow += DFS(s,inf);
        }
        return flow;
    }
};
Dicnic g;
int main()
{
    int T,w,h,d,f,b;
    scanf("%d",&T);
    while(T--){
        scanf("%d %d %d %d %d",&w,&h,&d,&f,&b);
        for(int i = 1;i<= h;i++){
            scanf("%s",mp[i]+1);
        }
        int ans = 0;
        for(int i = 1;i<=w;i++){
            if(mp[1][i]=='.'){
                ans += f;
                mp[1][i]='#';
            }
            if(mp[h][i]=='.'){
                ans += f;
                mp[h][i] = '#';
            }
        }
        for(int i= 2;i<=h-1;i++){
            if(mp[i][1]=='.'){
                ans += f;
                mp[i][1] = '#';
            }
            if(mp[i][w]=='.'){
                ans += f;
                mp[i][w] = '#';
            }
        }                                   ///啰里啰嗦一大段为的就是把四周的洞变成草，并且加上变成草的费用；
        g.init(w*h+1);
        for(int i=1;i<=h;i++){
            for(int j= 1;j<=w;j++){
                int id1 = w*(i-1) + j , id2;
                if(mp[i][j]=='#'&&i!=1&&i!=h&&j!=1&&j!=w)   g.AddEdge(0,id1,d);
                else if(mp[i][j]=='#'&&(i==1||i==h||j==1||j==w))   g.AddEdge(0,id1,inf);  ///这个判断十分重要；
                else                g.AddEdge(id1,w*h+1,f);
                for(int k = 0;k<4;k++){
                    int x = i+drx[k],y = j+dry[k];
                    if(x>=1&&x<=h&&y>=1&&y<=w){
                        id2 = w*(x-1)+y;
                        g.AddEdge(id1,id2,b);
                    }
                }
            }
        }
        ans += g.MaxFlow(0,w*h+1);
        printf("%d\n",ans);
    }
    return 0;
}