题目描述
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
首先,我们设选1为(1,1),选0为(1,-1)
目标就是(n+m,n-m)
总方案数为C(n+m,n),因为有n+m个位置,放n个1
然后要减去不合法的即线路通过y=-1的。将线路与y=-1交点的左边沿着y=-1做对称操作,则最后等价于从(0,-2)走到(n+m,n-m)的方案数
所以向上走n-m+2
则有x-y=n-m+2
x+y=n+m
x=n+1,y=m-1
所以不合法方案为C(n+m,n+1)
ans=C(n+m,n)-C(n+m,n+1)
求这些用模逆元,O(n)求解
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
long long A[],B[],ans1,ans2;
int n,m,Mod=;
int main()
{int i,j;
cin>>n>>m;
A[]=;
for (i=;i<=n+m;i++)
A[i]=((Mod-Mod/i)*A[Mod%i])%Mod;
for (i=;i<=n+m;i++)
A[i]=(A[i]*A[i-])%Mod;
B[]=;
for (i=;i<=n+m;i++)
B[i]=(B[i-]*i)%Mod;
ans1=(((B[n+m]*A[n])%Mod)*A[m])%Mod;
ans2=(((B[n+m]*A[n+])%Mod)*A[m-])%Mod;
cout<<(ans1-ans2+Mod)%Mod;
}