HDU2588-GCD(欧拉函数在gcd里的小运用)

题目

题意: 给出N,M,求出有多少整数X满足1<=X<=N and gcd(X,N)>=M.
我第一直觉:M可能不是N的因子,那就让M调大一点变成N的因子。然后我开始傻逼了。。gcd(X/M,N/M)>=1 ，那这个式子答案不就是X/M。。。妈的，这样求的话求出来的不是1-N中有多少是M的倍数?。。。唉。
既然‘>=’不会求。我们会求‘=’啊。‘=’就是gcd(X,N)=M时,X的个数就是phi(N/M).
把N的因子d都枚举出来,求sigma(phi(N/d))；
这个T<=100,2<=N<=1000000000, 1<=M<=N。完全不像这个一样需要打个表。

#include<cstdio>
#include<cmath>
#include<iostream>
#define en '\n'
using namespace std;
typedef long long ll;
int solve(int x)//计算phi(x)直接使用素数求欧拉值会比较快...
{
    if(x==1) return 1;
    int res=x;
    for(int i=2;(ll)i*i<=x;++i)
    {
        if(x%i==0)
        {
            res-=res/i;
            x/=i;
            while(x%i==0) x/=i;
        }
    }
    if(x>1) res-=res/x;
    return res;
}
int main()
{
    int T;scanf("%d",&T);
    while(T--)
    {
        ll n,m;scanf("%lld%lld",&n,&m);
        ll ans=0;
        for(int i=1;(ll)i*i<=n;++i)
        {
            if(n%i) continue;
            if(i>=m) ans+=solve(n/i);//枚举>=m的因子
            if(n/i>=m&&i*i!=n) ans+=solve(i);//枚举>=m的因子
        }
        printf("%lld\n",ans);
    }
}