CSU1320(详解卡特兰数)

题意:
zzy今天刚买了两个水瓢$A$和$B$，容量都是为$1$升。他有一个无限大的水缸，用$A$往里面加水，用$B$往外淘水，用$B$淘水时要保证水缸里至少还有$1$升水，就这样他用了$n$次$A$，也用了$n$次$B$。
给出$n$求出共有多少种淘水方式，输出的答案模$1000000007$

思路:
我们首先先转化下这个问题。我们可以看做有这样一个长度为$2n$的01串，共有$n$个0和$n$个1，我们要解出从左向右扫，在任意位置经过0不多于1的串的个数。
哈哈，这个问题是不是看着有些眼熟。嗯，确实许多问题都可以转化成这个模型，例如能够匹配的括号数、在方格中不经过对角线从左上走到右下的种数、n个数的出栈次序，甚至是。。。n个节点的二叉树种类个数。最后一个虽然看起来没什么关系，但是他们确实是一个问题。
实际上这些问题都是卡特兰(Catalan)数的问题，现在我们就来推导下它的通项。
要求出满足要求的串的种类，我们可以用总的个数减掉不合格的种类个数。首先总的个数很容易求，在$2n$个位置挑选出$n$个位置填$1$,剩下的补$0$，共有$\binom{2n}{n}$种选法。那么不合格的选法呢，我们考虑现在扫到了第$2m+1$位，此时前面有$m+1$个$0$和$m$个$1$，而对于后面的$n-m$个$1$和$n-m-1$个$0$，我们可以这些将$0$和$1$换下位置，就得到了一个含有$n+1$个$0$和$n-1$个$1$的01串，那么我们可以说每个不合格的序列一定会对应这样的一个01串。反过来，对于任意这样的串，我们也都可以找到一组不合格的序列(构造方法同上)。那么我们就得出了这种串和不合格序列一一对应的结论了，那么这样的串很明显有$\binom{2n}{n+1}$种。所以最后答案则是$\binom{2n}{n}-\binom{2n}{n+1}$种，即$\frac{\binom{2n}{n}}{n+1}$种。
那它又为什么和二叉树的形态个数等价呢？我们再来证明一下。首先我们设$f(n)$为$n$个结点的二叉树的形态个数，那么我们可以把它看成由根、$i$个结点的左子树和$n-i-1$个结点的右子树组成，其中$0\leq i \leq n-1$。那么即可得出

$$f(n) = \sum_{0}^{n-1}{f(i)f(n-i-1)}$$
这只是一个不太好算的递推式，我们还要求出它的通项。
我们这里使用生成函数法，考虑函数 $$F(x) = \sum_{k=0}^{\infty}{f(k)x^k}$$ 求平方，有 $$F^2(x) = \sum^{\infty}_{k=0}{(\sum^{k}_{i=0}{f(i)f(k-i)})x^k}$$ 由前式替换，有 $$F^2(x)=\sum^{\infty}_{k=0}{f(k+1)x^k}$$ 故可凑出 $$F(x) = xF^2(x)+1$$ 移项，有二次方程 $$xF^2(x)-F(x)+1=0$$ 解为 $$F(x)=\frac{1\pm\sqrt{(1-4x)}}{2x}$$ 由f(0)=1,故取 $$F(x)=\frac{1-\sqrt{(1-4x)}}{2x}$$ 由牛顿二项式展开定理， $$(1-4x)^{\frac{1}{2}}=\sum^{\infty}_{k=0}{\binom{\frac{1}{2}}{k}{(-4x)}^k}$$ 带回原式，整理得 $$F(x)=\sum_{k=0}^{\infty}{\binom{\frac{1}{2}}{k+1}{(-1)}^k2^{2k+1}x^k}$$ 所以有 $$f(n)=\binom{\frac{1}{2}}{n+1}{(-1)}^n2^{2n+1}=\frac{\frac{1}{2}(\frac{1}{2}-1)(\frac{1}{2}-2)···(\frac{1}{2}-n)}{(n+1)!}\\=\frac{1}{n+1}·\frac{(2n)!}{n!·n!}=\frac{\binom{2n}{n}}{n+1}$$ 然后我们就可以发现此式和上面的通项是一样的。
现在我们就得出了卡特兰数的通项了，也就可以直接解出该题的答案了。实际上，卡特兰数还有一个 $O(n)$的递推公式 $f(n+1)= \frac{4n-2}{n+1}f(n)$,这个式子也经常的使用。
现在我们就可以愉快的解出这道题了。
贴个代码:

#include <cstdio>

typedef long long ll;
const int MAXN = 10000 + 100;
const ll MOD = 1000000007;
ll ans[MAXN];

ll extgcd(ll a, ll b, ll &x, ll &y){
    ll d = a;
    if(b == 0LL){
        x = 1; y = 0;
    }else{
        d = extgcd(b, a % b, y, x);
        y -= (a / b) * x;
    }
    return d;
}

ll mod_inverse(ll a, ll MOD){
    ll x, y;
    extgcd(a, MOD, x ,y);
    return (x % MOD + MOD) % MOD;
}

void Init(){
    ans[0] = ans[1] = 1;
    for(int i = 2; i < MAXN; ++i){
        ans[i] = (((ans[i - 1] * (4 * i - 2)) % MOD) * mod_inverse(i + 1, MOD)) % MOD;
    }
    return;
}

int main(int argc, char const *argv[]){
    Init();
    int x;
    while(~scanf("%d", &x)){
        printf("%lld\n", ans[x]);
    }
    return 0;
}