FJUT1251 TMD 欧拉函数

LH Boy ，别看他平时TMD(挺萌的)，到处卖萌，疯疯癫癫的，然而，他却是一位常年混迹在各大Math论坛上，帮助无数误入迷途的Mather解决了各种难题的大神，然而，最近他遇到了一个很蛋疼的问题。这是一个来自论坛上Mather的难题：
　　设X∈[1~N]，存在多少个X使得GCD(X,N)>=M，统计符合要求X的个数。
Input
　　第一行输入一个整数 T( T<=100) ，表示有T组测试数据。
　　接下来有T组测试数据，每一组测试数据输入一行。每一行输入两个整数 N和 M (2<=N<=1000000000, 1<=M<=N)。
Output
　　对于每一组测试案例，在下一行中输出答案、
SampleInput
3
1 1
10 2
10000 72
SampleOutput
1
6
260

#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;
    }
}