Cod sursa(job #1637755)

Utilizator howsiweiHow Si Wei howsiwei Data 7 martie 2016 19:06:39
Problema Flux maxim Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.34 kb
#include <iostream>
#include <cstdio>
#include <vector>
#include <queue>
#include <functional>
using namespace std;
const int oo = 0x1f1f1f1f;

int main() {
#ifdef INFOARENA
	freopen("maxflow.in", "r", stdin);
	freopen("maxflow.out", "w", stdout);
#endif
	cin.sync_with_stdio(false);
	int n, m;
	cin >> n >> m;
	vector<vector<int>> al(n);
	vector<int> el;
	vector<int> cl;
	for (int i = 0; i < m; i++) {
		int u, v, c;
		cin >> u >> v >> c;
		u--, v--;
		al[u].push_back(el.size());
		el.push_back(v);
		cl.push_back(c);
		al[v].push_back(el.size());
		el.push_back(u);
		cl.push_back(0);
	}
	int s = 0, t = n-1;
	int mf = 0;
	vector<int> d;
	vector<vector<int>::iterator> cur(n);
	function<int(int,int)> dfs = [&] (int u, int c0) {
		if (u == t) return c0;
		int c = c0;
		for (; cur[u] != al[u].end(); ++cur[u]) {
			int i = *cur[u];
			if (!(d[u] > d[el[i]] && cl[i] > 0)) continue;
			int f = dfs(el[i], min(c, cl[i]));
			c -= f;
			cl[i] -= f;
			cl[i^1] += f;
			if (c == 0) break;
		}
		if (c > 0) d[u] = n;
		return c0-c;
	};
	while (true) {
		d.assign(n, n);
		d[t] = 0;
		queue<int> q;
		q.push(t);
		do {
			int v = q.front();
			q.pop();
			cur[v] = al[v].begin();
			for (int i: al[v]) if (d[el[i]] == n && cl[i^1] > 0) {
				d[el[i]] = d[v]+1;
				q.push(el[i]);
			}
		} while (!q.empty());
		if (d[s] == n) break;
		mf += dfs(s, oo);
	}
	printf("%d\n", mf);
}