BZOJ 1485 [HNOI2009]有趣的数列

卡特兰数

记把ai放在奇数项为0操作，放在偶数项为1操作

从a1->an的所有操作中，把0看作入栈，1看作出栈，则就是卡特兰数统计出栈次序（证明详见百度百科）。

公式：
h(n)=C(2n,n)/(n+1) (n=0,1,2,…)

顺便附上卡特兰数的前几个数，以便日后能及时反应过来：
1, 1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 208012, 742900, 2674440, 9694845, 35357670, 129644790, 477638700, 1767263190, 6564120420, 24466267020, 91482563640, 343059613650, 1289904147324, 4861946401452, …

#include<cstdio>
#define N 2000005
using namespace std;
int prime[N], minprime[N], cnt=0, cntprime[N];
bool notprime[N];
void init_prime(int lim)
{
    for(int i = 2; i <= lim; i++)
    {
        if(!notprime[i])prime[++cnt]=i,minprime[i]=cnt;
        for(int j = 1; j <= cnt && prime[j]*i<=lim; j++)
        {
            notprime[prime[j]*i]=1;
            minprime[prime[j]*i]=j;
            if(i%prime[j]==0)break;
        }
    }
}
void calc(int x, int k)
{
    while(x!=1)
    {
        cntprime[minprime[x]]+=k;
        x/=prime[minprime[x]];
    }
}
int main()
{
    int n, p;
    scanf("%d%d",&n,&p);
    init_prime(2*n);
    for(int i = n+1; i <= 2*n; i++)calc(i,1);
    for(int i = 1; i <= n; i++)calc(i,-1);
    calc(n+1,-1);
    int ans=1;
    for(int i = 1; i <= cnt; i++)
        while(cntprime[i])
        {
            ans=1ll*ans*prime[i]%p;
            cntprime[i]--;

        }
    printf("%d\n",ans);

}