題目描述
如題，給出一個網路圖，以及其源點和匯點，求出其網路最大流。

輸入輸出格式
輸入格式：
第一行包含四個正整數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

一堆註釋掉的中間輸出見證了腦殘錯誤的代價QAQ。。。
code：

//By Menteur_Hxy
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;

int rd() {
    int x=0,fla=1;
    char c=' ';
    while(c>'9'
 
|| c<'0') {if(c=='-') fla=-fla; c=getchar();}
    while(c<='9'&&c>='0') x=x*10+c-'0',c=getchar();
//  scanf("%d",&x);
    return x*fla;
}

const int MAX=10010;
const int INF=1e9;
int n,m,s,t,cnt;
int maxflow;
int head[MAX],deep[MAX],cur[MAX];

struct edges{
    int to,next,w;
};
edges edge[200010
 
];

void add(int x,int y,int z) {
     edge[++cnt].next=head[x];
     edge[cnt].to=y;
     edge[cnt].w=z;
     head[x]=cnt;
//     cout<<cnt<<" "<<edge[cnt].to<<" "<<edge[cnt].next<<" "<<edge[cnt].w<<endl;
}

queue <int> q;

bool bfs() {
    memset(deep,0,sizeof deep);
    deep[s]=1;
    for(int i=1;i<=n;i++) cur[i]=head[i];//!!!
    while(!q.empty()) q.pop();
    q.push(s);

    while(!q.empty()) {
        int u=q.front(),v; q.pop();
        for(int i=head[u];i!=-1;i=edge[i].next) 
            if(!deep[v=edge[i].to] && edge[i].w) {
                deep[v]=deep[u]+1;
                q.push(v);
            }
    }
//  cout<<deep[t]<<endl;
//  for(int i=1;i<=n;i++) cout<<deep[i]<<" ";
//  cout<<endl;
    return deep[t];
} 

int dfs(int now,int flow) {
    if(now==t) return flow;

    for(int& i=cur[now];i!=-1;i=edge[i].next) {
        int v=edge[i].to,di;
//      if(v==t) cout<<edge[i].w<<endl;
        if(deep[now]+1==deep[v] && edge[i].w 
            && (di=dfs(v,min(flow,edge[i].w))) ) {
//          if(v==t) cout<<edge[i].w<<endl;
            edge[i].w-=di;
            edge[i^1].w+=di;
            return di;
        }
    }
    return 0;
}

void dinic() {
    while(bfs())
        while(int di=dfs(s,INF))
            maxflow+=di;
}

int main() {
    n=rd(),m=rd(),s=rd(),t=rd();
    cnt=-1;
    memset(head,-1,sizeof head);
    for(int i=1;i<=m;i++) {
        int a=rd(),b=rd(),c=rd();
        add(a,b,c);
        add(b,a,0);
    }
    dinic();
    printf("%lld",maxflow);
    return 0;
}