网络流——洛谷P1345 [USACO5.4]奶牛的电信Telecowmunication

https://www.luogu.org/problem/show?pid=1345
就是建一个双向边的网络流，然后要用拆点把边权变成点权；
但是双向边的话就容易出错了；

 scanf("%d%d",&x,&y);
 init(x+n,y,1e9);init(y,x+n,0);
 init(y+n,x,1e9);init(x,y+n,0);

for(int i=1;i<=n;i++)init(i,i+n,1),init(i+n,i,0);

就是这两个啦；

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
using namespace std;
const int N=1e4+5;
struct cs{int to,nxt,v;}a[N*20];
int head[N],ll=1,cur[N];
int q[N],l,r;
int deep[N];
int n,m,S,E,x,y,z,ans;
void init(int x,int y,int z){a[++ll].to=y;a[ll].v=z;a[ll].nxt=head[x];head[x]=ll;}
bool bfs(){
    memset(deep,0,sizeof deep);//这个其实可以优化的，我不会 
    q[1]=S;l=r=1;deep[S]=1;
    for(;r>=l;l++){//当前弧优化赋值↓ 
        int x=q[l];cur[x]=head[x];if(x==E)return 1;//发现跟就可以退出了 
        for(int k=head[x];k;k=a[k].nxt)
            if(!deep[a[k].to]&&a[k].v!=0)
                deep[a[k].to]=deep[x]+1,q[++r]=a[k].to;
    }return 0;
}
int dfs(int x,int now){//now代表当前x的可用流 
    if(x==E)return now; 
    int ans=0; 
    for(int &k=cur[x];k;k=a[k].nxt)//只访问一次 
        if(deep[a[k].to]-1==deep[x]&&a[k].v){//这个优化a[k].v会快的 
            int temp=dfs(a[k].to,min(a[k].v,now));
            now-=temp; ans+=temp;
            a[k].v-=temp; a[k^1].v+=temp;//反向边 
            if(!now)break;//优化 
        }return ans;
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&S,&E);S+=n;
    for(int i=1;i<=m;i++){
        scanf("%d%d",&x,&y);
        init(x+n,y,1e9);init(y,x+n,0);
        init(y+n,x,1e9);init(x,y+n,0);
    }
    for(int i=1;i<=n;i++)init(i,i+n,1),init(i+n,i,0);//,init(i,i+n,0),init(i+n,i,1);
    while(bfs())ans+=dfs(S,1e9);
    printf("%d",ans);
}