hdu 2588 GCD---欧拉函数

题意：

求有多少 x ( 1 <= x <= n ), 使得 gcd ( x, n ) >= m。

其中(2<=N<=1000000000, 1<=M<=N)。

解析：

设 x = q * d， n = p * d ，且 gcd （ x， n） = d 。

则 q 和 p 必然是互素的，题目要求 d >= m，所以只要找出当 d >= m时：

p = （n / d ），小于p，且与p互素的数有多少个就行了，这就相当于求x有多少个了。

这个数就是欧拉函数...

然后求得时候肯定不是一遍爆过去，而是像上次多校赛那题一样。

现在换一种写法，以后这样写开根号：

for (int i = 1; i * i <= n; i++)

这样就把复杂度降下来了，并且找到了所有p (●'◡'●)

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <stack>
#include <vector>
#include <queue>
#include <map>
#include <climits>
#include <cassert>
#define LL long long

using namespace std;
const int inf = 0x3f3f3f3f;
const double eps = 1e-8;
const double pi = 4 * atan(1.0);
const double ee = exp(1.0);

int euler_phi(int n)
{
    int m = sqrt(n + 0.5);
    int res = n;
    for (int i = 2; i <= m; i++)
    {
        if (n % i == 0)
        {
            res = res / i * (i - 1);
            while (n % i == 0)
                n /= i;
        }
    }
    if (1 < n)
        res = res / n * (n - 1);
    return res;
}

LL gcd(LL a, LL b)
{
    return b ? gcd(b, a % b) : a;
}

int main()
{
#ifdef LOCAL
    freopen("in.txt", "r", stdin);
#endif // LOCAL
    int ncase;
    scanf("%d", &ncase);
    while (ncase--)
    {
        int n, m;
        scanf("%d%d", &n, &m);
        int ans = 0;
        for (int i = 1; i * i <= n; i++)
        {
            if (n % i == 0)
            {
                if (i >= m)
                    ans += euler_phi(n / i);
                if (n / i != i && n / i >= m)
                    ans += euler_phi(i);
            }
        }
        printf("%d\n", ans);
    }
    return 0;
}