【prufer序】树-dp-组合数学

今天本蒟蒻学了这么多东西，当然要写几篇水博文。

题意：

有n个点，第i个点的限制为度数不能超过ai。现在对于每一个s(1≤s≤n)(1≤s≤n)，问从这n个点中选出s个点组成有标号无根树的方案数mod1004535809

(n<=100)

 

题解：

先贴几篇博文：

当然要先看看morejarphone大佬的prufer讲解（贼好懂）：点击打开链接

还有这篇棒棒的博文：点击打开链接

关于这篇博文里一个地方，我有说明一下的必要：

 

2）一个有趣的推广是，n个节点的度依次为D1, D2, …, Dn的无根树共有   (n-2)! / [ (D1-1)!(D2-1)!..(Dn-1)! ]  个，因为此时Prüfer编码中的数字i恰好出现Di-1次。

即 n种元素，共n-2个，其中第i种元素有Di-1个，求排列数。

搜了一下，发现网上并没有对于这句话的解释，也可能是本蒟蒻太菜了，对组合数学一无所知，所以不能理解这个有趣的推广。

仔细想了后得出：

不考虑每个节点具体要求时填满n-2个位置是(n-2)!，设第i个树度数为Di,那么其中有Di-1个数均为i,而在最初的这任意Di-1位上，是任意取值的，所以有（Di-1）！种取法，现在变为了一种，所以此时总排列数是(n-2)!/（Di-1)!,以此类推，得出上式：

于是当序列中共有s个元素时,第i个数出现次数位Ci,得出下式：

S!/(Ci!*Ci+1!*...*Cn!)

但我们的输出是对于s=1 to n的，所以不能直接算，那么自然而然想到了dp。

设dp[i][j][k]表示前i个数，选择了j个数，序列中有k个元素时的种类。

但这道题妙的是求有标号无根树，所以还要在dp[n][s][s-2]基础上乘以s!

代码如下：

#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int mod=1004535809; 
typedef long long ll;
ll dp[105][105][105];
ll g[105],h[105];
int a[105],n;
inline int read()
{
    int X=0,w=0; char ch=0;
    while(ch<'0' || ch>'9') w|=ch=='-',ch=getchar();
    while(ch>='0' && ch<='9') X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
    return w?-X:X;
}

inline int min(int a,int b){
	return a>b? b:a;
}

inline ll fp(ll a,ll b)
{
	ll ret=1;
	while(b){
	  if(b&1) ret=ret*a%mod;
	  a=a*a%mod;
	  b>>=1; 
	} 
	return ret;
}

int main(){
	n=read();g[0]=1;dp[0][0][0]=1;h[0]=1;
	for(int i=1;i<=(int)n;i++) {
		a[i]=read();g[i]=(g[i-1]*i)%mod;
	}
	h[n]=fp(g[n],mod-2);
	for(int i=n-1;i>=1;i--) h[i]=h[i+1]*(i+1)%mod;
	for(int i=0;i<n;i++){
	  for(int j=0;j<=i;j++){
	  	for(int k=0;k<=n-2;k++){
		  if(dp[i][j][k]){
		  	dp[i+1][j][k]=(dp[i+1][j][k]+dp[i][j][k])%mod;
			int temp=min(n-2-k,a[i+1]-1); 
		  	for(int l=0;l<=temp;l++){
		  		dp[i+1][j+1][k+l]=(dp[i+1][j+1][k+l]+dp[i][j][k]*h[l]%mod)%mod;
		  	}
		  } 
	    }
      }
     }
     printf("%lld ",(ll)n);
	 for(int i=2;i<=n;i++) printf("%lld ",dp[n][i][i-2]*g[i-2]%mod);
	  return 0;
}