网络流--Dinic(自用，勿看)

注意：这是一篇个人学习笔记，如果有人因为某些原因点了进来并且要看一下，请一定谨慎地阅读，因为可能存在各种奇怪的错误，如果有人发现错误请指出谢谢！

---

https://www.luogu.org/problemnew/show/P3376

来自蓝书：

时间复杂度O(n^2*m)

所有容量均为1，可以证明时间复杂度O(min(n^(2/3),m^(1/2))*m)

除了源点和汇点之外，每个点要么只有一条入弧，且容量为1，要么只有一条出弧，且容量为1，其他弧的容量为任意整数，可以证明时间复杂度O(n^(1/2)*m)

 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<vector>
 #include<queue>
 using namespace std;
 #define fi first
 #define se second
 #define mp make_pair
 #define pb push_back
 typedef long long ll;
 typedef unsigned long long ull;
 typedef pair<int,int> pii;
 namespace F
 {

 struct E
 {
     int to,nxt,from,cap,flow;
 }e[];
 int f1[],ne=;
 int S,T,n;
 int d[];
 bool bfs()
 {
     int k,u;
     memset(d,,sizeof(int)*(n+));
     queue<int> q;
     q.push(S);d[S]=;
     while(!q.empty())
     {
         u=q.front();q.pop();
         for(k=f1[u];k;k=e[k].nxt)
             if(!d[e[k].to]&&e[k].cap>e[k].flow)
             {
                 d[e[k].to]=d[u]+;
                 //if(e[k].to==T)    return 1;
                 q.push(e[k].to);
             }
     }
     //return 0;
     return d[T];
 }
 int cur[];
 int dfs(int u,int x)
 {
     if(u==T||x==)    return x;
     int flow=,f;
     for(int &k=cur[u];k;k=e[k].nxt)
         if(e[k].cap>e[k].flow&&d[e[k].to]==d[u]+)
         {
             f=dfs(e[k].to,min(x-flow,e[k].cap-e[k].flow));
             e[k].flow+=f;e[k^].flow-=f;flow+=f;
             if(flow==x)    return flow;
         }
     return flow;
 }
 int solve()
 {
     int flow=;
     while(bfs())
     {
         memcpy(cur,f1,sizeof(int)*(n+));
         flow+=dfs(S,0x3f3f3f3f);
     }
     return flow;
 }
 void me(int a,int b,int c)
 {
     e[++ne].to=b;e[ne].nxt=f1[a];f1[a]=ne;
     e[ne].from=a;e[ne].cap=c;
     e[++ne].to=a;e[ne].nxt=f1[b];f1[b]=ne;
     e[ne].from=b;e[ne].cap=;
 }

 }
 int n,m;
 int main()
 {
     int i,a,b,c;
     scanf("%d%d%d%d",&n,&m,&F::S,&F::T);F::n=n;
     for(i=;i<=m;i++)
     {
         scanf("%d%d%d",&a,&b,&c);
         F::me(a,b,c);
     }
     printf("%d",F::solve());
     return ;
 }

（补上37行，42行改成41行，跑的要更快，但是感觉不太对所以没有写）