【bzoj1834】 network 网络扩容 【ZJOI2010】

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;
}