zoukankan      html  css  js  c++  java
  • 【刷题】UOJ #79 一般图最大匹配

    从前一个和谐的班级,所有人都是搞OI的。有 (n) 个是男生,有 (0) 个是女生。男生编号分别为 (1,…,n)

    现在老师想把他们分成若干个两人小组写动态仙人掌,一个人负责搬砖另一个人负责吐槽。每个人至多属于一个小组。

    有若干个这样的条件:第 (v) 个男生和第 (u) 个男生愿意组成小组。

    请问这个班级里最多产生多少个小组?

    输入格式

    第一行两个正整数,(n,m) 。保证 (n≥2)

    接下来 (m) 行,每行两个整数 (v,u) 表示第 (v) 个男生和第 (u) 个男生愿意组成小组。保证 (1≤v,u≤n) ,保证 (v≠u) ,保证同一个条件不会出现两次。

    输出格式

    第一行一个整数,表示最多产生多少个小组。

    接下来一行 (n) 个整数,描述一组最优方案。第 (v) 个整数表示 (v) 号男生所在小组的另一个男生的编号。如果 (v) 号男生没有小组请输出 (0)

    样例一

    input

    10 20
    9 2
    7 6
    10 8
    3 9
    1 10
    7 1
    10 9
    8 6
    8 2
    8 1
    3 1
    7 5
    4 7
    5 9
    7 8
    10 4
    9 1
    4 8
    6 3
    2 5
    

    output

    5
    9 5 6 10 2 3 8 7 1 4 
    

    样例二

    input

    5 4
    1 5
    4 2
    2 1
    4 3
    

    output

    2
    2 1 4 3 0 
    

    限制与约定

    (1≤n≤500)(1≤m≤124750)

    时间限制(1s)

    空间限制(256MB)

    题解

    带花树模板题,现在就是背代码。。。

    与二分图一样找增广路,如果找到奇交替链,那么返回 (1)

    如果期望匹配的点已经被匹配了,那么像匈牙利那样,看这个点匹配的那个点能不能再找一个,直接加入队列搜就好了

    如果找到奇环,那么就缩花了,暴力求 (lca) ,然后改 (fa) 数组就行

    #include<bits/stdc++.h>
    #define ui unsigned int
    #define ll long long
    #define db double
    #define ld long double
    #define ull unsigned long long
    const int MAXN=500+10,MAXM=130000+10,inf=0x3f3f3f3f;
    int n,m,e,fa[MAXN],beg[MAXN],nex[MAXM<<1],to[MAXM<<1],vis[MAXN],level[MAXN],clk,link[MAXN],pre[MAXN];
    std::queue<int> q;
    template<typename T> inline void read(T &x)
    {
    	T data=0,w=1;
    	char ch=0;
    	while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
    	if(ch=='-')w=-1,ch=getchar();
    	while(ch>='0'&&ch<='9')data=((T)data<<3)+((T)data<<1)+(ch^'0'),ch=getchar();
    	x=data*w;
    }
    template<typename T> inline void write(T x,char ch='')
    {
    	if(x<0)putchar('-'),x=-x;
    	if(x>9)write(x/10);
    	putchar(x%10+'0');
    	if(ch!='')putchar(ch);
    }
    template<typename T> inline void chkmin(T &x,T y){x=(y<x?y:x);}
    template<typename T> inline void chkmax(T &x,T y){x=(y>x?y:x);}
    template<typename T> inline T min(T x,T y){return x<y?x:y;}
    template<typename T> inline T max(T x,T y){return x>y?x:y;}
    inline int found(int x)
    {
    	if(fa[x]!=x)fa[x]=found(fa[x]);
    	return fa[x];
    }
    inline void insert(int x,int y)
    {
    	to[++e]=y;
    	nex[e]=beg[x];
    	beg[x]=e;
    }
    inline int LCA(int u,int v)
    {
    	for(++clk;;std::swap(u,v))
    		if(u)
    		{
    			u=found(u);
    			if(level[u]==clk)return u;
    			level[u]=clk,u=pre[link[u]];
    		}
    }
    inline void blossom(int u,int v,int lca)
    {
    	while(found(u)!=lca)
    	{
    		pre[u]=v,v=link[u];
    		if(vis[v]==2)vis[v]=1,q.push(v);
    		if(found(u)==u)fa[u]=lca;
    		if(found(v)==v)fa[v]=lca;
    		u=pre[v];
    	}
    }
    inline int bfs(int s)
    {
    	for(register int i=1;i<=n;++i)fa[i]=i;
    	memset(vis,0,sizeof(vis));
    	memset(pre,0,sizeof(pre));
    	while(!q.empty())q.pop();
    	vis[s]=1;
    	q.push(s);
    	while(!q.empty())
    	{
    		int x=q.front();
    		q.pop();
    		for(register int i=beg[x];i;i=nex[i])
    		{
    			if(found(x)==found(to[i])||vis[to[i]]==2)continue;
    			if(!vis[to[i]])
    			{
    				vis[to[i]]=2;
    				pre[to[i]]=x;
    				if(!link[to[i]])
    				{
    					for(register int p=to[i],las;p;p=las)las=link[pre[p]],link[p]=pre[p],link[pre[p]]=p;
    					return 1;
    				}
    				else vis[link[to[i]]]=1,q.push(link[to[i]]);
    			}
    			else
    			{
    				int lca=LCA(x,to[i]);
    				blossom(x,to[i],lca);blossom(to[i],x,lca);
    			}
    		}
    	}
    	return 0;
    }
    int main()
    {
    	static int ans=0;
    	read(n);read(m);
    	for(register int i=1;i<=m;++i)
    	{
    		int u,v;read(u);read(v);
    		insert(u,v);insert(v,u);
    	}
    	for(register int i=1;i<=n;++i)
    		if(!link[i])ans+=bfs(i);
    	write(ans,'
    ');
    	for(register int i=1;i<=n;++i)write(link[i],' ');
    	puts("");
    	return 0;
    }
    
  • 相关阅读:
    tornado-cookies+pycket 验证
    解决form表单通过ajax时,required失效问题
    如何判断一个文件是否存在
    django中@property装饰器的运用
    RestFul API接口设计风格介绍和核心功能(概念)
    python中各个response使用
    Python的命名规则
    break和continue的区别
    LOCK接口
    API
  • 原文地址:https://www.cnblogs.com/hongyj/p/9120185.html
Copyright © 2011-2022 走看看