描述
给定一张有向图,每条边都有一个容量C和一个扩容费用W。这里扩容费用是指将容量扩大1所需的费用。求:
1、 在不扩容的情况下,1到N的最大流;
2、 将1到N的最大流增加K所需的最小扩容费用。
输入
输入文件的第一行包含三个整数N,M,K,表示有向图的点数、边数以及所需要增加的流量。
接下来的M行每行包含四个整数u,v,C,W,表示一条从u到v,容量为C,扩容费用为W的边。
输出
输出文件一行包含两个整数,分别表示问题1和问题2的答案。
样例输入 [复制]
5 8 2
1 2 5 8
2 5 9 9
5 1 6 2
5 1 1 8
1 2 8 7
2 5 4 9
1 2 1 1
1 4 2 1
样例输出 [复制]
13 19
提示
30%的数据中,N<=100
100%的数据中,N<=1000,M<=5000,K<=10
思路:
显然 ”建“(改)图 跑最小费用流完事。。
只是注意”建图“常规思路:”按照功能区分不同的成分。“
此题跑完最大流后的残量网络中,有没有跑完的边,这些边容量跑完之后才”计费“
这时可以把 一条边的”计费“ 和 ”不计费“ 理解为两个不同的功能,”把一条边拆成两条“
在原来的图上建立”重边“,即
把原来的边 改为 容量无限,有费用 加入到原图中,
注:原来的边费用为0;
最后重新建立源点向1号连接一条 容量为 最大流+K 的边即可。(限制总流量)
#include<bits/stdc++.h>
using namespace std;
const int maxn=2e3,maxm=5e4,INF=2e9;
struct edge{
int v,nxt,c,f;
}e[maxm];
int head[maxn],cnt=1;
inline void _add(int u,int v,int c,int f){
e[++cnt]=(edge){v,head[u],c,f};
head[u]=cnt;
}
int s,t;
int dis[maxn];
bool vis[maxn];
inline bool spfa(int s,int t){
memset(vis,0,sizeof vis);
memset(dis,0x5f,sizeof dis);
dis[t]=0;
deque <int> q;
q.push_back(t);
while(!q.empty()){
int u=q.front();q.pop_front();vis[u]=0;
for(int i=head[u];i;i=e[i].nxt){
int v=e[i].v,ff=e[i].f;
if(e[i^1].c&&dis[v]>dis[u]-ff){
dis[v]=dis[u]-ff;
if(!vis[v]){
vis[v]=1;
if(!q.empty()&&dis[v]<dis[q.front()])q.push_front(v);
else q.push_back(v);
}
}
}
}
return dis[s]<0x5f5f5f5f;
}
int ans=0;
int dfs(int u,int f){
if(u==t){
vis[t]=1;
return f;
}
vis[u]=1;
int tot=0,w;
for(int i=head[u];i;i=e[i].nxt){
int v=e[i].v,&c=e[i].c,ff=e[i].f;;
if(!vis[v]&&c&&dis[v]==dis[u]-ff){
w=dfs(v,min(c,f-tot));
tot+=w,c-=w,e[i^1].c+=w,ans+=w*ff;
if(tot==f)break;
}
}
return tot;
}
int cflow(){
int ret=0;
while(spfa(s,t)){
vis[t]=1;
while(vis[t]){
memset(vis,0,sizeof vis);
ret+=dfs(s,INF);
}
}
return ret;
}
struct pi{
int a,b,c,d;
}p[maxm];
signed main(){
int n,m,K;scanf("%d%d%d",&n,&m,&K);
for(int i=1;i<=m;++i){
int a,b,c,d;scanf("%d%d%d%d",&a,&b,&c,&d);
p[i]=(pi){a,b,c,d};
_add(a,b,c,d),_add(b,a,0,-d);
}
s=1,t=n;
int ans1=cflow();
cnt=1;ans=0;
memset(head,0,sizeof head);
_add(0,1,ans1+K,0);_add(1,0,0,0);
for(int i=1;i<=m;++i){
int a=p[i].a,b=p[i].b,c=p[i].c,d=p[i].d;
_add(a,b,c,0),_add(b,a,0,0);
_add(a,b,INF,d),_add(b,a,0,-d);
}
s=0,t=n;
cflow();
printf("%d %d",ans1,ans);
return 0;
}