Pagini recente » Cod sursa (job #2046935) | Cod sursa (job #284502) | Istoria paginii runda/hard/clasament | Cod sursa (job #1952858) | Cod sursa (job #2889339)
#include <bits/stdc++.h>
using namespace std;
ifstream fin("cuplaj.in");
ofstream fout("cuplaj.out");
const int INF = 1e9;
const int NMAX = 20000;
int N, M, E, ans, s, d;
struct FlowEdge {
int v, u;
long long cap, flow = 0;
FlowEdge(int v, int u, long long cap) : v(v), u(u), cap(cap) {}
};
struct Dinic {
const long long flow_inf = 1e18;
vector<FlowEdge> edges;
vector<vector<int>> adj;
int n, m = 0;
int s, t;
vector<int> level, ptr;
queue<int> q;
Dinic(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, long long 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;
}
long long dfs(int v, long long 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;
long long 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;
}
long long flow() {
long long 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 (long long pushed = dfs(s, flow_inf)) {
f += pushed;
}
}
return f;
}
};
int main() {
fin >> N >> M >> E;
s = 0;
d = N + M + 1;
Dinic solver(N + M + 2, s, d);
for(int i = 1; i <= E; i++) {
int u, v;
fin >> v >> u;
u += N;
solver.add_edge(v, u, 1);
}
for(int node = 1; node <= N; node++) {
solver.add_edge(s, node, 1);
}
for(int node = N + 1; node <= N + M; node++) {
solver.add_edge(node, d, 1);
}
fout << solver.flow() << '\n';
for(auto it: solver.edges) {
if(it.cap == 0) {
continue;
}
if(it.v == s || it.u == d) {
continue;
}
if(it.flow == it.cap) {
fout << it.v << " " << it.u - N << '\n';
}
}
return 0;
}