埃式筛法尽管不错,但是确实做了许多无用功,某个数可能会被重复的筛好几次,欧拉筛解决了这个方法,下面为代码:
注意理解if(i%prim[j]==0) break;
大佬讲的不错的博客,我就不做复读机了。
点我传送
void ispirm(){ int cnt=0; memset(visited,true,sizeof(visited)); for(int i=2;i<=maxn;i++){ if(visited[i]) prim[cnt++]=i; for(int j=0;j<cnt&&prim[j]*i<=maxn;j++){ visited[prim[j]*i]=false; if(i%prim[j]==0) break; } } }
学会了欧拉筛法的话,下面有个经典例题,查询某段区间内质数的个数
题目链接:点我传送
这里用前缀和处理一下就可以再O(1)的时间复杂度情况查出了。
/*Keep on going Never give up*/ #pragma GCC optimize(3,"Ofast","inline") #include <bits/stdc++.h> const int maxn = 1e7; const int MaxN = 0x3f3f3f3f; const int MinN = 0xc0c0c00c; typedef long long ll; const int mod = 100000000; using namespace std; bool visited[maxn+10]; int prim[1000000+5]; int ans[maxn+10]; void ispirm(){ int cnt=0; memset(visited,true,sizeof(visited)); for(int i=2;i<=maxn;i++){ if(visited[i]) prim[cnt++]=i; for(int j=0;j<cnt&&prim[j]*i<=maxn;j++){ visited[prim[j]*i]=false; if(i%prim[j]==0) break; } } } int main() { ispirm(); int t; cin>>t; for(int i=2;i<=maxn;i++){ ans[i]=ans[i-1]+visited[i]; } while(t--){ int x,y; scanf("%d%d",&x,&y); printf("%d\n",ans[y]-ans[x-1]); } return 0; }