python使用集合实现筛选法求素数-python素数筛选法浅析
原理:
素数,指在一个大于1的自然数中,除了1和此整数自身外,不能被其他自然数整除的数。在加密应用中起重要的位置,比如广为人知的RSA算法中,就是基于大整数的因式分解难题,寻找两个超大的素数然后相乘作为密钥的。一个比较常见的求素数的办法是埃拉托斯特尼筛法(the Sieve of Eratosthenes),说简单一点就是画表格,然后删表格,如图所示:
从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
#include
#include
using namespace std;
#define N 100000001
unsigned f[(N>>5)+5];
int p[5761456],m;
void init()
{
int i,j;
for(i=4;i
f[i>>5]|=1<<(i&0x1F);
p[m++]=2;
for(i=3;i*i
if(!(f[i>>5]&(1<<(i&0x1F))))
{
p[m++]=i;
for(j=i*i;j
f[j>>5]|=1<<(j&0x1F);
}
for(;i
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 ",m,clock()-st);
/*while(~scanf("%d",&n),n)
{
i=lower_bound(p,p+m,n+1)-p;
printf("%d ",i);
}*/
return 0;
}
以上就是本文的全部内容,希望对大家的学习有所帮助,也希望大家多多支持脚本之家。
python使用集合实现筛选法求素数-python素数筛选法浅析相关推荐
- c语言筛选法求100之内素数,用筛选法求100之内的素数。
/*用筛选法求给出范围之内的素数.*/ #include #include int isPrime(int n){//判断素数 for(int i = 2; i < n; i ++){ if(n ...
- 【计算流体力学】Python实现加权余量法求微分方程数值解 比较伽辽金法(Galerkin法)、最小二乘法和矩法的求解精度 分析误差随n增大的变化情况
一.简介 微分方程的经典数值方法有Ritz法.加权余量法.有限元法.有限体积法等等.对于微分方程的边值问题,如果能找到与微分方程相对应的泛函,可以通过求取相应泛函的极小值,将微分方程转化为关于基函数待 ...
- 用欧几里得算法求最大公约数_欧几里得算法:GCD(最大公约数),用C ++和Java示例解释...
用欧几里得算法求最大公约数 For this topic you must know about Greatest Common Divisor (GCD) and the MOD operation ...
- python利用集合的无重复性_利用Python程序完成ABAQUS中的一些重复性操作
对于建立过大型的ABAQUS分析模型的同学们,肯定会面对繁琐的重复性建模过程.比如对上百个点施加集中力,比如建立几百个耦合约束.这些重复性操作,如果使用cae操作的话,不仅费时,而且极其容易出错. 利 ...
- c 语言 用矩形法求定积分,如何用矩形法(梯形法)求定积分
分析: 高中的时候,我们学习过,可以通过矩形法或者矩形法来求定积分. 思路就是将积分区间划分成n等份,然后将这n等份近似看成矩形(或梯形),然后对所有的矩形(或梯形)的面积进行求和. 简单的例子: 求 ...
- java筛选法求质数_用筛选法找出N内所有素数(质数)
//筛选法实现 // main.c #include #define ArraySize 100 int main(int argc, const char * argv[]) { int i ; i ...
- python已知两条直角边求斜边,Python实现“已知三角形两个直角边,求斜边”
用Python实现"已知三角形两个直角边,求斜边" 要求:用户输入两个直角边(数值为浮点类型),若非浮点类型,则提示用户,继续输入. 思路:伪代码描述下步骤 1.-input a ...
- c语言写一个用矩形法求,写一个用矩形法求定积分的通用函数
解: #include main() { float integral(float(8p)(float),float a,float b,int n); float a1,b1,a2,b2,a3,b3 ...
- python的集合和字典
python的集合和字典基础部分 因为python里面集合set以及字典dict既有联系,又相互区别,所以放在一起谈谈基础的部分 文章目录 python的集合和字典基础部分 前言 一.集合字典的定义 ...
- python中集合运算_python集合运算
广告关闭 腾讯云11.11云上盛惠 ,精选热门产品助力上云,云服务器首年88元起,买的越多返的越多,最高返5000元! 每个韦恩(venn)图中的红色部分是给定集合运算得到的结果.? python 集 ...
最新文章
- windows下设置计划任务自动执行PHP脚本
- chorme 下载文件 保留 放弃_西部数码使用指南:保留数据重装以后sqlserver数据库不能启动(sql2008,sql2012)...
- 【小白学习PyTorch教程】九、基于Pytorch训练第一个RNN模型
- 容量法和库仑法的异同点_快速搞懂「活性污泥法」,主要把握这4点!
- 【Ubuntu-Tensorflow】程序结束掉GPU显存没有释放的问题
- 别琢磨了,企业高效灵活运作的秘密拿走:企业邮箱5折起!分享会场抽取苹果手机和猫超卡!
- SAP 采购流程和销售流程
- 初识数据库 1006
- Egret 之 消除游戏 开发 PART 6 Egret elimination game development PART 6
- mysql字符串拼接返回
- 局域网QQ助手 扫描在线QQ号
- 国产高人气无线蓝牙耳机测评,抖音火爆的蓝牙耳机值得入手吗?
- 分布式服务架构-第三章 服务化系统容量评估和性能保障
- OpenHarmony 3.2 Release HDF的IDL文件初探(上)
- python结束子进程_如何清除python中的子进程
- 如何将mp4转换成gif?教你一招实现视频转gif高清转换
- Kafka工具--滴滴Logi-KafkaManager
- Unix/Linux头文件说明
- 开店三大分享:咖啡店装修风格怎么确定?门店氛围如何打造?
- 【java】Java中-> 是什么意思?