HDU6715 算术 [莫比乌斯函数]

$$算术$$

题目描述见链接 .

---

$正解部分$

现在要求解 $\sum _{i=1}^n\sum _{j=1}^m\mu (lcm(i,j))$,

因为 $\mu (lcm(i,j))=\mu (i)\mu (j)\mu (gcd(i,j))$,

$证$:

假设 $\mu (lcm(i,j))=0$, 说明 $\frac{ij}{gcd(i,j)}$ 分解质因数后含有幂数超过 $1$ 的质因子, 设其为 $p$,
若 $p|i$ 且 $p\nmid j$, 则 $\mu (i)=0$, 此时上式成立 .
若 $p|i$ 且 $p|j$, 设在 $i$ 中 $p$ 的幂数为 $a$, $j$ 中的幂数为 $b$, 则 $a+b-\min (a,b)>1$, 令 $a>b$, 则 $a>1$, 此时 $\mu (i)=0$, 此时上式成立 .

假设 $\mu (lcm(i,j))\ne 0$, 则幂数相加的奇偶性与幂数相减的奇偶性相同, 不会改变 $\mu$ 取值, 上式成立 .

$\therefore 综上所述,上式成立$

所以原式 $\sum _{i=1}^n\sum _{j=1}^m\mu (i)\mu (j)\mu (gcd(i,j))$,
枚举 $d=gcd(i,j)$ 得

$\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }\mu (id)\mu (jd)[gcd(i,j)==1]$

根据 这里 所述的 $\mu$ 性质 $2$,

$\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{d}\rfloor }\mu (id)\mu (jd)\sum _{k|gcd(i,j)}\mu (k)$

$\sum _{d=1}^{\min (n,m)}\mu (d)\sum _{k=1}^{\lfloor \frac{\min (n,m)}{d}\rfloor }\mu (k)\sum _{i=1}^{\lfloor \frac{n}{kd}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{kd}\rfloor }\mu (idk)\mu (jdk)$

设 $T=dk$,

$\sum _{T=1}^{\min (n,m)}\sum _{d\mid T}\mu (d)\mu (\frac{T}{d})\sum _{i=1}^{\lfloor \frac{n}{T}\rfloor }\mu (iT)\sum _{j=1}^{\lfloor \frac{m}{T}\rfloor }\mu (jT)$

---

$实现部分$

预处理括号内元素即可实现 $O(N\log N)$ 计算答案 .

$\sum _{T=1}^{\min (n,m)}\left(\sum _{d\mid T}\mu (d)\mu (\frac{T}{d})\right)\left(\sum _{i=1}^{\lfloor \frac{n}{T}\rfloor }\mu (iT)\right)\left(\sum _{j=1}^{\lfloor \frac{m}{T}\rfloor }\mu (jT)\right)$

#include<bits/stdc++.h>
#define reg register
typedef long long ll;

int read(){
        char c;
        int s = 0, flag = 1;
        while((c=getchar()) && !isdigit(c))
                if(c == '-'){ flag = -1, c = getchar(); break ; }
        while(isdigit(c)) s = s*10 + c-'0', c = getchar();
        return s * flag;
}

const int maxn = 1e6 + 5;
const int lim  = 1e6;

int N;
int M;
int pc;
int p[maxn];
int mu[maxn];
int s2[maxn];
int s3[maxn];
int vis[maxn];
int s1[maxn+1];

void Init(){
        mu[1] = 1; 
        for(reg int i = 2; i <= lim; i ++){
                if(!vis[i]) p[++ pc] = i, mu[i] = -1;
                for(reg int j = 1; j <= pc && p[j]*i <= lim; j ++){
                        vis[p[j]*i] = 1;
                        if(i % p[j] == 0){ mu[p[j]*i] = 0; break ; }
                        mu[p[j]*i] = -mu[i];
                }
        }
        for(reg int i = 1; i <= lim; i ++)
                for(reg int j = 1; j <= lim/i; j ++) s1[i*j] += mu[i] * mu[j];
}

void Work(){
        memset(s2, 0, sizeof s2); memset(s3, 0, sizeof s3);
        N = read(), M = read();
        for(reg int i = 1; i <= N; i ++)
                for(reg int j = 1; j <= N/i; j ++) s2[i] += mu[i * j];
        for(reg int i = 1; i <= M; i ++)
                for(reg int j = 1; j <= M/i; j ++) s3[i] += mu[i * j];
        N = std::min(N, M);
        ll Ans = 0;
        for(reg int i = 1; i <= N; i ++) Ans += 1ll*s1[i]*s2[i]*s3[i];
        printf("%lld\n", Ans);
}

int main(){
        Init();
        int test_cnt = read();
        while(test_cnt --) Work();
        return 0;
}