題目：Flow Problem

#include <iostream> #include <stdio.h> #include <string.h>using namespace std;const int N=2050; //N表示點的數量 const int M=2500000; //M表示邊的數量，注意應大于(N-1)*N，因為網絡是雙向的 const int INF=0x7fffffff;int e; int ver[M],next[M],flow[M]; //ver為邊指向的節點，flow為邊的容量，next為鏈表的下一條邊 int head[N],work[N],dis[N],q[N]; //head為節點的鏈表頭，work 用于算法中的臨時鏈表頭，dis 計算距離void init() {e=0;memset(head,-1,sizeof(head)); }void add(int u,int v,int c) //增加一條從u到v容量為c的邊 {ver[e]=v;flow[e]=c;next[e]=head[u];head[u]=e++;ver[e]=u;flow[e]=0;next[e]=head[v];head[v]=e++; }//廣搜計算出每個點與源點的最短距離，如果不能到達匯點說明算法結束 int bfs(int S,int T) {int rear=0;memset(dis,-1,sizeof(dis));dis[S]=0;q[rear++]=S;for(int i=0;i<rear;i++){for(int j=head[q[i]];j!=-1;j=next[j]){if(flow[j]&&dis[ver[j]]==-1){dis[ver[j]]=dis[q[i]]+1;q[rear++]=ver[j];if(ver[j]==T) return 1;}}}return 0; }//尋找可行流的增廣路算法，按節點的距離來找，加快速度 int dfs(int cur,int a,int T) {if(cur==T) return a;for(int &i=work[cur];i!=-1;i=next[i]) //work是臨時鏈表頭，這里用i引用它，這樣尋找過的邊不再尋找{if(flow[i]&&dis[ver[i]]==dis[cur]+1){int t=dfs(ver[i],min(a,flow[i]),T);if(t){flow[i]-=t;flow[i^1]+=t; //正反向邊容量改變return t;}}}return 0; }//求最大流，直到沒有可行流 long long Dinic(int S,int T) {long long ans=0;while(bfs(S,T)){memcpy(work,head,sizeof(head));while(int t=dfs(S,INF,T)) ans+=t;}return ans; }int main() {long long sum;int t,n,m,S,T,x,y,c,k=1;scanf("%d",&t);while(t--){scanf("%d%d",&n,&m);init();S=1;T=n;while(m--){scanf("%d%d%d",&x,&y,&c);add(x,y,c);}printf("Case %d: %I64d\n",k++,Dinic(S,T));}return 0; }

總結

以上是生活随笔為你收集整理的HDU3549(最大流算法的Dinic算法)的全部內容，希望文章能夠幫你解決所遇到的問題。

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