分解质因数-洛谷P3200 [HNOI2009]有趣的数列
https://www.luogu.org/problem/show?pid=3200
首先,我们不能保证要求的数的逆元和模域互质;
所以我们要用分解质因数来抵消除法;
其实逆元的话即使可行也会超时;
那么我转载了,实在没什么可以说的;
另外卡特兰数
http://baike.baidu.com/link?url=St3mmth0khr1jUoD9Vwdroupnfajo6hhTSgwvOkjAPrP0Htt12nZjsMue4T_5JhMopRqlhgAkCt2dDzd378Kg8xjsSYwGn3J_CMLgsvI4Psdhj3z0s4zTucxc1v6dRlP
话说这个方法好巧妙啊
http://blog.csdn.net/jiangshibiao/article/details/24009239
【转化】就是求卡特兰数。
【初始代码】
#include<cstdio>
using namespace std;
typedef long long ll;
ll prime[200005],a[200005];
bool f[2000005];
ll temp,n,p,i,j,cnt,mod;
ll pow(ll a,ll b)
{ ll ans; for (ans=1;b;b>>=1,a=a*a%mod) if (b&1) ans=ans*a%mod; return ans;
}
int main()
{ scanf("%lld%lld",&n,&mod); for (i=2;i<=n*2;i++) { if (!f[i]) prime[++cnt]=i; for (j=1;j<=cnt&&prime[j]*i<=n*2;j++) f[prime[j]*i]=true; } for (i=n+2;i<=n*2;i++) for (j=1,p=i;j<=cnt&&p;j++) while (p%prime[j]==0) a[j]++,p/=prime[j]; for (i=2;i<=n;i++) for (j=1,p=i;j<=cnt&&p;j++) while (p%prime[j]==0) a[j]--,p/=prime[j]; temp=1; for (j=1;j<=cnt;j++) if (a[j]) temp=(temp*pow(prime[j],a[j]))%mod; printf("%lld",temp);for (;;); return 0;
}
用欧拉筛法,O(n)的效率求出每个质数。然后枚举阶乘,像质数表一样把一个数给分解。但是效率很低。
【优化1】如果一个数是合数,我们可以把它的某个因子记下来。然后我们同样从开始枚举阶乘,而且是倒着枚举。对于每个数,如果它是合数,我就把它分解。比如,设f[n]为结果中含有n因子的个数。u是n的一个约数。那么我们可以f[u]+=f[n],f[n/u]+=f[n]。这样就不用多次用快速幂了。直到n是质数为止。
【优化2】开始可以把1–n的f[i]设为-1,把n+2–2*n(注意,最后要除n+1,所以从n+2开始)的f[i]设为1.这样只需1次循环。
【AC代码】
#include<cstdio>
using namespace std;
typedef long long ll;
ll prime[200005],a[2000005],come[2000005];
ll temp,n,p,i,j,cnt,mod;
ll pow(ll a,ll b)
{ ll ans; for (ans=1;b;b=b/2,a=a*a%mod) if (b&1) ans=ans*a%mod; return ans;
}
int main()
{ scanf("%lld%lld",&n,&mod); for (i=2;i<=n*2;i++) { if (!come[i]) prime[++cnt]=i; for (j=1;j<=cnt&&prime[j]*i<=n*2;j++) come[prime[j]*i]=i; } temp=1; for (i=2;i<=n;i++) a[i]=-1; for (i=n+2;i<=2*n;i++) a[i]=1; for (i=n*2;i>1;i--) if (come[i]) { a[come[i]]+=a[i]; a[i/come[i]]+=a[i]; } else temp=temp*pow(i,a[i])%mod; printf("%lld",temp); return 0;
}
转载于:https://www.cnblogs.com/largecube233/p/6797912.html
分解质因数-洛谷P3200 [HNOI2009]有趣的数列相关推荐
- P3200 [HNOI2009]有趣的数列 (巧妙转换,卡特兰数,分解质因数取模运算)
P3200 [HNOI2009]有趣的数列 分析: 这题是个思维题,通过打表也可以容易发现就是卡特兰数 讲一下,是如何巧妙的转换成卡特兰数的 由数列满足的三个条件,可以发现第 2i2i2i 位(偶数位 ...
- 【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, ...
- 信息学奥赛一本通 1937:【06NOIP普及组】数列 | 洛谷 P1062 [NOIP2006 普及组] 数列
[题目链接] ybt 1937:[06NOIP普及组]数列 洛谷 P1062 [NOIP2006 普及组] 数列 [题目考点] 1. 数制 [解题思路] 如果k为2,那么这个数列 第1项为202^02 ...
- [HNOI2009] 有趣的数列
题目描述 我们称一个长度为2n的数列是有趣的,当且仅当该数列满足以下三个条件: (1)它是从1到2n共2n个整数的一个排列{ai}: (2)所有的奇数项满足a1<a3<...<a2n ...
- [质因数分解]樱花 洛谷P1445
题目背景 又到了一年樱花盛开的时节.Vani 和妹子一起去看樱花的时候,找到了一棵大大的樱花树,上面开满了粉红色的樱花.Vani 粗略估计了一下,一共有足足 n! 片花瓣. Vani 轻柔地对她说:& ...
- 质因数分解(洛谷P1075题题解,Java语言描述)
题目要求 P1075题目链接 分析 其实很水,无非是从从小到大遍历,从2开始,每一次judge一下质数,如果是最小质数,那另一半就是答案. 一定要break,保证一个结果... 如果break不用的话 ...
- 洛谷 P2513 [HAOI2009]逆序对数列
题目描述 对于一个数列{ai},如果有i<j且ai>aj,那么我们称ai与aj为一对逆序对数.若对于任意一个由1~n自然数组成的数列,可以很容易求出有多少个逆序对数.那么逆序对数为k的这样 ...
最新文章
- pycharm中run pytest解决
- sqlmap 跑access_sqlmap注入Access
- Java并发编程实战笔记2:对象的组合
- 错误:pytube.exceptions.RegexMatchError
- python两数相加代码_一边学编程,一边学语数外,用python编程全排列10以内两数加法...
- centos 7.0 查看selinux状态|关闭|开启
- ICCV 2019《Zero-Shot Grounding of Objects from Natural Language Queries》论文笔记
- windows XP常见进程(个人总结)
- [MyBatis]DAO层只写接口,不用写实现类
- Fedora9中的gcc
- Dash的快速入门将使您在5分钟内进入“ Hello World”
- node.js编程错误记录集
- 作者:石乾新(1989-),男,贵州大学公共管理学院硕士生。
- python开发流程视频_自学Python的步骤和方法,Python入门全面视频无偿分享,使用Python开发的一些优点...
- c++ 原子操作 赋值_多线程操作可见性
- fiddler抓包第一课--手机数据抓包
- 医学生建议计算机系学生植发,00后都开始“秃”了?皮肤科专家:不想成为“地中海”,这件事情要记牢...
- 电脑桌面上文件夹图标右上角有双箭头
- IT人看十四主星代表
- 《房债》书中的精髓:现代金融系统让全世界的人同处一艘经济大船上,一些人搞坏了船,所有的人都跟着遭殃。