UOJ #79. 一般图最大匹配

http://uoj.ac/problem/79

有时间读读集训队论文理解一下。。。这里试用一下claris的板子

一般图最大匹配

#include<bits/stdc++.h>
using namespace std;

const int maxl=510;
const int N=510;

int n,m,ans;
int a[maxl];
vector<int> e[maxl];
struct Blossom{
    int mate[N],n,ret,nex[N],f[N],mark[N],vis[N],t;
    queue<int> Q;
    int F(int x){return x==f[x]?x:f[x]=F(f[x]); }
    void merge(int x,int y){ f[F(x)]=F(y); }
    int lca(int x,int y){
        for(t++;;swap(x,y)){
            if(~x){
                if(vis[x=F(x)]==t) return x; vis[x]=t;
                x=mate[x]!=-1?nex[mate[x]]:-1;
            }
        }
    }
    void group(int a,int p){
        for(int b,c;a!=p;merge(a,b),merge(b,c),a=c){
            b=mate[a],c=nex[b];
            if(F(c)!=p) nex[c]=b;
            if(mark[b]==2) mark[b]=1,Q.push(b);
            if(mark[c]==2) mark[c]=1,Q.push(c);
        }
    }
    void aug(int s,const vector<int> G[]){
        for(int i=0;i<n;++i)
            nex[i]=vis[i]=-1,f[i]=i,mark[i]=0;
        while(!Q.empty()) Q.pop();
        Q.push(s); mark[s]=1;
        while(mate[s]==-1&&!Q.empty()){
            int x=Q.front(); Q.pop();
            for(int i=0,y;i<(int)G[x].size();++i){
                if((y=G[x][i])!=mate[x]&&F(x)!=F(y)&&mark[y]!=2){
                    if(mark[y]==1){
                        int p=lca(x,y);
                        if(F(x)!=p) nex[x]=y;
                        if(F(y)!=p) nex[y]=x;
                        group(x,p); group(y,p);
                    }else if(mate[y]==-1){
                        nex[y]=x;
                        for(int j=y,k,l;~j;j=l){
                            k=nex[j]; l=mate[k];
                            mate[j]=k; mate[k]=j;
                        }
                        break;
                    }else{
                        nex[y]=x; Q.push(mate[y]);
                        mark[mate[y]]=1; mark[y]=2;
                    }
                }
            }
        }
    }
    void solve(int _n,const vector<int> G[]){
        n=_n; for(int i=0;i<n;++i) mate[i]=-1;
        for(int i=t=0;i<n;++i) if(mate[i]==-1) aug(i,G);
        for(int i=ret=0;i<n;++i) 
		{
			a[i]=mate[i];
			if(mate[i]>i) ans++;
		}
    }
}blo;

inline void prework()
{
	scanf("%d%d",&n,&m);
	int u,v;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&u,&v);
		--u;--v;
		e[u].push_back(v);
		e[v].push_back(u);
	}
}

inline void mainwork()
{
	blo.solve(n,e);
}

inline void print()
{
	printf("%d\n",ans);
	for(int i=0;i<n;i++)
		printf("%d%c",a[i]+1," \n"[i==n-1]);
}

int main()
{
	prework();
	mainwork();
	print();
	return 0;
}