RT,筛的过程貌似出了问题,其中vis是bitset类,p是int类型。解决后立删。
#include<cstdio>
#include<bitset>
#include<algorithm>
using namespace std;
const int M=100000000;
bitset<M+1> vis;
//……
void p_prime(){
for(int i=2;i<=M;i++){
if(!vis[i]) p[++cnt]=i;
for(int j=1;i*p[j]<=M;j++){
vis[i*p[j]]=true;
if(i%p[j]==0) break;
}
}
}