prufer序列学习小记

简介：

prufer序列是用来解决一系列有关有标号无根树计数问题的利器。
首先介绍一下相关的东西：
prufer序列的构造：
1.找到编号最小的度数为1的点。
2.删除该节点并在序列中添加与该节点相连的节点的编号。
3.重复$1,2$操作，直到整棵树只剩下两个节点。

将prufer序列对应回原树：
1.每次取出prufer序列中最前面的元素$u$。
2.在点集中找到编号最小的没有在prufer序列中出现的元素$v$。
3.给$u,v$连边然后分别删除。
4.最后在点集中剩下两个节点，给它们连边。
实际上也是把操作逆回去。其中第二个操作可以用数据结构来维护，复杂度可以做到$O(n\log n)$。

prufer序列一些性质：
prufer序列与无根树唯一对应。
点数为$n$的有标号无根树转化来的prufer数列长度为$n-2$，且每个数都$\in [1,n]$。
通过构造过程可以发现，一个在树中度数为$d_i$的点，在序列中出现了$d_i-1$次。

运用：

有如上性质，我们考虑一个$n$个点的，且每个点的度数确定的有标号无根树有多少个？
设每个点度数为$d_i$，则有$\frac{(n-2)!}{\prod (d_i-1)!}$这么多。
考虑每个数出现的次数然后用组合数相乘，可以化成如上形式。
运用到题目中是类似的。

有$n$个点，它们从$1$到$n$进行标号，第i个点的限制为度数不能超过$A_i$.
现在对于每个$s(1\le s\le n)$，问从这$n$个点中选出一些点组成大小为$s$的有标号无根树的方案数。

$n\le 100$
这是一道裸题。
考虑对于一个确定的$s$，用prufer序列计数，发现式子上面的$(n-2)!$已经确定，下面的不定，所以我们$dp$求出下面的乘积。
考虑设$f_{i,j,k}$表示做到第$i$个点，取了$j$个点，当前prufer序列长度为$k$。
转移只要考虑当前这一个选不选，度数为多少即可，非常简单。
复杂度$O(n^4)$

# include<cstdio>
# include<cstring>
# include<algorithm>
using namespace std;
const int N = 100 + 5;
const int mo = 1004535809;
typedef long long ll;
int fac_[N],fac[N],inv[N];
int f[N][N][N],a[N];
int n;
int main()
{
	freopen("tree.in","r",stdin);
	freopen("tree.out","w",stdout);
	scanf("%d",&n),fac[0] = fac_[0] = inv[1] = 1;
	for (int i = 1 ; i <= n ; ++i)
	{
		scanf("%d",&a[i]); if (i > 1) inv[i] = (ll)(mo - mo / i) * inv[mo % i] % mo;
		fac[i] = (ll)fac[i - 1] * i % mo,fac_[i] = (ll)fac_[i - 1] * inv[i] % mo;
	}
	f[0][0][0] = 1;
	for (int i = 0 ; i < n ; ++i)
		for (int j = 0 ; j <= i ; ++j)
			for (int k = 0 ; k <= n - 2 ; ++k)
			if (f[i][j][k])
			{
				f[i + 1][j][k] = (f[i + 1][j][k] + f[i][j][k]) % mo;
				for (int l = 0 ; l + k <= n - 2 && l < a[i + 1] ; ++l)
					f[i + 1][j + 1][k + l] = (f[i + 1][j + 1][k + l] + (ll)f[i][j][k] * fac_[l] % mo) % mo;
			}
	printf("%d ",n);
	for (int i = 2 ; i <= n ; ++i) printf("%d ",(ll)f[n][i][i - 2] * fac[i - 2] % mo);
	return 0;
}