题目要求小于n的并且最大因子为d的数的个数。
我们设m=kd,k是素数我们才能保证m的最大因子是d。
因为若k不是素数,我们就可以将k分解成ab,那么m的最大因子就可以是ad或者bd,反正不是d了。
所以我们考虑一下线性筛。
但是我们难道要找所有prime[i]*d<n的素数的个数吗?‘
显然不是。
以下字母变量和题目中的变量无关。
一个数m可以分解为m=a1p1a2p2…anpn,(a1,a2…an)为素数。他的最大因子d=a1p1-1a2p2…anpn
如果一个素数q。对于第一个m%q=0,那么这个q显然就是(a1,a2…an)中的最小的一个。那么如果下次遇到m%q2=0,q2就要比q大了,我们把q2乘到等式两边:q2*m=q2*(a1p1a2p2…anpn),这个q2现在是(a1,a2…an)中的一个,我们假设它是a2,(因为我们假设q是a1),那么现在就是q2*m=a1p1a2p2+1a3p3…anpn,那么对于m来说,它的最大因子现在就变成了a1p1-1a2p2+1a3p3…anpn,因为a2>a1,这显然比原本的d大。
所以我们一旦找到题目中的d%prime[i]==0,我们就不能继续往下找了。
下面是代码:
#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
const int maxn = 1e5+7;
int isprime[maxn];
int prime[maxn];
int cnt = 0;
void getprime()
{
memset(isprime,INF,sizeof(isprime));
isprime[0] = isprime[1] = 0;
for(int i=2;i<maxn;i++)
{
if(isprime[i])
{
prime[++cnt] = i;
}
for(int j=1;j<=cnt && i*prime[j]<maxn;j++)
{
isprime[i*prime[j]] = 0;
if(i%prime[j]==0) break;
}
}
}
int main()
{
int T;
scanf("%d",&T);
getprime();
int n,d;
while(T--)
{
int ans = 0;
scanf("%d%d",&n,&d);
for(int i=1;i<=cnt;i++)
{
if(d*prime[i]>=n) break;
ans++;
if(d%prime[i]==0) break;
}
printf("%d\n",ans);
}
return 0;
}