牛客网&HDOJ 卡特兰数相关题目

BU冰糖雪梨里的梨PT

已于 2023-01-23 22:29:36 修改

阅读量83

点赞数

分类专栏：刷刷算法文章标签：算法 Powered by 金山文档

于 2023-01-23 22:19:14 首次发布

本文链接：https://blog.csdn.net/weixin_45339670/article/details/128753351

版权

刷刷算法专栏收录该内容

18 篇文章 0 订阅

订阅专栏

文章探讨了卡特兰数在不同问题中的应用，如序列构造、质因数分解和高精度计算等。它指出卡特兰数的计算公式是C(2n,n)/(n+1)，并涉及模p的计算。文中提供的代码示例展示了如何进行高精度计算和快速幂运算来处理卡特兰数的模运算问题。

摘要由CSDN通过智能技术生成

对于n个0，m个1 的任意前缀0的个数不小于1的种数是 C(n+m, n)-C(n+m, m-1)。

卡特兰数是n=m 的特殊情况，通项公式为 C(2n, n) / (n+1)

NC15559-杯子

题解写在了这里，还挺详细的

杯子题解【更新】题解被删了，等恢复吧 ╮(﹀_﹀”)╭

【质因数分解】NC20075-[HNOI2009]有趣的数列

就是卡特兰数

依次考虑1~2n中的每个数，挨个放入序列中的某个位置（可以选奇数位置或者是偶数位置），要求中间的任何时候选的奇数位置的个数>=偶数位置的个数，走的方案数就是卡特兰数。卡特兰数计算公式C(2n, n)/ (n+1)

mod p的p是任给的，所以模数和除数不一定互质，不互质一定无法进行乘法逆元的运算

所以考察的是第n项卡特兰数模p

【高精度】NC19891-[AHOI2012]树屋阶梯

还是卡特兰数

n层阶梯n个角，一个角对应一个矩形，枚举填充第i个角的矩形拉到左下角，将整个图形一分为二

f(n) = f(0)f(n-1) + f(1)f(n-2) + ... + f(n-2)f(1) + f(n-1)f(0)

低精度的扩展卡特兰数（m！=n的情况）代码：

const int N = ■■■■■;
const int p = ■■■■■;
i64 fac[N + 10]; // 求i的阶乘mod p
i64 inv[N + 10]; // 求i的阶乘(mod p)的逆元

i64 pow(i64 a, i64 b) //快速幂
{
    i64 ans = 1;
    while (b) {
        if (b & 1) ans = ans * a % p;
        a = a * a % p; b >>= 1;
    }
    return ans;
}

void init(void)
{
    // i！
    fac[0] = 1;
    for (int i = 1; i <= N; i++) fac[i] = fac[i - 1] * i%p;

    // 费马小定理求i！的逆元
    for (int i = 0; i <= N; i++) inv[i] = pow(fac[i], p - 2);
}

i64 C(int n, int m) { return fac[n] * inv[m] % p*inv[n - m] % p; }

i64 Catalan(int n, int m) //卡特兰数的结论是：n个1和m个0组成序列，每个前缀里1个数>=0个数时，可行的序列个数是C(n+m,n)-C(n+m,m-1)
{
    return (C(n + m, n) - C(n + m, m - 1) + p) % p;
}