为什么我的代码没有正确生成素数

Question

小于10,000,000的素数是664,579，但我的代码只生成664,214。数字来源为https://primes.utm.edu/howmany.html

#include <iostream>
#include <bitset>
#include <vector>
using namespace std;

const int N = 10000001;
bitset<N>num;
vector<int>prime;
inline void sieve()
{
    num.flip();
    num[0] = num[1] = 0;
    for(int i=2;i<N;i++)
        if(num[i])
        {
            prime.push_back(i);
            for(long long unsigned j=i*i; j<N;j+=i)
                num[j] = 0;
        }
}

int main() {
    sieve();
    cout << prime.size() << endl;
    return 0;
}

Answer 1

计算i*i时会出现整数溢出。然后，您将结果分配给long long的事实不会使编译器在乘法之前提升类型。

如果我将i声明为long long unsigned int，那么您的程序会输出664579