HDU 3549 Flow Problem - 更新dinic模板..整成结构体...

我发现我以前写的dinic并没有完全做对...效率不高...这次更新了dinic的模板..整成一个结构体..用起来方便..效率不错....不过效率没有sap高..按道理来说dinic和sap的效率没有差别的...实在不想学sap了...先这样吧....

Program:

#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<queue>
#define MAXN 5005 
#define MAXM 50005
#define oo 2000000000
#define ll long long
using namespace std;
struct Dinic
{
      struct node 
      {
            int x,y,c,next;
      }line[MAXM];   
      int Lnum,_next[MAXN],dis[MAXN];
      void initial() { memset(_next,-1,sizeof(_next)),Lnum=-1;} 
      void addline(int x,int y,int c)
      {
             line[++Lnum].next=_next[x],_next[x]=Lnum;
             line[Lnum].x=x,line[Lnum].y=y,line[Lnum].c=c;
             line[++Lnum].next=_next[y],_next[y]=Lnum;
             line[Lnum].x=y,line[Lnum].y=x,line[Lnum].c=0;             
      }
      bool BFS(int s,int e)
      { 
             queue<int> Q;
             while (!Q.empty()) Q.pop();
             memset(dis,0,sizeof(dis));
             dis[s]=1;
             Q.push(s);
             while (!Q.empty())
             {
                   int h,k;
                   h=Q.front(),Q.pop();
                   if (h==e) return dis[e];
                   for (k=_next[h];k!=-1;k=line[k].next)
                      if (line[k].c && !dis[line[k].y])
                         dis[line[k].y]=dis[h]+1,Q.push(line[k].y);                 
             } 
             return false;
      }
      int dfs(int x,int flow,int e)
      {   
             if (x==e) return flow; 
             int temp,cost=0;
             for (int k=_next[x];k!=-1;k=line[k].next)
             if (line[k].c && dis[line[k].y]==dis[x]+1)
             {
                    temp=dfs(line[k].y,min(flow-cost,line[k].c),e); 
                    if (temp)
                    {
                           line[k].c-=temp,line[k^1].c+=temp;
                           cost+=temp;
                           if (flow==cost) return cost;
                    }else dis[line[k].y]=-1;
             }
             return cost;
      }
      int MaxFlow(int s,int e)
      {
             int MaxFlow=0;
             while (BFS(s,e)) 
                MaxFlow+=dfs(s,oo,e); 
             return MaxFlow;
      }
}T;
int main()
{      
       int i,x,y,c,n,m,cases,C; 
       scanf("%d",&C);
       for (cases=1;cases<=C;cases++)
       {
              scanf("%d%d",&n,&m);
              T.initial(); 
              while (m--)
              {
                     scanf("%d%d%d",&x,&y,&c);
                     T.addline(x,y,c);
              } 
              printf("Case %d: %d\n",cases,T.MaxFlow(1,n));
       }
       return 0;
}