PKU1459 Power Network - 最大流

模板题，不说了，整理一下模板：

 

1. /*
2. PKU1459 Power Network
3. */
5. #include <stdio.h>
6. #include <stdlib.h>
7. #include <memory.h>
9. #define clr(a)      memset(a,0,sizeof(a))
10. #define N 200
11. #define INF (1<<25)
13. int EdmondsKarp(int g[][N],int n,int s,int t,int f[][N])
14. {
15.     int i,j,k,c,head,tail,flow=0 ;
16.     int r[N][N];
17.     int prev[N],visit[N],q[N];
18.     for(i=0;i<n;i++)for(j=0;j<n;j++)
19.     {
20.         f[i][j]=0 ;
21.         r[i][j]=g[i][j];
22.     }
23.     //初始化流量网络和残留网络
24.     while(1)
25.     {
26.         //在残留网络中找到一条s到t的最短路径
27.         head=tail=0 ;
28.         memset(visit,0,sizeof(visit));
29.         q[tail++]=s ;
30.         prev[s]=-1 ;
31.         visit[s]=1 ;
32.         while(head<tail)
33.         {
34.             //宽度优先搜索从s到t的最短路径
35.             k=q[head++];
36.             for(i=0;i<n;i++)
37.             if(!visit[i]&&r[k][i]>0)
38.             {
39.                 visit[i]=1 ;
40.                 prev[i]=k ;
41.                 if(i==t)goto next ;
42.                 q[tail++]=i ;
43.             }
44.         }
45.         next :
46.         if(!visit[t])break ;
47.         //流量已达到最大
48.         for(c=INF,j=t;j!=s;j=i)
49.         {
50.             i=prev[j];
51.             if(c>r[i][j])c=r[i][j];
52.         }
53.         //下面改进流量
54.         for(j=t;j!=s;j=i)
55.         {
56.             i=prev[j];
57.             f[i][j]+=c ;
58.             f[j][i]=-f[i][j];
59.             r[i][j]=g[i][j]-f[i][j];
60.             r[j][i]=g[j][i]-f[j][i];
61.         }
62.         flow+=c ;
63.     }
64.     return flow ;
65. }
67. int g[N][N];
68. int f[N][N];
70. int main()
71. {
72.     int i,j,k,x,y;
73.     int n,np,nc,m;
74.     int s,t;
75.     
76.     while(scanf("%d%d%d%d",&n,&np,&nc,&m)!=EOF){
77.         //init
78.         clr(g); clr(f);
79.         //input
80.         for(i=0;i<m;i++){
81.             scanf(" (%d,%d)%d",&x,&y,&k); g[x][y]=k;
82.         }
83.         s=n++; t=n++;
84.         for(i=0;i<np;i++){
85.             scanf(" (%d)%d",&x,&k); g[s][x]=k;
86.         }
87.         for(i=0;i<nc;i++){
88.             scanf(" (%d)%d",&x,&k); g[x][t]=k;
89.         }
90.         //output
91.         printf("%d/n",EdmondsKarp(g,n,s,t,f));
92.     }
93.     
94.     return 0;
95. }