小x的质数(线性O(n)筛素数)

小x的质数

题目描述

小 X 是一位热爱数学的男孩子,在茫茫的数字中,他对质数更有一种独特的情感。小 X 认为,质数是一切自然数起源的地方。

在小 X 的认知里,质数是除了本身和 11 以外,没有其他因数的数字。

但由于小 X 对质数的热爱超乎寻常,所以小 X 同样喜欢那些虽然不是质数,但却是由两个质数相乘得来的数。

于是,我们定义,一个数是小 X 喜欢的数,当且仅当其是一个质数,或是两个质数的乘积。

而现在,小 X 想要知道,在 LL 到 RR 之间,有多少数是他喜欢的数呢?

输入格式

第一行输入一个正整数 QQ,表示询问的组数。

接下来 QQ 行。包含两个正整数 LL 和 RR。保证 L \le RL≤R。

输出格式

输出 QQ 行,每行一个整数,表示小 X 喜欢的数的个数。

数据范围与约定

样例

样例解释 1

66 以内的质数有 2,3,52,3,5,而 4=2 * 2,6 = 2 * 34=2∗2,6=2∗3。因此 2,3,4,5,62,3,4,5,6 都是小 X 喜 欢的数,而 1 不是。

样例输入1

1
1 6

样例输出1

5

样例输入2

10
282 491
31 178
645 856
227 367
267 487
474 697
219 468
582 792
315 612
249 307

样例输出2

97
78
92
65
102
98
114
90
133
29

样例输入3

10
20513 96703
15236 86198
23185 78205
40687 48854
42390 95450
63915 76000
36793 92543
35347 53901
44188 76922
82177 90900

样例输出3

24413
23001
17784
2669
16785
3833
17712
6028
10442
2734

code

 1 #include<cstdio>
 2 #include<cmath>
 3 const int MAXN = 10000100;
 4 bool lk[MAXN],noprime[MAXN];
 5 int prime[MAXN];
 6 int sum[MAXN];
 7
 8 int read() {
 9     int x = 0,f = 1;char ch = getchar();
10     for (; ch<‘0‘||ch>‘9‘; ch = getchar())
11         if (ch==‘-‘) f = -1;
12     for (; ch>=‘0‘&&ch<=‘9‘; ch = getchar())
13         x = x*10+ch-‘0‘;
14     return x*f;
15 }
16
17
18
19 int main() {
20
21     int tot = 0;
22
23     for (int i=2; i<=10000000; i++) {
24         if (!noprime[i]) {
25             prime[++tot]=i;
26             lk[i] = true;
27             for (int j=1; j<=tot && i*prime[j]<=10000000; j++) {
28                 noprime[i*prime[j]] = true;
29                 lk[i*prime[j]] = true;
30             }
31         }
32         else {
33             for (int j=1; j<=tot&&i*prime[j]<=10000000; j++) {
34                 noprime[i*prime[j]] = true;
35                 if(i % prime[j] == 0) break ;
36             }
37         }
38     }
39     for (int i=1; i<=10000000; ++i) {
40         if (lk[i]) sum[i] = sum[i-1]+1;
41         else sum[i] = sum[i-1];
42     }
43
44     int L,R,q = read();
45     while (q--) {
46         L = read();R = read();
47         printf("%d\n",sum[R]-sum[L-1]);
48     }
49
50     return 0;
51 }
时间: 2024-08-26 08:52:39

小x的质数(线性O(n)筛素数)的相关文章

常州模拟赛d2t1 小X的质数

题目背景 小 X 是一位热爱数学的男孩子,在茫茫的数字中,他对质数更有一种独特的 情感.小 X 认为,质数是一切自然数起源的地方. 题目描述 在小 X 的认知里,质数是除了本身和 1 以外,没有其他因数的数字. 但由于小 X 对质数的热爱超乎寻常,所以小 X 同样喜欢那些虽然不是质数, 但却是由两个质数相乘得来的数. 于是,我们定义,一个数是小 X 喜欢的数,当且仅当其是一个质数,或是两 个质数的乘积. 而现在,小 X 想要知道,在 L 到 R 之间,有多少数是他喜欢的数呢? 输入输出格式 输入

线性筛素数详细整理

