[BZOJ1856][Scoi2010]字符串（卡特兰数+组合数学）

题目描述

传送门

题解

答案 Cnn+m−Cn+1n+m ，模数是质数所以直接处理阶乘然后快速幂计算逆元就行
至于这个公式的推导可以参考卡特兰数的非常规分析
首先 Cnn+m 是总的方案数，从中减去不合法的方案
对于一个不合法的方案，假设从第2k+1位开始不合法，那么之前一定有k+1个0，k个1，如果将0,1互换就变成了一共有n+1个1，m-1个0
反过来，对于任意一个n+1个1，m-1个0组成的方案，任选一位将其前面的01互换，就变成了一个不合法的方案

代码

#include
#include
#include
#include
#include
using namespace std;
#define LL long long
#define N 2000005
#define Mod 20100403

int n,m;
LL mul[N],ans;

void calc()
{
    mul[0]=1LL;
    for (int i=1;i<=n+m;++i) mul[i]=mul[i-1]*(LL)i%Mod;
}
LL fast_pow(LL a,int p)
{
    LL ans=1;
    for (;p;p>>=1,a=a*a%Mod)
        if (p&1)
            ans=ans*a%Mod;
    return ans;
}
LL C(int n,int m)
{
    return mul[n]*fast_pow(mul[m]*mul[n-m]%Mod,Mod-2)%Mod;
}
int main()
{
    scanf("%d%d",&n,&m);
    if (n<m) {puts("0");return 0;}
    calc();
    ans=((C(n+m,n)-C(n+m,n+1))%Mod+Mod)%Mod;
    printf("%lld\n",ans);
}