【HDOJ】【BestCoder Round #84】1004 Dertouzos

问题描述

正整数$x$称为$n$的$positive$ $proper$ $divisor$, 当且仅当$x | n$并且$1 \le x < n$. 例如, $1, 2,$ 和$3$是$6$的$positive$ $proper$ $divisor$, 但是$6$不是.

$Peter$给你两个正整数$n$和$d$. 他想要知道有多少小于$n$的整数, 满足他们的最大$positive$ $proper$ $divisor$恰好是$d$.

输入描述

输入包含多组数据, 第一行包含一个整数$T (1 \le T \le 10^6)$表示测试数据组数. 对于每组数据:

第一行包含两个整数$n$和$d$ $(2 \le n, d \le 10^9)$.

输出描述

对于每组数据, 输出一个整数.

$Solution$

我们反过来考虑，给定一个$d$，它可以作为谁的$ppd$（原谅我懒）

以下为最终结论……

根据题目中给出的定义，我们可以知道，如果一个数$a*d$，它的$ppd$是$d$，那么需要满足的条件就是$a \le b$ ，$b$是$d$的非$1$最小因数，$a$为质数。
（如果$a$不为质数，它就可以拆出个最小非$1$因数$e$，则$a*d=e*(\frac ae*d)$ ，也就是说，它的$ppd$就是$(\frac ae*d)$，而不是$d$）
（类似的，我们可以知道如果$a>b$，也可以导出$d$不是$ppd$的结论）
那么问题就是统计有多少质数$a$，满足
$1、$ $a \le b$ （$b$为$d$最小非$1$约数）
$2、$ $a*d < n$
对于$d$较小时，采用第一种方式约束枚举次数，反之使用第二种方式
（然而$Claris$大大放的加强数据中有$4e5$个超大质数，第一版数据具体怎样也不知道，我只用了第二种方式约束就跑过了）
（自己简单估计，应该是$1e5$上下的质数一多，我的程序就 T  <script id="MathJax-Element-19509" type="math/tex">T</script>飞了）

#include<stdio.h>
#define lim 100007 

bool b[lim+1];
int p[lim],n,d,T,tot;

int main()
{
    for (int i=2;i<=lim;i++) if (!b[i]) for (int j=i<<1;j<=lim;j+=i) b[j]=1;
    for (int i=2;i<=lim;i++) if (!b[i]) p[tot++]=i;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d",&n,&d);
        int i;
        for (i=0;i<tot && d*p[i]<n;i++) if (d % p[i]==0) 
        {
            i++;
            break;
        }
        printf("%d\n",i);
    }
}