随笔 - 531  文章 - 0  评论 - 3  阅读 - 10215 

给定区间 [L,R]1≤R<(1<<30)  R−L≤1e6 ),请计算区间中素数的个数。

筛出 sqrt(R) 的质数p, 遍历 L~R 的数,看能否被p 约分,也就是合数,打个标记

 

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
#include <iostream>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std ;
  const int M=1e6+30;
 #define int long long
 int vis[M];
 int b[M],prime[M],tot,n;
  
 void init(int top){
     memset(b,1,sizeof b);
     b[1]=0;
      
     int i,j;
     for(i=2;i<=top;i++){
         if(b[i]) prime[++tot]=i;
          
        for(j=1;j<=tot&&i*prime[j]<=top;j++){
             b[i*prime[j]]=0;
             if(i%prime[j]==0) break;
        }
     }
 }
 signed main(){
    int L,R;
    cin>>L>>R; L= L==1?2:L;
    init(1e6);
    for(int i=1;i<=tot;i++){
        int t=prime[i];
        int k=(L+t-1)/t*t>2*t?(L+t-1)/t*t:2*t;
         
        for(int j=k;j<=R;j+=t) vis[j-L+1]=1;
    }
    int cnt=0;
    for(int i=L;i<=R;i++) if(vis[i-L+1]==0) cnt++;
    cout<<cnt;
 }

 

posted on   towboat  阅读(26)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 阿里最新开源QwQ-32B,效果媲美deepseek-r1满血版,部署成本又又又降低了!
· 单线程的Redis速度为什么快?
· SQL Server 2025 AI相关能力初探
· AI编程工具终极对决:字节Trae VS Cursor,谁才是开发者新宠?
· 展开说说关于C#中ORM框架的用法!
点击右上角即可分享
微信分享提示