用筛选法求质数的一个小程序

microzhao 2004-09-11 04:56:53

我碰到一个用筛选法求质数的题，程序如下：
#include <stdio.h>
#define N 50
#define LN 16

main()
{
unsigned int sieve[N],primes[N];
unsigned int j,w,p,c;
for(j=0;j<N;j++)
{
sieve[j]=0xFFFF;//筛子
primes[j]=0x00;//质数表
}
w=0;
j=0;
do{
while(((0x01<<j++))&&sieve[w]==0x00);
p=w*LN+j-1;
c=2*(p+1)+1;
primes[w]|=0x01<<(j-1);
do
{
sieve[p/LN]&=(~(0x01<<p%LN));//p％LN位置0
p+=c;
}while(p<N*LN-LN);
while((sieve[w]==0x00)&&(w<N-1))
{
w++;
j=0;
}

}while(sieve[w]);
printf("%5d",2);
for(w=0;w<N;w++){
for(j=0;j<LN;j++)
if((0x01<<j)&&primes[w])
printf("%10d",2*(w*LN+j+1)+1);
}
printf("\n");
}
可是结果不对，显示的是奇数，我看了半天也不知道怎么改，哪位帮忙调试调试，谢谢了

...全文

179 3 打赏收藏转发到动态举报

写回复

用AI写文章

3 条回复

切换为时间正序

请发表友善的回复…

发表回复

microzhao 2004-09-11

打赏
举报

非常感谢

shw014 2004-09-11

打赏
举报

....
while(((0x01<<j++)&sieve[w])==0x00);
应该用按位与& ,&优先级低与==
....
while(p<N*LN);
你最后一个sieve[49]漏了
....
if((0x01<<j)&primes[w])
应该用按位与&

shw014 2004-09-11

打赏
举报

#include <stdio.h>
#define N 50
#define LN 16

main()
{
unsigned int sieve[N],primes[N];
unsigned int j,w,p,c;
for(j=0;j<N;j++)
{
sieve[j]=0xFFFF;//筛子
primes[j]=0x00;//质数表
}
w=0;
j=0;
do{
while(((0x01<<j++)&sieve[w])==0x00);
// ^^ ^
p=w*LN+j-1;
c=2*(p+1)+1;
primes[w]|=0x01<<(j-1);
do
{
sieve[p/LN]&=(~(0x01<<p%LN));//p％LN位置0
p+=c;
}while(p<N*LN);
// ^^^^
while((sieve[w]==0x00)&&(w<N-1))
{
w++;
j=0;
}

}while(sieve[w]);
printf("%5d",2);
for(w=0;w<N;w++){
for(j=0;j<LN;j++)
if((0x01<<j)&primes[w])
// ^^
printf("%10d",2*(w*LN+j+1)+1);
}
printf("\n");
}

Matlab课程老师提出的计算素数比赛，在30秒内看谁能用机房里的电脑计算出最多的素数。兴趣使然，自己结合基本的素数定理和筛选法，成功编写这个小程序。目的无他，贴出以供交流和共同进步。

主要为大家详细介绍了python使用筛选法计算小于给定数字的所有素数，具有一定的参考价值，感兴趣的小伙伴们可以参考一下

课程概述本季度为uni-app实战项目第2季度，将实战商城类app和小程序，其中会包括发布到安卓端app，iOS端app、微信小程序，另外还会涉及Weex（nvue），Vuex模块化管理，商城支付模块等。课程大纲见以下图！

用python编写代码找出...prime=[] #用一个列表来存储素数 yes=True print(功能：找出所有小于等于1000的素数) for i in range(2,1000+1): #1不是素数，range“前闭后开” yes=True for j in range(2,i): if i%j==0:

我对其进行了革命性的数据结构改进，空间复杂度从2个O(n)降低到1/6个O(n)，程序的算法描述更加简洁，改用C++实现，我认为算法效率已经达到了素数算法的极限。注：创建的内存大小不要超过内存，否则效率下降