链接:
题意:
N * N的方格,从左上到右下画一条线。一个机器人从左上走到右下,只能向右或向下走。并要求只能在这条线的上面或下面走,不能穿越这条线,有多少种不同的走法?由于方法数量可能很大,只需要输出Mod 10007的结果。
输入样例
4
输入一个数N(2 <= N <= 10^9)。
输出样例
10
输出走法的数量 Mod 10007。
思路:
![](https://i-blog.csdnimg.cn/blog_migrate/a3c3c0493205d244896eccf597ec134f.png)
基本思路如上图,若没有对角线的限制,那么种类数就是C(2*n,n)
![](https://i-blog.csdnimg.cn/blog_migrate/d184a9c83cdae8d83a0ea20dcbebd8b5.png)
有了限制后,其实上图的思路,所有越过对角线的种类数是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;
}