带花树算法（一般图最大匹配）

我们首先介绍匈牙利算法，它可以处理不存在奇环的二分图图最大匹配问题，但是，当它处理含奇环的一般图时，要么无法保证复杂度，要么无法保证正确性。

因此，我们有了带花树算法

它的过程和匈牙利算法类似，只是加入了处理奇环（这里称为花）的情况。

首先，依次枚举每个点，找出未匹配的点s，将s染成蓝色，加入队列（使用宽搜），从s出发寻找增广路。

对于队列里的x，枚举它的邻居y，有以下几种情况

1.若y未匹配，则找到增广路

2.若y已匹配：

1）y未染色：将y染为红色，将y的匹配点染为蓝色，并加入队列

2）y染红色：说明遇到偶环，无需处理

3）y染蓝色：说明遇到奇环

 此时，需要用并查集进行缩点操作，将花上的点加入队列（此时花上所有点都可以进行扩展），并全部染成蓝色。

关于缩点的正确性，我们可以考虑奇环上的任意一个点，它与花柄共同将花分成一个奇链和偶链，而那个偶链必然可以作为增广路的一部分（可以自己画图看看），于是，我们不需要考虑花内点的走向，便可以将其缩为一个点。

然而，当我们找到一个增广路时，需要构造路径，于是我们用一个pre数组记录路径信息。

在奇环外，pre即为上一个到达它的点：（绿色即为pre）

 在奇环内，由于不知道增广路的走向，我们需要建立双向pre（可以自己举几个例子体会）：

可能有人会疑问，花柄的pre应该指向谁？对于上面那幅图，花柄pre的指向无所谓，因为到它时，必然会走匹配边。但是，对于下面这种花套花的情况，花柄的pre在之前的花中已经处理好，也不需要管它。

 

 4）若x，y已经缩成一点，则不必处理

自此，主要过程已经结束，可以看下具体代码：

寻找花柄（即搜索树上的lca）

x与y轮流向上跳，第一个重叠处即为lca，其中getfa操作目的是省去跳花内的边，保证每个边只会跳一次，保证复杂度

int lca(int x, int y) {
	++cnt; x=getfa(x); y=getfa(y); //getfa即找当前点属于的花的花柄
	while(v[x]!=cnt) {
		v[x]=cnt; 
		x=getfa(pre[match[x]]);
		if(y) swap(x, y);//y还没跳到根
	}
	return x;
}

建立双向pre

void modify(int x, int lc) {
	while(x!=lc) {
		int mx=match[x], p=pre[mx];
		q.push(mx); col[mx]=1;
		fa[x]=fa[mx]=lc;
		if(getfa(p)==lc) return;//不管花柄
		pre[p]=mx; //另一方向在bfs时已建立
		x=p;
	}
}

int lc=lca(x, y);
if(getfa(x)!=lc)/*不管花柄*/ pre[x]=y; if(getfa(y)!=lc)/*不管花柄*/ pre[y]=x;
modify(x, lc); modify(y, lc);

主过程bfs

int sol(int s) {
	while(!q.empty()) q.pop();
	cnt=0;
	for(int i=1; i<=n; ++i) fa[i]=i, pre[i]=col[i]=v[i]=vis[i]=0;
	q.push(s); col[s]=1; vis[s]=1;
	while(!q.empty()) {
		int x=q.front(); q.pop(); 
		for(int i=la[x]; i; i=g[i].nxt) {
			int y=g[i].y;
			if(getfa(x)==getfa(y)) continue;
			if(col[y]==0) {
				col[y]=2; pre[y]=x;
				if(match[y]==0) {//找到增广路
					while(x!=s) {
						int z=match[x];
						match[x]=y; match[y]=x;
						x=pre[z]; y=z;
					}
					match[x]=y; match[y]=x;
					return 1;
				} else if(vis[match[y]]==0) {col[match[y]]=1; vis[match[y]]=1; q.push(match[y]);}
			} else if(col[y]==1) {
				int lc=lca(x, y);
				if(getfa(x)!=lc) pre[x]=y; if(getfa(y)!=lc) pre[y]=x;
				modify(x, lc); modify(y, lc);
			}
		}
	}
	return 0;
}