[loj#101] 最大流 网络流模板

#101. 最大流

内存限制：512 MiB时间限制：5000 ms标准输入输出

题目类型：传统评测方式：文本比较

上传者： 匿名

提交提交记录统计讨论测试数据

 

题目描述

这是一道模板题。

给定 n nn 个点，m mm 条边，给定每条边的容量，求从点 s ss 到点 t tt 的最大流。

输入格式

第一行四个整数 n nn、m mm、s ss、t tt。
接下来的 m mm 行，每行三个整数 u uu、v vv、c cc，表示 u uu 到 v vv，流量为 c cc 的一条边。

输出格式

输出点 s ss 到点 t tt 的最大流。

样例

样例输入

7 14 1 7
1 2 5
1 3 6
1 4 5
2 3 2
2 5 3
3 2 2
3 4 3
3 5 3
3 6 7
4 6 5
5 6 1
6 5 1
5 7 8
6 7 7

样例输出

14

数据范围与提示

1≤n≤106,1≤m≤4×106,0≤c≤231−1 1 \leq n \leq 10 ^ 6, 1 \leq m \leq 4 \times 10 ^ 6, 0 \leq c \leq 2 ^ {31} - 11≤n≤10​6​​,1≤m≤4×10​6​​,0≤c≤2​31​​−1

 #include<iostream>
 #include<cstring>
 #include<cstdio>
 #include<cstdlib>
 #include<cmath>
 #include<algorithm>
 #define maxm 4000006
 #define maxn 1000006
 using namespace std;
 int read() {
     int x=,f=;char ch=getchar();
     while(!isdigit(ch)){ch=getchar();}
     while(isdigit(ch)){x=x*+ch-'';ch=getchar();}
     return x;
 }
 struct data {
     int from,to,next,w;
 }e[maxm*];
 int head[maxn],cnt;
 int cur[maxn];
 void add(int u,int v,int w){e[cnt].from=u;e[cnt].next=head[u];e[cnt].to=v;e[cnt].w=w;head[u]=cnt++;}
 int n,m,s,t;
 int q[maxn];
 bool vis[maxn];
 int dis[maxn];
 bool bfs() {
     memset(dis,-,sizeof(dis));
     int h=,tt=;
     q[h]=t;
     vis[t]=;
     dis[t]=;
     while(h!=tt) {
         int now=q[h];h++;vis[now]=;if(h==maxn) h=;
         for(int i=head[now];i>=;i=e[i].next) {
             int to=e[i].to;
             if(e[i^].w&&dis[to]<-) {
                 dis[to]=dis[now]-;
                 if(!vis[to]){
                     vis[to]=;
                     q[tt++]=to;if(tt==maxn) tt=;
                 }
             }
         }
     }
     return dis[s]>=-;
 }
 int dfs(int now,int a) {
     if(now==t||a==) return a;
     int flow=,f;
     for(int i=cur[now];i>=;i=e[i].next) {
         int to=e[i].to;
         if(dis[to]==dis[now]+&&e[i].w>&&(f=dfs(to,min(a,e[i].w)))) {
             e[i].w-=f;
             e[i^].w+=f;
             flow+=f;
             a-=f;
             if(a==) return flow;
         }
         cur[now]=i;
     }
     if(!flow) dis[now]=-;
     return flow;
 }
 int main() {
     memset(head,-,sizeof(head));
     n=read(),m=read(),s=read(),t=read();
     for(int i=;i<=m;i++) {
         int u=read(),v=read(),w=read();
         add(u,v,w);add(v,u,);
     }
     int ans=;
     while(bfs()){
         for(int i=;i<=n;i++) cur[i]=head[i];
         ans+=dfs(s,);
     }
     printf("%d",ans);
 }