目录
部分来自: 大佬的博客,加入了一些自己的理解。
1、 prufer编码
首先引入 prufer编码(其实是 Pru¨fer)
prufer数列,可以用来解一些关于无根树计数的问题。
prufer数列是一种无根树的编码表示,对于一 棵 n 个节点带编号的无根树,对应唯一一串长度为n-1的 prufer编码。
1)无根树转化为prufer序列
一棵无根树的 prufer编码的值运算如下:
首先定义无根树中度数为1的节点是叶子节点。 找到编号最小的叶子并删除,序列中添加与之相连的节点编号,重复执行直到只剩下2个节点。
如图:
它的 prufer编码就是4, 3, 3
先丢1,4加入序列,再丢2,3加入序列,再丢5,3加入序列,剩余两个,停止丢人,序列为4,3,3
显然,一棵有n个结点的无根树,它的 prufer 编码是唯一的,且有 n−2 个可能相同的元素。
具体实现可以用一个 set搞定,维护度数为 1 的节点。复杂度 O(nlogn)。
2)prufer序列转化为无根树。
设点集 V=(1,2,3,...,n),每次取出prufer序列中最前面的元素u,在V中找到编号最小的没有在 prufer序列中出现的元素v,给u,v连边然后分别删除,最后在V中剩下两个节点,给它们连边。最终得到的就是无根树。
具体实现也可以用一个 set,维护 prufer序列中没有出现的编号。复杂度 O(nlogn)。
很显然,每一个prufer序列与一棵无根树一一对应。
因此,对于一棵已知有n个结点的无根树,一定有一个n-2长度的序列,那么,我们枚举所有长度为n-2的序列,发现其与所有可能形态的无根树一一对应。而这种序列,根据乘法原理,有**n n-2**个可能的序列。
因此,对于一个已知的n,有n n-2 种不同的无根树。
最后有一个很重要的性质就是 prufer 序列中某个编号出现的次数就等于这个编号的节点在无根树中的度数-1。
(例如最上面的图中4,3,3,4的度数为2,3的度数为3。)
2、Cayley公式
Cayley公式是说,一个无向完全图有 nn−2 棵生成树,通俗的说就是n个节点的带编号的无根树有 nn−2个。
刚才Prufer有一个很重要的性质:序列与树是一一对应的
而Prufer序列有n-2项,序列中的每个数都在1到n的范围内。
所以我们可以直接推出n个点的无向完全图的生成树的计数: nn−2
一个有趣的推广是,n个节点的度依次为D1, D2, …, Dn的无根树共有 (n−2)!/[(D1−1)!(D2−1)!..(Dn−1)!] 个,因为此时 Pru¨fer编码中的数字i恰好出现 Di−1次(上面的结论)。
1)由Cayley公式得到四个推论
一棵n个节点的无根树唯一地对应了一个长度为n-2的数列,数列中的每个数都在1到n的范围内。
上面这句话比较重要。通过上面的定理,
1) 我们可以直接推出n个点的无向完全图的生成树的计数:n(n-2) 即n个点的有标号无根树的计数。
2) 一个有趣的推广是,n个节点的度依次为 D1,D2,…,Dn的无根树共有 (n−2)!/[(D1−1)!(D2−1)!..(Dn−1)!] 个,因为此时 Pru¨fer编码中的数字i恰好出现Di-1次。
即 n种元素,共 n−2个,其中第 i 种元素有 Di−1个,求排列数。
3) n个节点的度依次为 D1,D2,…,Dn,令有m个节点度数未知,求有多少种生成树?(BZOJ1005 明明的烦恼)
令每个已知度数的节点的度数为di,有n个节点,m个节点未知度数, left=(n−2)−(d1−1)−(d2−1)−...−(dk−1)
已知度数的节点可能的组合方式如下
(n−2)!/(d1−1)!/(d2−1)!/.../(dk−1)!/left!
剩余left个位置由未知度数的节点随意填补,方案数为mleft
于是最后有
ans=(n−2)!/(d1−1)!/(d2−1)!/.../(dk−1)!/left!∗mleft
例题1、P4981 父子
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<queue>
#define ls (p<<1)
#define rs (p<<1|1)
#define mid (l+r)/2
#define over(i,s,t) for(register int i=s;i<=t;++i)
#define lver(i,t,s) for(register int i=t;i>=s;--i)
using namespace std;
typedef long long ll;//全用ll可能会MLE,ll比int占的内存大
const ll N=1002;
const ll INF=1e9+9;
const ll mod=1e9+9;
const double EPS=1e-6;
ll n,m;
inline ll qpow(ll a,ll b)
{
ll res=1;
while(b)
{
if(b&1)res=(res*a)%mod;
a=(a*a)%mod;
b>>=1;
}
return res;
}
int main()
{
ll t;
scanf("%lld",&t);
while(t--)
{
scanf("%lld",&n);
printf("%lld\n",qpow(n,n-1));
}
return 0;
}
例二、P4430 小猴打架
P4430 小猴打架
由 Cayley定理,n个节点的带标号的形态不同的无根树有 nn−2n 个,这个 purfer序对应的树有 (n−1)条边, 如果还要考虑这棵树的生成顺序的话, 那就是 n−1条边全排列: (n−1)!
答案 (n−1)!∗nn−2mod9999991
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<queue>
#define ls (p<<1)
#define rs (p<<1|1)
#define mid (l+r)/2
#define over(i,s,t) for(register int i=s;i<=t;++i)
#define lver(i,t,s) for(register int i=t;i>=s;--i)
using namespace std;
typedef long long ll;//全用ll可能会MLE,ll比int占的内存大
const ll N=1002;
const ll INF=1e9+9;
const ll mod=9999991;
const double EPS=1e-6;
int main()
{
ll n,ans=1;
scanf("%lld",&n);
over(i,1,n-2)ans=(ans*n)%mod;
over(i,1,n-1)ans=(ans*i)%mod;
printf("%lld\n",ans);
return 0;
}
注:如果您通过本文,有(qi)用(guai)的知识增加了,请您点个赞再离开,如果不嫌弃的话,点个关注再走吧 ! 当然,也欢迎在讨论区指出此文的不足处,作者会及时对文章加以修正 !