Question

大家好：我正在使用pow（）函数来计算prime的功率。但是有些输出没有按预期显示。例如，当我输入200作为输入时，我希望打印出所有素数及其不超过200的幂。大多数数字都很好，但我注意到一些不同于正确的数字1（24应该是25,120应该是121,124应该是125,168应该是169）。我在这里遗漏了什么？我的平台是win7运行Coding :: Blocks GNUGCC编译。有什么建议吗？非常感谢......

60
2 3 4 5 7 8 9 11 13 16 17 19 23 24 27 29 31 32 37 41 43 47 49 53 59 61 64 67 71 73 79 81 83 89 97 101 103 107 109 113 120 124 127 128 131 137 139 149 151 157 163 167 168 173 179 181 191 193 197 199

#include <stdio.h>
#include <stdlib.h>
#include <stdbool.h>
#include <string.h>
#include <math.h>

bool is_primes[1000]= {true};
int primes[200]= {0};
int answer[1000]= {0};

void sieve_of_Eratosthenes (bool ary[],int);
void swap(int*,int*);
void quick_sort(int array[], int, int);

int main(int argc,char* argv[]) {
    int n =0,i = 0,j = 0,cnt = 0,temp = 1;

    memset(is_primes,true,sizeof(is_primes));
    memset(primes,0,sizeof(primes));
    memset(answer,0,sizeof(answer));

    sieve_of_Eratosthenes(is_primes,1000);

   //build primes table
   for(i=0; i<1000; ++i) {
       if(is_primes[i]) {
           primes[j++] = i;
       }
   }

   while(scanf("%d",&n)==1) {
       i = 0,j = 0,cnt = 0,temp = 1;

       //primes
       for(j=0,cnt=0;primes[j]!=0; ++j) {
           if(primes[j]<=n) {
            answer[cnt] = primes[j];
            ++cnt;
           }
       }

       //and power of primes
       for(j=0; primes[j]!=0; ++j) {
           for(temp=2;pow(primes[j],temp)<=n;) {
               answer[cnt] = pow(primes[j],temp);
               ++temp;
               ++cnt;
           }
       }

       //sort ascending
       quick_sort(answer,0,cnt-1);

       printf("%d\n",cnt);

       for(j=0; j<cnt; ++j) {
           printf("%d ",answer[j]);
       }

       printf("\n");
    }

    return 0;
}

void sieve_of_Eratosthenes (bool ary[],int n) {
    int i = 0,j = 0;
    is_primes[0] = is_primes[1] = false;
    for(i=2; i*i<=n; ++i) {
        if(is_primes[i]) {
            for(j=i*i; j<=n; j+=i) {
                is_primes[j] = false;
            }
        }
    }
}

void swap(int* a, int* b) {
    int temp = *a;
    *a = *b;
    *b = temp;
}

// sort interval [left, right]
void quick_sort(int array[], int left, int right) {
    if (left < right) {
    // divide (partition)
        int pivot = array[(left+right)/2];
        int i = left - 1, j = right + 1;
        while (i < j) {
            do ++i;
            while (array[i] < pivot);
            do --j;
            while (array[j] > pivot);
            if (i < j) swap(&array[i], &array[j]);
    }

    // then conquer
    quick_sort(array, left, i-1);
    quick_sort(array, j+1, right);

    // no need to combine sub-solutions
    }
}

Answer 1

如果您不了解pow，您可能会以不同方式编写该循环：

for(temp = primes[j]*primes[j];
    temp <= n;
    temp *= primes[j]) {
  answer[cnt++] = temp;
}

这不仅可以避免圆整问题;它也会快得多。（有可能提出一个比temp更好的变量名。）

c中的pow（）函数会产生截断（或舍入）错误吗？

1 个答案: