洛谷P3381 【模板】最小费用最大流(模板)

题目描述

如题,给出一个网络图,以及其源点和汇点,每条边已知其最大流量和单位流量费用,求出其网络最大流和在最大流情况下的最小费用。

输入输出格式

输入格式:

第一行包含四个正整数N、M、S、T,分别表示点的个数、有向边的个数、源点序号、汇点序号。

接下来M行每行包含四个正整数ui、vi、wi、fi,表示第i条有向边从ui出发,到达vi,边权为wi(即该边最大流量为wi),单位流量的费用为fi。

输出格式:

一行,包含两个整数,依次为最大流量和在最大流量情况下的最小费用。

输入输出样例

输入样例#1:  复制
4 5 4 3
4 2 30 2
4 3 20 3
2 3 20 1
2 1 30 9
1 3 40 5
输出样例#1:  复制
50 280

说明

时空限制:1000ms,128M

(BYX:最后两个点改成了1200ms)

数据规模:

对于30%的数据:N<=10,M<=10

对于70%的数据:N<=1000,M<=1000

对于100%的数据:N<=5000,M<=50000

样例说明:

如图,最优方案如下:

第一条流为4-->3,流量为20,费用为3*20=60。

第二条流为4-->2-->3,流量为20,费用为(2+1)*20=60。

第三条流为4-->2-->1-->3,流量为10,费用为(2+9+5)*10=160。

故最大流量为50,在此状况下最小费用为60+60+160=280。

故输出50 280。



题解:

最小费用最大流模板。

费用流就是在最大流的基础上用spfa算出最小费用就行了(因为跟最大流一样用了反向边,正确性可以保证)

抄了hzwer的模板(加了一些小小的优化,跑的还挺快):点击打开链接

为了更好的体现我对hzwer的崇拜,我还强行改了一波我的邻接表写法(改成了用结构体存,不然数组要开太多)


代码(洛谷1340ms):

#include<bits/stdc++.h>
using namespace std;
const int INF=500000000;
struct aaa{
	int from,to,v,c,next;
}e[200001];
int n,m,S,T,ans,sum,tot=2,g[10001],dis[10001],head[200001],q[1000001],from[200001];
void build(int u,int v,int t,int k){
	e[tot].from=u;e[tot].to=v;e[tot].v=t;e[tot].c=k;
	e[tot].next=head[u];head[u]=tot;
	tot++;
}
int spfa(int S,int T){
	int i,l=0,r=1,now;
	for(i=0;i<=n;i++)dis[i]=INF;
	dis[S]=0;q[0]=S;g[S]=1;
	while(l!=r){
		now=q[l++];
		if(l==1000001)l=0;
		for(i=head[now];i;i=e[i].next)
		if(e[i].v&&e[i].c+dis[now]<dis[e[i].to]){
		    dis[e[i].to]=dis[now]+e[i].c;
			from[e[i].to]=i;
			if(!g[e[i].to]){
				g[e[i].to]=1;
			if(dis[e[i].to]<dis[q[l]]){
				l--;if(l==-1)l=1000000;
				q[l]=e[i].to;
			}	
			 else{
			 	q[r++]=e[i].to;
			 	if(r==1000001)r=0;
			 }
		}
		}
		g[now]=0;
	}
	return dis[T]!=INF;
}
void MCMF(){
	while(spfa(S,T)){
		int i,x=INF;
		for(i=from[T];i;i=from[e[i].from])x=min(e[i].v,x);
		sum+=x;
		for(i=from[T];i;i=from[e[i].from]){
			ans+=x*e[i].c;
			e[i].v-=x;e[i^1].v+=x;
		}
	}
}
int main(){
	int i,u,v,t,k;
	scanf("%d%d%d%d",&n,&m,&S,&T);
	for(i=1;i<=m;i++){
		scanf("%d%d%d%d",&u,&v,&t,&k);
		build(u,v,t,k);build(v,u,0,-k);
	}
	MCMF();
	printf("%d %d\n",sum,ans);
}

猜你喜欢

转载自blog.csdn.net/qq_41510496/article/details/80960820