「LuoguP3376」 【模板】网络最大流

  1 /*
  2     qwerta
  3     P3376 【模板】网络最大流
  4     Accepted
  5     100
  6     代码 C++，1.68KB
  7     提交时间 2018-07-12 17:31:29
  8     耗时/内存
  9     140ms, 4847KB
 10 */
 11 #include<cmath>
 12 #include<queue>
 13 #include<cstdio>
 14 #include<cstring>
 15 #include<iostream>
 16 using namespace std;
 17 struct emm{
 18     int e,f,v;
 19 }a[200007];//邻接链表存边
 20 int h[10007],cur[10007];//cur:当前弧优化（没什么用
 21 int n,m,s,t;
 22 int tot=1;
 23 inline int read()//快读
 24 {
 25     char ch=getchar();
 26     int s=1,x=0;
 27     while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)s=-1;ch=getchar();}
 28     while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-48;ch=getchar();}
 29     return s*x;
 30 }
 31 inline void con(int l,int r,int w)//加边
 32 {
 33     a[++tot].f=h[l];
 34     h[l]=tot;
 35     //cur[l]=h[l];
 36     a[tot].e=r;
 37     a[tot].v=w;
 38     return;
 39 }
 40 inline void scan()//读入，建图
 41 {
 42     n=read(),m=read(),s=read(),t=read();
 43     for(int i=1;i<=m;++i)
 44     {
 45         int z=read(),y=read(),l=read();
 46         con(z,y,l);
 47         con(y,z,0);
 48     }
 49     return;
 50 }
 51 //Dinic
 52 queue<int>q;//bfs的queue
 53 int  d[100007];//分层图的深度
 54 inline bool bfs()
 55 {
 56     memset(d,0,sizeof(d));//初始化
 57     d[s]=1;//标记源点深度
 58     q.push(s);
 59     for(int i=1;i<=n;++i)cur[i]=h[i];//恢复cur数组
 60     while(!q.empty())
 61     {
 62         int now=q.front();
 63         q.pop();
 64         //扩张
 65         for(int i=h[now];i;i=a[i].f)
 66         if(!d[a[i].e]&&a[i].v)//若未被标记过并且该边在残量网络中
 67         {
 68             d[a[i].e]=d[now]+1;
 69             q.push(a[i].e);
 70         }
 71     }
 72     return d[t];//返回s,t是否联通
 73 }
 74 int dfs(int x,int al)
 75 {
 76     if(x==t||!al)return al;
 77     int tot=0;
 78     for(int i=cur[x];i;i=a[i].f)
 79     {
 80         cur[x]=i;//当前弧优化
 81         if(d[a[i].e]==d[x]+1&&a[i].v)
 82         {
 83             int f=dfs(a[i].e,min(al,a[i].v));//往下找
 84             if(f)//若非0
 85             {
 86                 a[i].v-=f;
 87                 a[i^1].v+=f;
 88                 tot+=f;
 89                 al-=f;
 90                 if(!al)break;
 91             }
 92         }
 93     }
 94     if(!tot)d[x]=-1;//最有用的优化！（敲黑板
 95     return tot;//返回流量值
 96 }
 97 inline void run()//运行
 98 {
 99     long long ans=0;
100     while(bfs())ans+=dfs(s,2147483647);
101     cout<<ans;
102     return;
103 }
104 int main()//超短主函数（当年码风真奇怪orz
105 {
106     scan();
107     run();
108     return 0;
109 }