---恢复内容开始---
题目描述
我们称一个长度为2n的数列是有趣的,当且仅当该数列满足以下三个条件:
(1)它是从1到2n共2n个整数的一个排列{ai};
(2)所有的奇数项满足a1<a3<...<a2n-1,所有的偶数项满足a2<a4<...<a2n;
(3)任意相邻的两项a2i-1与a2i(1<=i<=n)满足奇数项小于偶数项,即:a2i-1<a2i。
现在的任务是:对于给定的n,请求出有多少个不同的长度为2n的有趣的数列。因为最后的答案可能很大,所以只要求输出答案 mod P的值。
输入输出格式
输入格式:
输入文件只包含用空格隔开的两个整数n和P。输入数据保证,50%的数据满足n<=1000,100%的数据满足n<=1000000且P<=1000000000。
输出格式:
仅含一个整数,表示不同的长度为2n的有趣的数列个数mod P的值。
输入输出样例
输入样例#1: 复制
3 10
输出样例#1: 复制
5
对应的5个有趣的数列分别为(1,2,3,4,5,6),(1,2,3,5,4,6),(1,3,2,4,5,6),(1,3,2,5,4,6),(1,4,2,5,3,6)。 ------------------------------------------------------------------------------------------------- 做之前我知道这是一道
关于Catalan数的问题
但是
当我读完题目后
我是真的蒙了
它是怎么和Catalan数挂上关系的
于是
我再次被卑微了
————————————————————————————————————————————————————————————————————————————————————
考察这样相邻的两项:a_{2i-1}a2i−1与a_{2i}a2i,根据题目的第二条原则显然有a_{2i-1}<a_{2i}a2i−1<a2i。
而根据第一条原则又有奇数是递增的。
所以有a_1<a_3<...<a_{2i-1}<a_{2i}a1<a3<...<a2i−1<a2i。
这个时候可以联想到这道经典的题目。
我们可以将奇数项看为入栈,偶数项看为出栈。
发现和入栈次数必须大于出栈次数的条件恰好相符。
所以可以使用卡特兰数求解。
但是直接暴力用公式是会爆的
所以需要一些优化
#include <bits/stdc++.h>
using namespace std;
long long i,j,n,m,ans,sum,x,a[],b[],p[];
int main()
{
scanf("%lld%lld",&n,&m); for (i=; i<=*n; i++)
{
if (a[i]==)
{
sum++;
p[sum]=i;
a[i]=sum;
}//sum累加,p数组标记第sum个质数为i ,a[i]表示i这个数是由第sum个质数筛得的
for (j=; j<=sum; j++) //枚举当前已得到的所有质数,用这些质数筛选出在范围内的可以被这些质数筛得合数
if (p[j]*i<=*n)
a[p[j]*i]=j;
else
break;
//并把筛到的合数标记是被第几个质数 筛到的
} for (i=n+; i<=*n; i++)
{
x=i;
while (x>)
{
b[a[x]]++;
x=x/p[a[x]];
} //a[x]表示第x个数是由第a[x]个质数筛得的, 然后把表示第a[x]个质数个数的数组累加
//做完上述操作后将x除以已累加的质数
} for (i=; i<=n; i++)
{
x=i;
while (x>)
{
b[a[x]]--;
x=x/p[a[x]];
} //重复上述操作
}
ans=;
for (i=; i<=sum; i++)
for (j=; j<=b[i]; j++)
ans=ans*p[i]%m; //p[i]表示第i个质数是什么,b[i]表示有第i个质数共有几个
printf("%lld",ans);
return ;
}