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

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;
}