求质数表 素数筛求素数 效率比较 数论相关
2012-05-25 16:24
302 查看
转自:求质数表 素数筛求素数 效率比较 数论相关
第一个是平时最常用的基本写法(还有一种写法是for(j=i*2;j<N;j++) not_prime[i*j]=1; )
第二个是以i*i为起点开始筛的,要注意浮点溢出!
第三个不太好简洁地说清楚,详细请点击http://blog.csdn.net/dinosoft/article/details/5829550
程序中途分别输出每次计算了的prime[500000],从而为了验证计算结果的正确性的
三次时间分别为1.85691 1.30005 0.782549
第一个是平时最常用的基本写法(还有一种写法是for(j=i*2;j<N;j++) not_prime[i*j]=1; )
第二个是以i*i为起点开始筛的,要注意浮点溢出!
第三个不太好简洁地说清楚,详细请点击http://blog.csdn.net/dinosoft/article/details/5829550
程序中途分别输出每次计算了的prime[500000],从而为了验证计算结果的正确性的
#include<iostream> #include<time.h> using namespace std; const int N= 40000000; int prime ,cnt; bool not_prime ={1,1}; int main() { cin.sync_with_stdio(false); clock_t a,b,c,d,e,f; memset(prime,0,sizeof(prime)); memset(not_prime,0,sizeof(prime)); a=clock(); cnt=0; for(int i=2;i<N;i++) if(!not_prime[i]){ prime[cnt++]=i; for(int j=i+i;j<N;j+=i) not_prime[j]=1; } b=clock(); cout<<"第500000个素数:"<<prime[500000]<<endl; memset(prime,0,sizeof(prime)); memset(not_prime,0,sizeof(prime)); c=clock(); cnt=0; for(int i=2;i<N;i++) if(!not_prime[i]){ prime[cnt++]=i; if(i<10000) //防止浮点溢出 for(int j=i*i;j<N;j+=i) not_prime[j]=1; } d=clock(); cout<<"第500000个素数:"<<prime[500000]<<endl; memset(prime,0,sizeof(prime)); memset(not_prime,0,sizeof(prime)); e=clock(); cnt=0; for(int i=2;i<N;i++) { if(!not_prime[i]) prime[cnt++]=i; for(int j=0;j<cnt&&i*prime[j]<N;j++) { not_prime[i*prime[j]]=1; if(!(i%prime[j])) break; } } cout<<"第500000个素数:"<<prime[500000]<<endl; f=clock(); cout<<"time1:"<<double(b-a)/CLOCKS_PER_SEC<<endl; cout<<"time2:"<<double(d-c)/CLOCKS_PER_SEC<<endl; cout<<"time3:"<<double(f-e)/CLOCKS_PER_SEC<<endl; cout<<"speed:"<<CLOCKS_PER_SEC<<endl; return 0; }
三次时间分别为1.85691 1.30005 0.782549
相关文章推荐
- 求质数表 素数筛求素数 效率比较 数论相关
- hdu 2582(数论相关定理+素数筛选+整数分解)
- hihoCoder 九十二周 数论一·Miller-Rabin质数测试 (数论 筛素数)
- 求两个数字的最大公约数-Python实现,三种方法效率比较,包含质数打印质数的方法
- 数论基础_素数相关
- 质数的无穷性——从素数到数论
- 编程之美之质数相关--二分图最大匹配+素数筛选
- 【数论】【素数】素数相关基础——基础操作
- 小测C#与C++、Python、IronPython求100000以内素数效率之比较
- 【数论】【素数】素数相关基础——获得素数与判别素数
- 【数论】【素数】素数相关基础——欧拉函数与欧拉定理
- [数论]素数相关整理
- 【数论】【素数】素数相关基础——整数分解
- 20140915 【 数论 - Miller_Rabin素数测试 】 51nod 1106 . 质数检测
- 大素数统计1e11以内 【数论】
- 关于for循环与for-each的效率比较
- 为什么求模运算要用素数(质数)
- 对300万一张表数据,用游标进行循环,不同写法的效率比较
- java 素数(质数)筛法
- 筛选法求质数(素数)