bzoj4001 [TJOI2015]概率论（生成函数）

我以为我不会忘记的qaq
还是写一写推导吧qaq，方便后人
首先设hn表示n个节点的二叉树的形态数
那么我们有$h_n=\sum\limits_{i=1}^{n-1}h_i*h_{n-1-i},h_0=1,h_1=1$
写出$h_n$的生成函数：$f(x)=1+x+2x^2+5x^3+...$
$f^2(x)=1+2x+5x^2+...$
$f^2(x)=\frac{f(x)-1}{x}$
解方程得$f(x)=\frac{1±\sqrt{1-4x}}{2x}$
由于收敛性我们去掉一个根，得到
$f(x)=\frac{1-\sqrt{1-4x}}{2x}$
用广义二项式定理展开，可以得到$x^n$的系数为$\frac{C_{2n}^n}{n+1}$，也就是卡特兰数。
设sn表示n个节点的二叉树的叶子数
那么我们有$s_n=\sum\limits_{i=1}^{n-1}s_i*h_{n-1-i},s_0=0,s_1=1$
写出$s_n$的生成函数：$g(x)=0+x+2x^2+6x^3+10x^4...$
$f(x)*g(x)=x+3x^2+5x^3+...$
$f(x)*g(x)=\frac{g(x)-x}{2x}$
将$f(x)$带入，得到$g(x)=\frac{x}{\sqrt{1-4x}}$
同样用二项式展开可以得到$x^n$的系数为$C_{2n-2}^{n-1}$
因此答案就是$\frac{s_n}{h_n}=\frac{n(n+1)}{4n+2}$

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 500010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n;
int main(){
//  freopen("a.in","r",stdin);
    n=read();
    printf("%.9lf\n",n*1.0*(n+1)/((ll)4*n-2));
    return 0;
}