Pagini recente » Istoria paginii utilizator/buildersleagueunited | Istoria paginii runda/eusebiu_oji_2005_cls11-12/clasament | Cod sursa (job #2132546) | Cod sursa (job #2766678) | Cod sursa (job #2499336)
#include <iostream>
#include <cstdio>
#include <vector>
#include <queue>
using namespace std;
struct FlowEdge
{
int v, u;
int cap, flow = 0;
FlowEdge(int v, int u, int cap) : v(v), u(u), cap(cap) {}
};
struct Dinic
{
const int flow_inf = 1e18;
vector<FlowEdge> edges;
vector<vector<int>> adj;
int n, m = 0;
int s, t;
vector<int> level, ptr;
queue<int> q;
void init(int _n, int _s, int _t)
{
n = _n;
s = _s;
t = _t;
adj.resize(n);
level.resize(n);
ptr.resize(n);
}
void add_edge(int v, int u, int cap)
{
edges.emplace_back(v, u, cap);
edges.emplace_back(u, v, 0);
adj[v].push_back(m);
adj[u].push_back(m + 1);
m += 2;
}
bool bfs()
{
while (!q.empty())
{
int v = q.front();
q.pop();
for (int id : adj[v])
{
if (edges[id].cap - edges[id].flow < 1)
continue;
if (level[edges[id].u] != -1)
continue;
level[edges[id].u] = level[v] + 1;
q.push(edges[id].u);
}
}
return level[t] != -1;
}
int dfs(int v, int pushed)
{
if (pushed == 0)
return 0;
if (v == t)
return pushed;
for (int& cid = ptr[v]; cid < (int)adj[v].size(); cid++)
{
int id = adj[v][cid];
int u = edges[id].u;
if (level[v] + 1 != level[u] || edges[id].cap - edges[id].flow < 1)
continue;
int tr = dfs(u, min(pushed, edges[id].cap - edges[id].flow));
if (tr == 0)
continue;
edges[id].flow += tr;
edges[id ^ 1].flow -= tr;
return tr;
}
return 0;
}
int flow()
{
int f = 0;
while (true)
{
fill(level.begin(), level.end(), -1);
level[s] = 0;
q.push(s);
if (!bfs())
break;
fill(ptr.begin(), ptr.end(), 0);
while (int pushed = dfs(s, flow_inf))
{
f += pushed;
}
}
return f;
}
};
int main()
{
freopen ("maxflow.in", "r", stdin);
freopen ("maxflow.out", "w", stdout);
int n, m;
cin >> n >> m;
Dinic a;
a.init(n, 0, n - 1);
while (m--)
{
int x, y, z;
cin >> x >> y >> z;
x--;
y--;
a.add_edge(x, y, z);
}
cout << a.flow() << "\n";
}