首页 > 代码库 > 【CODEVS】1993 草地排水

【CODEVS】1993 草地排水

【算法】网络流-最大流(dinic)

【题解】网络流:http://m.blog.csdn.net/article/details?id=9401909

当前弧优化是因为DFS过程中访问x点时一旦流入量=流出量就退出,所以可以记录下此时正在考虑的弧,下次从此处继续考虑即可。

当前弧之前的弧,不能使流入量-流出量=0,那么一定该弧以及该弧之后的弧中有断裂,那么下次再考虑就没有意义了。

当前弧本身,使流入量-流出量=0,也就是使该点前面的弧中最小的一条断裂了,当前弧以及当前弧连出去之后的弧只是有可能断裂,那么下次就应该从这条当前弧开始考虑。

我的写法中直接用flow表示剩余流量,思路和程序结构参考紫书。

技术分享
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
const int maxn=210,inf=0x3f3f3f3f;
struct edge{int from,v,flow;}e[maxn*3];
int n,m,tot=2,S,T,d[maxn],q[510],first[maxn],cur[maxn];
void insert(int u,int v,int flow)
{tot++;e[tot].v=v;e[tot].flow=flow;e[tot].from=first[u];first[u]=tot;}
bool bfs()
{
    memset(d,-1,sizeof(d));
    int head=0,tail=1;q[head]=S;
    d[S]=0;
    while(head!=tail)
     {
         int x=q[head++];if(head>=501)head=0;
         for(int i=first[x];i;i=e[i].from)
         if(d[e[i].v]==-1&&e[i].flow>0)
          {
              q[tail++]=e[i].v;d[e[i].v]=d[x]+1;
              if(tail>=501)tail=0;
          }
     }
    if(d[T]==-1)return 0;else return 1;
}
int dfs(int x,int a)
{
    if(x==T||a==0)return a;
    int flow=0,f;
    for(int& i=cur[x];i;i=e[i].from)//当前弧优化 
     if(d[e[i].v]==d[x]+1&&(f=dfs(e[i].v,min(a,e[i].flow)))>0)
      {
          e[i].flow-=f;
          e[i^1].flow+=f;
          flow+=f;
          a-=f;
          if(a==0)break;
      }
    return flow;
}
int main()
{
    scanf("%d%d",&m,&n);
    S=1;T=n;
    for(int i=1;i<=m;i++)
     {
         int u,v,w;
         scanf("%d%d%d",&u,&v,&w);
         insert(u,v,w);
         insert(v,u,0);
     }
    int ans=0;
    while(bfs())
     {
         for(int i=1;i<=n;i++)cur[i]=first[i];
         ans+=dfs(S,inf);
     }
    printf("%d",ans);
    return 0;
}
View Code

 

【CODEVS】1993 草地排水