Description
正整数x称为n的positive proper divisor, 当且仅当x | n并且1≤x< n. 例如, 1, 2, 和3是6的positive proper divisor, 但是6不是
Peter给你两个正整数n和d. 他想要知道有多少小于n的整数, 满足他们的最大positive proper divisor恰好是d
Input
输入包含多组数据, 第一行包含一个整数T(1≤T≤10^6 )表示测试数据组数. 对于每组数据:包含两个整数n和d(2≤n,d≤10^9)
Output
对于每组数据, 输出一个整数
Sample Input
9
10 2
10 3
10 4
10 5
10 6
10 7
10 8
10 9
100 13
Sample Output
1
2
1
0
0
0
0
0
4
Solution
问题转化为找x的个数,使得d< d*x< n且d是d*x的最大因子,设d的最小素因子为m,简单推导知x只能取2~min((n-1)/d, m)中的素数,而min((n-1)/d, m)至多是sqrt(n)级别,所以只需要预处理出5*10^4中所有的素数即可
Code
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
#define maxn 55555
int mark[maxn],prime[maxn],num[maxn],res;
void get_prime(int n)
{
memset(mark,0,sizeof(mark));
res=0;
num[0]=num[1]=0;
for(int i=2;i<=n;i++)
{
if(!mark[i])
{
mark[i]=prime[res++]=i;
num[i]=num[i-1]+1;
}
else num[i]=num[i-1];
for(int j=0;j<res&&prime[j]*i<=n;j++)
{
mark[i*prime[j]]=prime[j];
if(i%prime[j]==0)
break;
}
}
}
int main()
{
get_prime(maxn-10);
int T,n,d;
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&d);
if(d>=(n+1)/2)printf("0\n");
else
{
int temp=(n-1)/d,flag=0;
for(int i=0;i<res&&prime[i]*prime[i]<=d;i++)
if(d%prime[i]==0)
{
flag=1;
temp=min(temp,prime[i]);
break;
}
if(!flag)temp=min(temp,d);
printf("%d\n",num[temp]);
}
}
return 0;
}