【ZJOI2010】#1762 网络扩容 【最小费用流+拆边】

描述

给定一张有向图，每条边都有一个容量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;
}