图论网络流 Dinic算法和模板

zz http://blog.sina.com.cn/s/blog_61533c9b0100iawn.html

Dinic算法是基于分层思想的网络流算法，它的的效率一般比E-K算法高
算法的思想:Dinic是在分阶段地在层次图中增广。它与最短路径增值算法不同之处是：在Dinic算法中，我们用一个dfs过程代替多次bfs来寻找阻塞流。下面给出其算法步骤：
1、初始化流量，计算出剩余图
    2、根据剩余图计算层次图。若汇点不在层次图内，则算法结束
    3、在层次图内用一次dfs过程增广
    4、转步骤2
 
算法模板：
//时间复杂度O（V^2E)
    #include<iostream>
#define Max 210
    int flow[Max][Max],d[Max]; //flow is the network
    int sta,end,N; //sta is the sourse ,end is the,N is the number of vector
    bool bfs(int s)
    {
                   int front=0,rear=0; int q[Max];
                   memset(d,-1,sizeof(d));  //d[] is the deep
                   q[rear++]=s;  d[s]=0;
                   while(front<rear)
                   {
                           int k=q[front++];
                           for(int i=0;i<=N;i++)
                              if(flow[k][i]>0&&d[i]==-1){
                                 d[i]=d[k]+1;
                                 q[rear++]=i;
                              }
                   }
             if(d[end]>=0)   return true;
                         return false;
            }
int dinic(int k,int sum)  //k is the sourse
   {
                  if(k==end)    return sum;
                  int os=sum;
                  for(int i=0;i<=N&&sum;i++)
                  if(d[i]==d[k]+1&&flow[k][i]>0)
                  {
                            int a=dinic(i,min(sum,flow[k][i])); //Deep to the end.
                            flow[k][i]-=a;
                            flow[i][k]+=a;
                            sum-=a;
                  }
                  return os-sum;
    }
int main()
   {
            int ret=0;
                while(bfs(sta))
                   ret+=dinic(sta,INT_MAX);
                printf("%d\n",ret);
                return 0;
    }