1. 程式人生 > >【模板】網路最大流 (Dinic)

【模板】網路最大流 (Dinic)

題目描述

給出一個網路圖,以及其源點和匯點,求出其網路最大流。

輸入輸出格式

輸入格式:

第一行包含四個正整數N、M、S、T,分別表示點的個數、有向邊的個數、源點序號、匯點序號。

接下來M行每行包含三個正整數ui、vi、wi,表示第i條有向邊從ui出發,到達vi,邊權為wi(即該邊最大流量為wi)

輸出格式:

一行,包含一個正整數,即為該網路的最大流。

輸入輸出樣例

輸入樣例#1:

4 5 4 3
4 2 30
4 3 20
2 3 20
2 1 30
1 3 40

輸出樣例#1:

50

說明

時空限制:1000ms,128M

資料規模:

對於30%的資料:N<=10,M<=25

對於70%的資料:N<=200,M<=1000

對於100%的資料:N<=10000,M<=100000

樣例說明:

這裡寫圖片描述

題目中存在3條路徑:

4–>2–>3,該路線可通過20的流量

4–>3,可通過20的流量

4–>2–>1–>3,可通過10的流量(邊4–>2之前已經耗費了20的流量)

故流量總計20+20+10=50。輸出50。

程式碼

裸的模板題就不多說了,EK也能過但是我要用Dinic

#include<cstdio>
#include<algorithm>
using namespace std;
struct node{
    int
u,v,w,next; }e[200005]; int n,m,s,t,cnt=1,ans; //cnt設為一是為了方便異或找邊 int h[10005],level[10005]; void add(int u,int v,int w) { e[++cnt].u=u; e[cnt].v=v; e[cnt].w=w; e[cnt].next=h[u]; h[u]=cnt; } bool bfs() { int head=0,tail=1,q[10005]={0}; bool vis[10005]={0}; q[1]=s,vis[s]=1; while(head<tail) { head++; for
(int i=h[q[head]];i;i=e[i].next) { int to=e[i].v; if(!vis[to]&&e[i].w>0) { vis[to]++; q[++tail]=to; level[to]=level[e[i].u]+1; } } } if(vis[t])return 1; return 0; } //BFS存層數 int dfs(int x,int minf) { if(x==0||minf==0) return 0; if(x==t) return minf; int zhi=minf; for(int i=h[x];i;i=e[i].next) { int v=e[i].v,u=e[i].u; if(level[v]==level[u]+1) { int f=dfs(v,min(minf,e[i].w)); e[i].w-=f; e[i^1].w+=f; minf-=f; } } return zhi-minf; } int main() { int x,y,z; scanf("%d%d%d%d",&n,&m,&s,&t); for(int i=1;i<=m;i++) { scanf("%d%d%d",&x,&y,&z); add(x,y,z),add(y,x,0);//記得建立反向弧 } while(bfs())//只要找的到就一直找下去 { for(int i=h[s];i;i=e[i].next) { int f=dfs(e[i].v,e[i].w); e[i].w-=f; e[i^1].w+=f; ans+=f; } } printf("%d",ans); return 0; }