素数
来源:互联网 发布:手机借款软件排行 编辑:程序博客网 时间:2024/06/14 10:43
源代码连接:http://blog.csdn.net/qq_38569113/article/details/71170279
素数:大于1,并且只能被1和自身整除的数。
方法一:简单判断
bool prime(int n){ if(n<2) return false; else { for(int i=2;i<sqrt(n);i++) { if(n%i==0) return false; } } return true;}方法二:埃氏筛法(时间复杂度为O(nlogn))
要得到自然数n以内的全部素数,必须把不大于
详细步骤:当n=25时
列出2以后的所有序列:
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
标出序列中的第一个素数,也就是2,划掉2的倍数,序列变成:
2 3 5 7 9 11 13 15 17 19 21 23 25
下一个素数是3,将序列中3的倍数划掉,序列变成:
2 3 5 7 11 13 17 19 23 25
下一个素数是5,同样将序列中5的倍数划掉,序列成了:
2 3 5 7 11 13 17 19 23
因为23小于5的平方,跳出循环.
结论:2到25之间的素数是:2 3 5 7 11 13 17 19 23。
代码实现:
void isprime(int n){ memset(prime,1,sizeof(prime)); prime[0]=prime[1]=0; for(int i=2;i<=sqrt(n);i++) { if(prime[i]) for(int j=i*i;j<n;j+=i) prime[j]=0; }}方法三:欧拉筛法
线性筛 效率 O(N)
可以先估算下N之内素数个数,num≈
如果N>10^7,就不适合全素数打表。
代码实现:
int visit[N];int prime[N];int cnt_prime(int n){ memset(visit,1,sizeof(visit)); visit[0]=visit[1]=0; int num=0; for(int i=2;i<n;i++) { if(visit[i]) prime[num++]=i; for(int j=1;j<=num&&i*prime[j]<=n;j++) { visit[i*prime[j]]=0; if(i%prime[j]==0) break; } } return num;//素数个数}因为 每一个合数都可以表示为若干个质数之积
if(i%prime[j]==0) break;
保证了每个合数都能被它最小的一个质因子筛掉,所以大大提高了效率。 例如x为偶数,只需要被x/2筛就好。
对于奇数,假设i为9时,prime素数表里有2、3、5、7,先筛掉2*9=18,然后筛掉3*9=27,判断9%3==0,跳出。这是因为9可以分成3*3,如果继续晒5*9,相当于筛3*3*5=3*15,所以在i=15的时候筛就好,而埃氏筛法就没有这个优化。
方法四:费马小定理判断大素数
互质数:公因数只有1的两个非零自然数。
费马小定理
对于任意正整数a,假如p是质数,且gcd(a,p)=1,那么ap−1 %p≡1。即:假如a是正整数,p是质数,且a,p互质(即两者只有一个公约数1),那么a的(p-1)次方除以p的余数恒等于1。
ap−1 %p=1是p为质数的必要条件而非充分条件!!
卡迈克尔数
对于任意正整数a,假如p是合数,且gcd(a,p)=1,那么ap−1 %p≡1。根据费马小定理和卡迈克尔数可以知道有些能满足费马小定理的数(p),不一定是质数,也就是说根据费马小定理判断素数并不一定可靠,但我们可以多测几个a的值增大正确率。
利用上述结论,对于给定的整数p,可以设计一个素数判定算法。
1. 随机选取整数a,2≤ a≤ n-1,计算d=ap−1 %p。
当d不等于1时,n是合数;
当d等于1时,n则很可能是素数,对于本次判断来说,判断错误的概率为1/2。
2. 如此重复多次,可以将判断错误的概率降低至期望值以下。
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- 素数
- java高并发程序设计总结三:JDK并发包之ReentrantLock重入锁
- 归并排序(逆序对的解决)
- VMware 黑屏解决方法
- cocos2d-x--使用doxygen+graphviz生成源码文档
- windows 环境下c++系统栈大小调整
- 素数
- HTML5 ul li 去前面的圆点 和 横排
- NGINX基于Cookie和Header的负载均衡会话保持
- Spring中标签属性name和id的区别
- 确定整数是否是回文?空间O(1)时间复杂度O(logn)
- VTK基础概念-圆柱体绘制
- 细说单点登录
- 无题
- clipToPadding 和clipChildren 的区别