题目描述
给定一张有向图,每条边都有一个容量C和一个扩容费用W。这里扩容费用是指将容量扩大1所需的费用。求: 1、 在不扩容的情况下,1到N的最大流; 2、 将1到N的最大流增加K所需的最小扩容费用。
输入格式
输入文件的第一行包含三个整数N,M,K,表示有向图的点数、边数以及所需要增加的流量。 接下来的M行每行包含四个整数u,v,C,W,表示一条从u到v,容量为C,扩容费用为W的边。
输出格式
输出文件一行包含两个整数,分别表示问题1和问题2的答案。
题解:
第一问直接最大流即可,第二问设置超级源点连向1,容量为K,费用为0,表示限制最大容量拓展了K,然后利用残余网络,在每条边之间再建一条边,容量为INF,费用为co[i]。
因为要跑最小费用最大流,而原图是无花费的最大流,所以只需要再残余网络上继续增广,就可以保证抛出最小费用了
AC代码:
#pragma GCC optimize(2)
#include<bits/stdc++.h>
#include<ext/rope>
using namespace std;
using namespace __gnu_cxx;
#define LL long long
#define pii pair<int,int>
#define mp(a,b) make_pair(a,b)
const int MAXN = 2e5+50;
const int MAXM = 2e6+50;
const int MOD = 1e9+7;
const int INF = 0x3f3f3f3f;
int n,m,k,s,t,tot=1,head[MAXN],to[MAXM],w[MAXM],co[MAXM],nxt[MAXM],h[MAXN];
int a[MAXN],b[MAXN],c[MAXN],d[MAXN],vis[MAXN],dis[MAXN],flow[MAXN],pre[MAXN];
inline void ade(int u,int v,int ww,int cost){
to[++tot]=v; w[tot]=ww; co[tot]=cost; nxt[tot]=head[u]; head[u]=tot;
}
inline void add(int u,int v,int w,int cost){ ade(u,v,w,cost); ade(v,u,0,-cost); }
inline int bfs(){
queue<int> que; que.push(s); memset(h,0,sizeof(h)); h[s]=1;
while(!que.empty()){
int u=que.front(); que.pop();
for(int i=head[u];i;i=nxt[i]){
if(w[i] && !h[to[i]]){
h[to[i]]=h[u]+1; que.push(to[i]);
}
}
}
return h[t];
}
inline int dfs(int x,int f){
if(x==t) return f; int fl=0;
for(int i=head[x];i;i=nxt[i]){
if(w[i] && h[to[i]]==h[x]+1){
int mi=dfs(to[i],min(f,w[i]));
w[i]-=mi; w[i^1]+=mi; fl+=mi; f-=mi;
}
}
if(!fl) h[x]=-1;
return fl;
}
inline int dinic(){
int res=0;
while(bfs()) res+=dfs(s,INF);
return res;
}
inline int spfa(){
for(int i=0;i<=t;i++) dis[i]=INF;
queue<int> que; que.push(s); dis[s]=0; vis[s]=1; flow[s]=INF;
while(!que.empty()){
int u=que.front(); que.pop(); vis[u]=0;
for(int i=head[u];i;i=nxt[i]){
if(w[i] && dis[to[i]]>dis[u]+co[i]){
dis[to[i]]=dis[u]+co[i]; pre[to[i]]=i;
flow[to[i]]=min(flow[u],w[i]);
if(!vis[to[i]]) vis[to[i]]=1,que.push(to[i]);
}
}
}
return dis[t]!=INF;
}
inline int EK(){
int res=0;
while(spfa()){
res+=flow[t]*dis[t]; int x=t;
while(x!=s){
int i=pre[x];
w[i]-=flow[t]; w[i^1]+=flow[t]; x=to[i^1];
}
}
return res;
}
signed main(){
#ifndef ONLINE_JUDGE
freopen("C:\\Users\\Administrator\\Desktop\\in.txt","r",stdin);
#endif // ONLINE_JUDGE
scanf("%d%d%d",&n,&m,&k); s=1; t=n;
for(int i=1;i<=m;i++) scanf("%d%d%d%d",&a[i],&b[i],&c[i],&d[i]);
for(int i=1;i<=m;i++) add(a[i],b[i],c[i],0);
int res1=dinic(); s=0; add(s,1,k,0);
for(int i=1;i<=m;i++) add(a[i],b[i],1e8,d[i]);
int res2=EK();
printf("%d %d\n",res1,res2);
return 0;
}