P3381 【模板】最小费用最大流（EK算法）

最新推荐文章于 2024-10-01 04:50:10 发布

summ1ts

最新推荐文章于 2024-10-01 04:50:10 发布

阅读量110

点赞数 4

分类专栏：一些模版文章标签：算法 c++ 图论费用流 EK算法网络流

本文链接：https://blog.csdn.net/summ1ts/article/details/142644695

版权

一些模版专栏收录该内容

13 篇文章 0 订阅

订阅专栏

费用流模版，用 EK 跑的。

实际上就是把 bfs 的过程换成了 spfa。

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N=1e5+10,INF=0x3f3f3f3f;

int read(){
	int x=0,f=1;char ch=getchar();
	while(!isdigit(ch)){
		if(ch=='-') f=-1;
		ch=getchar();
	}
	while(isdigit(ch)) x=x*10+ch-'0',ch=getchar();
	return x*f;
}

int n,m,S,T,head[N],tot=1,dist[N],vis[N],pre[N],mf[N],flow,cost;
struct node{
	int to,nxt,w,c;
}edge[N];
void add(int x,int y,int w,int c){
	edge[++tot].to=y;
	edge[tot].w=w,edge[tot].c=c;
	edge[tot].nxt=head[x];
	head[x]=tot;
}

bool spfa(){
	memset(dist,0x3f,sizeof(dist)),memset(mf,0,sizeof(mf));
	queue<int> q;q.push(S),vis[S]=1,dist[S]=0,mf[S]=INF;
	while(!q.empty()){
		int t=q.front();q.pop(),vis[t]=0;
		for(int i=head[t];i;i=edge[i].nxt){
			int y=edge[i].to;
			if(dist[y]>dist[t]+edge[i].w&&edge[i].c){
				dist[y]=dist[t]+edge[i].w;
				pre[y]=i,mf[y]=min(mf[t],edge[i].c);
				if(!vis[y]) q.push(y),vis[y]=1;
			}
		}
	}
	return mf[T]>0;
}

void EK(){
	while(spfa()){
		int v=T;
		while(v!=S){
			int i=pre[v];edge[i].c-=mf[T],edge[i^1].c+=mf[T],v=edge[i^1].to;
		}
		flow+=mf[T],cost+=mf[T]*dist[T];
	}
}

signed main(){
	
	n=read(),m=read(),S=read(),T=read();
	for(int i=1;i<=m;i++){
		int x=read(),y=read(),c=read(),w=read();
		add(x,y,w,c),add(y,x,-w,0);
	}
	EK();
	cout<<flow<<" "<<cost<<endl;
	
	return 0;
}