洛谷-3376 【模板】网络最大流

题目描述
如题，给出一个网络图，以及其源点和汇点，求出其网络最大流。
输入输出格式
输入格式：
第一行包含四个正整数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

#include<bits/stdc++.h>
using namespace std;
const int M=1e6+6,inf=0x3f3f3f3f;
inline int read(){
    int s=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) s=s*10+ch-'0',ch=getchar();
    return s;
}
int n,m;
int b[505][505],c[505],id[505][505];
struct edge{
    int to,w,rev;
};
int d[M],cur[M];vector<edge> G[M];
int s,t,tot=0;
queue<int> q;
inline int bfs(){
    memset(d,-1,sizeof(d));
    d[s]=0,q.push(s);
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=0;i<G[u].size();i++){
        	edge e=G[u][i];
            if(e.w>0&&d[e.to]<0){
                d[e.to]=d[u]+1;
                q.push(e.to);
            }
        }
    }
    return d[t]!=-1;
}
int dfs(int u,int f){
    if(u==t) return f;
    int w,tot=0;
    for(int &i=cur[u];i<G[u].size();i++){
        edge &e=G[u][i];
        if(e.w>0&&d[e.to]==d[u]+1&&(w=dfs(e.to,min(e.w,f)))){
            e.w-=w;G[e.to][e.rev].w+=w;
            f-=w,tot+=w;
            if(f<=0) break;
        } 
    }
    return tot;
}
inline void add_edge(int u,int v,int c){
    G[u].push_back({v,c,G[v].size()});
    G[v].push_back({u,0,G[u].size()-1});
}
inline int dinic(){
    int flow=0;
    while(bfs()){
        memset(cur,0,sizeof(cur));
        flow+=dfs(s,inf);
    }
    return flow;
}
int main(){
    n=read(),m=read();s=read();t=read();
    for(int i=1,x,y,z;i<=m;i++){
    	x=read();y=read();z=read();
    	add_edge(x,y,z);
	}
    printf("%d",dinic());
    return 0;
}