HDU1695 GCD(莫比烏斯反演)

傳送門
看了1個多小時，終於懂了一點了
題目大意：給n,m,k.求gcd(x,y) = k(1<=x<=n, 1<=y<=m)的個數
思路：令F(i)表示i|gcd(x,y)的(x,y)的對數，顯然F(x)=[nx]∗[mx] 。
設f(x)爲gcd(x,y)=x的對數。
因爲F(x)=∑i|xf(i) ，所以我們可以莫比烏斯反演它。
根據公式f(x)=∑x|dμ(d)F(d)
我們的目標就是f(1)（因爲n和m都可以除以k）
所以我們就可以在O(n) 的時間複雜度內求出答案了。

#include<cstdio>
#include<algorithm>
using namespace std;
#define MAXN 100005
int n, m, k, miu[MAXN], p[MAXN], cnt;
bool vis[MAXN];
void sieve() {
    miu[1] = 1;
    for(int i = 2; i < MAXN; ++ i) {
        if(!vis[i]) {
            miu[i] = -1;
            p[++ cnt] = i;
        }
        for(int j = 1; j <= cnt; ++ j) {
            if(p[j] * i > MAXN) break;
            vis[p[j]*i] = 1;
            if(i % p[j] == 0) {
                miu[i * p[j]] = 0;
                break;
            }
            else miu[i * p[j]] = - miu[i];
        }
    }
}
int main() {
    sieve(); int T, a, c; scanf("%d", &T);
    for(int C = 1; C <= T; ++ C) {
        printf("Case %d: ", C);
        scanf("%d%d%d%d%d", &a, &n, &c, &m, &k);
        if(!k) { puts("0"); continue; }
        n /= k; m /= k;
        if(n > m) swap(n, m);
        long long ans1 = 0, ans2 = 0;
        for(int i = 1; i <= n; ++ i)
            ans1 += 1LL * miu[i] * (n/i) * (m/i);
        for(int i = 1; i <= n; ++ i)
            ans2 += 1LL * miu[i] * (n/i) * (n/i);
        //因爲要求無序的數對，所以我們可以減去ans2/2
        ans1 -= ans2/2;
        printf("%I64d\n", ans1);
    }
}