素数

素数又称质数,是除了1和本身之外不能被其他数整除的一类数字。用数学语言描述即是,对给定的正整数n,如果对于任意的正整数a(1<a<n),都有n%a!=0成立,那么称n素数;否则称为合数
注:1 既不是素数也不是合数。

素数的判断

一个整数n如果不能被2,3,...,n-1整除,那么n为素数。使用遍历的方法判断,时间复杂度为O(n)。但在某些数据量大的题目中,该复杂度实际有点大,需要更快的方法。
注意到:如果在2~n-1中存在n的约数,设为k,即n%k==0,那么由k*(n/k)==n可知,n/k也是n的一个约数。且kn/k中一定满足其中一个小于等于sqrt(n),另一个大于等于sqrt(n)。因此只需判断n能否被2,3,...,sqrt(n)中的一个整除,即可判断是否为素数。该算法时间复杂度为O(sqrt(n))

代码如下:

bool is_prime(int n) {if (n <= 1) {return false;} else {for (int i = 2; i * i <= n; i++) {if (n % i == 0) return false;}return true;}
}

素数表的获取

通过判断是否为素数,可以直接得出1~n范围内的所有素数,即从1~n开始枚举,如果是素数就加入素数表。枚举时间复杂度为O(n),判断是否为素数时间复杂度为O(sqrt(n)),总的时间复杂度为O(n*sqrt(n))

求100以内的素数表
代码如下:


int prime[MAXN], p_num = 0; //prime存放所有的素数,p_num为素数的个数
bool p[MAXN] = {0};    //p[i]==ture即为素数
void find_prime() {for (int i = 1; i < MAXN; i++) {if (is_prime(i)) {prime[p_num++] = i;p[i] = true;}}
}

完整的代码如下:

#include <bits/stdc++.h>
using namespace std;const int MAXN = 101;
bool is_prime(int n) {if (n <= 1) {return false;} else {for (int i = 2; i * i <= n; i++) {if (n % i == 0) return false;}return true;}
}int prime[MAXN], p_num = 0;  // prime存放所有的素数,p_num为素数的个数
bool p[MAXN] = {0};          // p[i]==ture即为素数
void find_prime() {for (int i = 1; i < MAXN; i++) {if (is_prime(i) == true) {prime[p_num++] = i;p[i] = true;}}
}int main() {find_prime();for (int i = 0; i < p_num; i++) {printf("%d", prime[i]);}return 0;
}

埃氏筛法(Eratosthenes)

上述的算法对于10^5级别的都可以承受,但如果数字范围更大,那么将会崩溃。使用埃氏筛法时间复杂度可低至O(nloglogn)。该算法从小到大枚举所有的数,对每一个素数,筛去它的所有倍数,剩下的都是素数了。

举例子说明:求1~15所有的素数。
①2是素数(唯一需要事先确定的),因此筛去所有2的倍数,即{4,6,8,10,12,14},剩{2,3,5,7,9,11,13,15}。
②3没有被筛去,因此3是素数,筛去3所有的倍数,即{6,9,12,15},剩{2,3,5,7,11,13}。
③4已经被筛去,因此4不是素数。
④5没有被筛,因此5是素数,筛去所有5的倍数,即{10,15},剩{2,3,5,7,11,13}。
⑤6已经被筛去,因此6不是素数。
⑥7没有被筛去,因此7是素数,筛去所有7的倍数,即{14},剩{2,3,5,7,11,13}。

以此类推…
最终得到所有素数,即{2,3,5,7,11,13}。
由例子可以发现,如果到达某个数x时,它没有被前面的数筛掉,那么它就是素数。

100以内素数代码如下:

