Cod sursa(job #2027695)

Utilizator tziplea_stefanTiplea Stefan tziplea_stefan Data 26 septembrie 2017 16:34:59
Problema Cuplaj maxim in graf bipartit Scor 10
Compilator cpp Status done
Runda Arhiva educationala Marime 2.02 kb
#include <fstream>
#include <vector>
#include <queue>
#define VAL 2005
#define INF 1000000000
#define F first
#define S second

using namespace std;

ifstream fin("cuplaj.in");
ofstream fout("cuplaj.out");

int N, M, K, i, j, X, be, en;
int prec[VAL], cop, ANS, nod;
int flux[VAL][VAL];
bool viz[VAL];
pair <int, int> MC[VAL];
vector <int> G[VAL];
vector <int> :: iterator it;
queue <int> Q;

bool BFS()
{
    for (i=1; i<=en; i++)
      viz[i]=false;
    Q.push(be);
    viz[be]=true;
    while (Q.empty()==false)
    {
        nod=Q.front();
        Q.pop();
        if (nod==en)
          continue;
        for (i=0; i<G[nod].size(); i++)
        {
            if (viz[G[nod][i]]==false && flux[nod][G[nod][i]]==1)
            {
                viz[G[nod][i]]=true;
                prec[G[nod][i]]=nod;
                Q.push(G[nod][i]);
            }
        }
    }
    return viz[en]==true;
}

void MaxFlow()
{
    while (BFS()==true)
    {
        for (j=0; j<G[en].size(); j++)
        {
            prec[en]=G[en][j];
            X=INF;
            for (i=en; i!=be; i=prec[i])
              X=min(X, flux[prec[i]][i]);
            if (X==0)
              continue;
            for (i=en; i!=be; i=prec[i])
            {
                flux[prec[i]][i]-=X;
                flux[i][prec[i]]+=X;
            }
            ANS+=X;
        }
    }
}

int main()
{
    fin >> N >> M >> K;
    for (i=1; i<=K; i++)
    {
        fin >> MC[i].F >> MC[i].S;
        MC[i].S+=N;
        G[MC[i].F].push_back(MC[i].S);
        G[MC[i].S].push_back(MC[i].F);
        flux[MC[i].F][MC[i].S]=1;
    }
    be=N+M+1;
    en=N+M+2;
    for (i=1; i<=N; i++)
    {
        flux[be][i]=1;
        G[be].push_back(i);
        G[i].push_back(be);
    }
    for (i=N+1; i<=N+M; i++)
    {
        flux[i][en]=1;
        G[en].push_back(i);
        G[i].push_back(en);
    }
    MaxFlow();
    fout << ANS << '\n';
    for (i=1; i<=K; i++)
      if (flux[MC[i].F][MC[i].S]==0)
        fout << MC[i].F << " " << MC[i].S-N << '\n';
    fin.close();
    fout.close();
    return 0;
}