【BZOJ 1856】 [Scoi2010]字符串 折线法

方程ans=C(n+m,m)-C(n+m,m-1)类似卡特兰数的公式，卡特兰数通项公式：h[n]=c(2*n,n)-C(2*n,n-1)

证明用折线法，不过懒得画图了：

卡特兰数证明：http://blog.sina.com.cn/s/blog_6917f47301010cno.html

类似于上面我们每增加一个1就相当于斜向上走一步，0就是斜向下走一步，那么在没有限制的情况下方案数就是C(n+m,m)（从n+m中选m次向下走），现在减去不合法的情况，类似于上面，假定有一次越过了x轴那么剩下的能够向下走的就只剩下m-1次了，类似的，我们的方程也可以写成C(n+m,n)-C(n+m,n+1)均可

#include<cstdio>
#include<cstring>
#include<iostream>
#define LL long long
#define Mod 20100403
#define maxn 2000020
using namespace std;
LL n,m,fac[maxn];

void exgcd(LL a,LL b,LL& x,LL& y){
	if(b==0){x=1,y=0;}
	else{
		exgcd(b,a%b,y,x);
		y-=x*(a/b);
	}
}
LL C(LL a,LL b){
	LL x,y;
	LL c=fac[b]*fac[a-b]%Mod;
	exgcd(c,Mod,x,y);
	LL ans=fac[a]*x%Mod;
	while(ans<0)ans+=Mod;
	return ans;
}

int main(){
	fac[1]=1;
	for(LL i=2;i<=2000000;i++)fac[i]=fac[i-1]*i%Mod;
	scanf("%lld%lld",&n,&m);
	if(n<m)return puts("0"),0;
	n+=m;
	printf("%lld",(C(n,m)-C(n,m-1)+Mod)%Mod);
	return 0;
}