【51Nod四级题】机器人走方格 V3

链接：

https://www.51nod.com/Challenge/Problem.html#problemId=1120

题意：

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

输入样例

4

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

输出样例

10

输出走法的数量 Mod 10007。

思路：

基本思路如上图，若没有对角线的限制，那么种类数就是C(2*n,n)

有了限制后，其实上图的思路，所有越过对角线的种类数是C(2*n,n-1)，所以单个三角形的种类数是 C(2*n,n)- C(2*n,n-1)，通过它可以推导出等于C(2*n,n)/(n+1)。也就是说 C(2*n,n)- C(2*n,n-1)=C(2*n,n)/(n+1)。

也就是求卡特兰数的，Catalan[n] = C(2*n,n) / (n+1) = C(2*n,n) % mod * inv[n+1],inv[n+1]为n+1的逆元，根据费马小定理，可以通过快速幂快速求出。

因为n的数据范围较大，所以要用到卢卡斯定理：若p为素数，那么C(m,n)%p=Lucas(m,n)%p = Lucas(m/p,n/p)*C(m%p,n%p)%p，从而我们可以递归的可以求出Lucas(m,n)，当n==0，返回1。因为方格含有两个三角形，所以Catalan[n]*2 即是最终答案。当然也可以用递推的方式：

AC代码：

#include<bits/stdc++.h>
#define ll long long
#define endl "\n"
const int N=1e5+5;
const int INF=0x3f3f3f3f;
const int MOD=10007;
const double eps=1e-6;
using namespace std;
int n;
int quick(int a, int b){
    a%=MOD;
    int res=1;
    for(;b;b>>=1,a=a*a%MOD) if(b&1) res=res*a%MOD;
    return res;
}
int C(ll a, ll b){
    if(b>a) return 0;
    int ans=1;
    for(int i=1;i<=b;++i)
        ans=ans*(a-i+1)%MOD*quick(i,MOD-2)%MOD;
    return ans;
}
// 卢卡斯定理
int lucas(int a, int b){
    if(b==0) return 1;
    return lucas(a/MOD,b/MOD)%MOD*C(a%MOD,b%MOD)%MOD;
}
void run(){
    cin>>n;
    // 方格转化为点
    n--;
    int ans=lucas(2*n,n)%MOD*quick(n+1,MOD-2)%MOD;
    cout<<2*ans%MOD<<endl;
}
int main(){
    ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
    //freopen("../input.txt","r",stdin);
    //freopen("../output.txt","w",stdout);
    int _=1;
    //cin>>_;
    for(int i=0;i<_;++i) run();
    return 0;
}