题目描述
题解
答案是卡特兰数,因为…打表可知?
但是实际上这道题可以抽象成一个卡特兰数的经典问题——出栈序列问题。
也就是说,问题等价为1-n依次入栈有多少可能的出栈序列
题目要求满足两个条件:奇数项和偶数项分别单增,并且每两个相邻的奇数项和偶数项满足奇数项小于偶数项
可以把这个问题看做:将1-n填到2n个空格中,填到奇数项时必须满足上一个奇数项已经有数,这时相当于是一个元素入栈;填到偶数项时必须满足上一个偶数项有数,并且它的前一个奇数项有数,这时相当于是一个元素出栈
这是由于栈的性质是:入栈时新元素放入栈的顶端;出栈是必须满足它的上面元素全部出栈,并且栈不为空。这两个性质分别和上面两个条件对应
代码
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define N 2000000
#define LL long long
int 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);
}