HDU 5750 Dertouzos（数论）

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;
}