从Prufer编码到各种树的计数

从$Prufer$编码到各种树的计数

1. $Prufer$编码

1. 将无根树进行编码的一种方式，无根树就是边没有方向的树，即以树上的每一个点做根，只算一种方案。

2. 无根树转序列做法

   ​ 每次找到编号最小的叶子，输出它父亲的编号，然后删除，直到只剩2个节点。

   3. 序列转无根树做法

   ​ 设集合 S = { 1 , 2 , ⋯   , n } S=\left \{ 1,2,\cdots,n \right \} S={1,2,⋯,n}

   ​ 第$i$次，找到在$S$中没有在$prufer$编码第$i$到$n-2$位出现的最小的数，将它与$prufer$编码的第$i$位连边，并把它从S中删除。

   ​ 最后剩下的两个数，它们两个连一条边即可。

   ​ 代码实现的时候，可以将每个点的度数初始化为$prufer$序列中出现的次数+1，每次找度数为1最小的那一个与$prufer$序列当前位连边。

2. 有标号无根树的计数($Cayley$公式)

1. 由上面的做法显然可以得出，每一个长度为n-2的序列，唯一的对应了一棵无根树，且一颗无根树唯一的对应了一个序列。

2. 因此可以得到n个点的有标号无根树总共有$n^{n-2}$个。

3. 性质与应用

   1. 在$prufer$编码中，每个点出现的次数为$该点度数-1$。

   2. 如果限制n个点的读数为$d_1,d_2,d_3,\cdots ,d_n$，那么这样的无根树共有
      $$\frac{(n-2)!}{(d_1-1)!(d_2-1)!(d_3-1)!\cdots (d_n-1)!}个$$

4. 如果只限制了k个呢，$prufer$编码剩了m个位置呢？那么假装限制了，再乘以随便填的方案数,即
   $$\frac{(n-2)!}{(d_1-1)!(d_2-1)!\cdots (d_k-1)!(m)!}\ast (n-k{)}^m$$

3. 有标号有根树的计数

1. $prufer$编码是无根树，但是树形均不同，那么每个生成的树以每个点做根均可，因此共有
   $$n^{n-2}\ast n=n^{n-1}个$$

4.无标号有根树的计数

1. 方法1朴素形式，设$a_n$为树大小为n时的答案。设对于根，大小为k的子树共有$c_k$个，而每个子树都可以选择$1\sim a_k$的一种方案，即${\sum }_{k=1}^{a_k}{x}_k=c_k$,该方程非负整数解个数为${\prod }_{k>0}{C}_{a_k+c_k-1}^{c_k}$($a_1=1$)。
   $$a_n=\sum _{c_1+2c_2+\cdots +n{c}_n=n-1}\prod _{k>0}{C}_{a_k+c_k-1}^{c_k}$$

2. 方法1高级形式经过复杂的母函数推导，我们得到
   $$\begin{gathered} 令s_{n,j}=\sum _{1\le i\le \frac{n}{j}}{a}_{n+1-i\ast j},则有s_{n,j}=s_{n-j,j}+a_{n+1-j} \\ 则有a_{n+1}=\frac{{\sum }_{1\le j\le n}j{a}_j{s}_{n,j}}{n} \end{gathered}$$

3. 方法2朴素形式如果我们把定义改为$f_{n,j}$表示树大小为n，根节点度数为j的方案数，那么$a_n={\sum }_{j=1}^n{f}_{n,j}$,就有
   $$f_{i,j}=\sum _{c_1+2c_2+\cdots +j{c}_j=i-1}\prod _{k>0}^j{C}_{a_k+c_k-1}^{c_k}$$

4. 方法2高级形式我们考虑对$f_{i,j}$进行$dp$，可以通过限制最大子树进行背包，即从小到大枚举$maxson$，则有
   f i , j = ∑ k = 1 m a x { j , ⌊ i m a x s o n ⌋ } C a m a x s o n + k − 1 k f i − k ∗ m a x s o n , j − k f_{i,j}=\sum_{k=1}^{max\left\{j,\lfloor\frac i {maxson\rfloor}\right\}}C_{a_{maxson}+k-1}^k f_{i-k*maxson,j-k} fi,j​=k=1∑max{j,⌊maxson⌋i​}​Camaxson​+k−1k​fi−k∗maxson,j−k​

5. 这样有一个好处就是可以限制最大度数，如构造烷烃。

5.无标号无根树的计数

1. 利用一下上题的$a_n$

2. 方法1，设$b_n$为本题的答案，根据质心的关系，可以得到
   $$\begin{gathered} n为奇数时 \\ {b}_n=a_n-\sum _{0<i\le \lfloor \frac{n}{2}\rfloor }{a}_i{a}_{n-i} \\ n为偶数时 \\ {b}_n=a_n-\sum _{0<i\le n}{a}_i{a}_{n-i}+\frac{1}{2}{a}_{\frac{n}{2}}(a_{\frac{n}{2}}+1) \end{gathered}$$

3. 方法2,可以通过枚举重心的方法来解决，强制$maxson<\lceil \frac{n}{2}\rceil$，则保证是是以重心为根的无标号有根树计数，但如果n为偶数，可能会有两个重心，这两个重心子树大小都是$\frac{n}{2}$，因此最终结果为
   $$ans=\sum _{j=0}^m{f}_{n,j}+[nmod2=0]C_{a_{\frac{n}{2}+1}}^2$$

4. 烷烃例题 [BZOJ4271]chemistry，方法2(没有加高精度)代码如下

   //无高精度版本 
   #include<iostream>
   #include<stdio.h>
   #include<string.h>
   #include<algorithm>
   #define maxn 505
   using namespace std;
   typedef long long ll;
   int n,m;
   ll dp[maxn][maxn];
   ll c(ll n,ll m)
   {
   	ll ans=1;
   	for(int i=m;i>=1;i--)
   	{
   		ans=ans*n;
   		n--;
   	}
   	for(int i=m;i>1;i--)
   	{
   		ans/=i;
   	}
   	return ans;
   }
   int main()
   {
   	dp[1][0]=1;
   	scanf("%d",&n); m=4;
   	for(int maxson=1;maxson<((n+1)/2);maxson++)
   	{
   		ll a=0;
   		for(int j=0;j<m;j++) 
   			a+=dp[maxson][j];
   		for(int i=n;i>maxson;i--)
   		{
   			for(int j=1;j<=m;j++)
   			{
   				for(int k=1;k<=j&&maxson*k<i;k++)
   				{
   					dp[i][j]+=dp[i-maxson*k][j-k]*c(a+k-1,k);
   				}
   			}
   		}
   	}
   	ll ans=0;
   	for(int j=0;j<=m;j++)
   		ans+=dp[n][j];
   	if(!(n&1))
   	{
   		ll a=0;
   		for(int j=0;j<m;j++)
   			a+=dp[n/2][j];
   		ans+=c(a+1,2);
   	}
   	printf("%lld\n",ans);
   }
   

参考资料

1. 顾誉州WC交流https://wenku.baidu.com/view/3fbabbcc647d27284a735113.html
2. 华师一附中赵爽论文 https://max.book118.com/html/2018/0403/159821895.shtm
3. Flashhu https://www.cnblogs.com/flashhu/p/9457830.html