Cod sursa(job #299180)

Utilizator hasegandaniHasegan Daniel hasegandani Data 6 aprilie 2009 16:53:22
Problema Cuplaj maxim in graf bipartit Scor 40
Compilator cpp Status done
Runda Arhiva educationala Marime 2.55 kb
#include<stdio.h>
#include<string.h>

#define Nmax 20010

struct nod
{
    int drum,cost;
    nod *urm;
} *l[Nmax];
int n,n2,e,v[Nmax],w[Nmax],viz[Nmax],poz,in,sf;

void add(nod *&,int,int);
void addb(nod *&,int,int);
void ford();
void afis();
int DF(int);

int main()
{
    int i,a,b;
    freopen("cuplaj.in","r",stdin);
    freopen("cuplaj.out","w",stdout);
    scanf("%d%d%d",&n,&n2,&e);
    for(i=1;i<=e;++i)
        {
            scanf("%d%d",&a,&b);
            add(l[a],b+n,1);
            add(l[b+n],a,0);
        }
    in=n+n2+1;
    sf=in+1;
    for(i=1;i<=n;++i)
        {
        add(l[in],i,1);
        add(l[i],in,0);
        }
    for(i=1;i<=n2;++i)
        {
        add(l[i+n],sf,1);
        add(l[sf],i+n,0);
        }
    
    ford();
    afis();
    return 0;
}

int DF(int k)
{
    nod *t=l[k];
    while (t)
        {
            if (t->cost && !viz[t->drum])
                {
                if (t->drum==sf)
                    {
                    v[++poz]=sf;
                    return 1;
                    }
                    
                ++viz[t->drum];
                v[++poz]=t->drum;
                
                int a=DF(t->drum);
                if (a)
                    return 1;

                --viz[t->drum];
                --poz;
                }
            t=t->urm;
        }
    return 0;
}

void ford()
{
    v[1]=in;
    viz[in]=1;
    poz=1;
    while ( DF(in) )
        {
            for(int i=2;i<=poz;++i)
                {
                addb(l[v[i-1]],v[i],-1);
                addb(l[v[i]],v[i-1],1);
                }
        poz=1;
        memset(viz,0,Nmax);
        viz[in]=1;
        }
}

void afis()
{
    int i;
    poz=0;
    for(i=1;i<=n;++i)
        {
            nod *t=l[i];
            while (t)
                {
                    if (t->cost==0 && t->drum!=in)
                        {
                            v[++poz]=i;
                            w[poz]=t->drum;
                            break;
                        }
                    
                    t=t->urm;
                }
        }
        
    printf("%d\n",poz);
    for(i=1;i<=poz;++i)
        printf("%d %d\n",v[i],w[i]-n);
}

void add(nod *&inc,int info1,int info2)
{
    nod *p=new nod;
    p->drum=info1;
    p->cost=info2;
    p->urm=inc;
    inc=p;
}

void addb(nod *&inc,int info1,int info2)
{
    nod *t=inc;
    while (t->drum!=info1 && t)
        t=t->urm;
    if (t)
        t->cost+=info2;
}