读书人

c语言求素数的有关问题

发布时间: 2012-04-27 11:57:44 作者: rapoo

c语言求素数的问题
用下面这个程序可以得出正确的结果:
#include <stdio.h>
#include <math.h>

int main(void)
{
int k;
int n, i;

printf("请输入您要求得的数:\n");
scanf("%d",&n);

k = sqrt(n + 1);

for (i = 2; i <= k; i++)
{
if (n%i == 0) break;
}
if (i>=k+1) printf("%d 是一个素数!\n",n);
else printf("%d 不是一个素数!\n",n);

return 0;
}

但是当k定义为double的时候,结果是错误的,比如输入17,就会显示不是素数,这是什么问题呢?


[解决办法]
看错了,这个问题是发生在计算sqrt和比较时的错误。
1、计算sqrt时,由于k是double,则k的值可能会含有小数,比如5.00001。
2、i > k后,不一定i>k+1。因为k是double,在比较时i被提升为double。比如5>5.00001,但是5<5.00001+1。
[解决办法]
如果n=17,则k=4.24
所以在出循环时i=5
则i<k+1所以显示为不是素数。

而4楼的数
如果n=17,则k=4.24
所以在出循环时i=5
则i>k所以显示为是素数。

所以i在等于17时不论k为sqrt(17)还是sqrt(18)民,n的值都在k与k+1之间。
[解决办法]
本来在判断素数的时候,从2判断到sqrt(n)就行了,你判断到了sqrt(n+1),然后最后在比较的时候,你又判断i≥k+1,两次都超出一点,最后完全有可能走到else部分。

C/C++ code
/* * 打印100以内的素数,每行显示5个。 */#include <math.h>#include <stdio.h>int isPrime(int i);int main(void){    int c = 0, i;    for (i=1; i<=100; i++)        if (isPrime(i))        {            printf("%3d", i);            c++;            if (c%5 == 0)                printf("\n");/*            if (c == 5)            {                printf("\n");                c = 0;            }*/        }    return 0;}int isPrime(int i){    int x;    int y = (int)ceil(sqrt(i));  /* int y = (int)(sqrt(i + 1)); */    if (i == 2)        return 1;    if (i%2==0 || i<2)  /* 过滤所有偶数及既不是素数又不是合数的情况 */        return 0;    for (x=3; x<=y; x++)        if (i%x == 0)            return 0;    return 1;} 

读书人网 >C++

热点推荐