费用流_最小费用最大流模板_洛谷3381

题目：
题目描述
如题，给出一个网络图，以及其源点和汇点，每条边已知其最大流量和单位流量费用，求出其网络最大流和在最大流情况下的最小费用。

输入格式
第一行包含四个正整数N、M、S、T，分别表示点的个数、有向边的个数、源点序号、汇点序号。

接下来M行每行包含四个正整数ui、vi、wi、fi，表示第i条有向边从ui出发，到达vi，边权为wi（即该边最大流量为wi），单位流量的费用为fi。

输出格式
一行，包含两个整数，依次为最大流量和在最大流量情况下的最小费用。

输入输出样例
输入 #1 复制
4 5 4 3
4 2 30 2
4 3 20 3
2 3 20 1
2 1 30 9
1 3 40 5
输出 #1 复制
50 280
说明/提示
时空限制：1000ms,128M

（BYX：最后两个点改成了1200ms）

数据规模：

对于30%的数据：N<=10，M<=10

对于70%的数据：N<=1000，M<=1000

对于100%的数据：N<=5000，M<=50000

样例说明：

如图，最优方案如下：

第一条流为4–>3，流量为20，费用为3*20=60。

第二条流为4–>2–>3，流量为20，费用为（2+1）*20=60。

第三条流为4–>2–>1–>3，流量为10，费用为（2+9+5）*10=160。

故最大流量为50，在此状况下最小费用为60+60+160=280。

故输出50 280。

解题思路：
直接构图，然后跑最大流，模板

AC代码:

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>

using namespace std;

const int N = 5005, M = 1e5 + 5;
const int inf = 0x3f3f3f3f;

int n, m, s, t, ans, maxflow;
int h[N], e[M], w[M], c[M], ne[M], idx;
int d[N], pre[N], incf[N];
bool vis[N];

inline void add(int a, int b, int we, int co) {
	e[idx] = b, w[idx] = we, c[idx] = co, ne[idx] = h[a], h[a] = idx ++;
	e[idx] = a, w[idx] = 0, c[idx] = -co, ne[idx] = h[b], h[b] = idx ++;
}

inline bool spfa(void) {
	queue<int> q;
	while(q.size()) q.pop();
	memset(d, 0x3f3f3f3f, sizeof d);
	memset(vis, false, sizeof vis);
	
	incf[s] = inf; d[s] = 0; q.push(s); vis[s] = true;
	
	while(q.size()) {
		int u = q.front(); q.pop(); vis[u] = false;
		
		for(int i = h[u]; i != -1; i = ne[i]) {
			int v = e[i];
			if(w[i]) {
				if(d[v] > d[u] + c[i]) {
					d[v] = d[u] + c[i];
					pre[v] = i;
					incf[v] = min(incf[u], w[i]);
					if(!vis[v]) q.push(v), vis[v] = true;
				}
			} 
		}
	}
	
	if(d[t] == inf) return false;
	return true;
}

inline void update(void) {
	int x = t;
	while(x != s) {
		int i = pre[x];
		w[i] -= incf[t];
		w[i ^ 1] += incf[t];
		x = e[i ^ 1];
	}
	
	maxflow += incf[t];
	ans += d[t] * incf[t];
}

int main(void) {
	scanf("%d%d%d%d", &n, &m, &s, &t);
	memset(h, -1, sizeof h);
	for(int i = 1; i <= m; i ++) {
		int a, b, we, co;
		scanf("%d%d%d%d", &a, &b, &we, &co);
		add(a, b, we, co);
	}
	
	while(spfa()) update();
	
	printf("%d %d\n", maxflow, ans);
	
	return 0;
}