[BZOJ1485][HNOI2009]有趣的数列（卡特兰数+组合数学）

题目描述

传送门

题解

答案是卡特兰数，因为…打表可知？
但是实际上这道题可以抽象成一个卡特兰数的经典问题——出栈序列问题。
也就是说，问题等价为1-n依次入栈有多少可能的出栈序列
题目要求满足两个条件：奇数项和偶数项分别单增，并且每两个相邻的奇数项和偶数项满足奇数项小于偶数项
可以把这个问题看做：将1-n填到2n个空格中，填到奇数项时必须满足上一个奇数项已经有数，这时相当于是一个元素入栈；填到偶数项时必须满足上一个偶数项有数，并且它的前一个奇数项有数，这时相当于是一个元素出栈
这是由于栈的性质是：入栈时新元素放入栈的顶端；出栈是必须满足它的上面元素全部出栈，并且栈不为空。这两个性质分别和上面两个条件对应

代码

#include<algorithm>#include<iostream>#include<cstring>#include<cstdio>#include<cmath>using namespace std;#define N 2000000#define LL long longint n,Mod;int p[N+1],prime[N+1],cnt[N+1],num[N+1];LL ans;void get(){    for (int i=2;i<=N;++i)    {        if (!p[i]) prime[++prime[0]]=i,num[i]=prime[0];        for (int j=1;j<=prime[0]&&i*prime[j]<=N;++j)        {            p[i*prime[j]]=1;            if (i%prime[j]==0) break;        }    }}int main(){    get();    scanf("%d%d",&n,&Mod);    for (int i=n+1;i<=2*n;++i)    {        int x=i;        for (int j=1;x>1&&j<=prime[0]&&prime[j]*prime[j]<=x;++j)            while (x%prime[j]==0) ++cnt[j],x/=prime[j];        if (x>1) ++cnt[num[x]];    }    for (int i=1;i<=n+1;++i)    {        int x=i;        for (int j=1;x>1&&j<=prime[0]&&prime[j]*prime[j]<=x;++j)            while (x%prime[j]==0) --cnt[j],x/=prime[j];        if (x>1) --cnt[num[x]];    }    ans=1LL;    for (int i=1;i<=prime[0];++i)        while (cnt[i])        {            ans=ans*(LL)prime[i]%Mod;            --cnt[i];        }    printf("%lld\n",ans);}