hdu 2588 欧拉函数

hdu 2588
http://acm.hdu.edu.cn/showproblem.php?pid=2588

gcd(x,n)=p≥m
gcd(x,n)=p -> gcd(x/p,n/p)=1
∵ 1≤x≤n
∴求x的个数 等价 若n的因数pi≥m Σφ(n/pi) (1≤i≤n)
φ(n/p)为欧拉函数

#include<bits/stdc++.h>
typedef long long ll;
using namespace std;
ll hsy(ll n)
{
    ll s=n;
    for(ll i=2; i*i<=n; i++)
    {
        if(n%i==0)
        {
            s=s/i*(i-1);
            while(n%i==0)
                n/=i;
        }
    }
    if(n>1)
        s=s/n*(n-1);
    return s;
}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        ll n,m,ans=0;
        cin>>n>>m;
        for(ll i=1; i*i<=n; i++)
        {
            if(n%i==0)
            {
                if(i*i<n)
                {
                    if(i>=m)
                        ans+=hsy(n/i);
                    if(n/i>=m)
                        ans+=hsy(i);
                }
                else if(i*i==n&&i>=m)
                    ans+=hsy(i);
            }
        }
        cout<<ans<<endl;
    }
}