导读:
筛法的应用
一.用筛法求指定区间上的所有素数,并统计素数的个数.
筛法: 这是一种效率较高的算法.
求素数的筛法是公元前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
筛法的应用
一.用筛法求指定区间上的所有素数,并统计素数的个数.
筛法: 这是一种效率较高的算法.
求素数的筛法是公元前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