luogu2604网络扩容

题目

　　https://www.luogu.org/problem/show?pid=2604

题解

　　建边的时候先建正常的边，容量为c，权值为0。再另外建一条容量为+∞，权值为扩增权值的边。先以1为S，N为T跑一边正常的最大流（此时要忽略那些有权值的边）。然后令T=N+1，连一条从N到T的边，容量为K，权值为0，此时跑最小费用最大流就是第二问的答案。（这一步要在残余网络上进行）

代码

//最小费用最大流
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <queue>
#define inf 0x3f3f3f3f
#define maxn 100000
using namespace std;
int tot=1, head[maxn], to[maxn], w[maxn], c[maxn], d[maxn], in[maxn], cost, N, M, S, T,
	pre[maxn], K, next[maxn];
queue<int> q;
void adde(int a, int b, int cc, int ww)
{to[++tot]=b;c[tot]=cc;w[tot]=ww;next[tot]=head[a];head[a]=tot;}
void adde2(int a, int b, int cc, int ww){adde(a,b,cc,ww);adde(b,a,0,-ww);}
bool spfa(int flag)
{
	int x, p, i;
	for(i=1;i<=T;i++)d[i]=inf;
	in[S]=1, d[S]=0, q.push(S);
	while(!q.empty())
	{
		x=q.front();q.pop();
		for(p=head[x];p;p=next[p])
			if(c[p] and d[to[p]]>d[x]+w[p])
			{
				if(flag==1 and w[p])continue;
				d[to[p]]=d[x]+w[p];
				pre[to[p]]=p;
				if(!in[to[p]])in[to[p]]=1,q.push(to[p]);
			}
		in[x]=0;
	}
	return d[T]!=inf;
}
int augment()
{
	int x, f=inf;
	for(x=T;x!=S;x=to[pre[x] xor 1])f=min(f,c[pre[x]]);
	for(x=T;x!=S;x=to[pre[x] xor 1])c[pre[x]]-=f, c[pre[x] xor 1]+=f;
	cost+=f*d[T];
	return f;
}
void input()
{
	int i, a, b, cc, ww;
	scanf("%d%d%d",&N,&M,&K);
	for(i=1;i<=M;i++)
	{
		scanf("%d%d%d%d",&a,&b,&cc,&ww);
		adde2(a,b,cc,0),adde2(a,b,inf,ww);
	}
	S=1, T=N;
}
void solve()
{
	int flow;
	flow=cost=0;
	while(spfa(1))flow+=augment();
	adde2(N,T=N+1,K,0);
	while(spfa(2))augment();
	printf("%d %d",flow,cost);
}
int main()
{
	input();
	solve();
	return 0;
}