C:使用阵列筛选Eratosthenes

时间:2012-03-19 05:34:49

标签: c sieve-of-eratosthenes

我正在进行C编程分配以实现Eratosthenes的Sieve而不使用C的平方根函数。下面是我的输出和我的教授输出,我不确定我的代码中是什么导致它出错。有任何想法吗?

这是预期的输出

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; 
}

3 个答案:

答案 0 :(得分:3)

您的zap函数将始终消除输入值。例如,当您使用2的除数调用zap时,它会检查2%2,找到0并将其删除,即使2是素数。

要解决此问题,您可以让它开始在divisor+1进行切换。

但是,我注意到它根本不是在做Sieve。 zap不需要执行任何模数,只需按divisor的步长运行数组。仔细检查Eratosthenes的Sieve究竟是什么。

答案 1 :(得分:2)

这不是真正的Eratosthenes筛子。这个算法的重点是避免测试可分性(%盲目(即没有任何计算)排除除2之外的每个第2个数字,然后每3个除3之外,然后每个4除了4等等。

您需要修复zap功能:首先,如果数字等于divisor,请不要删除数字,并且不检查余数,只需删除数字。

答案 2 :(得分:1)

你可以使用这样的东西:

(初始化筛子是一个足够大的bool数组,每个条目都设置为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;
}