c++ - 埃拉托色尼筛法 C++

Question

目前我正在做一个项目，我想计算所有素数。当我编译 (MINGW Windows Comp.) 时，程序崩溃并返回一个随机错误号。这是我写的代码:

http://pastebin.com/4vVnAM2v

/*
    Sieb des Eratosthenes
*/

#include <iostream>
#include <math.h>
using namespace std;

main()
{
    //variablendeklaration
    unsigned int nmax=100;
    unsigned int i,j,erg;
    bool prim[nmax];

    //Initialisieren
    prim[0]=false;
    prim[1]=false;

    //array prim[i] erstellen
    for(i=2;i<nmax;i++)
    {
        prim[i]=true;
    }




    for(i=2;i<nmax;i++) //alle Primzahlen durchlaufen
    {
        if(prim[i] == true) //auf Prim prüfen
        {
            for(j=2;j<nmax;j++) //multiplizieren und wegstreichen
            {
                erg = j * i;
                prim[erg] = false;
            }
        }
    }

    for(i=2;i<nmax;i++)
    {
        cout << prim[i] << endl;
    }


}

Answer 1

此时:

            erg = j * i;
            prim[erg] = false;

您最终将访问 prim 的范围之外，因为 i 和 j 的值最多为 nmax - 1，因此 erg 的最大值为 (nmax - 1) * (nmax - 1)。您需要检查此条件并在 erg >= nmax 时中断，例如

            erg = j * i;
            if (erg < nmax)          // if i * j still within bounds
                prim[erg] = false;   // set prim[i * j] to false
            else                     // else i * j now >= nmax
                break;               // break out of loop and try next i value