P3200 [HNOI2009]有趣的数列
题目描述
我们称一个长度为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的值。
输入输出样例
3 10
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)。
Solution:
本题zyys。
1.观察下列几种简单情况:
$(1)n=1:(1,2);$
$(2)n=2:(1,2,3,4),(1,3,2,4);$
$(3)n=3:(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);$
可以发现每组中$1$一定在第一个位置,$2n$一定在最后一个位置,由数列的性质可以证明;
每组数列都可:增加方案数为n-1;移动上一次$2n$的位置,增加方案数为$1$;在此基础上添加$2n-1$,可以发现$2n-1$允许插入的范围为$n+1,n+2,...,2n-1$,由乘法原理知,总方案数为$\frac{C(2n,n)}{n+1}$;
2.所以本题化简为求解模$p$剩余系下的卡特兰数,那么通过卡特兰数通项公式化简知$c[n]=2n\times (2n-1).....\times (n+2)/n!$,易证分子是可以整除分母的,那么统计约分后各个因子个数即可;
3.用线性筛法求出$[1,2n]$的$mindiv$(最小质因数),将分母分子分解质因数;
4.计算各质因数的幂取模相乘即可;
代码:
#include<bits/stdc++.h> #define il inline #define ll long long #define For(i,a,b) for(int (i)=(a);(i)<=(b);(i)++) #define Bor(i,a,b) for(int (i)=(b);(i)>=(a);(i)--) #define Max(a,b) ((a)>(b)?(a):(b)) #define Min(a,b) ((a)>(b)?(b):(a)) using namespace std; const int N=2000005; ll n,mod,ans,cnt,prime[N],isprime[N],minn[N],tot[N];il ll fast(ll s,ll k){ll ans=1;while(k){if(k&1)ans=ans*s%mod;k>>=1;s=s*s%mod;}return ans; }int main(){ios::sync_with_stdio(0);cin>>n>>mod;For(i,2,n<<1){if(!isprime[i])prime[++cnt]=i;for(int j=1;j<=cnt&&prime[j]*i<=(n<<1);j++){isprime[prime[j]*i]=1;minn[prime[j]*i]=prime[j];if(i%prime[j]==0)break;}}For(i,n+2,n<<1) {int p=i;while(minn[p])tot[minn[p]]++,p=p/minn[p];tot[p]++;}For(i,1,n) {int p=i;while(minn[p])tot[minn[p]]--,p=p/minn[p];tot[p]--;}ll ans=1;For(i,1,cnt) ans=ans*fast(prime[i],tot[prime[i]])%mod;cout<<ans;return 0; }
转载于:https://www.cnblogs.com/five20/p/9215107.html
P3200 [HNOI2009]有趣的数列相关推荐
- P3200 [HNOI2009]有趣的数列 (巧妙转换,卡特兰数,分解质因数取模运算)
P3200 [HNOI2009]有趣的数列 分析: 这题是个思维题,通过打表也可以容易发现就是卡特兰数 讲一下,是如何巧妙的转换成卡特兰数的 由数列满足的三个条件,可以发现第 2i2i2i 位(偶数位 ...
- 分解质因数-洛谷P3200 [HNOI2009]有趣的数列
https://www.luogu.org/problem/show?pid=3200 首先,我们不能保证要求的数的逆元和模域互质: 所以我们要用分解质因数来抵消除法: 其实逆元的话即使可行也会超时: ...
- 【BZOJ1485】[HNOI2009]有趣的数列(组合数学)
[BZOJ1485][HNOI2009]有趣的数列(组合数学) 题面 BZOJ 洛谷 题解 从小往大填数,要么填在最小的奇数位置,要么填在最小的偶数位置. 偶数位置填的数的个数不能超过奇数位置填的数的 ...
- 「BZOJ1485」[HNOI2009] 有趣的数列 (卡特兰数列)
「BZOJ1485」[HNOI2009] 有趣的数列 Description 我们称一个长度为2n的数列是有趣的,当且仅当该数列满足以下三个条件: (1)它是从1到2n共2n个整数的一个排列{ai ...
- [HNOI2009]有趣的数列
[HNOI2009]有趣的数列 有一个长度为2n的1~2n的全排列,保证其奇数项递增,偶数项递增,并且相邻的奇数项和偶数项,后面的偶数项大于奇数项的方案数\(mod\ p,n<=1000000, ...
- [HNOI2009] 有趣的数列
题目描述 我们称一个长度为2n的数列是有趣的,当且仅当该数列满足以下三个条件: (1)它是从1到2n共2n个整数的一个排列{ai}: (2)所有的奇数项满足a1<a3<...<a2n ...
- BZOJ1485: [HNOI2009]有趣的数列
题解:考虑按顺序从小到大,加入数字,将加入奇数位置看作入栈,加入偶数位置看作出栈.为什么可以?考虑我们要保证相邻奇数小于偶数,所以必须先填上一个奇数的位置才能填偶数的位置,既时刻保证奇数>=偶数 ...
- [Catalan数三连]网格有趣的数列树屋阶梯
如何让孩子爱上打表 Catalan数 Catalan数是组合数学中一个常出现在各种计数问题中的数列. 以比利时的数学家欧仁·查理·卡塔兰 (1814–1894)的名字来命名. 先丢个公式(设第n项为$ ...
- P3200-[HNOI2009]有趣的数列【卡特兰数】
正题 题目链接:https://www.luogu.com.cn/problem/P3200 题目大意 求一个长度为2∗n2*n2∗n的排列要求 奇数位和偶数位分别递增 相邻的偶数位大于奇数位 解题思 ...
- 《信息学奥赛一本通提高篇》第6章 组合数学
例1 计算系数(NOIP2011提高) 信息学奥赛一本通(C++版)在线评测系统 NOIP2011计算系数_nanhan27的博客-CSDN博客 「NOIP2011」 计算系数 - 组合数_TbYan ...
最新文章
- 经济学人:清华大学或将迅速成为全球科研领军者
- 大型网站核心架构要素--可用性
- 雷观(十五):提高生产力和程序员价值的2种方法
- oracle回收ddl权限,oracle禁止指定用户DDL操作触发器
- python制作界面怎么触发事件_python模拟事件触发机制详解
- linux捕捉信号sigint失败,为shell布置陷阱:trap捕捉信号方法论
- OpenShift 4 - 配置OpenShift可使用的外部Image Registry和Mirror Registry
- mysql中名词解析
- 拓端tecdat|R语言中的生存分析Survival analysis晚期肺癌患者4例
- MySQL入门(详细总结)
- Spring Security(二) UserDetailsService 和 PasswordEncoder 密码解析器 详解
- Redis 命令 - 在线参考
- 2021-4-25科学计算器开高次根号
- 春运抢火车票攻略汇总(持续更新中)
- 【字符集五】c++标准库<locale>
- Pycharm:解决点击重运行(Rerun)没有反应
- 有梦想的闲鱼之团队介绍
- 《卡耐基三部曲》(Yanlz+VR云游戏+Unity+SteamVR+云技术+5G+AI+人性的弱点+人性的优点+语言的突破+术业有专攻+世界观+人生观+价值观+志同道合+不卑不亢+立钻哥哥++==)
- tplink无线桥接WDS小结
- 【记忆化搜索】UVA 10118 free candy