bzoj 1485: [HNOI2009]有趣的数列

卡特兰数（貌似和卡常数有什么py交易）

黄学长是用DP，输出看的。。然而本蒟蒻只知道有卡特兰数这个名字。。（其实，，这个sb  DP能不能写出来还不一定。。）

公式好多，百度一下。。（看起来最初的定义式特别像卷积的样子，然而并不是。。）

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<map>
#include<queue>
#define N 1000005
#define inf 1000000000
#define LL long long
using namespace std;
const int mod=1000000007;
LL ans=1;
int n,P,cnt;
int prime[1000005],mn[2000005],num[2000005];
bool del[2000005];
void get_prime()
{
    for (int i=2; i<=2*n; i++)
    {
        if (!del[i]) prime[++cnt]=i,mn[i]=cnt;
        for (int j=1; prime[j]*i<=2*n && j<=cnt; j++)
        {
            del[prime[j]*i]=1; mn[prime[j]*i]=j;
            if (i%prime[j]==0) break;
        }
    }
}
void add(int x, int f)
{
    while (x!=1)
    {
        num[mn[x]]+=f;
        x/=prime[mn[x]];
    }
}
int main()
{
    scanf("%d %d",&n,&P); get_prime();
    for (int i=2*n; i>n+1; i--) add(i,1);
    for (int i=1; i<=n; i++) add(i,-1);
    for (int i=1; i<=cnt; i++)
        while (num[i]--) ans=(ans*prime[i])%P;
    printf("%d
",ans);
    return 0;
}