hdu 5750 Dertouzos

##代码

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#define LL long long
using namespace std;
bool isprime[100010];
int prime[100010],cnt;
void init(){
    memset(isprime,1,sizeof isprime);
    isprime[0]=isprime[1]=false;
    cnt=0;
    for(int i=2;i<=100000;++i){
        if(isprime[i]){
            prime[cnt++]=i;
            for(LL j=(LL)i*i;j<=100000;j+=i)
                isprime[j]=false;    
        }    
    }    
}
int main(){
    int t;
    init();
    scanf("%d",&t);
    while(t--){
        int n,d;
        int ans=0;
        scanf("%d%d",&n,&d);
        n=(n-1)/d;
        for(int i=0;i<cnt&&prime[i]<=n;++i){
            ++ans;
            if(d%prime[i]==0)
                break;    
        }
        printf("%d\n",ans);
    }
    //system("pause");
    return 0;    
}

##总结

1. 设p为d的最小质因子，a为2到p的所有质数，以d为最大非平凡（真）因子的数为a*d。

###证明：

设$a\ast d=n$,d的最小的质因子为p，$d=p\ast p^{{}^{\prime }}$,$p^{{}^{\prime }}$最小为一

假设$a>p$

$n=a\ast p\ast p^{{}^{\prime }}=p\ast (a\ast p^{{}^{\prime }})$

$\because a>p$

$\therefore a\ast p^{{}^{\prime }}>p\ast p^{{}^{\prime }}=d$

$\therefore a>p时，d不是n的最大因子$

$\therefore d时n的最大因子时，a<=p$

$\because d为n的最大因子时，a一定为n的最小因子，而一个数的最小因子一定为质数$

至此结论得证