最小费用最大流模板WA最后四个点求调
查看原帖
最小费用最大流模板WA最后四个点求调
347664
菲斯斯夫斯基楼主2024/1/27 22:06

rt,谢谢。

#include<bits/stdc++.h>
using namespace std;
const int N=5e4+10;
int n,m,s,t,cnt=1;
int head[N<<1],in[N],val[N];
struct edge
{
	int to,c,v,next;
}a[N<<1];
struct ccf
{
	int to,x;
}f[N<<1];
void add(int x,int y,int z,int w)
{
	a[++cnt]={y,z,w,head[x]};
	head[x]=cnt;
}
bool spfa()
{
	queue<int>q;
	memset(val,127,sizeof(val));
	memset(in,0,sizeof(in));
	memset(f,-1,sizeof(f));
	q.push(s);
	val[s]=0;
	while(!q.empty())
	{
		int x=q.front();
		q.pop();
		in[x]=1;
		for(int i=head[x];i;i=a[i].next)
		{
			int to=a[i].to;
			if(!a[i].c||val[to]<=val[x]+a[i].v)continue;
			val[to]=val[x]+a[i].v;
			f[to]={x,i};
			if(!in[to])q.push(to);
			in[to]=1;
		}
	}
	return val[t]!=val[N-1];
}
void EK()
{
	int ans=0,sum=0;
	while(spfa())
	{
		int mi=1e9;
		for(int i=t;i!=s;i=f[i].to)
			mi=min(mi,a[f[i].x].c);
		for(int i=t;i!=s;i=f[i].to)
		{
			a[f[i].x].c-=mi;
			a[f[i].x^1].c+=mi;
		}
		ans+=mi,sum+=mi*val[t];
	}
	cout<<ans<<' '<<sum;
}
int main()
{
//	freopen("P3381_8.in","r",stdin);
	cin>>n>>m>>s>>t;
	for(int i=1;i<=m;i++)
	{
		int x,y,z,w;
		scanf("%d%d%d%d",&x,&y,&z,&w);
		add(x,y,z,w);
		add(y,x,0,-w); 
	}
	EK();
	return 0;
}

2024/1/27 22:06
加载中...