hdu 2588 欧拉函数的应用

/*
 * hdu2588/linux.c
 * Created on: 2011-8-2
 * Author    : ben
 */
#include <stdio.h>
#include <math.h>

int phi(int x) {
    int i, res = x, temp = (int) sqrt(x);
    for (i = 2; i < temp + 1; i++) {
        if (x % i == 0) {
            res = res / i * (i - 1);
            while (x % i == 0) {
                x /= i;
            }
        }
    }
    if (x > 1) {
        res = res / x * (x - 1);
    }
    return res;
}

void work() {
    int T, N, M, i, sqrtn, ans;
    scanf("%d", &T);
    while (T--) {
        scanf("%d%d", &N, &M);
        sqrtn = sqrt(N);
        ans = 0;
        for (i = 1; i < sqrtn; i++) {
            if (N % i != 0) {
                continue;
            }
            if (i >= M) {
                ans += phi(N / i);
            }
            if (N / i >= M) {
                ans += phi(i);
            }
        }
        if (sqrtn * sqrtn == N && sqrtn >= M) {
            ans += phi(sqrtn);
        }
        printf("%d\n", ans);
    }
}

int main() {
#ifndef ONLINE_JUDGE
    freopen("data.in", "r", stdin);
#endif
    work();
    return 0;
}