卡特兰数
记把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);
}