eratosthenes筛不工作c++

Sieve of eratosthenes not working C++

本文关键字：c++ 工作 eratosthenes      更新时间：2023-10-16

所以我试图制作一个c++程序来生成所有素数，直到某一点，但由于某种原因，它将2之后的所有数字打印为非素数。

int A[1000000];
void sieve(int till)
{
    for(int i = 2; i < till; i++)
    {
        if(A[i] == 0)
        {
            for(int j = i*i; j < till; j+=i)
            {
                A[j - 1] = 1;
                //printf("%i is NOT primen", ij);
            }
        }
    }
}

和在主I do:

int N;
scanf("%i", &N);
sieve(N);

但是当我试着调试的时候它说:

NOT PRIME 0
NOT PRIME 1
NOT PRIME 2
PRIME 3
NOT PRIME 4
PRIME 5
NOT PRIME 6
PRIME 7
NOT PRIME 8
PRIME 9
NOT PRIME 10
PRIME 11
NOT PRIME 12
NOT PRIME 13

有谁知道我做错了什么吗?

非常简单。我们有几个简单的例子:

  1. 大于两个因数，所以不是素数
  2. 1只有一个因数，所以它不是素数
  3. 2是唯一为素数的偶数

所以，基本上我们不应该计算偶数。我们只对奇数感兴趣。

我们将所有非素数标记为1。

int A[1000000];
void sieve(int till)
{
    // mark all even numbers (except 2) as not prime
    for(int i = 4; i < till; i += 2)
    {
        A[i] = 1;
    }
    // now calculate not primes for only odd numbers
    for(int i = 3; i < till; i += 2 )
    {
        if(A[i] == 0)
        {
            for(int j = i*i; j < till; j+=(i+i))
            {
                A[j] = 1;
            }
        }
    }
}
int main()
{
    int N;
    scanf("%i", &N);
    sieve(N);
    for(int i = 2; i < N; i++)
    {
        if(A[i] == 1) printf("%d is not primen", i);
        else printf("%d is primen", i);
    }
    return 0;
}

这一行:

A[j - 1] = 1;
应:

A[j] = 1;

所以，现在应该运行良好:

#include<stdio.h>
int A[1000000];
void sieve(int till)
{
    for(int i = 2; i < till; i++)
    {
        if(A[i] == 0)
        {
            for(int j = i*i; j < till; j+=i)
            {
                A[j] = 1;
            }
        }
    }
}
int main()
{
    int N,i;
    scanf("%i", &N);
    for (i = 0;i < N;i++)
        A[i] = 0;
    sieve(N);
    for (i = 2;i < N;i++)
        if (A[i])
            printf("%i is not primen",i);
        else
            printf("%i is primen",i);
}

你的代码有很多问题。

  1. int A[1000000];

    没有初始化数组。您应该初始化A[0]=1, A[1]=1和所有其他索引为0。

  2. for(int j = i*i; j < till; j+=i)

    i*i可能导致整数溢出。您可能应该使用j=i+i来减少溢出的可能性。

  3. A[j - 1] = 1;

    不确定是否使用了正确的数组索引。

相关文章：