实战演练 test16 T1

题目描述

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

      在小 X 的认知里，质数是除了本身和 1 以外，没有其他因数的数字。
      但由于小 X 对质数的热爱超乎寻常，所以小 X 同样喜欢那些虽然不是质数，但却是由两个质数相乘得来的数。
      于是，我们定义，一个数是小 X 喜欢的数，当且仅当其是一个质数，或是两个质数的乘积。
      而现在，小 X 想要知道，在 L 到 R 之间，有多少数是他喜欢的数呢？

输入

      第一行输入一个正整数 Q，表示询问的组数。
      接下来 Q 行，包含两个正整数 L 和 R，保证 L≤R。

输出

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

样例输入

1
1 6

样例输出

5

提示

      样例 1 解释：6 以内的质数有 2、3、5，而 4 = 2 * 2,6 = 2 * 3，因此，2,3,4,5,6 都是小 X 喜欢的数，而 1 不是。

      更多样例：下载 

       

题解（参考自xyx）：

【测试点 1-4、11-12】  考虑枚举 L 到 R 之间所有的数，判断这个数数是否是小 X 喜欢的数。  首先，如果数 A 是质数，显然它是小 X 喜欢的数，然后，我们考虑枚举每一 个 A 的因数 B，然后判断 B 和 A/B 是否是质数，如果是，那么 A 是小 X 喜欢的 数。 
由于一个数至多有 O（√N）个因数，因此， 
该算法的时间复杂度是 O（Q * N * √N * √N）的。 
【测试点 1-7、11-12】  那么，如何更快地判断一个数是否是小 X 喜欢的数呢？  我们可以采用分解质因数的做法，统计一个数的质因数个数，然后判断这个 数的质因数个数是否是 1 或 2 即可。 
如果分解质因数的过程采用 O（√N）的做法的话， 
该算法的时间复杂度是 O（Q * N * √N）的。 
【测试点 1-14】 考虑换一种思路，通过质数来生成小 X 喜欢的数。 我们可以通过线性筛的方式得出数据范围内所有的质数，同时，我们也得到 了每一个数的最小质因数。那么，判断一个数是否是小 X 喜欢的数只需将其最小 质因数除去后，判断剩下的部分是否为质数或 1 即可。 该算法的时间复杂度是 O（Q * N）的。 【测试点 1-7、11-14、17-18】  事实上，L 到 R 之间小 X 喜欢的数的个数，等于 1 到 R 之间小 X 喜欢的数的 个数减去 1 到 L - 1 之间小 X 喜欢的数的个数。因此，我们只需要知道 1 到某数 之间小 X 喜欢的数的个数，就可以 O（1）回答询问。  那么，我们采用分解质因数的做法在询问前预处理 1 到某数之间小 X 喜欢的 数的个数，然后再回答询问即可。 
该算法的时间复杂度是 O（N * √N + Q）的。 
【测试点 1-20】  有了上述思路，将分解质因数的做法替换为线性筛即可。  该算法的时间复杂度是 O（N + Q）的。 

#include<bits/stdc++.h>
#define maxn 10000005
using namespace std;
template <typename T> void read(T &x){
	x=0;int f=1;char ch=getchar();
	for(;!isdigit(ch);ch=getchar())if(ch=='-')f=-1;
	for(;isdigit(ch);ch=getchar())x=(x<<1)+(x<<3)+ch-'0';
	x*=f;
}

int f[maxn],prime[maxn],sum[maxn],tot; 

void work(){
	for(int i=2;i<=maxn;++i){
		if(f[i]==0){
			prime[++tot]=i;
			f[i]=i;
		}
		if(f[i]==i||f[i/f[i]]==i/f[i]) sum[i]=1;
		sum[i]+=sum[i-1];
		for(int j=1;j<=tot&&prime[j]<=f[i];++j){
			if(i*prime[j]>maxn) break;
			f[i*prime[j]]=prime[j];
		}
	}
}

int main(){
	work();
	int p,l,r;
	read(p);
	for(int ti=1;ti<=p;++ti){
		read(l),read(r);
		cout<<sum[r]-sum[l-1]<<endl;
	}
 	return 0;
}