洛谷P3381 最小费用最大流模板

记得把数组开大一点,不然就RE了。。。

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

int n,m,s,t,ans,d[N],maxflow;
int tot=1,adj[N],nex[M],to[M],cap[M],cost[M];
bool walk[N],vis[N];

void add(int u,int v,int w,int c){
	nex[++tot]=adj[u];
	adj[u]=tot;
	to[tot]=v;
	cap[tot]=w;
	cost[tot]=c;
}

bool SPFA(){
	queue<int> q;
	for(int i=1;i<=n;i++)
		vis[i]=false,walk[i]=false,d[i]=INF;
	q.push(s),d[s]=0,vis[s]=true;
	while(!q.empty()){
		int u=q.front();q.pop();
		vis[u]=false;
		for(int i=adj[u];i;i=nex[i]){
			int v=to[i];
			if(cap[i]>0&&d[u]+cost[i]<d[v]){
				d[v]=d[u]+cost[i];
				if(!vis[v]){
					vis[v]=true;
					q.push(v);
				}
			}
		}
	}
	return d[t]<INF;
}

int dinic(int u,int fw){
	if(u==t){
		ans+=fw*d[t];
		return fw;
	}
	walk[u]=true;
	int rest=fw;
	for(int i=adj[u];i&&rest;i=nex[i]){
		int v=to[i];
		if(!walk[v]&&cap[i]>0&&d[u]+cost[i]==d[v]){
			int k=dinic(v,min(cap[i],rest));
			if(k){
				cap[i]-=k;cap[i^1]+=k;
				rest-=k;
			}
		}
	}
	return fw-rest;
}

void solve(){
	maxflow=ans=0;
	while(SPFA()) maxflow+=dinic(s,INF);
}

signed main(){
	cin>>n>>m>>s>>t;
	for(int i=1;i<=m;i++){
		int u,v,w,c;
		cin>>u>>v>>w>>c;
		add(u,v,w,c);
		add(v,u,0,-c);
	}
	solve();
	cout<<maxflow<<" "<<ans;
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值