c – 200万以下所有素数的总和

我做了一个程序,返回200万以下所有素数的总和.我真的不知道这个是怎么回事,当得到正确的答案是142913828922时,我得到了142891895587作为我的答案.看起来它缺少了几个素数.我很确定getPrime函数可以正常工作.我曾经使用它几次并且工作正常.代码如下：

vector<int> getPrimes(int number);

int main()
{

    unsigned long int sum = 0;
    vector<int> primes = getPrimes(2000000);

    for(int i = 0; i < primes.size(); i++)
    {
        sum += primes[i];
    }

    cout << sum;

    return 0;
}


vector<int> getPrimes(int number)
{

    vector<bool> sieve(number+1,false);
    vector<int> primes;
    sieve[0] = true;
    sieve[1] = true;

    for(int i = 2; i <= number; i++)
    {
        if(sieve[i]==false)
        {
            primes.push_back(i);
            unsigned long int temp = i*i;
            while(temp <= number)
            {
                sieve[temp] = true;
                temp = temp + i;
            }
        }
    }
    return primes;
}

解决方法

表达式i * i溢出,因为我是一个int.它在被分配给temp之前被截断.为避免溢出,请将其强制转换为：static_cast< unsigned long>(i)* i.

更好的是,在该条件发生之前终止迭代：for(int i = 2; i * i< = number; i). 测试固定. 顺便说一句,你有点(不)幸运的是,这不会产生额外的素数以及缺少一些素数：int值是有符号的,并且在溢出时可能是负数,并且通过我对§4.7/ 2的读数,这将导致要跳过的内循环.