Uva1515 Pool construction

Time Limit: 3000MS64bit IO Format: %lld & %llu

网络流 最小割

心生绝望，用了好久的网络流模板居然是错的。

↑居然之前还侥幸能过一堆(并不)题。

______

利用边容量来限制花费：

　　边缘的坑必须填上，每格花费为f；

　　如果一个地方id是坑，将其填上，从S到id连边，容量为f；

　　如果一个地方id是草，在这挖坑，从id到T连边，容量为d；

　　建立栅栏：每格向四周连边，容量为b；

↑如此建边起到了限制最小花费的作用，如果挖坑的代价比种草小，挖坑的边先流满，反之同理，求出的最小割就是达成要求的最小花费。

求最小割即可。

第53行 if(!f)dep[u]=0; 作用是当前路径无法扩展时，切断路径，避免重复无意义DFS。

↑以前写dinic一直不知道要加这个，居然没被卡……

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<queue>
 using namespace std;
 const int INF=<<;
 const int mxn=;
 const int mx[]={,,,-,};
 const int my[]={,,,,-};
 struct edge{
     int v,nxt,f;
 }e[];
 int hd[mxn],mct=;
 int n,m,k,ans;
 int id[][];
 int S,T;
 void add_edge(int u,int v,int c){
     e[++mct].v=v;e[mct].nxt=hd[u];e[mct].f=c;hd[u]=mct;return;
 }
 int dep[];
 bool BFS(){
     queue<int>q;
     memset(dep,,sizeof dep);
     dep[S]=;
     q.push(S);
     while(!q.empty()){
         int u=q.front();q.pop();
         for(int i=hd[u];i;i=e[i].nxt){
             int v=e[i].v;
             if(!dep[v] && e[i].f){
                 dep[v]=dep[u]+;
                 q.push(v);
             }
         }
     }
     return dep[T];
 }
 int DFS(int u,int lim){
     if(u==T) return lim;
     int tmp,f=;
     for(int i=hd[u];i;i=e[i].nxt){
         int v=e[i].v;
         if(dep[v]==dep[u]+ && e[i].f){
             tmp=DFS(v,min(lim,e[i].f));
             lim-=tmp;
             f+=tmp;
             e[i].f-=tmp;
             e[i^].f+=tmp;
             if(!lim)return f;
         }
     }
     if(!f)dep[u]=;
     return f;
 }
 int Dinic(){
     int res=;
     while(BFS())res+=DFS(S,1e9);
     return res;
 }
 inline void init(){
     memset(hd,,sizeof hd);
     for(int i=;i<=n;i++)
      for(int j=;j<=m;j++)
          id[i][j]=(i-)*m+j;
     mct=;
     ans=;
     return;
 }
 int d,f,b;
 char mp[][];
 int main()
 {
     int cas;
     scanf("%d",&cas);
     int i,j;
     while(cas--){
         scanf("%d%d",&m,&n);
         scanf("%d%d%d",&d,&f,&b);
         init();
         S=;T=n*m+;
         for(i=;i<=n;i++)scanf("%s",mp[i]+);
         for(i=;i<=n;i++)
             for(j=;j<=m;j++){
                 if(i== || i==n || j== || j==m){
                     if(mp[i][j]=='.')ans+=f;
                     add_edge(S,id[i][j],INF);
                     add_edge(id[i][j],S,);
                 }
                 else if(mp[i][j]=='#'){
                     add_edge(S,id[i][j],d);
                     add_edge(id[i][j],S,);
                 }
                 else{
                     add_edge(id[i][j],T,f);
                     add_edge(T,id[i][j],);
                 }
                 if(i>){
                     add_edge(id[i][j],id[i-][j],b);add_edge(id[i-][j],id[i][j],);
                     add_edge(id[i][j],id[i-][j],);add_edge(id[i-][j],id[i][j],b);
                 }
                 if(j>){
                     add_edge(id[i][j],id[i][j-],b);add_edge(id[i][j-],id[i][j],);
                     add_edge(id[i][j-],id[i][j],b);add_edge(id[i][j],id[i][j-],);
                 }
             }
         ans+=Dinic();
         printf("%d\n",ans);
     }
     return ;
 }