【HDU-2588】GCD

题目链接

题目描述

输入 N 和 M (2<=N<=1000000000, 1<=M<=N), 找出所有的X满足1<=X<=N 且 gcd(X,N)>=M.

题解

1.首先若m==1,则答案为N
2.m≠1时，初步想法，将N进行质因数分解，则答案为其质因数的组合的积中大于m的组合个数，但实现起来不现实。
不妨设$gcd(X,N)=Y >=M$
则: $令X=K*Y，其中一定有gcd(K,N)=1 且 Y|N，即Y是N的约数$
于是有以下做法:先求出N的所有约数Y，$Φ(N/Y)$之和(即符合条件的K的个数)
即为答案
为什么是$Φ(N/Y)$呢？ 因为这样才能保证 $K*Y<=N$
实际实现时可通过折半枚举来降低复杂度,即若x是N的约数,那么N/x也是。

#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<algorithm>
inline int phi(int a)
{
    int ans=a;
    for(register int i=2;i<=sqrt(a);i++)
    {
        if(a%i==0)
        {
            ans=ans/i*(i-1);
            while(a%i==0) a/=i;
        }
    }
    if(a>1) ans=ans-ans/a;
    return ans;
}
int main()
{
    int T;register int n,m;
    scanf("%d",&T);
    while(T--){
        scanf("%d %d",&n,&m);
        register int ans=0;
        if(m==1) {printf("%d\n",n);continue;}
        for(register int i=1;i*i<=n;i++){
            if(n%i==0){
                if(i>=m)ans+=phi(n/i);
                if(n/i>=m&&(i*i!=n)) ans+=phi(i);
            }
        }
        printf("%d\n",ans);
    }
}