[复习]prufer 序

对于一棵 n n n 个点无根树,其 prufer 序长度为 n − 2 n-2 n2,具体来说,每次选择一个编号最小的叶子节点并且删除它,记录下它的父亲,显然 n − 2 n-2 n2 次操作后还剩下 2 2 2 个点。

首先显然任何一个无根树都能转化为 prufer 序。

然后需要说明任意 prufer 序都能唯一对应一棵无根树,度数为 x x x 的点在 prufer 序中出现 x − 1 x-1 x1 次,那么考虑按照度数来构造,从左到右扫一遍,对于每个点和当前最小的叶子连边即可。

这里再说一下构造的复杂度,用堆 O ( n log ⁡ n ) O(n\log n) O(nlogn) 比较显然,事实上,通过指针扫一遍复杂度可以变成 O ( n ) O(n) O(n),这个不难实现。

由于两者是一一对应的,我们可以得到: n n n 个点的无根树的数量为 n n − 2 n^{n-2} nn2 次方。

拓展结论:

1.求完全二分图的生成树的个数(文艺计算机)

考虑对这个二分图生成树构造 prufer 序,最后一定剩下一个黑点一个白点。上面两个性质同理可证。所以答案就是 n m − 1 m n − 1 n^{m-1}m^{n-1} nm1mn1

2.求有多少棵以 1 1 1 为根的有根树满足位于奇数层的点数恰好为 m m m(雅礼集训)

考虑先钦定哪些点在奇数,方案数是 C n − 1 m − 1 C^{m-1}_{n-1} Cn1m1,然后点固定下来之后,原问题还是可以看成一个二分图生成树问题,故套用文艺计算机的结论即可。

由于机智的发文助手认为我太短了,那就只好支配树镇楼

#include<bits/stdc++.h>
using namespace std;
const int maxn=2e6+5;
int n,m;
struct Edge{
	int v,nxt;
}e[maxn];
int cnt,h[maxn];
void add(int u,int v){
	e[++cnt].v=v;
	e[cnt].nxt=h[u];
	h[u]=cnt;
}
int dfn[maxn],idom[maxn],sdom[maxn],dt,bl[maxn],fath[maxn],mn[maxn];
void dfs(int u,int fa){
	dfn[u]=++dt;
	bl[dt]=u;
	fath[u]=fa;
	for(int i=h[u];i;i=e[i].nxt){
		int v=e[i].v;
		if(!(i&1))continue;
		if(!dfn[v]){
			dfs(v,u);
		}
	}
}
int fa[maxn];
int gf(int u){
	if(u==fa[u])return u;
	int f=gf(fa[u]);
	if(dfn[sdom[mn[u]]]>dfn[sdom[mn[fa[u]]]])mn[u]=mn[fa[u]];
	return fa[u]=f; 
}
vector<int>G[maxn];
void solve(int S){
	dfs(S,0);
	for(int i=1;i<=n;i++)sdom[i]=fa[i]=mn[i]=i;
	for(int i=n;i>=2;i--){
		int u=bl[i];
		for(int i=h[u];i;i=e[i].nxt){
			int v=e[i].v;
			if(i&1)continue;
			gf(v);
			if(dfn[sdom[u]]>dfn[sdom[mn[v]]])sdom[u]=sdom[mn[v]];
		}
		fa[u]=fath[u];
		G[sdom[u]].push_back(u);
		u=fath[u];
		for(int i=0;i<(int)G[u].size();i++){
			int v=G[u][i];
			gf(v);
			if(u==sdom[mn[v]])idom[v]=u;
			else idom[v]=mn[v];
		}              
		G[u].clear();
	}	
	for(int i=2;i<=n;++i){ 
		int u=bl[i];
		if(idom[u]^sdom[u]){
			idom[u]=idom[idom[u]];
		} 
	}
} 
int ans[maxn];
int main(){
	scanf("%d%d",&n,&m);
	for(int i=1,u,v;i<=m;i++){
		scanf("%d%d",&u,&v);
		add(u,v);
		add(v,u);
	}
	solve(1);
	for(int i=n;i>=2;--i)ans[idom[bl[i]]]+=(++ans[bl[i]]);
	++ans[1];
	for(int i=1;i<=n;++i)printf("%d ",ans[i]);
	return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值