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

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
,ll=1,cur
;
int q
,l,r;
int deep
;
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);
}