bzoj 1856: [Scoi2010]字符串（卡特兰数）

1856: [Scoi2010]字符串

Time Limit: 5 Sec   Memory Limit: 64 MB
Submit: 1604   Solved: 904
[ Submit][ Status][ Discuss]

Description

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

Input

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

Output

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

Sample Input

2 2

Sample Output

2

题目可以转化为从(0, 0)走到(n+m, n)不穿过y=-1的方案数，数字为1表示往上走，数字为0表示往右走

很像卡特兰数：

给出一个n*n的正方形棋盘，求从左下角到右上角的不经过对角线的所有走法刚好是F[n]

即C(2n, n)-C(2n, n-1)

那么类比推理，这题答案就是C(n+m, m)-C(n+m, m-1)

#include<stdio.h>
#define LL long long
#define mod 20100403
LL Pow(LL a, LL b);
LL C(LL m, LL n);
LL Lucas(LL m, LL n);
long long jc[2000005] = {1};
int main(void)
{
	int i;
	LL m, n;
	for(i=1;i<=2000000;i++)
		jc[i] = (jc[i-1]*i)%mod;
	while(scanf("%lld%lld", &n, &m)!=EOF)
		printf("%lld\n", (Lucas(n+m, m)-Lucas(n+m, m-1)+mod)%mod);
	return 0;
}

LL Pow(LL a, LL b)
{
	LL ans;
	ans = 1;
	while(b)
	{
		if(b%2==1)
			ans = (ans*a)%mod;
		a = (a*a)%mod;
		b /= 2;
	}
	return ans;
}

LL C(LL n, LL m)    
{
	LL ans;
	if(n<m)
		return 0;
	ans = (jc[n]*Pow((jc[m]*jc[n-m])%mod, mod-2)%mod)%mod;
	return ans;
}

LL Lucas(LL n, LL m) 
{
	if(m==0)
		return 1;
	return (Lucas(n/mod, m/mod)*C(n%mod, m%mod))%mod;
}