HDU-2588（水题1）

来一道数论水题

题目大意呢，就是给你一个数n,m,问你在1~n中有多少个x，使得gcd（x,n)>=m。

这个题啊，我是一点也不觉得哪里算水，毕竟我太菜了。

题解：我们可以先假设找出了这个公约数，设为s，那么x和n提出这个s后，我们分别设为a和b，即（x，y）->(s x a，s x b)，那么重点来了，注意集中精力啊，可以得出的a和b就是两个互质的数！因为x<=n，所以a<=b，所以我们这么来说，a是所有与b互质的数的其中一个！那么与b互质的并且小于b的这样的a有多少个呢？那么这个题就出来了——欧拉函数！
我们枚举所有的s (1~n)，（为什么s不是（m ~n）？见代码 (●’◡’●)），也就是for（1-n）,如果n%i==0，并且i>=m，那么这个i就可以当作一个s，n/i就可以当作一个b，然后phi（b）就可以当作a的可取的个数，ans+=phi（b)就可以了。然后i可以当作s，我们再看看n/i能不能当作s，这样枚举一边就ok了，是不是很妙？（我觉得妙，因为不会）

接下来奉上我孱弱的代码

//#pragma GCC optimize(3,"Ofast","inline")
//#include<unordered_map>
//#include<unordered_set>
#include<cstdio>
#include<iostream>
#include<cmath>
#include<functional>
#include<cstring>
#include<string>
#include<cstdlib>
#include<queue>
#include<map>
#include<algorithm>
#include<set>
#include<stack>
#include<vector>
using namespace std;
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=1e7+10;
const ll mod=1e9+7;
ll phi(ll x)
{
    ll ans=x;
    for(ll i=2;i*i<=x;i++)
    {
        if(x%i==0)ans=ans/i*(i-1);
        while(x%i==0)x/=i;
    }
    if(x>1)ans=ans/x*(x-1);
    return ans;
}
int main()
{
    ll x,n,m,t,ans;
    scanf("%lld",&t);
    while(t--)
    {
        ans=0;
        scanf("%lld%lld",&n,&m);
        for(ll i=1;i*i<=n;i++)//这里为啥不从m开始？因为等于是是求因数嘛，边界是i*i<=n,只有这个方法才能保证不会遗漏并且只枚举一遍！
        {
            if(n%i==0)
            {
                if(i>=m)ans+=phi(n/i);
                if(i*i!=n&&n/i>=m)ans+=phi(i);
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}