施罗德数

问题引入

施罗德数（又称超级卡特兰数）用来描述表格中从$(0,0)$到$(n,n)$的格路中，只能使用$(1,0),(0,1),(1,1)$三种移动方式，且始终在对角线下方而不能超过对角线的路径数

（图片来自百度百科）

通项公式

施罗德数前几项为$1,2,6,22,90,394...$

递推公式为$F_n=F_{n-1}+{\sum }_{i=0}^{n-1}{F}_i\ast F_{n-i-1}$

计算

通过通项公式计算的时间复杂度为$O(n^2)$，显然不能接受，这时需要引入一个新的递推序列：

$1,1,3,11,45,197,...$

不难发现该序列和施罗德数除了第一项外，只需要乘以二就能得到施罗德数，该序列的递推公式为：

$(n+1)F_n=(6n-3)\ast F_{n-1}-(n-2)\ast F_{n-2}$

递推代码

ll f[maxn], g[maxn];

ll qkp(ll x, ll n, ll p) {
    ll ans = 1;
    x %= p;
    while (n) {
        if (n & 1) ans = ans * x % Mod;
        x = x * x % Mod;
        n >>= 1;
    }
    return ans;
}

ll inv(ll a, ll p) {
    return qkp(a, p - 2, p);
}

//(n+1)f[n]=(6n-3)*f[n-1]-(n-2)*f[n-2]
void init(int n, ll p) {
    f[0] = f[1] = 1;
    g[0] = 1, g[1] = 2;
    for (int i = 2; i <= n; i++) {
        f[i] = ((6LL * i - 3) * f[i - 1] - (1LL * i - 2) * f[i - 2] + p) % p * inv(i + 1, p) % p;
        g[i] = f[i] * 2 % p;
    }
}