Description
给定一张有向图,每条边都有一个容量C和一个扩容费用W。这里扩容费用是指将容量扩大1所需的费用。求: 1、 在不扩容的情况下,1到N的最大流; 2、 将1到N的最大流增加K所需的最小扩容费用。
Input
输入文件的第一行包含三个整数N,M,K,表示有向图的点数、边数以及所需要增加的流量。 接下来的M行每行包含四个整数u,v,C,W,表示一条从u到v,容量为C,扩容费用为W的边。
Output
输出文件一行包含两个整数,分别表示问题1和问题2的答案。
Sample Input
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
Sample Output
13 19
30%的数据中,N<=100
100%的数据中,N<=1000,M<=5000,K<=10
这道题的第一问是最大流,用Dinic可以解决,第二问只要在第一问的残量网络再次加入原先的边,并将这些边的容量设为INF(实际上k就可以了),价格设为为c,同时将残量网络中的价格设为0,之后设置一个源点与1号相连,最后求一遍0-n的最小费用最大流就可以了,下面是程序:
#include<stdio.h>
#include<string.h>
#include<iostream>
#define min(a,b) a<b?a:b
using namespace std;
const int N=50005,INF=0x3f3f3f3f;
struct queue{
int l,r,a[N];
void clear(){
l=0,r=1;
}
bool empty(){
return l+1==r;
}
int front(){
return a[(l+1)%N];
}
void pop(){
l=(l+1)%N;
}
void push(int x){
a[r]=x;
r=(r+1)%N;
}
}q;
struct edge{
int u,v,w,c,next;
}e[N];
int head[N],last[N],d[N],k,s,t;
bool vis[N];
void add(int u,int v,int w,int c){
e[k]=(edge){u,v,w,c,head[u]};
head[u]=k++;
e[k]=(edge){v,u,0,-c,head[v]};
head[v]=k++;
}
bool BFS(){
memset(d,0,sizeof(d));
q.clear();
q.push(s);
d[s]=1;
int u,i,v;
while(!q.empty()){
u=q.front();
q.pop();
for(i=head[u];i!=-1;i=e[i].next){
if(!d[e[i].v]&&e[i].w){
d[e[i].v]=d[u]+1;
q.push(e[i].v);
}
}
}
return d[t];
}
int dfs(int u,int f){
if(u==t||f==0){
return f;
}
int use=0,i,t;
for(i=head[u];i!=-1;i=e[i].next){
if(d[u]+1==d[e[i].v]&&(t=dfs(e[i].v,min(f-use,e[i].w)))){
e[i].w-=t;
e[i^1].w+=t;
use+=t;
if(use==f){
return f;
}
}
}
return use;
}
bool SPFA(){
memset(d,0x3f,sizeof(d));
q.clear();
q.push(s);
vis[s]=1;
d[s]=0;
int u,v,i;
while(!q.empty()){
u=q.front();
q.pop();
vis[u]=0;
for(i=head[u];i!=-1;i=e[i].next){
v=e[i].v;
if(d[u]+e[i].c<d[v]&&e[i].w){
d[v]=d[u]+e[i].c;
last[v]=i;
if(!vis[v]){
q.push(v);
vis[v]=1;
}
}
}
}
return d[t]!=INF;
}
int cost(){
int i,ans=0,mn=INF;
for(i=t;i!=s;i=e[last[i]].u){
mn=min(e[last[i]].w,mn);
}
for(i=t;i!=s;i=e[last[i]].u){
ans+=mn*e[last[i]].c;
e[last[i]].w-=mn;
e[last[i]^1].w+=mn;
}
return ans;
}
int dinic(){
int ans=0;
while(BFS()){
ans+=dfs(s,INF);
}
return ans;
}
int mcf(){
int ans=0;
while(SPFA()){
ans+=cost();
}
return ans;
}
int main(){
int n,m,q,u,v,w,c;
scanf("%d%d%d",&n,&m,&q);
s=1,t=n;
memset(head,-1,sizeof(head));
while(m--){
scanf("%d%d%d%d",&u,&v,&w,&c);
add(u,v,w,c);
}
printf("%d ",dinic());
add(0,1,q,0);
c=k-3;
for(u=0;u<=c;u+=2){
add(e[u].u,e[u].v,INF,e[u].c);
e[u].c=0;
e[u^1].c=0;
}
s=0;
printf("%d\n",mcf());
return 0;
}