有源汇有上下界最小流

link

洛谷上的板子被合并了。不是很妙。

直接说整体流程吧。第一步统计每个点出流量和入流量的差值，如果出的多就从源点向出点连边，否则就向汇点连边。然后原图上每条边都重新赋值为上下流量之差，跑最大流，如果源汇点都能流满说明存在可行流，需要注意的是要额外连一条从原汇点到原源点的无限边以维护流量的平衡。最后删除无限边，跑一条从原汇点到原源点的最大流，用无限边上的流量减去这个最大流就可以得到有源汇上下界最小流了。如果要跑最大流那么从原源点向原汇点跑个最大流然后加起来即可。当时调了很久甚至锤了门一下，就很悲惨。

代码比较好写。

#include<bits/stdc++.h>
//#define feyn
#define int long long
const int N=50010;
const int M=5e6;
const int maxn=1e15;
using namespace std;
inline void read(int &wh){
	wh=0;int f=1;char w=getchar();
	while(w<'0'||w>'9'){if(w=='-')f=-1;w=getchar();}
	while(w>='0'&&w<='9'){wh=wh*10+w-'0';w=getchar();}
	wh*=f;return;
}
inline int min(int s1,int s2){
	return s1<s2?s1:s2;
}

int ss1,tt1,sum=0;

struct edge{
	int t,v,next;
}e[M];
int head[N],esum=1;
inline void adde(int fr,int to,int val){
	e[++esum]=(edge){to,val,head[fr]};head[fr]=esum;
}
inline void add(int fr,int to,int val){
	adde(fr,to,val);adde(to,fr,0);
}

int m,ss,tt,s0,t0,a[N],xx;

int q[N],ll,rr,d[N],t[N],nt;
inline bool check(){
	q[ll=rr=1]=ss;t[ss]=++nt;d[ss]=1;
	while(ll<=rr){
		int wh=q[ll++];//if(tt==ss1)printf("%lld\n",wh);
		for(int i=head[wh],th;i;i=e[i].next){
			if(e[i].v==0||t[th=e[i].t]==nt)continue;
			t[th]=nt,d[th]=d[wh]+1,q[++rr]=th;
		}
	}
	//printf("%lld %lld\n",ss,tt);
	return t[tt]==nt;
}
inline int dinic(int wh,int val){
	if(wh==tt)return val;
	int used=0;
	for(int i=head[wh],th;i;i=e[i].next){
		if(e[i].v==0||d[th=e[i].t]!=d[wh]+1)continue;
		int now=dinic(th,min(val,e[i].v));
		if(now){
			//if(wh<=m&&th<=m)printf("%lld %lld\n",th,wh);
			used+=now,val-=now,e[i].v-=now,e[i^1].v+=now;
		}
		if(val==0)break;
	}
	if(val)d[wh]=-1;return used;
}


signed main(){
	
	#ifdef feyn
	freopen("in.txt","r",stdin);
	//freopen("out.txt","w",stdout);
	#endif
	
	read(m);read(ss1);read(ss1);read(tt1);
	s0=m+1,t0=m+2;
	int s1,s2,s3,s4;
	while(~scanf("%lld%lld%lld%lld",&s1,&s2,&s3,&s4)){
		a[s1]+=s3,a[s2]-=s3;
		add(s1,s2,s4-s3);
	}
	add(tt1,ss1,maxn);xx=esum;//printf("xx:%lld %lld\n",xx,e[xx].t);
	for(int i=1;i<=m;i++){
		if(a[i]<0)add(s0,i,-a[i]);
		else add(i,t0,a[i]),sum+=a[i];
	}
	int ans=0;
	ss=s0,tt=t0;
	while(check())ans+=dinic(ss,maxn);
	ss=tt1,tt=ss1;
	if(ans!=sum){
		printf("please go home to sleep\n");
		return 0;
	}
	ans=e[xx].v; 
	e[xx].v=0;
	e[xx^1].v=0;
	//printf("%lld %lld\n",sum,ans);
	while(check());ans-=dinic(ss,maxn);
	printf("%lld",ans);
	
	return 0;
}