const int MAXN = 101;
int prime[MAXN], p_num = 0;
bool p[MAXN] = {0};  // i is a prime,return false
void find_prime() {for (int i = 2; i < MAXN; i++) {if (p[i] == false) {  // if i is a primeprime[p_num++] = i;for (int j = i + i; j < MAXN; j += i) {p[j] = true;}}}
}

总结

  1. 1不是素数,也不是合数;
  2. 涉及到素数问题,通常打表是一个不错的方法;
  3. 素数的表长最少要比n大于1;
  4. 注意时间复杂度。

质因子分解

一个正整数可以写成一个或者多个质数的乘积形式,可以先把素数表打印出来。前面说过,如果一个正整数n存在除了1和本身之外的因子,那么一定在sqrt(n)左右成对出现。当应用到质因子上来,会有一个强化结论

对于一个正整数n,如果存在[2,n]之间的质因子,那么这些质因子要不全都小于等于sqrt(n),要不只存在一个大于等于sqrt(n)的质因子,其余的质因子都小于等于sqrt(n)

①首先定义一个结构体,用来存放质因子和对应的个数:

struct factor{int x,cnt;
}fac[10];

fac数组只需要到10个就足够了,因为十个素数相乘已经超出了int范围。
②枚举1~sqrt(n)范围内的所有质因子p,判断p是否为n的因子。

  • 如果p是n的因子,那么给fac数组增加质因子p,并初始化其个数为0。如果p还是n的质因子,就让n不断的除以p,每次令p的个数加一,直到p不再是n的因子为止。
if(n%prime[i]==0){fac[num].x=prime[i];fac[num].cnt=0;while(n%prime[i]==0){fac[num].cnt++;n/=prime[i];}num++;
}
  • 如果p不是n的因子就直接跳过。

③如果在上述步骤结束后,n依旧大于1,说明有且仅有一个大于sqrt(n)的质因子,有可能是n本身,这时将该质因子加入fac数组,并令其个数为1。

if(n!=1){fac[num].x=n;fac[num++].cnt=1;
}

最后,fac数组里存放的就是质因子分解的结果。

素数问题与质因子分解相关推荐

  1. Lightoj-1356 Prime Independence(质因子分解)(Hopcroft-Karp优化的最大匹配)

    题意: 找出一个集合中的最大独立集,任意两数字之间不能是素数倍数的关系. 思路: 最大独立集,必然是二分图. 最大数字50w,考虑对每个数质因子分解,然后枚举所有除去一个质因子后的数是否存在,存在则建 ...

  2. 2021年中国大学生程序设计竞赛 女生专场 - 热身赛 Problem C. 口算训练(质因子分解)

    题目分析:判断 xxx 是 yyy 的倍数,等价于质因子分解后, yyy 中的每个质因子的出现次数都小于等于其在 xxx 中的出现次数. 那么对于每次询问 [l,r,d][l,r,d][l,r,d], ...

  3. CodeForces - 1419E Decryption(质因子分解+构造)

    题目链接:点击查看 题目大意:给出一个数 n ,现在需要将其所有大于 1 的因子重新排列成一个首尾相接的环,规定每次操作可以在相邻两个位置中加入这两个数的 lcm,问最少需要操作多少次,才能使得环上所 ...

  4. 牛客多校9 - Groundhog Chasing Death(质因子分解+思维)

    题目链接:点击查看 题目大意:给出 a , b , c , d , x , y ,求 题目分析:因为涉及到了 gcd 的乘积运算,那么易知不同质因子的贡献是相互独立的,首先我们就可以先将 x 和 y ...

  5. 洛谷——P2043 质因子分解

    P2043 质因子分解 题目描述 对N!进行质因子分解. 输入格式 输入数据仅有一行包含一个正整数N,N<=10000. 输出格式 输出数据包含若干行,每行两个正整数p,a,中间用一个空格隔开. ...

  6. python质因子分解_质因子分解_个人文章 - SegmentFault 思否

    质因子分解的问题就是给定一个n使得n能够分解为多个因子的乘积形式,并且相同因子用指数形式表示: 例如180=2^23^25; 对于这个问题,很好理解,我们的目的就是寻找其因子,通常的方法也就是从0开始 ...

  7. 质因子分解 Python

    质因子分解 PythonPythonPython 要做质因子分解,首先需要明白什么是质数,以及如何快速判断质数. 质数 质数,也称素数,是只能被1和其本身整除的数,规定1不是质数. 质数的判断可以详见 ...

  8. 质因子个数————质因子分解

    文章目录 1 题目 2 解析 3 参考代码 1 题目 质因子个数 时间限制 1000 ms 内存限制 32768 KB 代码长度限制 100 KB 判断程序 Standard (来自 小小) 题目描述 ...

  9. E Groundhog Chasing Death(2020牛客暑期多校训练营(第九场))(思维+费马小定理+质因子分解)

    E Groundhog Chasing Death(2020牛客暑期多校训练营(第九场))(思维+费马小定理+质因子分解) 链接:https://ac.nowcoder.com/acm/contest ...

最新文章

  1. php转换文字Unicode,php实现将中文转为unicode的方法
  2. 移远NB-IOT模块的BC28与M26与BC26对比
  3. 基本数据类型转换 || 自动类型转换与强制类型转换
  4. [上海]LinkCoder第四期活动——Jeffrey Richter:Win 8应用开发与.NET4.5
  5. C++连接mysql的两种方式(ADO连接和mysql api连接)
  6. win10安装部署网络测试工具 NetCat 用于Spark跑测试
  7. 史上最丧心病狂的商品定价套路:如何从数学角度,榨干你身上的每一分钱
  8. 基于PHPEnv的本地环境搭建—PHP第一个项目:HelloWorld(从安装到运行)
  9. linux进程泄露命令明文参数,sshpass 使Linux可以明文参数输入SSH密码(示例代码)
  10. Android Testing学习02 HelloTesting 项目建立与执行
  11. 三维数据入库发布流程之3dMAX数据
  12. Xpose 建立新工程
  13. JavaScript就这么回事(好收藏,哪天忘了可以查一查)
  14. 计算机上的放大快捷键,电脑放大镜热键 怎么取消电脑放大镜快捷键?
  15. 《IBM SPSS Modeler数据与文本挖掘实战》之文本挖掘算法
  16. 【虾皮选品分享】马来西亚/菲律宾市场时尚美妆热销趋势
  17. 2022年中职网络空间安全国赛竞赛题解析仅代表自己的建议——2022年中职网络安全国赛竞赛试题9解析
  18. 洞悉物联网发展1000问之为什么新零售不仅仅是无人零售?
  19. 自行车测速器怎么调_山地车变速器的调整技巧
  20. JSOI 2008-help

热门文章

  1. 推荐三款实用的录音转文字软件给你
  2. SP_Flash_Tool对应MTK CPU型号速查表汇总分析
  3. 全球及中国商业航天产业应用模式探析及投资布局态势分析报告2021-2027年版
  4. R语言使用epiDisplay包的logistic.display函数获取二分类logistic回归模型的汇总统计信息(自变量初始和调整后的优势比及置信区间,回归系数的Wald检验的p值、简化汇总结果
  5. php 倒计时 考试,php实现倒计时
  6. 利用四位共阳数码管显示小数
  7. QT 多线程异步获取htpp信息 网络通信
  8. 2018.11.4 东华杯(骇极) REVERSE What's it wp
  9. python连接数据库mysql失败_mysql数据库连接失败是什么原因
  10. 怎样才能压缩PDF文件大小