#include <iostream>
#include<math.h>
#include<stdio.h>
using namespace std;/*取得前n个素数*/
void get_primes(int *primeArray,int n){
primeArray[0]=2;
primeArray[1]=3;
int count=1;
int i;
int j; for(i=3;count<n;i+=2){
int sqrt_i=sqrt(i);
for( j=1;primeArray[j]<=sqrt_i && i%primeArray[j]!=0; j++); if(primeArray[j]>sqrt_i){
primeArray[count]=i;
count++;
} }
}
/*输出前n个素数*/
void display(int* array,int length){ for(int i=0;i<length;i++){
if(i!=0&&i%8==0){
cout<<endl;
}
printf("% 4d",array[i]); }
} /*初始化筛选数组
除了2以外的偶数不可能是素数,所以为false。*/
void initArray(bool * array,int length){
int i;
for( i=2;i<length;i+=2){
array[i]=false;
array[i-1]=true;
} if(i==length){
array[length-1]=true;
}
array[2]=true;
} /*输出小于等于n的素数。*/
void primesLessThanN(bool* primeArray,int n){int length=n+1;
initArray(primeArray,length);
int sqrt_n=sqrt(n);
/*获取已知每个素数,然后将
每个素数的奇数倍设置为合数,即元素值为false.
(因为偶数倍在初始化的时候已经设置为false 了。)
*/
for(int prime=3;prime<=sqrt_n;prime+=2){ if(primeArray[prime]==true){
for(int j=3*prime;j<n;j+=2*prime)
{
primeArray[j]=false;
} }
}
}
/*输出每个素数,元素值为true的元素的下标就是素数。*/
void displayPrimesLessThanN(bool* primesArray,int length){
cout<<" "<<2;
int count=1; for(int i=3;i<length;i+=2){
if(primesArray[i]==true)
{
if(count%8==0){
cout<<endl;
}
printf("% 4d",i);
count++;
} }
} int main()
{
cout<<"get the former n primes "<<endl;
int n=101;
int * primeArray=new int[n];
get_primes(primeArray,n);
display(primeArray,n); cout<<endl<<"get the primes less than N"<<endl;
bool* primesArray=new bool[n+1];
primesLessThanN(primesArray,n);
displayPrimesLessThanN(primesArray,n+1);
cout << endl<<"Hello world!" << endl;
return 0;
}