标签:
输入一个数N(2 <= N <= 10^9)。
输出走法的数量 Mod 10007。
4
10
思路:这个在对角线的上方,就可以转换为,火车进站的问题,火车进站出来的肯定要小于等于进站的,那么也就相当于y>=x;那么就是卡特兰数,然后因为严格上方,所以我们要求的是F(n-1)卡特兰数,然后我们用lucas来求组合数取模。
1 #include<stdio.h> 2 #include<algorithm> 3 #include<iostream> 4 #include<string.h> 5 #include<queue> 6 #include<math.h> 7 using namespace std; 8 const int mod = 10007; 9 typedef long long LL; 10 LL N[10008]; 11 LL quick(LL n,LL m); 12 LL lucas(LL n,LL m); 13 int main(void) 14 { 15 N[0] = 1; 16 int i,j; 17 for(i = 1; i <= 10007; i++) 18 N[i] = N[i-1]*(LL)i%mod; 19 LL n; 20 scanf("%lld",&n);n--; 21 if(n==0)printf("1\n"); 22 else 23 {LL ask = lucas(2*n,n) - lucas(2*n,n-1); 24 ask = ask%mod + mod; 25 printf("%lld\n",2*ask%mod);} 26 return 0; 27 } 28 LL lucas(LL n,LL m) 29 { 30 if(m == 0)return 1; 31 LL x = n%mod; 32 LL y = m%mod; 33 if(y > x) 34 return 0; 35 LL ni = N[x-y]*N[y]%mod; 36 ni =quick(ni,mod-2); 37 ni = ni*N[x]%mod; 38 return ni*lucas(n/mod,m/mod); 39 } 40 LL quick(LL n,LL m) 41 { 42 LL ans = 1; 43 n %= mod; 44 while(m) 45 { 46 if(m&1) 47 ans = ans*n%mod; 48 n = n*n%mod; 49 m>>=1; 50 } 51 return ans; 52 }
标签:
原文地址:http://www.cnblogs.com/zzuli2sjy/p/5814025.html