Cod sursa(job #1953464)

Utilizator sebi110Ciobanu Sebastian sebi110 Data 4 aprilie 2017 20:39:29
Problema Cuplaj maxim in graf bipartit Scor 8
Compilator cpp Status done
Runda Arhiva educationala Marime 1.89 kb
#include <fstream>
#include <vector>
#include <queue>
#include <string.h>
#define INF 2000000000
using namespace std;
ifstream fin("cuplaj.in");
ofstream fout("cuplaj.out");
int c[1001][1001],f[1001][1001],p[1001];
bool viz[1001];
vector <int> vec[1001];
void getAP(int x)
{
    queue <int> q;
    q.push(x);
    viz[x]=1;
    while(!q.empty())
    {
        x=q.front();q.pop();
        int vecin;
        for(int i=0;i<vec[x].size();i++)
        {
            vecin=vec[x][i];
            if(viz[vecin]==0 && c[x][vecin]-f[x][vecin]>0)
            {
                p[vecin]=x;
                q.push(vecin);
                viz[vecin]=1;
            }
        }
    }
}
int main()
{
    int n,m,x,y,flux,i,maxf,muchii,S,T;
    fin>>n>>m>>muchii;
    for(i=1;i<=muchii;i++)
    {
        fin>>x>>y;
        vec[x].push_back(y+n+1);
        vec[y+n+1].push_back(x);
        c[x][y+n+1]=1;
    }
    S=0;
    for(i=1;i<=n;i++){
        vec[0].push_back(i);
        vec[i].push_back(0);
        c[0][i]=1;
    }
    T=n+m+1;
    for(i=1;i<=m;i++){
        vec[i+n+1].push_back(T);
        c[i+n+1][T]=1;
    }
    int F=0;
    while(1)
    {
        memset(viz,0,sizeof(viz));
        memset(p,-1,sizeof(p));
        getAP(S);
        if(p[T]==-1)
            break;
        maxf=INF;
        i=T;
        while(i!=S)
        {
            maxf=min(maxf,c[p[i]][i]-f[p[i]][i]);
            i=p[i];
        }
        i=T;
        while(i!=S)
        {
            f[p[i]][i]+=maxf;
            f[i][p[i]]-=maxf;
            i=p[i];
        }
        F=F+maxf;
    }
    fout<<F<<'\n';
    for(i=1;i<=n;i++)
    {
        int vecin,j;
        for(j=0;j<vec[i].size();j++)
        {
            vecin=vec[i][j];
            if(f[i][vecin]==1)
            {
                fout<<i<<' '<<vecin-n-1<<'\n';
            }
        }
    }
    return 0;
}