题意:给你数n和d,让你求小于n的数中,以d作为最大因子的数有多少个。
我们简单分析可得,这个数范围必定是 d + 1到n - 1。
本题刚开始就会被数据范围所限制,开1e9的数组显然不行,而这题作为一道练习题想到了线性筛。
但是不管怎样开1e9的数组是不合理的,所以就分析题意,果然,d作为数的最大真因子有 d * x = num(<n)
且x一定是个大于1的质数,即x一定是数的最小质因子,那么我们可以得到x的范围
即x <= d 且 x <= (n - 1) / d 即xmax = min(d, (n - 1) / d);
那么我们只需枚举小于等于xmax的素数即可。
我们通过分析得知,xmax最大与d相同,两者之积小于n,故我们只需求1e9开平方的范围内素数即可,1e5足够了。
但是在枚举过程需要注意的是,我们所枚举的2~xmax范围内的所有素数,如果d是当前素数的倍数,即d = pri[i] * x', 那么到下一个素数的时候,因为d' = pri[i + 1] * x' > d,d就不能保证是最大的了,所以当d可以被某个素数整除的时候, 就要退出循环了
比如我们求小于8的最大真因子为2的数有几个,这里即n = 16, d = 4; xmax = min(4, 3) = 3; pri[] = {2, 3};
枚举素数即 pri[0] * d = 8, 但是 d % pri[0] == 0所以要退出循环。
当枚举pri[1]时,pri[1] * d = 3 * 4 = 12,我们很容易得到12的最大真因子为6
而可以判断出的退出循环的原因就是:d = pri[0] * x',当枚举到pri[1]时,d * pri[1] = pri[0] * x' * pri[1],我们知道pri[1] > pri[0],所以d' = pri[1] * x' > d,所以这时候d就不是最大真因子了。那么当d可以被某个素数整除时,下一次循环d就必定不是最大真因子。
代码在此:
/*
由于这题是用最小质因子反推
最次的情况就是最小质因子和最大因子相同
所以只需枚举10w的素数即可
*/
#include<cstdio>
#include<algorithm>
using namespace std;
const int N = 100010;
int pri[N];
bool isp[N] = {true, true};
int cnt = 0;
void xs(){
for(int i = 2; i < N; i++){
if(!isp[i])
pri[cnt++] = i;
for(int j = 0; j < cnt && 1ll * i * pri[j] < N; j++){
isp[i * pri[j]] = true;
if(i % pri[j] == 0)
break;
}
}
}
int main()
{
xs();
int t;
scanf("%d", &t);
while(t--){
int n, d;
scanf("%d%d", &n, &d);
int ans = 0;
/*
思想:已知d,我们要求出从d+1到n-1之间以d为最大真因子的数
设该数为x 则有 x <= d 且 d * x <= (n - 1) 得:x <= d 且 x <= (n - 1) / d
即Xmax = min(d, (n-1)/d)
*/
int xmax = min(d, (n - 1) / d);
for(int i = 0; pri[i] <= xmax; i++){
ans++;
//如果d是当前素数的倍数,即d = pri[i] * x',
//那么到下一个素数的时候,因为d' = pri[i + 1] * x' > d
//d就不能保证是最大的了,所以当d可以被某个素数整除的时候,
//就要退出循环了
if(d % pri[i] == 0)
break;
}
printf("%d\n", ans);
}
return 0;
}