Pagini recente » Cod sursa (job #67172) | Cod sursa (job #2119328) | Cod sursa (job #343995) | Cod sursa (job #677812) | Cod sursa (job #2961542)
#include <bits/stdc++.h>
#define MAX 131072
using namespace std;
const int NMAX = 20300;
struct obj {
int node;
int capacity, flow;
obj *revNode;
};
struct obj_1 {
int node, idx;
};
int N1, N2, M, S, D, ansFlow;
bool wasSeen[NMAX];
obj_1 dad[NMAX];
vector <obj> edges[NMAX];
vector <int> directedEdges[NMAX];
// flow.first = capacity
// flow.second = flow
void read(){
scanf("%d%d%d", &N1, &N2, &M);
int x, y;
S = 0;
D = N1 + N2 + 1;
for(int i = 1; i <= M; i++){
scanf("%d%d", &x, &y);
directedEdges[x].push_back(N1 + y);
edges[x].push_back({N1 + y, 1, 0, nullptr});
edges[N1 + y].push_back({x, 1, 0, nullptr});
edges[x].back().revNode = &edges[N1 + y].back();
edges[N1 + y].back().revNode = &edges[x].back();
}
for(int node = 1; node <= N1; node++){
directedEdges[S].push_back(node);
edges[S].push_back({node, 1, 0, nullptr});
edges[node].push_back({S, 1, 0, nullptr});
edges[S].back().revNode = &edges[node].back();
edges[node].back().revNode = &edges[S].back();
}
for(int node = N1 + 1; node <= N1 + N2; node++){
directedEdges[node].push_back(D);
edges[D].push_back({node, 1, 0, nullptr});
edges[node].push_back({D, 1, 0, nullptr});
edges[D].back().revNode = &edges[node].back();
edges[node].back().revNode = &edges[D].back();
}
}
// rezolvare cerinta a
// aplic algoritmul de flux, pornesc un BFS din nodul 1 si incerc sa ajung in nodul N pe muchiile care inca nu
// si-au atins capacitatea maxima, iar pentru fiecare nod retin nodul de unde am venit
// pentru toti vecinii lui N, calculez fluxul cu care pot sa ajung intr-un vecin de a lui, uitandu-ma la tatii acestuia
// insumez toate fluxurile gasite si obtin rezultatul
void printMinimumCut(int node){
for(int i = 0; i <= D; i++)
wasSeen[i] = false;
vector <pair <int, int> > cutEdges;
queue <int> Q;
Q.push(node);
wasSeen[node] = true;
while(!Q.empty()){
node = Q.front();
Q.pop();
}
for(auto it : cutEdges)
printf("%d %d\n", it.first, it.second);
}
bool generateFlow(int node){
for(int i = 0; i <= D; i++){
wasSeen[i] = false;
dad[i] = {0, 0};
}
queue <int> Q;
Q.push(node);
wasSeen[node] = true;
while(!Q.empty()){
node = Q.front();
Q.pop();
for(int i = 0; i < edges[node].size(); i++){
auto it = edges[node][i];
if(!wasSeen[it.node] && it.capacity - it.flow > 0){
wasSeen[it.node] = true;
dad[it.node] = {node, i};
Q.push(it.node);
}
}
}
if(!dad[D].node)
return false;
node = D;
for(auto it : edges[node]){
int MinFlow = it.revNode->capacity - it.revNode->flow;
int PosFlow = 0;
if(MinFlow > 0){
for(int j = it.node; j != 0; j = dad[j].node){
auto it = edges[dad[j].node][dad[j].idx];
PosFlow = it.capacity - it.flow;
if(PosFlow < MinFlow)
MinFlow = PosFlow;
}
it.revNode->flow += MinFlow;
it.flow -= MinFlow;
for(int j = it.node; j != 0; j = dad[j].node){
auto it = edges[dad[j].node][dad[j].idx];
it.flow += MinFlow;
it.revNode->flow -= MinFlow;
}
ansFlow += MinFlow;
}
}
return true;
}
int main() {
freopen("cuplaj.in", "r", stdin);
freopen("cuplaj.out", "w", stdout);
read();
bool repeat = true;
while(repeat)
repeat = generateFlow(0);
printf("%d\n", ansFlow);
printMinimumCut(0);
return 0;
}