前n个素数 && 小于等于n的素数_legend



#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;
}