求质数表 素数筛求素数 效率比较 数论相关
来源:互联网 发布:最会写文案的淘宝店 编辑:程序博客网 时间:2024/06/06 07:21
第一个是平时最常用的基本写法(还有一种写法是for(j=i*2;j<N;j++) not_prime[i*j]=1; )
第二个是以i*i为起点开始筛的,要注意浮点溢出!
第三个不太好简洁地说清楚
程序中途分别输出每次计算了的prime[500000],从而为了验证计算结果的正确性的
#include<iostream>#include<time.h>using namespace std;const int N= 40000000;int prime[N],cnt;bool not_prime[N]={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;}
下面是在XCODE下运行这段代码结果的截图,三次时间分别为1.85691 1.30005 0.782549 效率差距有点小大啊!
- 求质数表 素数筛求素数 效率比较 数论相关
- 求质数表 素数筛求素数 效率比较 数论相关
- 求质数(素数)
- 求素数(质数)的算法
- 求1-1000质数(素数)
- 求质数(素数)算法思想
- 求质数(素数)算法思想
- 求素数(质数)算法
- 求范围内质数(素数)
- 求1~1000的质数(素数)
- 筛选法求质数(素数)
- 求100到200的质数(素数)
- 求100内的所有素数(质数)
- python3练习题--求质数(素数)
- Python求1000以内质数(素数)
- 求质数(素数)算法及其改进
- 求质数(素数)算法和代码
- python 求100内的素数/质数
- AsyncTask的用法
- dip,px,pt,sp 的区别
- WPF DynamicDataDisplay 实现动态曲线图 双击获得数据点
- 2011中国云计算调研:海量数据引发IT架构变革
- 分享Silverlight/WPF/Windows Phone/HTML5一周学习导读(1月9日-1月15日)
- 求质数表 素数筛求素数 效率比较 数论相关
- 树的遍历
- OpenID Summary1 (Facebook)
- OpenId Summary2 (Google Account)
- First Version Image service for Mobile App(Odata.svc)
- lsmod, modprobe
- 怎么在网站中插入天气预报
- Asp.net上传图片同时生成文字水印图,图片水印图,缩略图
- SQL语句的基本操作