2

我正在从事 C 编程任务以在不使用 C 的平方根函数的情况下实现 Eratosthenes 的 Sieve。下面是我的输出和我的教授输出，我不确定我的代码中是什么导致它出错。有任何想法吗？

这是预期的输出

Program initiated
   1   2   3   5   7  11  13  17  19  23  29  31
  37  41  43  47  53  59  61  67  71  73  79  83
  89  97 101 103 107 109 113 127 131 137 139 149
 151 157 163 167 173 179 181 191 193 197 199 211
 223 227 229 233 239 241 251 257 263 269 271 277
 281 283 293 307 311 313 317 331 337 347 349 353
 359 367 373 379 383 389 397 401 409 419 421 431
 433 439 443 449 457 461 463 467 479 487 491 499
 503 509 521 523 541 547 557 563 569 571 577 587
 593 599 601 607 613 617 619 631 641 643 647 653
 659 661 673 677 683 691 701 709 719 727 733 739
 743 751 757 761 769 773 787 797 809 811 821 823
 827 829 839 853 857 859 863 877 881 883 887
Program terminated

这是我的输出： 

Program initiated
   1  37  41  43  47  53  59  61  67  71  73  79
  83  89  97 101 103 107 109 113 127 131 137 139
 149 151 157 163 167 173 179 181 191 193 197 199
 211 223 227 229 233 239 241 251 257 263 269 271
 277 281 283 293 307 311 313 317 331 337 347 349
 353 359 367 373 379 383 389 397 401 409 419 421
 431 433 439 443 449 457 461 463 467 479 487 491
 499 503 509 521 523 541 547 557 563 569 571 577
 587 593 599 601 607 613 617 619 631 641 643 647
 653 659 661 673 677 683 691 701 709 719 727 733
 739 743 751 757 761 769 773 787 797 809 811 821
 823 827 829 839 853 857 859 863 877 881 883 887
Program terminated

这是我的代码：

#include <stdio.h>

void zap(int data[], int divisor)
{
    for(int i=0;i<900;i++)
    {
        if(data[i]%divisor==0) // if mod is not 0, 0 out the index. 
        {
            data[i] = 0; 
        }
    }
}
// the display method 
void display(int data[])
{
    int count = 0; // init counter on the out side
    for(int i=0;i<900;i++)
    {
        if(data[i]>0)// don't print 0s
        {
            printf("%4d",data[i]);// print the data in a column 

            count++;// increment count 

            if(count==12) // print rows and columns 
            {
                count=0; // reset count
                printf("\n"); // print new line 
            }
        }
    }
    if(count<12)// we terminate loop and we now need print a new line 
    {
        printf("\n"); 
    }
}

int main()
{   
    // start the program, with a message
    printf("Program initiated\n"); 

    // needs to be 900 long 
    int primes[900];

    // populate the array 
    for(int i=1; i <= 900; i++)
    {
            primes[i] = i; 
    }

    // elminate the bad numbers 
    for(int i=2; i < 35; i++)
    {
        zap(primes,i); 
    }

    // display the array. 
    display(primes);

    // print the end message    
    printf("Program terminated\n"); 

    return 0; 
}
4

3 回答 3

3

您的zap函数将始终删除输入值。例如，当您zap使用除数 2调用时，即使 2 是素数，它也会检查2%2、查找和删除它。0

要解决此问题，您可以让它在divisor+1.

但是，我注意到它实际上根本没有做筛子。zap不需要做任何模数，只需按divisor. 仔细检查埃拉托色尼筛实际上是什么。

于 2012-03-19T05:42:05.667 回答
2

这不是真正的埃拉托色尼筛子。该算法的要点是%通过盲目（即不进行任何计算）排除除 2 之外的每 2 个数字，然后除 3 之外的每 3 个数字，然后除 4 之外的每 4 个数字，以此类推，完全避免测试可除性（用 ）。

您需要修复zap功能：首先，如果它等于 ，则不要删除数字divisor，并且不要检查余数，只需删除数字。

于 2012-03-19T08:30:08.693 回答
1

那么你可以使用这样的东西：

（将 sieve 初始化为足够大的布尔数组，每个条目都设置为 true - 因为我想保持简单设置sieve[0] = false; sieve[1] = false;

for(int i = 2; i < endOfNumbers; i++)
{
   if (sieve[i] == false) continue;
   for (int m = 2*i; m < endOfNumbers; m += i)
      sieve[m]=false;
}
于 2012-03-19T05:44:13.593 回答