【BZOJ】1101 [POI2007]Zap（莫比乌斯反演）

题目

传送门：QWQ

分析

莫比乌斯反演。

还不是很熟练qwq

代码

//bzoj1101
//给出a,b,d,询问有多少对二元组(x,y)满足gcd(x,y)=d.x<=a,y<=b

#include <bits/stdc++.h>
using namespace std;
const int maxn = ;
int n, p[maxn];
int mu(int m) {
    int tmp=, k=m;
    for(int i=;i*i<=k;i++) {
        if(!(m%i)) {
            tmp++; m/=i;
            if(!(m%i)) return ;
        }
    }
    if(m>) tmp++;
    return (tmp&)?-:;
}

int main() {
    for(int i=;i<=maxn-;i++)
        p[i] = p[i-] + mu(i);
    scanf("%d",&n);
    for(int i=;i<=n;i++) {
        int a, b, d;
        scanf("%d%d%d",&a,&b,&d);
        a/=d; b/=d;
        if(a>b) swap(a,b);
        int c=,ans=;
        for(int j=;j<=a;j=c+) {
            c=min(a/(a/j), b/(b/j));
            ans+=(p[c]-p[j-])*(a/j)*(b/j);
        }
        printf("%d\n",ans);
    }
    return ;
}