线性筛+最大因子【HDU5750】

传送门藏起来啦

题目要求小于n的并且最大因子为d的数的个数。
我们设m=kd，k是素数我们才能保证m的最大因子是d。
因为若k不是素数，我们就可以将k分解成ab，那么m的最大因子就可以是ad或者bd，反正不是d了。

所以我们考虑一下线性筛。
但是我们难道要找所有prime[i]*d<n的素数的个数吗？‘
显然不是。

以下字母变量和题目中的变量无关。

一个数m可以分解为m=a1^p1a2^p2…an^pn，(a1,a2…an)为素数。他的最大因子d=a1^p1-1a2^p2…an^pn
如果一个素数q。对于第一个m%q=0，那么这个q显然就是(a1,a2…an)中的最小的一个。那么如果下次遇到m%q2=0，q2就要比q大了，我们把q2乘到等式两边:q2*m=q2*(a1^p1a2^p2…an^pn)，这个q2现在是(a1,a2…an)中的一个，我们假设它是a2，(因为我们假设q是a1)，那么现在就是q2*m=a1^p1a2^p2+1a3^p3…an^pn，那么对于m来说，它的最大因子现在就变成了a1^p1-1a2^p2+1a3^p3…an^pn，因为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;
}