Pagini recente » Cod sursa (job #1815658) | Cod sursa (job #1942213) | Cod sursa (job #1527282) | Cod sursa (job #1055141) | Cod sursa (job #3258113)
// Hopcroft Karp
#include <iostream>
#include <fstream>
#include <algorithm>
#include <cmath>
#include <map>
#include <unordered_map>
#include <set>
#include <unordered_set>
#include <queue>
#include <stack>
#include <deque>
#include <iomanip>
#include <vector>
#include <cassert>
#pragma GCC optimize("O3")
#define fi first
#define se second
#define pb push_back
#define pf push_front
using namespace std;
ifstream fin ("cuplaj.in");
ofstream fout ("cuplaj.out");
typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
const ll Nmax=1e4+5, inf=1e9+5;
int n, m, e;
// avem nevoie doar de drumurile de la un nod stang la un nod drept,
// deorece ne dorim sa augmentam drumuir ce incep de la un nod stang
vector<int> ltor[Nmax];
// nodurile din dreapta cu care s-au cuplat nodurile din stanga
int pl[Nmax];
// nodurile din stanga cu care s-au cuplat nodurile din dreapta
int pr[Nmax];
// ce noduri au fost vizitate in dfs (si au fost sau nu augmentate)
bool vis[Nmax];
// cuplajul maxim
vector<pii> cuplaj;
bool augment(int nod){ // nod stang
for (auto it:ltor[nod]) // iteram prin vecinii din dreapta
if (!vis[it]){ // daca un vecin a fost deja vizitat, atunci ori nu se poate augmenta, ori a augmentat deja alt nod
vis[it]=1;
if (!pr[it] || augment(pr[it])){ // daca nodul din dreapta nu are pereche, atunci il putem folosi la augmentare
pl[nod]=it; // update pentru cuplaj
pr[it]=nod;
return 1;
}
}
return 0;
}
int main()
{
fin>>n>>m>>e;
int a, b;
for (int i=0; i<e; i++){
fin>>a>>b;
ltor[a].pb(b);
}
bool ok=1;
while (ok){
ok=0;
for (int i=1; i<=n; i++)
vis[i]=0;
for (int i=1; i<=n; i++)
if (!pl[i] && augment(i)) // a mers augmentarea
ok=1;
}
for (int i=1; i<=n; i++)
if (pl[i])
cuplaj.pb({i, pl[i]});
fout<<cuplaj.size()<<'\n';
for (auto it:cuplaj)
fout<<it.fi<<' '<<it.se<<'\n';
return 0;
}