HDU6715 算术 [莫比乌斯函数]

$算术$

题目描述见链接 .

---

$color{red}{正解部分}$

现在要求解 $sumlimits_{i=1}^nsumlimits_{j=1}^mmu(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 ot| 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)) ot= 0$, 则幂数相加的奇偶性与幂数相减的奇偶性相同, 不会改变 $mu$ 取值, 上式成立 .

$herefore 综上所述, 上式成立$

所以原式 $sumlimits_{i=1}^nsumlimits_{j=1}^m mu(i)mu(j)mu(gcd(i,j))$,
枚举 $d = gcd(i, j)$ 得

$sumlimits_{d=1}^{min(n,m)}mu(d)sumlimits_{i=1}^{lfloorfrac{n}{d} floor}sumlimits_{j=1}^{lfloorfrac{m}{d} floor}mu(id)mu(jd)[gcd(i,j)==1]$

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

$sumlimits_{d=1}^{min(n,m)}mu(d)sumlimits_{i=1}^{lfloorfrac{n}{d} floor}sumlimits_{j=1}^{lfloorfrac{m}{d} floor}mu(id)mu(jd)sumlimits_{k|gcd(i,j)} mu(k)$

$sumlimits_{d=1}^{min(n, m)}mu(d) sumlimits_{k=1}^{lfloor frac{min(n, m)}{d} floor } mu(k) sumlimits_{i=1}^{lfloor frac{n}{kd} floor} sumlimits_{j=1}^{lfloor frac{m}{kd} floor} mu(idk)mu(jdk)$

设 $T = dk$,

$sumlimits_{T=1}^{min(n,m)}sumlimits_{dmid T}mu(d)mu(frac{T}{d})sumlimits_{i=1}^{lfloorfrac{n}{T} floor}mu(iT)sumlimits_{j=1}^{lfloorfrac{m}{T} floor}mu(jT)$

---

$color{red}{实现部分}$

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

$sumlimits_{T=1}^{min(n,m)} left(sumlimits_{dmid T}mu(d)mu(frac{T}{d}) ight) left(sumlimits_{i=1}^{lfloorfrac{n}{T} floor}mu(iT) ight) left(sumlimits_{j=1}^{lfloorfrac{m}{T} floor}mu(jT) ight)$

#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
", Ans);
}

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