最快素数打表,比欧拉筛快一倍。

版本一:1e7内比欧拉筛子快一倍,2e7持平,之后略慢 不论N多大整体计算次数,都是欧拉筛子的1/3,求大神解答1e8之后为什么变慢
思路:相较于欧拉筛考虑1-n所有数,这里基于孪生素数,只需要考虑孪生素数,n/3;
版本二:解决了1e8的问题,多加了一个sqrt(N)的处理,缺点代码更繁琐。


版本一代码:

点击查看代码
#include<bits/stdc++.h>
using namespace std;
const int N=1e7+9;
bool a[N+2];
int b[N+2];
int idx=0;
int w=0,s;
clock_t startTime;
int getCurrentTime() {
	return 1000.0 * (clock() - startTime) / CLOCKS_PER_SEC;
}
void init(int N){//孪生素数
	a[3]=a[2]=true; 
	for(int i=5;i<N;i+=6) a[i+2]=a[i]=true;
	for(int i=5;i<=N/5;i+=6){
		b[idx++]=i;
		b[idx++]=i+2;
	}
	for(int i=0;i<idx;i++)
		for(int j=0;j<=i&&b[j]*b[i]<N;j++)
			a[b[i]*b[j]]=false;
	
}
bool prime[N];
int p[N], tot;
void init2(int N){//欧拉 
    for(int i = 2; i < N; i ++) prime[i] = true;//w++;
    for(int i = 2; i < N; i++){
        if(prime[i]) p[tot ++] = i;
        for(int j = 0; j < tot && i * p[j] < N; j++){
            prime[i * p[j]] = false;// w++;
            if(i % p[j] == 0) break; //w++;
        }
    }
}
int main()
{
	int n;
	
	init(N);	
	cout<<"my   :"<<w<<endl;
	cerr << "Execution time: " << getCurrentTime() << " ms." << endl;	
	init2(N);
	cout<<"欧拉 :"<<w<<endl;
	cerr << "Execution time: " << getCurrentTime() << " ms." << endl;
	
	
	
	return 0;
 } 
版本二代码 更新1e8完全持平
点击查看代码
#include<bits/stdc++.h>
using namespace std;
const int N=1e8+9;
bool a[N+200000];
int b[N];
int c[N];
int idx=0;
int w=0,s;
bool prime[N];

int p[N], tot;
clock_t startTime;
int getCurrentTime() {
	return 1000.0 * (clock() - startTime) / CLOCKS_PER_SEC;
}
void init(int N){//孪生素数
	a[3]=a[2]=true; 
	for(int i=5;i<N;i+=6) a[i+2]=a[i]=true;
	for(int i=5;i<sqrt(N);i+=6){
		b[idx++]=i;
		b[idx++]=i+2;
	}
	for(int i=0;i<idx;i++)
		for(int j=0;j<=i;j++)
			a[b[i]*b[j]]=false;
	int w=idx;
	for(int i=b[idx-1]+4;i<N/5;i+=6){
		b[idx++]=i;
		b[idx++]=i+2;
	}
	for(int i = w; i < idx; i++){
        	for(int j=0;j<=i&&b[i]*b[j]<N;j++){
        			if(a[b[j]])
           			a[b[i]*b[j]]=false;// w++;
		}
    }
}

void init2(int N){//欧拉 
    for(int i = 2; i < N; i ++) prime[i] = true;//w++;
    for(int i = 2; i < N; i++){
        if(prime[i]) p[tot ++] = i;
        for(int j = 0; j < tot && i * p[j] < N; j++){
            prime[i * p[j]] = false;// w++;
            if(i % p[j] == 0) break; //w++;
        }
    }
}
int main()
{
	int n;
		init(N);	
	cout<<"my   :"<<w<<endl;
	cerr << "Execution time: " << getCurrentTime() << " ms." << endl;
	
	init2(N);
	cout<<"欧拉 :"<<w<<endl;
	cerr << "Execution time: " << getCurrentTime() << " ms." << endl;
	
	
	int ans=0;
	for(int i=1;i<N-9;i++)
	{
		if(a[i]) //cout<<i<<endl;
		ans++;
	}
	cout<<ans<<endl;
	cout<<endl;
	 ans=0;
	for(int i=1;i<N-9;i++)
	{
		if(prime[i]) //cout<<i<<endl,
		ans++;
	}
	cout<<ans<<endl;
	return 0;
 } 
posted @ 2023-02-20 21:52  xxj112  阅读(70)  评论(0)    收藏  举报