洛谷3200 bzoj1485 HNOI2009 有趣的数列 卡特兰数

题目链接
题解：
手玩了一下小数据发现好像和卡特兰数有关。具体的原因是，我们把奇数和偶数位置上的数分别拿出来形成两个数列，那么要求就是这两个数列均为单调递增的并且偶数数列在对应位置上大于奇数数列。那么我们就从小到大往两个数列里填数，那么只要在对应位置上奇数数列先填数，偶数数列后填数即可得到一个满足题意的情况。而换句话说，就是必须奇数数列的该位置有数了，偶数数列的该位置才可能被填。那么我们把这个与卡特兰数经典的括号匹配问题做一个类比，在奇数数列填数就相当于在括号序列末尾加一个左括号，在偶数数列填数就相当于在括号末尾加一个右括号，所以这题的答案也是卡特兰数。
但是还有一个要解决的问题是$n<=1e6,p<=1e9$，而且$p$不一定是质数。由于算组合数的时候要做除法，但是不一定有逆元，所以我们要想一个解决办法。这个数据范围估计扩展卢卡斯的方法和扩展CRT的方法是会TLE的，但是我没试过。我们要求$\frac{1}{n+1}C_{2n}^n$，也就是$\frac{(2n)!}{n!(n+1)!}$，那么我们用一种类似扩展卢卡斯的办法。我们利用线性筛找到每一个质因数，我们算出每一个质因数在最后的答案中出现在多少次，然后用快速幂计算即可，这样只需要算乘法，就避免了算除法却可能没有逆元的问题。
代码：

#include <bits/stdc++.h>
using namespace std;

int n,book[2000010],p[2000010],cnt;
long long mod,ans=1;
inline long long ksm(long long x,long long y,long long mod)
{
    long long res=1;
    while(y)
    {
        if(y&1)
        res=(res*x)%mod;
        x=(x*x)%mod;
        y>>=1;
    }
    return res;
}
inline long long solve(long long x)
{
    long long res=0,y=x;
    while(y<=2*n)
    {
        res+=(2*n)/y-n/y-(n+1)/y;
        y*=x;
    }
    res=ksm(x,res,mod);
    return res;
}
int main()
{
    scanf("%d%lld",&n,&mod);
    for(int i=2;i<=2*n;++i)
    {
        if(!book[i])
        {
            p[++cnt]=i;
            ans=(ans*solve(i))%mod;
        }
        for(int j=1;j<=cnt&&p[j]*i<=2*n;++j)
        {
            book[p[j]*i]=1;
            if(i%p[j]==0)
            break;
        }
    }
    printf("%lld\n",ans);
    return 0;
}