[BZOJ1856][SCOI2010]字符串

最新推荐文章于 2021-05-17 03:45:19 发布

原创最新推荐文章于 2021-05-17 03:45:19 发布 · 532 阅读

0 ·

CC 4.0 BY-SA版权

数学专栏收录该内容

23 篇文章

订阅专栏

本文介绍了一种基于Catalan数的推广——类Catalan数的计算方法，并给出了一段C++代码实现。该算法利用组合数学原理，通过快速幂运算提高计算效率。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

类Catalan数定义
按照推倒Catalan数的方法推倒这题的公式
答案即为

C m n + m - C m - 1 n + m

$C_{n+m}^m-C_{n+m}^{m-1}$

/**************************************************************
    Problem: 1856
    User: di4CoveRy
    Language: C++
    Result: Accepted
    Time:800 ms
    Memory:1288 kb
****************************************************************/

#include <iostream>
#include <cstdio>
#define mod 20100403LL

using namespace std;
typedef long long LL;

LL n,m;
LL qpow(LL a,LL b) {
    if (b == 0) return 0;
    if (b == 1) return a;
    if (b % 2 == 0) {
        LL tmp = qpow(a,b/2);
        return (tmp * tmp) % mod;
    } else {
        LL tmp = qpow(a,b-1);
        return (tmp * a) % mod;
    }
}

LL C(LL n,LL m) {
    LL JCN = 1LL , JCM = 1LL , JCNM = 1LL;
    for (int i=1;i<=n;i++) JCN = JCN * i % mod;
    for (int i=1;i<=m;i++) JCM = JCM * i % mod;
    for (int i=1;i<=n-m;i++) JCNM = JCNM * i % mod;
    return (JCN * qpow( (JCM * JCNM) % mod , mod - 2 )) % mod;
}

int main(){
    cin >> n >> m;
    if (m > n) {
        cout << 0 << endl;
        return 0;
    }
    LL ans = (C(n+m,n) - C(n+m,m-1) + mod) % mod;
    cout << ans << endl;
    return 0;
}