bzoj1485: [HNOI2009]有趣的数列 卡特兰数

50%的dp应该很好想。一看题解是卡特兰数，我就震惊了。

我们来思考为什么是卡特兰数。。。

每个数只能走一次，设奇数位的为x 偶数位的为y 那么到达(N,N)一直不能穿过直线y=x

由于路径只能往右和往上走，所以这就是个卡特兰数。

如果还不清楚，那么可以理解为如果向上到(a,b)就是把下一个奇数位放a+b，如果向右到(a,b)就是把下一个偶数数位放a+b

#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
using namespace std;
#define maxn 2000010
bool vis[maxn];
int pri[maxn],cnt,a[maxn];
int n,p;
void init()
{
    for(int i=2;i<=2*n;i++)
    {
        if(!vis[i])
        {
            cnt++;pri[cnt]=i;
            a[i]=cnt;
        }
        for(int j=1;j<=cnt&&i*pri[j]<=2*n;j++)
        {
            vis[i*pri[j]]=1;a[i*pri[j]]=j;
            if(i%pri[j]==0)
            {
                break;
            }
        }
    }
}
int num[maxn];
void add(int now,int f)
{
    while(now!=1)
	{
		num[a[now]]+=f;
		now/=pri[a[now]];
	}
}
int main()
{
    scanf("%d%d",&n,&p);
    init();
    long long ans=1;
    for(int i=n+1;i<=n*2;i++) add(i,1);
    for(int i=1;i<=n;i++)add(i,-1);
    add(n+1,-1);
    for(int i=1;i<=cnt;i++)
		while(num[i]--)ans=(ans*pri[i])%p;
    printf("%lld\n",ans);
    return 0;
}