素数打印函数

/*a[]： 用于返回素数的数组 n： 产生n以内的素数，按升序放入a[]中 返回值： n以内素数的个数*/#include<stdio.h>#include<stdlib.h>#include<iostream>using namespace std;const int MAX=250000;int prime(int a[],int n){ int i,j,k,x,num,*b; n++; n/=2; b=(int *)malloc(sizeof(int)*(n+1)*2); a[0]=2;a[1]=3;num=2; for(i=1;i<=2*n;i++) b[i]=0; for(i=3;i<=n;i+=3)//这边可能不太好理解你可以写出来跟踪下 for(j=0;j<2;j++){ x=2*(i+j)-1; if(b[x]==0){ a[num++]=x; for(k=x;k<=2*n;k+=x)//把x的倍数标记掉 b[k]=1;}}return num;}int main(){int n,a[MAX],i;cout<<"输入n"<<endl;cin>>n;int count=prime(a,n);for(i=0;i<count;i++)cout<<a[i]<<" ";cout<<endl<<"总共"<<count<<"个素数"<<endl;return 0;}

 

上面的方法其实其步长为6，这个规律可以多写一些数出来，然后首先把偶数和3的倍数

排去。然后i从6开始，步长为6。代码如下：

 

int prime(int a[],int n){ int i,j,k,x,num,*b; n++; n/=2; b=(int *)malloc(sizeof(int)*(n+1)*2); a[0]=2;a[1]=3;num=2; for(i=1;i<=2*n;i++) b[i]=0; for(i=6;i<=n;i+=6) { x=i-1; while(b[x]==0) { a[num++]=x; for(k=x;k<=2*n;k+=x) b[k]=1; }x=i+1; while(b[x]==0) { a[num++]=x; for(k=x;k<=2*n;k+=x) b[k]=1; } } return num;}