python素数筛选法浅析


Posted in Python onMarch 19, 2018

原理:

素数,指在一个大于1的自然数中,除了1和此整数自身外,不能被其他自然数整除的数。在加密应用中起重要的位置,比如广为人知的RSA算法中,就是基于大整数的因式分解难题,寻找两个超大的素数然后相乘作为密钥的。一个比较常见的求素数的办法是埃拉托斯特尼筛法(the Sieve of Eratosthenes) ,说简单一点就是画表格,然后删表格,如图所示:

python素数筛选法浅析

从2开始依次往后面数,如果当前数字一个素数,那么就将所有其倍数的数从表中删除或者标记,然后最终得到所有的素数。

有一个优化:

标记2和3的倍数的时候,6被标记了两次。所以从i的平方开始标记,减少很多时间。

比如3的倍数从9开始标记,而不是6,并且每次加6。

除了2以外,所有素数都是奇数。奇数的平方还是奇数,如果再加上奇数就变成了偶数一定不会是素数,所以加偶数(2倍素数)。

预先处理了所有偶数。

注意:1既不是素数也不是合数,这里没有处理1。

#! prime.py 
import time 
 
def primes(n): 
 P = [] 
 f = [] 
 for i in range(n+1): 
  if i > 2 and i%2 == 0: 
   f.append(1) 
  else: 
   f.append(0) 
 
 i = 3 
 while i*i <= n: 
  if f[i] == 0: 
   j = i*i 
   while j <= n: 
    f[j] = 1 
    j += i+i 
  i += 2 
 
 P.append(2) 
 for i in range(3,n,2): 
  if f[i] == 0: 
   P.append(i) 
 
 return P 
 
def isPrime(n): 
 if n > 2 and n%2 == 0: 
  return 0 
 
 i = 3 
 while i*i <= n: 
  if n%i == 0: 
   return 0 
  i += 2 
 
 return 1 
 
def primeCnt(n): 
 cnt = 0 
 for i in range(2,n): 
  if isPrime(i): 
   cnt += 1 
 return cnt 
 
