【网络流】- LA5905-Pool construction

链接：我是链接啊啊啊哈哈哈

题意：

 泳池最外一圈必须全搞成#

对其他的点有三种处理方式

1.放着不动

2.花费 d 把 # 变成 .

3.花费 f 把 . 变成 #

最后 # 和 . 不能相邻，要在他们之间修围墙

围墙单位边造价是b

问最小总花费

图最后变成什么样都行,只要能让花费最小 

思路：

一看是修围墙，那就网络流呗，最小割啥的，好吧，不会建图

网络流难就难在建图。想了一阵没想出来

看答案去了。

这样建图

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

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

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

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

嘛。。也许见多了就会建图了。。学习着吧

代码：

#include<cstdio>#include<cstring>#include<iostream>using namespace std;const int maxn = 3000;const int maxm = 60000;const int INF = 0x7FFFFFF;#define mem(name,value) memset((name),(value),sizeof(name))struct Side{    int to,next,c;}side[maxm];int top,node[maxn];void add_side(int u,int v,int c,int rc){    //cout<<u<<" --> "<<v<<" need "<<c<<endl;    side[top]=(Side){v,node[u],c};node[u]=top++;    side[top]=(Side){u,node[v],rc};node[v]=top++;}int start,end,cnt,dis[maxn],gap[maxn];int get_flow(int u,int flow){    //printf("%d %d\n",u,flow);    if(u==end)return flow;    int ans=0;    for(int i=node[u];i!=-1;i=side[i].next){        int v=side[i].to,c=side[i].c;        if(dis[u]>dis[v]&&c){            int f=get_flow(v,min(flow-ans,c));            ans+=f;            side[i].c-=f;            side[i^1].c+=f;            if(ans==flow)return ans;        }    }    if(!(--gap[dis[u]]))dis[start]=cnt+2;    gap[++dis[u]]++;    return ans;}char grid[55][55];int fans,f,d,b,cnn;int n,m;int getid(int x,int y){    return x*n+y+1;}void deal_ini(){    //cout<<"****"<<endl;    int i,j;    int id;    //cout<<"n="<<n<<endl;    for(i=0;i<n;i++)    {        //cout<<"****"<<endl;        if(grid[0][i] == '.') fans += f;        id = getid(0,i);        //cout<<0<<" --> "<<id<<" need "<<INF<<endl;        add_side(0,id,INF,0);        if(grid[m-1][i] == '.') fans += f;        id = getid(m-1,i);        add_side(0,id,INF,0);    }    for(i=1;i<m-1;i++)    {        if(grid[i][0] == '.') fans += f;        id = getid(i,0);        add_side(0,id,INF,0);        if(grid[i][n-1] == '.') fans += f;        id = getid(i,n-1);        add_side(0,id,INF,0);    }}void makemap(){    int i,j,id;    for(i=1;i<m-1;i++)    {        for(j=1;j<n-1;j++)        {            id = getid(i,j);            if(grid[i][j] == '#')                add_side(0,id,d,0);            else                add_side(id,m*n+1,f,0);        }    }    for(i=1;i<m;i++)    {        for(j=0;j<n;j++)        {            int id1 = getid(i,j);            int id2 = getid(i-1,j);            add_side(id1,id2,b,0);            add_side(id2,id1,b,0);        }    }    for(j=1;j<n;j++)    {        for(i=0;i<m;i++)        {            int id1 = getid(i,j);            int id2 = getid(i,j-1);            add_side(id1,id2,b,0);            add_side(id2,id1,b,0);        }    }}int main(){    //freopen("input.txt","r",stdin);    int T;    scanf("%d",&T);    while(T--)    {        fans=0;        cnn=1;        memset(grid,0,sizeof(grid));        scanf("%d%d",&n,&m);        scanf("%d%d%d",&d,&f,&b);        for(int i=0;i<m;i++)            scanf("%s",grid[i]);        /*for(int i=0;i<m;i++)            cout<<grid[i]<<endl;*/        top=0;        mem(node,-1);        mem(gap,0);        mem(dis,0);        deal_ini();        makemap();        int ans=0;        start=0;        end=n*m+1;        cnt=n*m+2;        gap[0]=cnt;        while(dis[start]<cnt)ans+=get_flow(start,INF);        printf("%d\n",ans+fans);    }    return 0;}