數據輸入格式：首先輸入頂點個數n和弧數m，然后輸入每條弧的數據。規定源點為頂點0，匯點為頂點n-1.每條弧的數據格式為：u，v，w，分別表示這條弧的起點，終點，容量。頂點序號從0開始。

代碼：

?

1 #include <iostream> 2 #include <cstdio> 3 #include <cstring> 4 #include <algorithm> 5 #include <cmath> 6 #include <string> 7 #include <map> 8 #include <stack> 9 #include <vector> 10 #include <set> 11 #include <queue> 12 #pragma comment (linker,"/STACK:102400000,102400000") 13 #define maxn 1005 14 #define MAXN 2005 15 #define mod 1000000009 16 #define INF 0x3f3f3f3f 17 #define pi acos(-1.0) 18 #define eps 1e-6 19 typedef long long ll; 20 using namespace std; 21 22 int d[maxn]; //標號 23 int mp[maxn][maxn]; //殘留網絡，初始為原圖 24 int num[maxn]; //num[i]表示標號為i的頂點數有多少 25 int pre[maxn]; //記錄前驅 26 int n,m,s,t; //n個頂點，m條邊，源點s，匯點t 27 28 void init() //bfs計算標號，匯點t的標號為零 29 { 30 int k; 31 queue<int>Q; 32 memset(d,INF,sizeof(d)); //初始化d數組為無窮大 33 memset(num,0,sizeof(num)); //初始化num數組num 34 Q.push(t); //匯點t入隊列 35 d[t]=0; //匯點標號為零 36 num[0]=1; //標號為0的頂點數為1 37 while (!Q.empty()) 38 { 39 k=Q.front(); 40 Q.pop(); 41 for (int i=0;i<n;i++) 42 { 43 if (d[i]>=n&&mp[i][k]>0) 44 { 45 d[i]=d[k]+1; 46 Q.push(i); 47 num[d[i]]++; 48 } 49 } 50 } 51 } 52 53 int findAlowArc(int i) //從i出發尋找允許弧 54 { 55 int j; 56 for (j=0;j<n;j++) 57 if (mp[i][j]>0&&d[i]==d[j]+1) 58 return j; 59 return -1; 60 } 61 62 int reLable(int i) //重新標號 63 { 64 int mm=INF; 65 for (int j=0;j<n;j++) 66 if (mp[i][j]>0) 67 mm=min(mm,d[j]+1); 68 return mm==INF?m:mm; 69 } 70 71 int maxFlow(int s,int t) //求從源點s出發的最大流 72 { 73 int flow=0,i=s,j; 74 int delta; //增量 75 memset(pre,-1,sizeof(pre)); 76 while (d[s]<n) 77 { 78 j=findAlowArc(i); 79 if (j>=0) 80 { 81 pre[j]=i; 82 i=j; 83 if (i==t) //更新殘留網絡 84 { 85 delta=INF; 86 for (i=t;i!=s;i=pre[i]) 87 delta=min(delta,mp[pre[i]][i]); 88 for (i=t;i!=s;i=pre[i]) 89 mp[pre[i]][i]-=delta,mp[i][pre[i]]+=delta; 90 flow+=delta; 91 } 92 } 93 else 94 { 95 int x=reLable(i); //重新標號 96 num[x]++; 97 num[d[i]]--; 98 if (num[d[i]]==0) //間隙優化 99 return flow; 100 d[i]=x; 101 if (i!=s) 102 i=pre[i]; 103 } 104 } 105 return flow; 106 } 107 108 int main() 109 { 110 int u,v,c; //弧的起點，終點，容量 111 while (scanf("%d%d",&n,&m)!=EOF) //讀入頂點個數n和弧數m 112 { 113 memset(mp,0,sizeof(mp)); 114 for (int i=0;i<m;i++) 115 { 116 scanf("%d%d%d",&u,&v,&c); 117 mp[u][v]=c; 118 } 119 s=0,t=n-1; 120 init(); 121 printf("%d\n",maxFlow(0,n-1)); 122 } 123 return 0; 124 } 125 /* 126 6 10 127 0 1 8 128 0 2 4 129 1 3 2 130 1 4 2 131 2 1 4 132 2 3 1 133 2 4 4 134 3 4 6 135 3 5 9 136 4 5 7 137 */ View Code

?

轉載于:https://www.cnblogs.com/i8888/p/4048928.html

總結

以上是生活随笔為你收集整理的网络流之 最短增广路算法模板（SAP）的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。