if __name__ == '__main__': 
 start = time.clock() 
 n = 10000000 
 P = primes(n); 
 print("There are %d primes less than %d"%(len(P),n)) 
 #for i in range(10): 
 # print(P[i]) 
 print("Time: %f"%(time.clock()-start)) 
 #for n in range(2,100000): 
 # if isPrime(n): 
 #  print("%d is prime"%n) 
  #print("%d is "%n + ("prime" if isPrime(n) else "not prime")) 
 
 start = time.clock() 
 n = 1000000 
 print("There are %d primes less than %d"%(primeCnt(n),n)) 
 print("Time: %f"%(time.clock()-start)

用素数筛选法求1千万以内的素数用了5.767s,

普通素数判断法求1百万以内的素数用了9.642s,

用C++素数筛选法求1亿以内的素数用了0.948s,

用C++普通素数判断法求1千万以内的素数用了3.965s,

可见解释语言确实比编译语言慢很多。

附C++程序,用了位压缩优化空间

#include <iostream> 
#include <cstdio> 
#include <algorithm> 
using namespace std; 
#define N 100000001 
 
unsigned f[(N>>5)+5]; 
int p[5761456],m; 
void init() 
{ 
  int i,j; 
  for(i=4;i<N;i+=2) 
    f[i>>5]|=1<<(i&0x1F); 
  p[m++]=2; 
  for(i=3;i*i<N;i+=2) 
    if(!(f[i>>5]&(1<<(i&0x1F)))) 
    { 
      p[m++]=i; 
      for(j=i*i;j<N;j+=i+i) 
        f[j>>5]|=1<<(j&0x1F); 
    } 
  for(;i<N;i+=2) 
    if(!(f[i>>5]&(1<<(i&0x1F)))) 
      p[m++]=i; 
} 
int is_prime(int n) 
{ 
  int i; 
  for(i=0;p[i]*p[i]<=n;i++) 
    if(n%p[i]==0) 
      return 0; 
  return 1; 
} 
int isPrime(int n) 
{ 
  if(n>2 && n%2==0) 
    return 0; 
  int i=3; 
  while(i*i<=n) 
  { 
    if(n%i==0) 
      return 0; 
    i+=2; 
  } 
  return 1; 
} 
int main() 
{ 
  int n=0,i; 
  clock_t st=clock(); 
  init(); 
  /*for(i=2;i<10000000;i++) 
    if(isPrime(i)) 
      n++;*/ 
  printf("%d %dms\n",m,clock()-st); 
  /*while(~scanf("%d",&n),n) 
  { 
    i=lower_bound(p,p+m,n+1)-p; 
    printf("%d\n",i); 
  }*/ 
  return 0; 
}

以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持三水点靠木。

Python 相关文章推荐
学习python的几条建议分享
Feb 10 Python
在Python中使用poplib模块收取邮件的教程
Apr 29 Python
在Python中处理时间之clock()方法的使用
May 22 Python
python面向对象_详谈类的继承与方法的重载
Jun 07 Python
python opencv旋转图像(保持图像不被裁减)
Jul 26 Python
win10 64bit下python NLTK安装教程
Sep 19 Python
Python基于滑动平均思想实现缺失数据填充的方法
Feb 21 Python
关于 Python opencv 使用中的 ValueError: too many values to unpack
Jun 28 Python
pandas如何处理缺失值
Jul 31 Python
python将图片转base64,实现前端显示
Jan 09 Python
pytorch 计算Parameter和FLOP的操作
Mar 04 Python
Python中Permission denied的解决方案
Apr 02 Python
python实现堆和索引堆的代码示例
Mar 19 #Python
python实现一个简单的并查集的示例代码
Mar 19 #Python
python使用筛选法计算小于给定数字的所有素数
Mar 19 #Python
python将每个单词按空格分开并保存到文件中
Mar 19 #Python
python将文本分每两行一组并保存到文件
Mar 19 #Python
python: line=f.readlines()消除line中\n的方法
Mar 19 #Python
Python File readlines() 使用方法
Mar 19 #Python
You might like
不用数据库的多用户文件自由上传投票系统(2)
2006/10/09 PHP
php单链表实现代码分享
2016/07/04 PHP
thinkPHP模板算术运算相关函数用法分析
2016/07/12 PHP
Nginx下ThinkPHP5的配置方法详解
2017/08/01 PHP
JS 建立对象的方法
2007/04/21 Javascript
javascript实现二分查找法实现代码
2007/11/12 Javascript
javascript 原型模式实现OOP的再研究
2009/04/09 Javascript
jQuery的实现原理的模拟代码 -1 核心部分
2010/08/01 Javascript
jQuery判断元素是否是隐藏的代码
2011/04/24 Javascript
js实现单一html页面两套css切换代码
2013/04/11 Javascript
IE的事件传递-event.cancelBubble示例介绍
2014/01/12 Javascript
jQuery实现简易的天天爱消除小游戏
2015/10/16 Javascript
Eclipse编辑jsp、js文件时卡死现象的解决办法汇总
2016/02/02 Javascript
实例讲解jQuery EasyUI tree中state属性慎用
2016/04/01 Javascript
Bootstrap基本组件学习笔记之列表组(11)
2016/12/07 Javascript
非常实用的vue导航钩子
2017/03/20 Javascript
vue2.0的contextmenu右键弹出菜单的实例代码
2017/07/24 Javascript
vue中SPA单页面应用程序详解
2017/11/07 Javascript
详解性能更优越的小程序图片懒加载方式
2018/07/18 Javascript
小程序实现多选框功能
2018/10/30 Javascript
Vue.js 父子组件通信的十种方式
2018/10/30 Javascript
swiper.js插件实现pc端文本上下滑动功能示例
2018/12/03 Javascript
vue-cli3单页构建大型项目方案
2020/04/07 Javascript
Vue作用域插槽实现方法及作用详解
2020/07/08 Javascript
理解JavaScript中的对象
2020/08/25 Javascript
vue el-upload上传文件的示例代码
2020/12/21 Vue.js
[48:44]2014 DOTA2国际邀请赛中国区预选赛5.21 TongFu VS HGT
2014/05/22 DOTA
[02:33]2018 DOTA2亚洲邀请赛回顾视频 再次拾起那些美妙的时刻
2018/04/10 DOTA
Python实现将Excel转换成xml的方法示例
2018/08/25 Python
python-Web-flask-视图内容和模板知识点西宁街
2019/08/23 Python
sklearn和keras的数据切分与交叉验证的实例详解
2020/06/19 Python
Fossil美国官网:化石手表、手袋、首饰及配饰
2019/02/17 全球购物
师范毕业生个人求职信
2013/12/09 职场文书
2014年党员承诺书范文
2014/05/20 职场文书
党员评议表自我评价范文
2014/10/20 职场文书
文化大革命观后感
2015/06/17 职场文书