51nod1120 卡特兰数+Lucas定理

最新推荐文章于 2022-02-11 13:58:53 发布

Dave_L

最新推荐文章于 2022-02-11 13:58:53 发布

阅读量574

点赞数

分类专栏：暑期集训刷题算法复习（新手）数学-公式数论

本文为博主原创文章，未经博主允许也可以转载。

本文链接：https://blog.csdn.net/FrankAx/article/details/76390282

版权

暑期集训刷题算法复习（新手）同时被 3 个专栏收录

123 篇文章 0 订阅

订阅专栏

数论

45 篇文章 0 订阅

订阅专栏

数学-公式

43 篇文章 1 订阅

订阅专栏

1120机器人走方格 V3

基准时间限制： 1 秒空间限制： 131072 KB 分值: 80 难度：5级算法题

关注

取消关注

N * N的方格，从左上到右下画一条线。一个机器人从左上走到右下，只能向右或向下走。并要求只能在这条线的上面或下面走，不能穿越这条线，有多少种不同的走法？由于方法数量可能很大，只需要输出Mod 10007的结果。

Input

输入一个数N(2 <= N <= 10^9)。

Output

输出走法的数量 Mod 10007。

Input示例

Output示例

   
   
    
    #include <bits/stdc++.h>
#define LL long long
using namespace std;
const int MOD = 10007;
const int AX = 1e6+66;
LL f[AX];
void init(int mod){                //f[n]为n的阶乘
 f[0] = 1;
 for(int i = 1 ; i <= mod ;i ++ ){
  f[i] = f[i-1] * i % mod;
 }
}
void ex_gcd(LL a,LL b,LL& d,LL& x,LL& y)
{
    if(!b) { d = a; x = 1; y = 0; }
    else{ ex_gcd(b, a%b, d, y, x); y -= x*(a/b); }
}
LL inv(LL a,LL m)           //欧几里得求逆元
{
    LL d, x, y;
    ex_gcd(a, m, d, x, y);
    return d == 1 ? (x+m)%m : -1;
} 
LL Lucas(LL m , LL n , LL p){             //Lucas定理
 LL  res = 1;
 while(n && m){
  LL n1 = n % p ;
  LL m1 = m % p ;
  res = res * f[n1] * inv(f[n1-m1],p) * inv(f[m1],p)%p;
  n /= p;
  m /= p;
 }
 return ( res % p + p ) % p;
}
int main(){
 int n;
 init(MOD);
 scanf("%d",&n);
 n--;
 LL ans = Lucas(n,2*n,MOD)*inv(n+1,MOD)%MOD;
 printf("%lld\n",2*ans%MOD);
 return 0;
}