大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。
Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺
给定一个包含 n 个点 m 条边的有向图,并给定每条边的容量和费用,边的容量非负。
图中可能存在重边和自环,保证费用不会存在负环。
求从 S 到 T 的最大流,以及在流量最大时的最小费用。
输入格式
第一行包含四个整数 n,m,S,T。
接下来 m 行,每行三个整数 u,v,c,w,表示从点 u 到点 v 存在一条有向边,容量为 c,费用为 w。
点的编号从 1 到 n。
输出格式
输出点 S 到点 T 的最大流和流量最大时的最小费用。
如果从点 S 无法到达点 T 则输出 0 0。
数据范围
2≤n≤5000,
1≤m≤50000,
0≤c≤100,
−100≤w≤100
S≠T
输入样例:
5 5 1 5
1 4 10 5
4 5 5 10
4 2 12 5
2 5 10 15
1 5 10 10
输出样例:
20 300
#include<bits/stdc++.h>
using namespace std;
const int N = 5010, M = 2 * 100010;
const int INF = 1e8;
struct Edge{
int v,next,w,f;
}edge[M];
int n,m,s,e;
int head[N],cnt;
void add(int u,int v,int f,int w){
edge[cnt].v = v;
edge[cnt].f = f;
edge[cnt].w = w;
edge[cnt].next = head[u];
head[u] = cnt ++;
}
int d[N],q[N],hh = 0,tt = 0,pre[N],curf[N],st[N];
bool spfa(){
//最大流最大费用的话就求最长路即可
memset(d,0x3f,sizeof d);
memset(curf,0,sizeof curf);
memset(st,0,sizeof st);
hh = tt = 0;
d[s] = 0,q[tt ++] = s,curf[s] = INF;
st[s] = true;
while(hh != tt){
int t = q[hh ++];
if(hh == N)hh = 0;
st[t] = false;
for(int i = head[t];~i;i = edge[i].next){
int v = edge[i].v,w = edge[i].w;
if(edge[i].f && d[v] > d[t] + w){
d[v] = d[t] + w;
pre[v] = i;
curf[v] = min(edge[i].f,curf[t]);
if(!st[v]){
q[tt ++] = v;
if(tt == N)tt = 0;
st[v] = true;
}
}
}
}
return curf[e] > 0;
}
void EK(int &flow,int &cost){
flow = cost = 0;
while(spfa()){
int t = curf[e];
flow += t;cost += t * d[e];
for(int i = e;i != s;i = edge[pre[i] ^ 1].v){
edge[pre[i]].f -= t,edge[pre[i] ^ 1].f += t;
}
}
}
int main(){
cin>>n>>m>>s>>e;
int x,y,f,w;
memset(head,-1,sizeof head);
for(int i = 0;i < m;i ++){
scanf("%d %d %d %d",&x,&y,&f,&w);
add(x,y,f,w),add(y,x,0,-w);
}
int flow = 0,cost = 0;
EK(flow,cost);
printf("%d %d",flow,cost);
return 0;
}
发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/168567.html原文链接:https://javaforall.cn
【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛
【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...