[SCOI2010]生成字符串 卡特兰数

题目描述
lxhgww最近接到了一个生成字符串的任务，任务需要他把n个1和m个0组成字符串，但是任务还要求在组成的字符串中，在任意的前k个字符中，1的个数不能少于0的个数。现在lxhgww想要知道满足要求的字符串共有多少个，聪明的程序员们，你们能帮助他吗？
输入输出格式
输入格式：

输入数据是一行，包括2个数字n和m

输出格式：

输出数据是一行，包括1个数字，表示满足要求的字符串数目，这个数可能会很大，只需输出这个数除以20100403的余数

输入输出样例
输入样例#1： 复制
2 2
输出样例#1： 复制
2
云端评测
标签
各省省选
2010
四川
高性能
难度
省选/NOI-
时空限制
1000ms / 128MB

---

思路：
考虑从(0,0)出发的路径；
选1时 x+1,y+1;
选0时 x+1,y-1;
那么满足约束条件时应该是不能走到
y=-1的线下面；即求上述的方案数；
直接求肯定很困难，考虑间接求解；总的方案数为C(n+m,n)，那么不合法的方案数是多少呢？显然------不合法自然就是穿过 y = -1 的方案数，那么试着将 y=-1 的翻转下去，可以发现就是num_1 +1 ,num_0 -1 ，此时的方案数就是 C(n+m,n+1)

---

答案就是 C(n+m,n)-C(n+m,n+1) ；
lucas处理即可；

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;

#define mod 20100403
ll mod_pow(ll x, ll n, ll p){
    ll res = 1;
    while(n){
        if(n & 1) res =res * x % p;
        x = x * x % p;
        n >>= 1;
    }
    return res;
}
 
ll comb(ll n, ll m, ll p){
    if(m > n) return 0;
    ll ret = 1;
    m = min(n - m, m);
    for(int i = 1; i <= m; i ++){
        ll a = (n + i - m) % p;
        ll b = i % p;
        ret = ret * (a * mod_pow(b, p - 2, p) % p) % p;
    }
    return ret;
}
 
ll Lucas(ll n, ll m, ll p){
    if(m == 0) return 1;
    return comb(n % p, m % p, p) * Lucas(n / p, m / p, p) % p;
}
 
int main(){
    int T;
    ll n, m, p;
  //  scanf("%d", &T);
  //  while(T--){
        scanf("%lld%lld", &n, &m);
        printf("%lld\n", (Lucas(n+m, n, mod)-Lucas(n+m,n+1,mod)+mod)%mod);
  //  }
    return 0;
}