传送门
题意大意是:给定的 n, d, 求不超过 n 且最大因子是 d 的数有多少个
则 m*d = n 中 d 是最大因子
显然如果 m 是非素数的话,那么一定可以分解出一个数,这个数会和 d 构成更大的 n 的因子。
此外,如果 d 是非素数的话,那么可能 d 分解出的数乘上某个小于 d 的素数的乘积比 d 本身要大。所以这题中,如果第一次发现 d 是某个素数的整数倍,那么就应当到此为止。为什么呢?
我们知道,一个数可以分解成素数的乘积,即:
d = p1^k1*p2^k2*p3^k3*p4^k4*p5^k5...
它乘上某个素数 pm 的结果是
m = pm*p1^k1*p2^k2*p3^k3*p4^k4*p5^k5...
它的最大的非自身的因数是
p1^k1*p2^k2*p3^k3*p4^k4*p5^k5...
那你如果找到 d 的第一个素数因子 p1 没有结束的话,枚举下一个素数 pt 乘上去,那么新的数就是
m = p1^k1*pt*p2^k2*p3^k3*p4^k4*p5^k5...
它的最大的非自身的因数是
p1^(k1-1)*pt*p2^k2*p3^k3*p4^k4*p5^k5...
就不再是 d 了,所以所枚举的素数最大最大是 d 的最小因子
所以先打个素数表然后一个个找就可以了
#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e5+5;
bool not_prime[maxn];
int prime[maxn];
void init() {
memset(not_prime, false,sizeof(not_prime));
int k = 0;
for (int i = 2; i < maxn; ++i) {
if (!not_prime[i]) {
prime[k++] = i;;
}
for (int j = 0; j<k && i*prime[j] < maxn; ++j) {
not_prime[i*prime[j]] = true;
if (i%prime[j] == 0)
break;
}
}
}
int main() {
std::ios::sync_with_stdio(false);
init();
int T;
cin >> T;
int n, d;
while (T--) {
cin >> n >> d;
int ans = 0;
for (int i = 0;;++i) {
if (d*prime[i] >= n)
break;
ans++;
if (d%prime[i] == 0)
break;
}
cout << ans << endl;
}
return 0;
}