[校内模拟][树形DP] 201031 NOIP T1 BST

题意

已知一个二叉搜索树，求依次插入而形成这颗树的序列数量

正解

建树，手玩发现答案与子树大小有关，可以DP
记答案（以$u$为根的子树的序列数）为$f_u$
手玩发现答案与子树大小有关

考虑转移
发现两子树合并其实就是子树的序列合并，两棵子树的序列不能变
子树的序列有$f_s$种，子树的序列合起来要填到$si{z}_u-1$个空位上（$-1$是因为根节点一定在第一个）
考虑左子树，应该选$si{z}_s$个位置放，组合数
所以转移：
$$f_u=f_{lch}{f}_{rch}\left(\genfrac{}{}{0ex}{}{si{z}_{lch}+si{z}_{rch}}{si{z}_{lch}}\right)$$

#define FILE(x) freopen(x".in","r",stdin);freopen(x".out","w",stdout);
#include<bits/stdc++.h>
using namespace std;
#define in Read()
int in{
	int i=0,f=1;char ch=0;
	while(!isdigit(ch)&&ch!='-') ch=getchar();
	if(ch=='-') ch=getchar();
	while(isdigit(ch)) i=(i<<1)+(i<<3)+ch-48,ch=getchar();
	return i*f;
}

const int N=1e3+5,mod=1e9+7;
int n,lch[N],rch[N];
int siz[N],f[N],root;
int fac[N],ifac[N];

int add(int a,int b){return a+b>=mod?a+b-mod:a+b;}
int dec(int a,int b){return a<b?a+mod-b:a-b;}
int mul(int a,int b){return 1ll*a*b%mod;}
int qpw(int a,int b){
	int res=1;
	while(b){
		if(b&1) res=mul(res,a);
		a=mul(a,a);
		b>>=1;
	}
	return res;
}
int C(int a,int b){return mul(fac[a],mul(ifac[b],ifac[a-b]));}

void print_tree(int rt){
	printf("%d %d %d\n",rt,lch[rt],rch[rt]);
	if(lch[rt]) print_tree(lch[rt]);
	if(rch[rt]) print_tree(rch[rt]);
	return;
}

void insert(int x,int rt){
	if(!root){
		root=x;
		return;
	}
	if(x>rt){
		if(!rch[rt]){
			rch[rt]=x;
			return;
		}else insert(x,rch[rt]);
	}else{
		if(!lch[rt]){
			lch[rt]=x;
			return;
		}else insert(x,lch[rt]);
	}
	return;
}

void DFS(int u){
	f[u]=1;
	siz[u]=1;
	if(lch[u]) DFS(lch[u]);else f[lch[u]]=1;
	if(rch[u]) DFS(rch[u]);else f[rch[u]]=1;
	
	siz[u]+=siz[lch[u]]+siz[rch[u]];
	f[u]=mul(mul(f[lch[u]],f[rch[u]]),C(siz[u]-1,siz[lch[u]]));
	return;
}

void solve(){
	memset(lch,0,sizeof lch);
	memset(rch,0,sizeof rch);
	memset(siz,0,sizeof siz);
	memset(f,0,sizeof f);
	root=0;
	
	n=in;
	for(int i=1;i<=n;++i) insert(in,root);
	
	DFS(root);
	printf("%d\n",f[root]-1);
	return;
}

int main(){
	FILE("bst");
	
	fac[0]=1;
	for(int i=1;i<=N-3;++i) fac[i]=mul(fac[i-1],i);
	ifac[N-3]=qpw(fac[N-3],mod-2);
	for(int i=N-4;i>=0;--i) ifac[i]=mul(ifac[i+1],i+1);
	
	int T=in;
	while(T--) solve();
	return 0;
}