P3200-[HNOI2009]有趣的数列【卡特兰数】

QuantAsk

于 2020-11-30 07:00:27 发布

阅读量253

点赞数

分类专栏：数论and数学文章标签： luogu HNOI 卡特兰数

本文链接：https://blog.csdn.net/Mr_wuyongcong/article/details/110355265

版权

数论and数学专栏收录该内容

261 篇文章

订阅专栏

博客主要解析了一道洛谷上的编程题目，涉及到序列填充和卡特兰数的概念。题目要求构造一个特定的2n阶排列，其中奇数位置和偶数位置分别递增，并且偶数位置的数值始终大于相邻的奇数位置。解题策略转换为计算第n个卡特兰数，并通过质因数分解处理除数，以达到O(nlogn)的时间复杂度。代码实现中包含了质数判断、质因数分解和计数数组等算法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

正题

题目链接:https://www.luogu.com.cn/problem/P3200

题目大意

求一个长度为 $2 * n$ 的排列要求

奇数位和偶数位分别递增
相邻的偶数位大于奇数位

解题思路

可以看做是一个 $2 * n$ 的序列按顺序填进奇数和偶数位，然后因为第二个要求所以奇数位在任何时候都得比偶数位的要多，就转换为了求第 $n$ 个卡特兰数了。

然后因为 $p$ 不是质数所以要质因数分解来除数，时间复杂度 $O(n\log n)$

$c o d e$

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=2e6+10;
int n,p,tot,pri[N],cnt[N];
bool v[N];
int main()
{
	scanf("%d%d",&n,&p);
	for(int i=2;i<=2*n;i++){
		if(!v[i])pri[++tot]=i;
		for(int j=1;j<=tot&&i*pri[j]<=n;j++){
			v[i*pri[j]]=1;
			if(i%pri[j]==0)break;
		}
	}
	for(int i=n+2;i<=2*n;i++){
		int x=i;
		for(int j=1;pri[j]*pri[j]<=x&&j<=tot;j++)
			while(x%pri[j]==0)x/=pri[j],cnt[j]++;
		if(x!=1)cnt[lower_bound(pri+1,pri+1+tot,x)-pri]++;
	}
	for(int i=1;i<=n;i++){
		int x=i;
		for(int j=1;pri[j]*pri[j]<=x&&j<=tot;j++)
			while(x%pri[j]==0)x/=pri[j],cnt[j]--;
		if(x!=1)cnt[lower_bound(pri+1,pri+1+tot,x)-pri]--;
	}
	long long ans=1;
	for(int i=1;i<=tot;i++)
		while(cnt[i])ans=ans*pri[i]%p,cnt[i]--;
	printf("%d\n",ans);
}