带花树算法

算法思想:

①找增广路

②过奇环中任意一点有增广路,则过奇环亦有增广路,所以可以将奇环缩点

算法流程:

以一个当前未匹配上的点为s点,染为白色,bfs染色,只许白点进队。

当两个白色点有边时,说明遇到奇环,缩点。

当遇到一个未染色且未匹配点时,找到增广路,修改匹配后退出。


代码:

#include<bits/stdc++.h>

#define N 200020

using namespace std;

inline int read()

{

	int x=0,f=1;char ch=getchar();

	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

	while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}

	return x*f;

}

queue<int>Q;

int n,m,head[N],pos;

int tot,match[N],fri[N],top[N],ring[N];

bool odd[N],vis[N];

struct edge{int to,next;}e[N<<1];

inline void add(int a,int b)

{

	pos++;e[pos].to=b,e[pos].next=head[a],head[a]=pos;

}

inline void insert(int u,int v){add(u,v);add(v,u);}

int fa[N];

inline int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}

int clo;

inline int lca(int x,int y)

{

	clo++;

	for(;x;x=find(top[x]))ring[x]=clo;

	for(x=y;ring[x]!=clo;x=find(top[x]));

	return x;

}

inline void blossom(int x,int y,int p)

{

	for(;find(x)!=find(p);x=fri[y])

	{

		fri[x]=y,y=match[x];

		fa[find(x)]=fa[find(y)]=p;

		Q.push(y),odd[y]=0;

	}

}

inline bool dfs(int s)

{

	clo=0;

	for(int i=1;i<=n;i++)odd[i]=vis[i]=ring[i]=top[i]=fri[i]=0,fa[i]=i;

	while(!Q.empty())Q.pop();vis[s]=1,Q.push(s);

	while(!Q.empty())

	{

		int u=Q.front();Q.pop();

		for(int i=head[u];i;i=e[i].next)

		{

			int v=e[i].to;

			if(!vis[v])

			{

				vis[v]=odd[v]=1,fri[v]=top[v]=u;

				if(!match[v])

				{

					for(int x,y,j=v;j;j=y)

					{

						x=fri[j],y=match[x];

						match[j]=x,match[x]=j;

					}return 1;

				}

				vis[match[v]]=1,top[match[v]]=v,Q.push(match[v]);

			}

			else if(!odd[v]&&find(v)!=find(u))

			{

				int p=lca(u,v);

				blossom(u,v,p);

				blossom(v,u,p);

			}

		}

	}return 0;

}

int main()

{

	n=read(),m=read();

	for(int i=1;i<=m;i++)

	{

		int x=read(),y=read();

		insert(x,y);

	}

	int ans=0;

	for(int i=1;i<=n;i++)if(!match[i])ans+=dfs(i);

	printf("%d\n",ans);

	for(int i=1;i<=n;i++)printf("%d ",match[i]);

}


  • 2
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值