ISAP

[置顶]【整合】网络流ISAP算法模板

网络流 图论

[cpp] view plain copy print ?

1. 注释应该比较清楚了吧= = 

注释应该比较清楚了吧= =

[cpp] view plain copy print ?

1. #include<iostream> 
2. #include<cstdio> 
3. #include<cstring> 
4. #include<cmath> 
5. #include<algorithm> 
6. #define MAXN 2048 
7. #define MAXINT 0x7fffffff 
8. using namespace std; 
9. int n,m;//点的数目和边的数目 
10. int dis[MAXN+1];//距离标号  
11. int num[MAXN+1]; 
12. int src,des;//src源点，des汇点  
13. int u,v,w; 
14. struct edge 
15. { 
16.     int f;//某条弧当前流量  
17.     int c;//容量 
18.     int ver;//弧的起点 
19.     edge *rev;//反弧  
20.     edge *next;//弧的终点  
21.     edge(){}; 
22.     edge(int v,int cap,edge *nex):ver(v),c(cap),next(nex),rev(NULL),f(0){}; 
23.     void* operator new(size_t, void *p) 
24.     { 
25.         return p; 
26.     } 
27. }*s[MAXN+1]; 
28. void init() 
29. { 
30.     int queue[MAXN+1],head=0,tail=0; 
31.     for (int i=1;i<=n;i++) 
32.         dis[i]=MAXN,num[i]=0; 
33.     queue[tail++]=des; 
34.     dis[des]=0; 
35.     num[0]=1; 
36.     while (head!=tail) 
37.     { 
38.         int v=queue[head++]; 
39.         edge *e=s[v]; 
40.         while (e) 
41.         { 
42.             if ((e->rev)&&(e->rev->c==0)||dis[e->ver]<MAXN) {}//若改点后的边不可通过 
43.             //流或者已经更新就无需更新 
44.             else  
45.             { 
46.                 dis[e->ver]=dis[v]+1; 
47.                 ++num[dis[e->ver]]; 
48.                 queue[tail++]=e->ver; 
49.             } 
50.             e=e->next; 
51.         } 
52.     } 
53. } 
54. int maxflow()//ISAP算法实现过程  
55. { 
56.     int st=src,ret=0; 
57.     edge *E[MAXN],*rep[MAXN]; 
58.     for (int i=1;i<=n;i++) E[i]=s[i]; 
59.     while (dis[src]<n)//不存在增广路径就停止  
60.     { 
61.         if (st==des)//找到了增广路径  
62.         { 
63.             int delta=MAXINT; 
64.             for (int i=src;i!=des;i=E[i]->ver) 
65.                 delta=min(delta,E[i]->c); 
66.             for (int i=src;i!=des;i=E[i]->ver) 
67.             { 
68.                 E[i]->c-=delta; 
69.                 E[i]->f+=delta; 
70.                 E[i]->rev->c+=delta; 
71.                 E[i]->rev->f-=delta; 
72.             } 
73.             ret+=delta; 
74.             st=src; 
75.         } 
76.         edge *e; 
77.         for (e=E[st];e;e=e->next) 
78.             if (e->c>0&&dis[st]==dis[e->ver]+1) break; 
79.         if (e)//存在允许弧  
80.         { 
81.             E[st]=e; 
82.             rep[e->ver]=e->rev; 
83.             st=e->ver; 
84.         } 
85.         else 
86.         { 
87.             if ((--num[dis[st]])==0) break;//GAP优化  
88.             E[st]=s[st]; 
89.             int mind=n; 
90.             for (edge *t=s[st];t;t=t->next) 
91.                 if (t->c>0) 
92.                     mind=min(mind,dis[t->ver]); 
93.             dis[st]=mind+1; 
94.             ++num[dis[st]]; 
95.             if (st!=src) st=rep[st]->ver; 
96.         } 
97.     } 
98.     return ret; 
99. } 
100. int main() 
101. { 
102.     freopen("ditch.in","r",stdin); 
103.     freopen("ditch.out","w",stdout); 
104.     scanf("%d%d",&m,&n); 
105.     src=1,des=n; 
106.     edge *buffer=new edge[2*m]; 
107.     edge *data=buffer; 
108.     while (m--) 
109.     { 
110.         scanf("%d%d%d",&u,&v,&w); 
111.         s[u]=new((void*) data++) edge(v,w,s[u]); 
112.         s[v]=new((void*) data++) edge(u,0,s[v]); 
113.         s[u]->rev=s[v]; 
114.         s[v]->rev=s[u]; 
115.     } 
116.     init(); 
117.     printf("%d",maxflow()); 
118. }