洛谷-1641 [SCOI2010]生成字符串

题目描述
lxhgww最近接到了一个生成字符串的任务，任务需要他把n个1和m个0组成字符串，但是任务还要求在组成的字符串中，在任意的前k个字符中，1的个数不能少于0的个数。现在lxhgww想要知道满足要求的字符串共有多少个，聪明的程序员们，你们能帮助他吗？
输入格式
输入数据是一行，包括2个数字n和m
输出格式
输出数据是一行，包括1个数字，表示满足要求的字符串数目，这个数可能会很大，只需输出这个数除以20100403的余数

输入输出样例
输入 #1
2 2

输出 #1
2

说明/提示
limitation
每点2秒
对于30%的数据，保证1<=m<=n<=1000
对于100%的数据，保证1<=m<=n<=1000000
来源：SCOI 2010

解释：把它抽象成在一个二维坐标系中，在$（0,0）$点到$（n+m,m）$点的个数，其中不经过$y=-1$这条直线，首先$（0,0）$点到$（n+m,m）$总共为$C(n+m,m)$个，那么经过$y=-1$的点，我们把它对称为从$（0,-2）$出发到达$（n+m,m）$的个数为$C(n+m,m-1)$个，所以
$ans=C(n+m,m)-C(n+m,m-1)$

#include<iostream>
#define mod 20100403
using namespace std;
long long fact[2000003]={1};
long long pow(long long a,long long b){
    long long ret=1%mod;
    while(b){
        if(b&1) ret=ret*a%mod;
        a=a*a%mod;b>>=1;
    }
    return ret;
}
long long C(long long n,long long m){
    return fact[n]*pow(fact[m],mod-2)%mod*pow(fact[n-m],mod-2)%mod;
}
int main(){
    for(int i=1;i<2000003;i++) fact[i]=fact[i-1]*i%mod;
    long long n=0,m=0;
    cin>>n>>m;
    long long ret=C(n+m,m)-C(n+m,m-1);
    ret%=mod;ret+=mod;ret%=mod;
    cout<<ret<<endl;
    return 0;
}