UOJ 171 【WC2016】挑战NPC

　　一开始还真没想到是一般图匹配这种模型(毕竟才会的带花树)

　　把每一个盒子拆成3个，每一个可以放置进它的小球分别向这三个点连边，然后这三个点在连成一个三元环，最终答案就是小球数目-匹配数。

　　由于是一般图，所以套一个带花树就可以了。

　　

　　NOTICE：寻找增广路时，应该从球先找起，这样子才保证了每个球有地方放置。

　　

---

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<cstring>
using namespace std;
#define maxn 10010
#define llg long long 
#define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
llg n,m,sett,father[maxn],pre[maxn],match[maxn],vis[maxn],id[maxn],dl[maxn*10],head,tail,ans,T,totn,e;
llg c[1000][1000];

vector<llg> a[maxn];

llg find(llg x){if (father[x]!=x) father[x]=find(father[x]); return father[x];}

llg lca(llg x,llg y)
{
    sett++;
    while (vis[x]!=sett)
    {
        if (x)
        {
            x=find(x);
            if (vis[x]==sett) return x;
            vis[x]=sett;
            if (match[x]!=0) x=find(pre[match[x]]);
            else x=0;
        }
        
            swap(x,y);
    }
    return x;
}

void change(llg x,llg y,llg fa)
{
    llg z;
    while (find(x)!=fa)
    {
        pre[x]=y; z=match[x];
        if (id[z]==1) {id[z]=0,dl[++tail]=z;}
        if (find(z)==z) father[z]=fa;
        if (find(x)==x) father[x]=fa;
        y=z; x=pre[y];
    }
}

bool bfs(llg p)
{
    llg x,w,v;
    for (llg i=1;i<=n;i++) id[i]=-1,father[i]=i;
    head=0,tail=1,dl[1]=p; id[p]=0;
    do
    {
        x=dl[++head],w=a[x].size();
        for (llg i=0;i<w;i++)
        {
            v=a[x][i];
            if (id[v]==-1)
            {
                pre[v]=x;
                id[v]=1;
                if (!match[v])
                {
                    llg last,t,now=v;
                    while (now!=0)
                    {
                        t=pre[now],last=match[t];
                        match[t]=now; match[now]=t;
                        now=last;// v=t;
                    }//把原lai的匹配和非匹配bian取反
                    return 1;
                }
                id[match[v]]=0,dl[++tail]=match[v];
            }
            else
                if (id[v]==0 && find(v)!=find(x))
                {
                    llg dad=lca(x,v);
                    change(x,v,dad);
                    change(v,x,dad);
                }
        }
    }while (head!=tail);
    return 0;

}

void link(llg x,llg y) {a[x].push_back(y),a[y].push_back(x);}

void init()
{
    for (llg i=1;i<=n;i++) a[i].clear(),vis[i]=match[i]=pre[i]=father[i]=id[i]=dl[i]=0;
    //for (llg i=1;i<=n;i++) for (llg j=1;j<=n;j++) c[i][j]=c[j][i]=0;
    scanf("%lld%lld%lld",&n,&m,&e);
    totn=n;
    llg x,y;
    for (llg i=1;i<=e;i++)
    {
        scanf("%lld%lld",&x,&y);
        
        llg po=m*3+x;
        link(po,y*3-2),link(po,y*3-1),link(po,y*3);
        //a[x].push_back(y),a[y].push_back(x);
    }
    n=n+m*3;
    for (llg i=1;i<=m;i++)
    {
        link(i*3-2,i*3-1),link(i*3-1,i*3),link(i*3,i*3-2);
    }
}

void oupt()
{
    cout<<ans-totn<<endl;
    for (llg i=1;i<=totn;i++)
    {
        printf("%lld ",(match[3*m+i]-1)/3+1);
    }
    printf("
");
}

int main()
{
    //yyj("a");
    cin>>T;
    while (T--)
    {
        ans=0;
        init();
        for (llg i=n;i>1;i--) 
            if (!match[i] && bfs(i)) ans++;
        oupt();
    }    
    return 0;
}