Uva-1515-Pool construction

       这个题刚开始确实没想到是用最小割来做，建图也非常的巧妙，

       1、设定源点S,和汇点T.

       2、把周围的洞全部变为草，然后对于周围的每个草建立S到t容量为无穷大的边

       3、对于其他的草，建立S到该点容量为d的边，对于洞建立该点到汇点容量为f的边

       4、对于图中的每个点与它周围建立容量为b的双向边。

      根据最小割最大流定理，跑一遍网络流问题应该就可以解决了。

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
using namespace std;
const int inf=1<<29;
const int maxn=100;
const int maxm=1e5+100;
char map[maxn][maxn];
int w,h,d,f,b,ans,e,st,des,head[3000],nxt[maxm],flow[maxm],pnt[maxm];
int level[3000];
bool vis[3000];
queue<int> q;
void AddEdge(int u,int v,int f)
{
    pnt[e]=v;nxt[e]=head[u];flow[e]=f;head[u]=e++;
    pnt[e]=u;nxt[e]=head[v];flow[e]=0;head[v]=e++;
}
void Init()
{
    for(int i=0;i<w;i++)
    {
        if(map[0][i]=='.')
        {
             ans+=f;
             map[0][i]='#';
        }
        if(map[h-1][i]=='.')
        {
            ans+=f;
            map[h-1][i]='#';
        }
        AddEdge(st,i,inf);
        AddEdge(st,(h-1)*w+i,inf);
    }
    for(int i=0;i<h;i++)
    {
        if(map[i][0]=='.')
        {
            ans+=f;
            map[i][0]='#';
        }
        if(map[i][w-1]=='.')
        {
            ans+=f;
            map[i][w-1]='#';
        }
        AddEdge(st,i*w,inf);
        AddEdge(st,i*w+w-1,inf);
    }
    for(int i=1;i<h-1;i++)
        for(int j=1;j<w-1;j++)
        {
            if(map[i][j]=='#')
                AddEdge(st,i*w+j,d);
            else
                AddEdge(i*w+j,des,f);
        }
    for(int i=0;i<h;i++)
        for(int j=0;j<w;j++)
        {
            if(j>0)
            {
                AddEdge(i*w+j,i*w+j-1,b);
                AddEdge(i*w+j-1,i*w+j,b);
            }
            if(i>0)
            {
                AddEdge(i*w+j,(i-1)*w+j,b);
                AddEdge((i-1)*w+j,i*w+j,b);
            }
        }
}
bool BFS(int st,int des)
{
    memset(level,0,sizeof(level));
    level[st]=1;
    q.push(st);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int i=head[u];i!=-1;i=nxt[i])
        {
            int v=pnt[i];
            if(flow[i]&&!level[v])
            {
                level[v]=level[u]+1;
                q.push(v);
            }
        }
    }
    return level[des];
}
int DFS(int u,int maxf)
{
    if(u==des||!maxf)
        return maxf;
    int ans=0;
    for(int i=head[u],t;i!=-1;i=nxt[i])
        if(level[pnt[i]]==level[u]+1&&(t=DFS(pnt[i],min(flow[i],maxf))))
        {
            flow[i]-=t;
            flow[i^1]+=t;
            return t;
        }
    return level[u]=0;
}
int maxflow()
{
    int ans=0;
    while(BFS(st,des))
        while(1)
        {
            int val=DFS(st,inf);
            if(!val)
                break;
            ans+=val;
        }
    return ans;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        ans=0;e=0;
        memset(head,-1,sizeof(head));
        scanf("%d%d",&w,&h);
        scanf("%d%d%d",&d,&f,&b);
        for(int i=0;i<h;i++)
            scanf("%s",map[i]);
        st=w*h+1,des=w*h+2;
        Init();
        printf("%d\n",maxflow()+ans);
    }
    return 0;
}