1856: [Scoi2010]字符串
Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 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;
}