BZOJ1856:[SCOI2010]字符串

Description

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

Input

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

Output

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

Sample Input

2 2

Sample Output

HINT

对于30%的数据，保证1<=m<=n<=1000
对于100%的数据，保证1<=m<=n<=1000000

题解：

从坐标(0,x)出发，每次向右上或右下移动√2个单位，移动过程中触碰过x轴，直到移动到(z,y)，其方案集合为F1；从坐标(0,-x)出发，每次向右上或右下移动√2个单位，直到移动到(z,y)，其方案集合为F2。若x,y,z>0，则有|F1|=|F2|。

证明：对于F1中每一种走法，对其进行这样一种变换：设其与x轴第一个交点为(w,0)，将横坐标0~w的路线以x轴对称，则唯一地得到了F2中的一种方案。对于F2的一种方案，显然在F1中有且只有一种方案可以变换为成它（该变换是自己的逆变换），所以F1与F2中元素一一对应，集合大小相等。

在该题中，考虑将方案变为坐标轴中的路线，每加入一个字符即向右移动一次，加入1为右上，加入0为右下，则问题变成了：从(0,1)移动到(n+m,n-m+1)，不触碰x轴的方案数。

ANS=C(n+m,n-m)-C(n+m,m-1)

代码：

 1 const
 2   mo:int64=20100403;
 3 var
 4   i,j,k,l,n,m:longint;
 5   a:array[0..2000001]of int64;
 6   ans:int64;
 7 function ksm(x,y:int64):int64;
 8 var z:int64;
 9 begin
10   z:=1;
11   while y>0 do
12   begin
13     if y mod 2=1 then z:=(z*x)mod mo;
14     y:=y div 2; x:=(x*x)mod mo;
15   end;
16   exit(z);
17 end;
18 function c(x,y:int64):int64;
19 begin
20   if x=y then exit(1);
21   exit(((a[x]*ksm(a[y],mo-2))mod mo*ksm(a[x-y],mo-2))mod mo);
22 end;
23 begin
24   readln(n,m);
25   a[0]:=1;
26   for i:=1 to n+m do a[i]:=a[i-1]*i mod mo;
27   ans:=c(n+m,n-m)-c(n+m,m-1);
28   ans:=(ans mod mo+mo)mod mo;
29   writeln(ans);
30 end.

View Code

转载于:https://www.cnblogs.com/GhostReach/p/6253758.html