艾拉斯托尼素数筛选法

 #include<iostream.h>
#include<math.h>
int sieve(int prime[], int n, int k)
{ 
 int i;
 for(i=k+1;i<n; i++)
  if(prime[i] && prime[i]%prime[k]==0)
   prime[i]=0;
  k++;                //求下一个可用的筛下标
  while(prime[k]==0)
   k++;
  return k;
}
void erat_sieve(int prime[],int n)
{ 
 int i,m,k;
 prime[0]=2;
 n=n/2;
 for(i=1;i<n;i++)
  prime[i]=2*i+1;
 m=sqrt(n);
 k=1;
 while(prime[k]<m)
  k=sieve(prime,n,k); //用prime[k]作筛子，对数组prime作筛选
}
void main()
{
 int prime[10000];
 int i,n,cnt;
 cout<<"输入一个正整数作为区间最大值"<<endl;
 cin>>n;
 erat_sieve(prime,n);
 cnt=0;    //cnt用于统计prime总的素数个数
 for(i=0;i<n/2;i++)
  if(prime[i])
  {
   cnt++;
   if(cnt%8==0)
   {                  //每行输出8个素数
    cout<<prime[i]<<" ";
    cout<<endl;
   }
   else
    cout<<prime[i]<<" ";
  }
  cout<<endl;
  cout<<"1到"<<n<<"之间的素数个数="<<cnt;
  cout<<endl;
}