bzoj4001概率论 卡特兰数+生成函数+导数

#include <iostream>
#include <cstdio>
using namespace std;

int main(){
    double n;
    cin >> n;
    printf("%.9lf\n", n*(n+1)/2/(2*n-1));
    return 0;
}

题解:
首先节点为$n$的二叉树个数, 卡特兰数$C_i={\sum }_{j=0}^{i-1}{C}_j\ast C_{i-j-1}$, 其个数便是$C_n$.
令其生成函数为F(x)，由其递推式可以列出方程：$F(x)=x{F}^2(x)+1$，

$F(x)=\frac{1-\sqrt{1-4x}}{2x}$.
接下来算期望,假设节点为$i$个的二叉树期望是$p_i$,那么对于节点数为$n$,左子树节点数为$l$,右子树节点数为$r$的树,总期望 * 总形态数 = $\sum$(左形态数 * 左期望 * 右形态 + 右形态数 *右期望 * 左形态数)
即为:
$p_i\ast c_i={\sum }_{j=0}^{i-1}(c_j\ast c_{i-j-1}\ast (p_j+p_{i-j-1}))$
我们设$t_i=p_i\ast c_i$
那么有$t_i={\sum }_{j=0}^{i-1}(c_i\ast t_{i-j-1}+t_i\ast c_{i-j-1})=2{\sum }_{j=0}^{i-1}{c}_i\ast t_{i-j-1}$

回想求卡特兰数列生成函数的递推方程的思维方式,由于卡特兰数列是多个项乘积的和,我们可以令生成函数相乘来相消:
$F^2(x)=(c_0+c_{1}x+c_2{x}^2+......)\ast (c_0+c_{1}x+c_2{x}^2+......)$
$=c_0{c}_0+(c_0{c}_1+c_1{c}_0)x+(c_0{c}_2+c_1{c}_1+c_2{c}_0)x^2+.......$
$=c_1+c_{2}x+c_3{x}^2+........$
$=\frac{F(x)-c_0}{x}=\frac{F(x)-1}{x}$
于是$F(x)=x{F}^2(x)+1$
我们类比,这次不是两个$c$相乘而是$c\ast t$, 设$t_i$的生成函数是$G(x)$:
$F(x)\ast G(x)=(c_0+c_{1}x+c_2{x}^2+....)\ast (t_0+t_{1}x+t_2{x}^2+....)$
$=c_0{t}_0+(c_0{t}_1+c_1{t}_0)x+(c_0{t}_2+c_1{t}_1+c_2{t}_0)x^2+......$

$=\frac{t_1+t_{2}x+t_3{x}^2+....}{2}$

$=......$
于是$G(x)=2xF(x)G(x)+x$

$G(x)=\frac{x}{\sqrt{1-4x}}$

接下来开始玄学,注意到:

$(xF(x){)}^{\prime }=\frac{1}{\sqrt{1-4x}}=\frac{G(x)}{x}$

对两者之间有$c_i{x}^{i+1}$求导后$c_i(i+1)x^i$,$\frac{t_{i+1}{x}^{i+1}}{x}=t_{i+1}{x}^i$

说明$c_i(i+1)=t_{i+1}=c_{i+1}{p}_{i+1}$,所以$p_i=\frac{i\ast c_{i-1}}{c_i}$

由卡特兰数公式$\frac{c_i}{c_{i-1}}=\frac{4i-2}{i+1}$

得到$p_i=\frac{i\ast (i+1)}{4\ast i-2}$