如果你在1个月前让我判断素数,我一定会猛敲出以下代码: bool check( int num ) { int tmp =sqrt( num); for(int i= 2;i <=tmp; i++) if(num %i== 0) return 0 ; return 1 ; //实在是太慢了! } $ $ 下面给大家带来3种筛选素数和一种直接判断素数 $ $ $ $ 什么是线性筛? 对于求多个质数时与其一个个判断不如用排除法,用空间换取大量时间. $ $ $ $ $ $ 一般筛法(埃拉托斯特尼筛法

洛谷 P3383 【模板】线性筛素数

P3383 [模板]线性筛素数 题目描述 如题,给定一个范围N,你需要处理M个某数字是否为质数的询问(每个数字均在范围1-N内) 输入输出格式 输入格式: 第一行包含两个正整数N.M,分别表示查询的范围和查询的个数. 接下来M行每行包含一个不小于1且不大于N的整数,即询问概数是否为质数. 输出格式: 输出包含M行,每行为Yes或No,即依次为每一个询问的结果. 输入输出样例 输入样例#1: 100 5 2 3 4 91 97 输出样例#1: Yes Yes No No Yes 说明 时空限制:5

睡前数学一小时之线性筛素数:

睡前数学一小时之线性筛素数:1,朴素的筛素数算法:埃拉托斯特尼筛法.这是个简单再简单不过的一个素数的筛法.只是名字很拉风.这就告诉我们,往往东西不好这没什么,名字很拉风.别人也不会记住.hhhhh.这个的思路就是.每一个数都是由一个质数与和数(质数也可以)的积组成.这也是质数与和数的定义.而这个它这个筛发,就是当遇到一个质数的时候开始枚举,枚举[1,n]中间关于这个质数的倍数.每次都枚举,每次都将算出的这个数打上标记.而最后整个区间内的质数枚举完后,整个区间内的质数也就筛选出来了.这个很简单.时

【算法学习】线性筛素数

声明:本文所涉及部分内容可能并非原创.如发现本文侵犯了您的权益,可申请博主删除. 嘛……好久没有写博客了,今天无聊来写一篇~ 为什么要写这个呢?因为这个算法让我知道了我什么都不会啊…… 闲话少说,切入正题. 1. 筛素数 素数即质数,定义就不说啦 那么我们经典的筛素数的方法呢都很简单,开一个bool数组,从1到$\sqrt n$,如果是素数就把这个数的倍数给筛掉. 时间复杂度显然:$O(nlog_n)$ 代码也很简单: void prime(int n) { for(int i = 2; i *

线性筛素数(欧拉筛)+前缀和优化

关于素数的定义:在大于1的自然数中,除了1和它本身以外不再有其他因数. 判断一个数是否是素数: 1 int x; // 要求的数 2 for(int i=2;i<=sqrt(x);++i) 3 { 4 if(x%i==0) 5 { 6 cout << "这不是素数" << endl; 7 break; 8 } 9 } 埃氏筛法(时间复杂度:$O(NloglogN)$): 1 int num_prime = 0; // 素数的数量 2 int prime[5

Goldbach&#39;s Conjecture POJ - 2262 线性欧拉筛水题 哥德巴赫猜想

题意 哥德巴赫猜想:任一大于2的数都可以分为两个质数之和 给一个n 分成两个质数之和 线行筛打表即可 可以拿一个数组当桶标记一下a[i]  i这个数是不是素数  在线性筛后面加个装桶循环即可 #include<cstdio> #include<cstring> using namespace std; bool Is_Primes[1000005]; int Primes[1000005]; int cnt; void Prime(int n){ cnt=0; memset(Is_

[线性筛][筛素数/筛约数个数]

线性筛中每个数只会被最小的素因子筛一次 1)线性筛素数 1 void getprime(int siz){ 2 memset(isprime,1,sizeof(isprime)); 3 isprime[1]=0; 4 for(int i=2;i<=siz;i++){ 5 if(isprime[i])prime[++tot]=i; 6 for(int j=1;j<=tot&&i*prime[j]<=siz;j++){ 7 isprime[i*prime[j]]=0; 8 i

线性筛素数模板

传送门:线性筛素数 Prime: 1 #include<cstdio> 2 3 const int MAXN = 10000100; 4 int Prime[MAXN],n,m,Size; 5 bool Vis[MAXN]={1,1}; 6 7 int main() 8 { 9 scanf("%d%d",&n,&m); 10 for(int i=2;i<n;i++) 11 { 12 if(!Vis[i]) 13 Prime[++Size]=i; 14