带花树模板（一般图最大匹配）

wc出了一道带花树的题 然后暴力滚粗了…………
然后回家后就学了带花树…………
大体学习过程是这样的：
先花了10分钟看了一下范浩强的博客
http://fanhq666.blog.163.com/blog/static/8194342620120304463580
然后感觉一脸茫然……就去uoj找模板，找到了那个最短的代码开始看，然后1.2k的代码盯着看了一下午……对你没有看错，就是盯着看了一下午…………
感觉整个人都成带花树了……
下面进入正题——————————————
带花树是解决一般图的最大匹配的算法
具体就是找增广路，然后要对环进行处理，偶数环直接无视，奇数环就缩环，因为奇数环中可以看成所有的点都是等价的，只要有一个点连出边去，那么，整个环就可以匹配！
下面是uoj模板 #79. 一般图最大匹配

#include<bits/stdc++.h>
#define T 505
using namespace std;
int sc()
{
    int i=0;char c=getchar();
    while(c<'0'||c>'9')c=getchar();
    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
    return i;
}
int head[T],lst[T*T],nxt[T*T];
int tot,n,m;
void insert(int x,int y)
{
    lst[++tot]=y;
    nxt[tot]=head[x];
    head[x]=tot;
}
struct Work
{
    int ma[T],st[T],pr[T],fa[T],q[T],v[T];
    int ans,TI,u,t;
    int lca(int x,int y)
    {
        for(TI++;;swap(x,y))if(x)
        {
            if(v[x]==TI)return x;
            v[x]=TI;
            x=fa[pr[ma[x]]];
        }
    }
    void up(int x,int y,int f)
    {
        while(fa[x]!=f)
        {
            pr[x]=y;
            if(st[ma[x]]>0)st[q[++t]=ma[x]]=0;
            if(fa[x]==x)fa[x]=f;
            if(fa[ma[x]]==ma[x])fa[ma[x]]=f;
            x=pr[y=ma[x]];
        }
    }
    int match(int x)
    {
        for(int i=1;i<=n;i++)fa[i]=i,st[i]=-1;
        st[q[t=1]=x]=0;
        for(int l=1;l<=t;l++)for(int i=head[q[l]];i;i=nxt[i])if(st[lst[i]]<0)
        {
            st[lst[i]]=1;
            pr[lst[i]]=q[l];
            if(!ma[lst[i]])
            {
                for(int j=q[l],k=lst[i];j;j=pr[k=u])
                {
                    u=ma[j];
                    ma[j]=k;
                    ma[k]=j;
                }
                return 1;
            }
            st[q[++t]=ma[lst[i]]]=0;
        }
        else if(fa[lst[i]]!=fa[q[l]]&&!st[lst[i]])
        {
            int f=lca(lst[i],q[l]);
            up(q[l],lst[i],f);
            up(lst[i],q[l],f);
            for(int j=1;j<=n;j++)fa[j]=fa[fa[j]];
        }
        return 0;
    }
    void solve()
    {
        for(int i=1;i<=n;i++)ans+=!ma[i]&&match(i);
        cout<<ans<<endl;
        for(int i=1;i<=n;i++)printf("%d ",ma[i]);
    }
}mp;
int main()
{
    n=sc(),m=sc();
    for(int i=1;i<=m;i++)
    {
        int x=sc(),y=sc();
        insert(x,y);
        insert(y,x);
    }
    mp.solve();
    return 0;
}