我写了一个程序来打印所有素数到给定的数字

时间:2018-08-10 05:29:23

标签: c primes

#include<stdio.h>

float abso(float a)
{
    if(a<0.0)
        return(-1*a);
    else
        return a;
}

int sqert(int x)
{
    float x1 = (float)x;
    float g1, g2;
    g1 = x1/2.0;
    g2 = x1/g1;
    double e=0.000000001;
    int sr;
    while(abso(g1-g2)>=e)
    {
        g2 = (g1+x1/g1)/2.0;
        g1 = x1/g2;
    }
    sr = (int)g1;
    return(sr);
}

int main()
{
    int num;
    //num = 0;
    printf("Enter a num between 1 and 1000 ");
    scanf("%d",&num);
    //printf("hello");
    int flag;
    flag = 0;
    int i;
    i = 2;
    int m, j;
    m = j = 0;

    while(i<=num)
    {
        flag = 1;
        if(i==2)
        {
            printf("%d ",i);
            i++;
        }
        else if(i==3)
        {
            printf("%d ",i);
            i= i+2;
        }
        else if(i>3)
        {
            m = sqert(i);
            for(j=2;j<=m;j++)
            {
                if((i%j) == 0)
                {
                    flag = 0;
                    break;
                }
            }

            if(flag == 0)
            {
                i = i + 2;
                continue;
            }    

            if(j==m+1)
                printf("%d ",i);

            i = i + 2;
        }
    }
    printf("\n");
    return(0);
}

这里的sqert函数取输入值的平方根,abso函数取给定值的绝对值。

使用的逻辑很简单,我们迭代到该数字的平方根,然后检查每个数字是否可被除自身以外的其他任何数字整除,那么它不是质数,否则它是质数。

但是该程序仅适用于5以下的值。此后,它不打印任何内容。它也不打印输入的号码。

任何帮助将不胜感激。

1 个答案:

答案 0 :(得分:1)

在您的代码中,您尝试将floatdouble进行比较。

while(abso(g1-g2)>=e)

这是不正确的,因为floatdouble之间的精度有所不同。

double的精度是float的2倍。

float是32位IEEE 754单精度浮点数
(符号为1位,指数为8位,值为23 *),
即float的精度为7位小数。

double是64位IEEE 754双精度浮点数
(1位代表符号,11位代表指数,52 *位代表值),
即double的精度为15位小数。

e声明为float,并确保为其分配的值适合float的精度。
像这样:

float e = 0.00001;

然后您的代码将按预期工作。

@EricPostpischil在评论中说,sqert函数可能并不总是返回正确的近似值。这是对的。 因此最好使用sqrt中的math.h,而不是sqert