筛选法的应用-淡淡的冰-搜狐博客

导读： 
           筛法的应用
一.用筛法求指定区间上的所有素数,并统计素数的个数.
   筛法: 这是一种效率较高的算法.
   求素数的筛法是公元前3世纪的厄拉多塞(Eratosthenes)提出来的:对于一个大整数X,只要知道不超过   sqrt(X)的所有素数,划去所有p的倍数2p,3p,...剩下的整数就是不超过X的全部素数.
   应用筛法求素数,为了方便实施"划去"操作,应设置数组.每一数组元素对应一个待判别的奇数,并赋初值0.如果该奇数为p的倍数则应划去,对应元素加一个划去标记,通常给该元素赋值-1.最后,打印元素值不是-1(即没有划去)的元素对应的奇数即所求素数.
   在实际应用筛法的过程中,p通常不限于取不超过sqrt(x)的素数,而是适当放宽取不超过sqrt(x)的奇数(从3开始),这样做尽管多了一些重复划去操作,但程序实现要简便些,而且只用考虑所取p的奇数倍(p的偶数倍肯定不是素数).
   应用筛法求素数的程序设计
   在指定区间[c,d](约定c为奇数)上所有奇数表示为 j=c+2k,(k=0,1,2,3,...e,这里e=(d-c)/2.).于是k=(j-c)/2是奇数j在数组中的序号(下标).如果j为奇数的倍数时,对应数组元素作划去标记,即a[(j-c)/2]=-1.最后,凡数组元素a[k]不等于-1,对应的奇数j=c+2k则为素数.
 
 /*筛法求指定区间上的素数程序*/
    #include
#include
using namespace std;
int main()
{
  long int c,d,i,j,k;
  static long int a[11000];
  int e,n;
  cout<<"求区间[c,d]上的素数"<
  cout<<"请输入c,d(c>2):"<
  cin>>c>>d;
  if(c%2==0)
   c++;
  e=(d-c)/2;
  i=1;
  while(i<=sqrt(d))  /*在[c,d]中筛选素数*/
  {
   i+=2;
   j=i;
   while(j<=d)
   {
    j+=2*i;  /*计算3*i,5*i,7*i.... */
    if(j>=c)
      a[(j-c)/2]=-1; /*筛去标记-1*/
   }
  }
   for(n=0,k=0;k<=e;k++)
    if(a[k]!=-1)
    {
     n++;
     cout< <<" ";
     if(n%10==0)
       cout<
    }   
  cout< <<"共有素数"< <

本文转自
http://yxdark.blog.sohu.com/83198402.html