Cod sursa(job #2741255)

Utilizator smoc_georgemarianSmoc George-Marian smoc_georgemarian Data 15 aprilie 2021 19:37:31
Problema Flux maxim Scor 100
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.46 kb
#include <bits/stdc++.h>
#define NMAX 1009
#define INF 9999999
using namespace std;
ifstream fin("maxflow.in");
ofstream fout("maxflow.out");
vector<int>g[NMAX];
int c[NMAX][NMAX];
int f[NMAX][NMAX];
bool uz[NMAX];
int tata[NMAX];
int n, m;
void citire();
bool bfs();
void solve();
int main()
{
	citire();
	solve();
}
void citire()
{
	int i;
	int x, y,cap;
	fin >> n >> m;
	for (i = 1; i <= m; i++)
	{
		fin >> x >> y>>cap;
		g[x].push_back(y);
		g[y].push_back(x);
		c[x][y] += cap;

	}
	

}
void solve()
{
	int i, j;
	int rez = 0;
	while (bfs()) //am gasit ddrumul
	{
		for (i = 0; i < g[n].size(); i++)
		{
			int vec = g[n][i];
			int minflow = INF;
			if (c[vec][n] != f[vec][n] && uz[n]) /// daca aceasta e calea gasit
			{
				tata[n] = vec;
				for (j = n; j != 1; j = tata[j])
					minflow = min(minflow, c[tata[j]][j] - f[tata[j]][j]);
				if (!minflow)
					continue;
				for (j = n; j != 1; j = tata[j])
				{
					f[tata[j]][j] += minflow;
					f[j][tata[j]] -= minflow;

				}
				rez += minflow;
			}
		}
	}
	fout << rez;
}

bool bfs()
{
	queue<int>H;
	H.push(1);
	memset(uz, 0, sizeof(uz));
	uz[1] = 1;
	while (!H.empty())
	{
		int act = H.front(); H.pop();
		if (act == n)
			continue;

		for (int i = 0; i < g[act].size(); i++)
		{
			int vec = g[act][i];
			if (!uz[vec] && c[act][vec] > f[act][vec])
			{
				uz[vec] = 1;
				H.push(vec);
				tata[vec] = act;
			}
		}
	}
	return uz[n];
}