莫比乌斯反演例题

欢迎访问本菜鸡的独立博客：Codecho

UPDATED on $\text{Apr}.22,2019$.

0. 套路

$(1)$ $\left[n=1\right]={\sum }_{d|n}\mu \left(d\right)$

$(2)$ ${\sum }_{i=1}^a{\sum }_{j=1}^b\left[\text{gcd}\left(i,j\right)=d\right]={\sum }_{i=1}^{\lfloor \frac{a}{d}\rfloor }{\sum }_{j=1}^{\lfloor \frac{b}{d}\rfloor }\left[\text{gcd}\left(i,j\right)=1\right]$

$(3)$ $\left[ d|\text{gcd}\left( i,j \right) \right] =\left[ d|i \right] \left[ d|j \right] $

$(4)$ ${\sum }_{i,j,...}f\left(i,j,...\right)={\sum }_{k=1}^{maxvalue}{\sum }_{i,j,...}\left[f\left(i,j,...\right)=k\right]\cdot f\left(i,j,...\right)$

$(5)$ 若 $l=kd$，则 ${\sum }_{k=1}^n{\sum }_{d=1}^{\lfloor \frac{n}{k}\rfloor }f={\sum }_{l=1}^n{\sum }_{k|l}f={\sum }_{l=1}^n{\sum }_{d|l}f$，之后一般可推出一个积性函数乘上一些式子

$(6)$ ${\sum }_{d|n}1={\sum }_{d=1}^n\left[d|n\right]$ 表示 $[1,n]$ 范围内 $d$ 的倍数的个数。

$(7)$ ${\sum }_{d|n}d={\sum }_{d=1}^n\left[d|n\right]\cdot d$ 表示 $[1,n]$ 范围内 $d$ 的倍数之和。

1. BZOJ - 1101

1.1 题面

求 ${\sum }_{i=1}^a{\sum }_{j=1}^b\left[\text{gcd}\left(i,j\right)=d\right]$ 的值。

数据范围：$1\le d\le a,b\le 50000$，共 $T(1\le T\le 50000)$ 组数据。

1.2 解题过程

$$\begin{array}{rl}& \sum _{i=1}^a\sum _{j=1}^b\left[\text{gcd}\left(i,j\right)=d\right]\\ & =\sum _{i=1}^{\lfloor \frac{a}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{b}{d}\rfloor }\left[\text{gcd}\left(i,j\right)=1\right]\\ & =\sum _{i=1}^{\lfloor \frac{a}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{b}{d}\rfloor }\sum _{d^{\prime }|\text{gcd}\left(i,j\right)}\mu \left(d^{\prime }\right)\\ & =\sum _{i=1}^{\lfloor \frac{a}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{b}{d}\rfloor }\sum _{d^{\prime }=1}^{\min \left(\lfloor \frac{a}{d}\rfloor ,\lfloor \frac{b}{d}\rfloor \right)}\left[d^{\prime }|\text{gcd}\left(i,j\right)\right]\mu \left(d\right)\\ & =\sum _{i=1}^{\lfloor \frac{a}{d}\rfloor }\sum _{j=1}^{\lfloor \frac{b}{d}\rfloor }\sum _{d^{\prime }=1}^{\min \left(\lfloor \frac{a}{d}\rfloor ,\lfloor \frac{b}{d}\rfloor \right)}\left[d^{\prime }|i\right]\left[d^{\prime }|j\right]\mu \left(d\right)\\ & =\sum _{d^{\prime }=1}^{\min \left(\lfloor \frac{a}{d}\rfloor ,\lfloor \frac{b}{d}\rfloor \right)}\mu \left(d\right)\cdot \left(\sum _{i=1}^{\lfloor \frac{a}{d}\rfloor }\left[d^{\prime }|i\right]\right)\cdot \left(\sum _{j=1}^{\lfloor \frac{b}{d}\rfloor }\left[d^{\prime }|j\right]\right)\\ & =\sum _{d^{\prime }=1}^{\min \left(\lfloor \frac{a}{d}\rfloor ,\lfloor \frac{b}{d}\rfloor \right)}\mu \left(d\right)\lfloor \frac{\lfloor \frac{a}{d}\rfloor }{d^{\prime }}\rfloor \lfloor \frac{\lfloor \frac{b}{d}\rfloor }{d^{\prime }}\rfloor \end{array}$$

1.3 优化技巧

本题需要分块，以达到 $O(T\cdot (\sqrt{a}+\sqrt{b}))$ 的时间复杂度。

分块代码如下：

int answer = 0;
int pos = 0, last = 0;
int p = a / d, q = b / d;
int minvalue = min(p, q);
while(pos + 1 <= minvalue)
{
    pos = min(p / (p / (pos + 1)), q / (q / (pos + 1)));
    answer += (prefix[pos] - prefix[last]) * (p / pos) * (q / pos);
    last = pos;
}

注意使用 long long 可能会被卡常数。

2. BZOJ - 2301

2.1 题面

求 ${\sum }_{i=a}^b{\sum }_{j=c}^d[\text{gcd}(i,j)=k]$ 的值。

数据范围：$1\le n,k\le 5\cdot 10^4$，$1\le a\le b\le 5\cdot 10^4$，$1\le c\le d\le 5\cdot 10^4$。

2.2 解题过程

设 $sum(a,b)={\sum }_{i=1}^a{\sum }_{j=1}^b[\text{gcd}(i,j)=k]$ ，则答案为二维前缀和：
$$sum(b,d)-sum(a-1,d)-sum(b,c-1)+sum(a-1,c-1)$$
每个 $sum$ 都可以用与 BZOJ - 1101 相同的思路进行求解。

3. HDOJ - 1695

3.1 题面

给定 $b,d,k$ 求 $i\in [1,b],j\in [1,d]$ 范围内满足 $\text{gcd}(i,j)=k$ 的 $(i,j)$ 对数。

规定 $(i,j)$ 和 $(j,i)$ 算相同的一对。

共有 $T(1\le T\le 3000)$ 组数据。

数据范围：$1\le c,d\le 10^5$。

3.2 解题过程

设 $sum(a,b)={\sum }_{i=1}^a{\sum }_{j=1}^b[\text{gcd}(i,j)=k]$ ，可以发现重复的数对一定出现在 $[1,\min (b,d)]$ 的范围内，则答案为
$$sum(b,d)-\frac{sum(\min (b,d))}{2}$$
。

4. Luogu - P2303

4.1 题面

给定 $n$，求 ${\sum }_{i=1}^n\text{gcd}\left(i,n\right)$ 的值。

数据范围：$0<n\le 2^{32}$。

4.2 解题过程

$$\begin{array}{rl}& \sum _{i=1}^n\text{gcd}\left(i,n\right)\\ =& \sum _{k|n}\sum _{i=1}^n\left[\text{gcd}\left(i,n\right)=k\right]\cdot k,\left(k|\text{gcd}\left(i,n\right)\right)\\ =& \sum _{k|n}k\sum _{i=1}^{\frac{n}{k}}\left[\text{gcd}\left(i,n\right)=1\right]\\ =& \sum _{k|n}k\cdot \phi \left(\frac{n}{k}\right)\end{array}$$

我们可以枚举 $n$ 的因子，将每个因子的贡献加到答案中。

时间复杂度：$O(\sqrt{n\sqrt{n}}=n^{\frac{3}{4}})$。

4.3 代码

ll getphi(ll n)
{
    ll res = n;
    for(ll i = 2; i * i <= n; i++)
    {
        if(n % i == 0)
        {
            res = res - res / i;
            while(n % i == 0)
                n /= i;
        }
    }
    if(n > 1)
        res = res - res / n;
    return res;
}

int main()
{
    ll n;
    scanf("%lld", &n);

    ll answer = 0;
    for(ll i = 1; i * i <= n; i++)
    {
        if(n % i)   continue;
        answer += i * getphi(n / i);
        if(n / i != i)  answer += n / i * getphi(i);
    } 
    printf("%lld\n", answer);
    return 0;
}

5. BZOJ - 2154

5.1 题面

求 ${\sum }_{i=1}^n{\sum }_{j=1}^m\text{lcm}\left(i,j\right)$ 的值，答案对 $20101009$ 取模。

数据范围：$1\le n,m\le 10^7$。

5.2 解题过程

$$\begin{array}{rl}& \sum _{i=1}^n\sum _{j=1}^m\text{lcm}\left(i,j\right)\\ & =\sum _{k=1}^{\min \left(n,m\right)}\sum _{i=1}^n\begin{array}{c}{\sum }_{j=1}^m\frac{i\cdot j}{k}\left[\text{gcd}\left(i,j\right)=k\right]\end{array}\\ & =\sum _{k=1}^{\min \left(n,m\right)}\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }ijk\left[\text{gcd}\left(i,j\right)=1\right],\left(i\leftarrow k\cdot i,j\leftarrow k\cdot j\right)\\ & =\sum _{k=1}^{\min \left(n,m\right)}k\sum _{d=1}^{\min \left(\lfloor \frac{n}{k}\rfloor ,\lfloor \frac{m}{k}\rfloor \right)}\mu \left(d\right)\left(\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\left[d|i\right]\cdot i\right)\left(\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }\left[d|j\right]\cdot j\right)\\ & =\sum _{k=1}^{\min \left(n,m\right)}k\sum _{d=1}^{\min \left(\lfloor \frac{n}{k}\rfloor ,\lfloor \frac{m}{k}\rfloor \right)}\mu \left(d\right)\left(\sum _{p=1}^{\lfloor \frac{n}{kd}\rfloor }dp\right)\left(\sum _{q=1}^{\lfloor \frac{m}{kd}\rfloor }dq\right),\left(i\leftarrow dp,j\leftarrow dq\right)\\ & =\sum _{k=1}^{\min \left(n,m\right)}k\sum _{d=1}^{\min \left(\lfloor \frac{n}{k}\rfloor ,\lfloor \frac{m}{k}\rfloor \right)}{d}^2\cdot \mu \left(d\right)\left(\sum _{p=1}^{\lfloor \frac{n}{kd}\rfloor }p\right)\left(\sum _{q=1}^{\lfloor \frac{m}{kd}\rfloor }q\right)\\ & =\sum _{k=1}^{\min \left(n,m\right)}k\sum _{d=1}^{\min \left(\lfloor \frac{n}{k}\rfloor ,\lfloor \frac{m}{k}\rfloor \right)}{d}^2\cdot \mu \left(d\right)\left(\frac{\left(1+\lfloor \frac{n}{kd}\rfloor \right)\cdot \lfloor \frac{n}{kd}\rfloor }{2}\right)\left(\frac{\left(1+\lfloor \frac{m}{kd}\rfloor \right)\cdot \lfloor \frac{m}{kd}\rfloor }{2}\right)\\ & =\sum _{l=1}^{\min \left(n,m\right)}\left(\frac{\left(1+\lfloor \frac{n}{l}\rfloor \right)\cdot \lfloor \frac{n}{l}\rfloor }{2}\right)\left(\frac{\left(1+\lfloor \frac{m}{l}\rfloor \right)\cdot \lfloor \frac{m}{l}\rfloor }{2}\right)\sum _{d|l}\mu \left(d\right)\cdot ld\end{array}$$

5.3 优化技巧

在推导的过程中，我们得到了一个积性函数，因此优化的方法是 $O(n)$ 地通过线性筛预处理出积性函数的所有值。

这样，我们便可在 $O(1)$ 的时间内查询到函数的值。

6. BZOJ - 2693

6.1 题面

求 ${\sum }_{i=1}^n{\sum }_{j=1}^m\text{lcm}\left(i,j\right)$ 的值，共有 $T$ 组询问，答案对 $10^8+9$ 取模。

数据范围：$1\le n,m\le 10^7$，$1\le T\le 10^4$。

6.2 解题过程

分块完成 $5.2$ 节的计算即可，注意取模。

时间复杂度： O ( min ⁡ { n , m } ) ​ O(\sqrt {\min\left\{n,m\right\}})​ O(min{n,m} ​)​。

7. CCPC-Wannafly Winter Camp Day-3 Div.2 F 小清新数论

7.1 题面

求 ${\sum }_{i=1}^n{\sum }_{j=1}^n\mu \left(\text{gcd}\left(i,j\right)\right)$ 的值。

数据范围：$1\le n\le 10^7$。

7.2 解题过程

$$\begin{array}{rl}& \sum _{i=1}^n\sum _{j=1}^n\mu \left(\text{gcd}\left(i,j\right)\right)\\ & =\sum _{k=1}^n\sum _{i=1}^n\sum _{j=1}^n\mu \left(k\right)\left[\text{gcd}\left(i,j\right)=k\right]\\ & =\sum _{k=1}^n\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{k}\rfloor }\mu \left(k\right)\left[\text{gcd}\left(i,j\right)=1\right]\\ & =\sum _{k=1}^n\mu \left(k\right)\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{d|\text{gcd}\left(i,j\right)}\mu \left(d\right)\\ & =\sum _{k=1}^n\mu \left(k\right)\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\left[d|i\right]\left[d|j\right]\mu \left(d\right)\\ & =\sum _{k=1}^n\mu \left(k\right)\sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\mu \left(d\right)\left(\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\left[d|i\right]\right)\left(\sum _{j=1}^{\lfloor \frac{n}{k}\rfloor }\left[d|j\right]\right)\\ & =\sum _{k=1}^n\mu \left(k\right)\sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }\mu \left(d\right)\lfloor \frac{n}{kd}{\rfloor }^2\\ & =\sum _{l=1}^n\lfloor \frac{n}{l}{\rfloor }^2\sum _{k|l}\mu \left(k\right)\mu \left(\frac{l}{k}\right),l=kd\end{array}$$

我们设
$$f\left(l\right)=\sum _{k|l}\mu \left(k\right)\mu \left(\frac{l}{k}\right)$$
则很显然，$f(l)$ 为狄利克雷卷积的形式，又因为 $\mu (x)$ 为积性函数，因此 $f(l)$ 也为积性函数。

因此，我们可以通过 $O(n)$ 的线性筛求出 $f(l)$ 的前 $n$ 项。

由 $f(l)$ 的定义，我们知道，

$f(1)=1,f(p)=-2,f(p^k)=0,k\ge 2,p$ 为素数。

因此 $f(l)$ 的求解代码如下：

bool isnotprime[maxn];
ll prime[maxn];
int func[maxn];
int cnt = 0;
void functionsieve(ll range)
{
    func[1] = 1;
    for(ll i = 2; i <= range; i++)
    {
        if(!isnotprime[i])
        {
            prime[++cnt] = i;
            func[i] = -2;
        }
        for(int j = 1; j <= cnt; j++)
        {
            if(i * prime[j] > range)    break;
            isnotprime[i * prime[j]] = true;
            if(i % prime[j] == 0)
            {
                if((i / prime[j]) % prime[j] == 0)  func[i * prime[j]] = 0;
                else                                func[i * prime[j]] = func[i / prime[j]];
                break;
            }
            else func[i * prime[j]] = -2 * func[i];
        }
    }
}

7.3 优化技巧

在推导的过程中，我们得到了一个积性函数，因此优化的方法是 $O(n)$ 地通过线性筛预处理出积性函数的所有值。

这样，我们便可在 $O(1)$ 的时间内查询到函数的值。

总的时间复杂度为：$O(n)$。

8. BZOJ - 2818

8.1 题面

求 ${\sum }_{i=1}^n{\sum }_{j=1}^n\left[\text{gcd}\left(i,j\right)\text{is a prime}\right]$ 的值。

数据范围：$1\le n\le 10^7$。

8.2 解题过程

$$\begin{array}{rl}& \sum _{i=1}^n\sum _{j=1}^n\left[\text{gcd}\left(i,j\right)\text{is a prime}\right]\\ & =\sum _{p=2}^{maxprime}\sum _{i=1}^n\sum _{j=1}^n\left[\text{gcd}\left(i,j\right)=p\right]\\ & =\sum _{p=2}^{maxprime}\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{p}\rfloor }\left[\text{gcd}\left(i,j\right)=1\right]\\ & =\sum _{p=2}^{maxprime}\sum _{i=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{p}\rfloor }\sum _{k=1}^{\lfloor \frac{n}{p}\rfloor }\left[k|i\right]\left[k|j\right]\cdot \mu \left(k\right)\\ & =\sum _{p=2}^{maxprime}\sum _{k=1}^{\lfloor \frac{n}{p}\rfloor }\mu \left(k\right)\cdot \lfloor \frac{n}{pk}{\rfloor }^2\end{array}$$

注意到，$10^7$ 以内的素数只有 $664579$ 个，而
$$\underset{n\to \infty }{\lim }\sum _{i=1}^n\frac{n}{i}=\log n$$
因此，总的时间复杂度为 $O(664579\cdot \log n)$。

9. BZOJ - 4407

9.1 题面

给定 $n,m,k$，求 ${\sum }_{i=1}^n{\sum }_{j=1}^m\text{gcd}{\left(i,j\right)}^k$ 对 $10^9+7$ 取模之后的值。

共有 $T(1\le T\le 2000)$ 组数据。

数据范围：$1\le n,m,k\le 5\cdot 10^6$。

9.2 解题过程

$$\begin{array}{rl}& \sum _{i=1}^n\sum _{j=1}^m\text{gcd}{\left(i,j\right)}^k\\ =& \sum _{q=1}^{\min \left(n,m\right)}\sum _{i=1}^n\sum _{j=1}^m\left[\text{gcd}\left(i,j\right)=q\right]\cdot q^k\\ =& \sum _{q=1}^{\min \left(n,m\right)}{q}^k\sum _{i=1}^{\lfloor \frac{n}{q}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{q}\rfloor }\left[\text{gcd}\left(i,j\right)=1\right]\\ =& \sum _{q=1}^{\min \left(n,m\right)}{q}^k\sum _{d=1}^{\min \left(\lfloor \frac{n}{q}\rfloor ,\lfloor \frac{m}{q}\rfloor \right)}\mu \left(d\right)\lfloor \frac{n}{qd}\rfloor \lfloor \frac{m}{qd}\rfloor \\ =& \sum _{l=1}^{\min \left(n,m\right)}\lfloor \frac{n}{l}\rfloor \lfloor \frac{m}{l}\rfloor \sum _{q|l}{q}^k\mu \left(\frac{l}{q}\right),l=qd\end{array}$$

设 $f\left(l\right)={\sum }_{q|l}{q}^k\mu \left(\frac{l}{q}\right)$，显然 $f(l)$ 为积性函数。

可以发现，$f(1)=1$，$f(p)=p^k-1$，$f(p^l)=p^{kl}-p^{k(l-1)}=p^{k(l-1)}(p^k-1)$。

因此我们可以利用线性筛求出 $f(l)$ 的前缀和。

之后通过根号分块求解最终答案。

9.3 代码

bool isnotprime[maxn];
int prime[maxn];
ll po[maxn], f[maxn], prefix[maxn];
int cnt = 0;

ll pow_mod(ll a, ll b)
{
    ll res = 1;
    a = a % Mod;
    while(b > 0)
    {
        if(b % 2 == 1)  res = (res * a) % Mod;
        b = b / 2;
        a = (a * a) % Mod;
    }
    return res;
}

void linearsieve(int range)
{
    f[1] = 1;
    prefix[1] = 1;
    for(int i = 2; i <= range; i++)
    {
        if(!isnotprime[i])
        {
            prime[++cnt] = i;
            f[i] = (po[i] - 1 + Mod) % Mod;
        }
        for(int j = 1; j <= cnt; j++)
        {
            if(i * prime[j] > range)    break;
            isnotprime[i * prime[j]] = true;
            if(i % prime[j] == 0)
            {
                f[i * prime[j]] = f[i] * po[prime[j]] % Mod;
                break;
            }
            else f[i * prime[j]] = f[i] * f[prime[j]] % Mod;
        }
        prefix[i] = (prefix[i - 1] + f[i]) % Mod;
    }
}

int main()
{
    int t, k, n, m;
    scanf("%d%d", &t, &k);

    for(int i = 1; i <= 5000000; i++)   po[i] = pow_mod(1LL * i, k);
    linearsieve(5000000);

    while(t--)
    {
        scanf("%d%d", &n, &m);
        ll answer = 0;
        int pos = 0, last = 0;
        int minvalue = min(m, n);
        while(pos + 1 <= minvalue)
        {
            pos = min(m / (m / (pos + 1)), n / (n / (pos + 1)));
            answer = (answer + (prefix[pos] - prefix[last] + Mod) % Mod * (m / pos) % Mod * (n / pos) % Mod) % Mod;
            last = pos;
        }
        printf("%lld\n", answer);
    }
    return 0;
}

10. BZOJ - 2820

10.1 题面

求 ${\sum }_{i=1}^n{\sum }_{j=1}^m\left[\text{gcd}\left(i,j\right)\text{is a prime}\right]$ 的值，有 $T$ 组询问。

数据范围：$1\le n,m\le 10^7$，$1\le T\le 10^4$。

10.2 解题过程

本题为 BZOJ - 2818 的多组询问版本。
∑ i = 1 n ∑ j = 1 m [ gcd ( i , j ) &ThinSpace;&ThinSpace; is a prime ] = ∑ p = 2 m a x p r i m e ∑ i = 1 n ∑ j = 1 m [ gcd ( i , j ) = p ] = ∑ p = 2 m a x p r i m e ∑ i = 1 ⌊ n p ⌋ ∑ j = 1 ⌊ m p ⌋ [ gcd ( i , j ) = 1 ] = ∑ p = 2 m a x p r i m e ∑ i = 1 ⌊ n p ⌋ ∑ j = 1 ⌊ m p ⌋ ∑ k = 1 min ⁡ { ⌊ n p ⌋ , ⌊ m p ⌋ } [ k ∣ i ] [ k ∣ j ] ⋅ μ ( k ) = ∑ p = 2 m a x p r i m e ∑ k = 1 min ⁡ { ⌊ n p ⌋ , ⌊ m p ⌋ } μ ( k ) ⋅ ⌊ n p k ⌋ ⌊ m p k ⌋ = ∑ l = 1 min ⁡ { n , m } ⌊ n l ⌋ ⌊ m l ⌋ ∑ p ∣ l μ ( l p ) , p &ThinSpace;&ThinSpace; is prime \begin{aligned} &amp;\sum_{i=1}^n{\sum_{j=1}^m{\left[ \text{gcd}\left( i,j \right) \,\,\text{is a prime} \right]}} \\ &amp;=\sum_{p=2}^{maxprime}{\sum_{i=1}^n{\sum_{j=1}^m{\left[ \text{gcd}\left( i,j \right) =p \right]}}} \\ &amp;=\sum_{p=2}^{maxprime}{\sum_{i=1}^{\lfloor \frac{n}{p} \rfloor}{\sum_{j=1}^{\lfloor \frac{m}{p} \rfloor}{\left[ \text{gcd}\left( i,j \right) =1 \right]}}} \\ &amp;=\sum_{p=2}^{maxprime}{\sum_{i=1}^{\lfloor \frac{n}{p} \rfloor}{\sum_{j=1}^{\lfloor \frac{m}{p} \rfloor}{\sum_{k=1}^{{\min\left\{\lfloor \frac{n}{p}\rfloor, \lfloor \frac{m}{p}\rfloor\right\}} }{\left[ k|i \right] \left[ k|j \right] \cdot \mu \left( k \right)}}}} \\ &amp;=\sum_{p=2}^{maxprime}{\sum_{k=1}^{\min \left\{ \lfloor \frac{n}{p} \rfloor ,\lfloor \frac{m}{p} \rfloor \right\}}{\mu \left( k \right) \cdot \lfloor \frac{n}{pk} \rfloor \lfloor \frac{m}{pk} \rfloor}} \\ &amp;=\sum_{l=1}^{\min \left\{ n,m \right\}}{\lfloor \frac{n}{l} \rfloor \lfloor \frac{m}{l} \rfloor \sum_{p|l}{\mu \left( \frac{l}{p} \right) ,p\,\,\text{is prime}}} \\ \end{aligned} ​i=1∑n​j=1∑m​[gcd(i,j)is a prime]=p=2∑maxprime​i=1∑n​j=1∑m​[gcd(i,j)=p]=p=2∑maxprime​i=1∑⌊pn​⌋​j=1∑⌊pm​⌋​[gcd(i,j)=1]=p=2∑maxprime​i=1∑⌊pn​⌋​j=1∑⌊pm​⌋​k=1∑min{⌊pn​⌋,⌊pm​⌋}​[k∣i][k∣j]⋅μ(k)=p=2∑maxprime​k=1∑min{⌊pn​⌋,⌊pm​⌋}​μ(k)⋅⌊pkn​⌋⌊pkm​⌋=l=1∑min{n,m}​⌊ln​⌋⌊lm​⌋p∣l∑​μ(pl​),pis prime​
设 $f\left(l\right)={\sum }_{p|l}\mu \left(\frac{l}{p}\right)$，其中 $p$ 为素数。

尽管 $f(l)$ 不是积性函数，我们仍然可以通过线性筛将其前缀和预处理出来。

注意到，$f(1)=0$，$f(p)=1$，$f(p_1{p}_2...p_{n-1}{p}_n^2)=\mu (p_1{p}_2...p_{n-1}{p}_n)$，

设 $l=p_1^{k_1}{p}_2^{k_2}{p}_3^{k_3}...p_n^{k_n}$，若 max ⁡ { k 1 , k 2 , k 3 , . . . , k n } ≥ 3 \max\left\{k_1,k_2,k_3,...,k_n\right\} \ge 3 max{k1​,k2​,k3​,...,kn​}≥3 或者 { k 1 , k 2 , k 3 , . . . , k n } \left\{k_1,k_2,k_3,...,k_n\right\} {k1​,k2​,k3​,...,kn​} 中出现了大于等于 $2$ 个的 $2$，则 $f(l)=0$。

除此之外，若 $a$ 和 $b$ 互质，且 $b$ 为质数，则对于 $l=ab$ 来说，

$(1)$ $b$ 的引入会使 $f(a)$ 在原基础上全部取负值（因为自变量乘上一个原来未出现过的质因子，$\mu (x)$ 会变为原来的相反数）；

$(2)$ $b$ 这一因子本身对 $f(l)$ 的贡献为 $\mu (a)$，因为 $\frac{l}{b}=a$。

因此，在此情况下，$f(l)=-f(a)+\mu (a)$。

所以，我们便可利用线性筛，在 $O(n)$ 的时间复杂度内，求出 $f(l)$ 的前缀和。

之后，我们便可 O ( min ⁡ { n , m } ) O(\sqrt {\min\left\{n,m\right\}}) O(min{n,m} ​) 地求出每一次询问的答案。

总的时间复杂度为 O ( T ⋅ min ⁡ { n , m } ) O(T \cdot \sqrt {\min\left\{n,m\right\}}) O(T⋅min{n,m} ​)。

10.3 线性筛求 $f(l)$ 的前缀和

bool isnotprime[maxn];
ll prime[maxn];
ll f[maxn], prefix[maxn], mu[maxn];
int cnt = 0;
void mobiussieve(ll range)
{
    f[1] = 0;
    mu[1] = 1;
    prefix[1] = 0;
    for(ll i = 2; i <= range; i++)
    {
        if(!isnotprime[i])
        {
            prime[++cnt] = i;
            f[i] = 1;
            mu[i] = -1;
        }
        for(int j = 1; j <= cnt; j++)
        {
            if(i * prime[j] > range)    break;
            isnotprime[i * prime[j]] = true;
            if(i % prime[j] == 0)
            {
                if(i / prime[j] % prime[j] != 0)    f[i * prime[j]] = mu[i];
                else f[i * prime[j]] = 0;
                mu[i * prime[j]] = 0;
                break;
            }
            else
            {
                f[i * prime[j]] = -f[i] + mu[i];
                mu[i * prime[j]] = -mu[i];
            }
        }
        prefix[i] = prefix[i - 1] + f[i];
    }
}

11. Luogu - P3935

11.1 题面

求 ${\sum }_{i=l}^{r}d\left(i\right)$ 对 $998244353$ 取模之后的值。

数据范围：$1\le l\le 10^{14}$，$1\le r\le 1.6\cdot 10^{14}$。

11.2 解题过程

$$\begin{array}{rl}& \sum _{i=1}^{n}d\left(i\right)\\ =& \sum _{i=1}^n\sum _{d|i}1\\ =& \sum _{k=1}^n\sum _{d=1}^{\lfloor \frac{n}{k}\rfloor }1,i=kd\\ =& \sum _{k=1}^n\lfloor \frac{n}{k}\rfloor \end{array}$$

可以通过根号分块得到答案。

12. BZOJ - 3309

12.1 题面

求 ${\sum }_{i=1}^a{\sum }_{j=1}^{b}f(\text{gcd}(i,j))$ 的值。

其中 $f(n)$ 表示 $n$ 所含质因子的最大幂指数。

有 $T$ 组输入，满足 $T\le 10^4$。

数据范围：$1\le a,b\le 10^7$。

12.2 解题过程

$$\begin{array}{rl}& \sum _{i=1}^a\sum _{j=1}^{b}f\left[\text{gcd}\left(i,j\right)\right]\\ =& \sum _{k=1}^{\min \left(a,b\right)}\sum _{i=1}^a\sum _{j=1}^{b}f\left(k\right)\cdot \left[\text{gcd}\left(i,j\right)=k\right]\\ =& \sum _{k=1}^{\min \left(a,b\right)}f\left(k\right)\sum _{i=1}^{\lfloor \frac{a}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{b}{k}\rfloor }\left[\text{gcd}\left(i,j\right)=1\right]\\ =& \sum _{k=1}^{\min \left(a,b\right)}f\left(k\right)\sum _{i=1}^{\lfloor \frac{a}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{b}{k}\rfloor }\sum _{d|\text{gcd}\left(i,j\right)}\mu \left(d\right)\\ =& \sum _{k=1}^{\min \left(a,b\right)}f\left(k\right)\sum _{i=1}^{\lfloor \frac{a}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{b}{k}\rfloor }\sum _{d=1}^{\min \left(\lfloor \frac{a}{k}\rfloor ,\lfloor \frac{b}{k}\rfloor \right)}\mu \left(d\right)\left[d|i\right]\left[d|j\right]\\ =& \sum _{k=1}^{\min \left(a,b\right)}f\left(k\right)\sum _{d=1}^{\min \left(\lfloor \frac{a}{k}\rfloor ,\lfloor \frac{b}{k}\rfloor \right)}\mu \left(d\right)\lfloor \frac{a}{kd}\rfloor \lfloor \frac{b}{kd}\rfloor \\ =& \sum _{l=1}^{\min \left(a,b\right)}\lfloor \frac{a}{l}\rfloor \lfloor \frac{b}{l}\rfloor \sum _{k|l}f\left(k\right)\mu \left(\frac{l}{k}\right)\end{array}$$

设 $g(l)={\sum }_{k|l}f\left(k\right)\mu \left(\frac{l}{k}\right)$，则我们的任务就是将 $g(l)$ 的前缀和筛出来。

此时，我们不妨先打表，求出 $g(l)$ 前几百项的值。

可以发现，$f(1)=0$，$f(p_1^k{p}_2^k...p_l^k)=(-1{)}^{l+1}$, $f(others)=0$。

于是，我们可以先用线性筛筛出 $f(l)=-\mu (l)$（$f(1)=0$）。

之后，枚举所有满足 $f(d)̸=0$ 的 $d$，暴力更新 $f(d^k)=f(d)$。

然后，就可以线性地求出 $f(l)$ 的前缀和。

预处理的时间复杂度为 $O(n+n\log n)$。

最后，我们可以 $O(\sqrt{n})$ 地分块，求解出最终的值。

12.3 代码

bool isnotprime[maxn];
int prime[maxn], mu[maxn];
int prefix[maxn];
int cnt = 0;

void linearsieve(int range)
{
    mu[1] = 0;
    for(int i = 2; i <= range; i++)
    {
        if(!isnotprime[i])
        {
            prime[++cnt] = i;
            mu[i] = 1;
        }
        for(int j = 1; j <= cnt; j++)
        {
            if(i * prime[j] > range)    break;
            isnotprime[i * prime[j]] = true;
            if(i % prime[j] == 0)
            {
                mu[i * prime[j]] = 0;
                break;
            }
            else
            {
                mu[i * prime[j]] = -mu[i];
            }
        }
    }
}

int main()
{
    int t, a, b;

    linearsieve(10000000);
    int tot = 0;
    for(int i = 1; i <= 10000000; i++)
    {
        if(!mu[i])  continue;
        tot++;
        int factor = i;
        while(true)
        {
            mu[factor] = mu[i];
            if(10000000 / factor < i) break;
            factor *= i;
        }
    }
    for(int i = 1; i <= 10000000; i++) prefix[i] = prefix[i - 1] + mu[i];

    scanf("%d", &t);

    while(t--)
    {
        scanf("%d%d", &a, &b);
        ll answer = 0;
        int pos = 0, last = 0;
        int minvalue = min(a, b);
        while(pos + 1 <= minvalue)
        {
            pos = min(a / (a / (pos + 1)), b / (b / (pos + 1)));
            answer += 1LL * (prefix[pos] - prefix[last]) * (a / pos) * (b / pos);
            last = pos;
        }
        printf("%lld\n", answer);
    }
    return 0;
}

13. BZOJ - 3561

13.1 题面

给定 $n,m$，求${\sum }_{i=1}^n{\sum }_{j=1}^m\text{lcm}{\left(i,j\right)}^{\text{gcd}\left(i,j\right)}$ 对 $10^9+7$ 之后的值。

共有 $T(1\le T\le 3)$ 组数据。

数据范围：$1\le n,m\le 5\cdot 10^5$。

13.2 解题过程

$$\begin{array}{rl}& \sum _{i=1}^n\sum _{j=1}^m\text{lcm}{\left(i,j\right)}^{\text{gcd}\left(i,j\right)}\\ =& \sum _{i=1}^n\sum _{j=1}^m{\left(\frac{ij}{\text{gcd}\left(i,j\right)}\right)}^{\text{gcd}\left(i,j\right)}\\ =& \sum _{k=1}^{\min \left(n,m\right)}\sum _{i=1}^n\sum _{j=1}^m{\left(\frac{ij}{k}\right)}^k\left[\text{gcd}\left(i,j\right)=k\right]\\ =& \sum _{k=1}^{\min \left(n,m\right)}\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }{\left(ijk\right)}^k\left[\text{gcd}\left(i,j\right)=1\right]\\ =& \sum _{k=1}^{\min \left(n,m\right)}\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }{\left(ijk\right)}^k\sum _{d=1}^{\min \left(\lfloor \frac{n}{k}\rfloor ,\lfloor \frac{m}{k}\rfloor \right)}\left[d|i\right]\left[d|j\right]\mu \left(d\right)\\ =& \sum _{k=1}^{\min \left(n,m\right)}{k}^k\sum _{d=1}^{\min \left(\lfloor \frac{n}{k}\rfloor ,\lfloor \frac{m}{k}\rfloor \right)}\mu \left(d\right)\sum _{i=1}^{\lfloor \frac{n}{k}\rfloor }{i}^k\left[d|i\right]\sum _{j=1}^{\lfloor \frac{m}{k}\rfloor }{j}^k\left[d|j\right]\\ =& \sum _{k=1}^{\min \left(n,m\right)}{k}^k\sum _{d=1}^{\min \left(\lfloor \frac{n}{k}\rfloor ,\lfloor \frac{m}{k}\rfloor \right)}\mu \left(d\right)\sum _{i=1}^{\lfloor \frac{n}{kd}\rfloor }{\left(di\right)}^k\sum _{j=1}^{\lfloor \frac{m}{kd}\rfloor }{\left(dj\right)}^k\\ =& \sum _{k=1}^{\min \left(n,m\right)}{k}^k\sum _{d=1}^{\min \left(\lfloor \frac{n}{k}\rfloor ,\lfloor \frac{m}{k}\rfloor \right)}\mu \left(d\right)\cdot d^{2k}\sum _{i=1}^{\lfloor \frac{n}{kd}\rfloor }{i}^k\sum _{j=1}^{\lfloor \frac{m}{kd}\rfloor }{j}^k\end{array}$$

我们可以预处理出所有的 $k^k$，$d^{2k}$ 以及 $i^k$ 的前缀和。

注意预处理时应该使用递推的思想，不能使用快速幂，否则会多一个 $\log$ 的复杂度。

计算答案时，第一维枚举 $k$，第二维枚举 $d$，总的时间复杂度为 $O(n\log n)$。

处理方法参见下面代码。

13.3 代码

vector <ll> po[maxn], po2[maxn];
ll pk[maxn];
bool isnotprime[maxn];
ll prime[maxn];
int mu[maxn];
int cnt = 0;

ll pow_mod(ll a, ll b)
{
    ll res = 1;
    a = a % Mod;
    while(b)
    {
        if(b & 1)  res = (res * a) % Mod;
        b = b >> 1;
        a = (a * a) % Mod;
    }
    return res;
}

void mobiussieve(ll range)
{
    mu[1] = 1;
    for(ll i = 2; i <= range; i++)
    {
        if(!isnotprime[i])
        {
            prime[++cnt] = i;
            mu[i] = -1 + Mod;
        }
        for(int j = 1; j <= cnt; j++)
        {
            if(i * prime[j] > range)    break;
            isnotprime[i * prime[j]] = true;
            if(i % prime[j] == 0)
            {
                mu[i * prime[j]] = 0;
                break;
            }
            else mu[i * prime[j]] = (-mu[i] + Mod) % Mod;
        }
    }
}

int main()
{
    int n, m;
    scanf("%d%d", &n, &m);
    int maxvalue = max(n, m);
    mobiussieve(maxvalue);
    for(int i = 1; i <= maxvalue; i++)
    {
        ll sum = 0;
        po[i].pb(0);
        po2[i].pb(0);
        for(int j = 1; j <= maxvalue / i; j++)
        {
            ll temp;
            if(i == 1)  temp = j;
            else        temp = po2[i - 1][j] * 1LL * j % Mod;
            sum = (sum + temp) % Mod;
            po[i].pb(sum);
            po2[i].pb(temp);
        }
        pk[i] = pow_mod(1LL * i, 1LL * i);
    }
    int range = min(n, m);
    ll answer = 0;
    for(int k = 1; k <= range; k++)
    {
        int range2 = min(n / k, m / k);
        ll sum = 0;
        for(int d = 1; d <= range2; d++)
        {
            ll temp1 = 1LL * mu[d] * po2[k][d] % Mod * po2[k][d] % Mod;
            ll temp2 = po[k][n / k / d] * po[k][m / k / d] % Mod;
            sum = (sum + temp1 * temp2 % Mod) % Mod;
        }
        sum = sum * pk[k] % Mod; 
        answer = (answer + sum) % Mod;
    }
    printf("%lld\n", answer);
    return 0;
}

14. BZOJ - 3529

14.1 题面

求 ${\sum }_{i=1}^n{\sum }_{j=1}^m\left[\sigma \left[\text{gcd}\left(i,j\right)\right]\le a\right]\cdot \sigma \left[\text{gcd}\left(i,j\right)\right]$ 对 $2^{31}$ 取模之后的值。

对于每组数据，一共有 $q$ 组询问。

数据范围：$1\le n,m\le 10^5$, $1\le q\le 2\cdot 10^4$,$|a|\le 10^9$。

14.2 解题过程

我们先考虑没有 $a$ 的限制的情况，即求 ${\sum }_{i=1}^n{\sum }_{j=1}^m\sigma \left[\text{gcd}\left(i,j\right)\right]$ 的值。
∑ i = 1 n ∑ j = 1 m σ [ gcd ( i , j ) ] = ∑ k = 1 min ⁡ { n , m } ∑ i = 1 n ∑ j = 1 m σ ( k ) [ gcd ( i , j ) = k ] = ∑ k = 1 min ⁡ { n , m } ∑ p = 1 min ⁡ { ⌊ n k ⌋ , ⌊ m k ⌋ } μ ( p ) ⋅ ⌊ n p k ⌋ ⌊ m p k ⌋ = ∑ l = 1 min ⁡ { m , n } ⌊ n l ⌋ ⌊ m l ⌋ ∑ k ∣ l σ ( k ) μ ( k ) , l = p k \begin{aligned} &amp;\sum_{i=1}^n{\sum_{j=1}^m{\sigma \left[ \text{gcd}\left( i,j \right) \right]}} \\ &amp;=\sum_{k=1}^{\min \left\{ n,m \right\}}{\sum_{i=1}^n{\sum_{j=1}^m{\sigma \left( k \right) \left[ \text{gcd}\left( i,j \right) =k \right]}}} \\ &amp;=\sum_{k=1}^{\min \left\{ n,m \right\}}{\sum_{p=1}^{\min \left\{ \lfloor \frac{n}{k} \rfloor ,\lfloor \frac{m}{k} \rfloor \right\}}{\mu \left( p \right) \cdot \lfloor \frac{n}{pk} \rfloor \lfloor \frac{m}{pk} \rfloor}} \\ &amp;=\sum_{l=1}^{\min \left\{ m,n \right\}}{\lfloor \frac{n}{l} \rfloor \lfloor \frac{m}{l} \rfloor \sum_{k|l}{\sigma \left( k \right) \mu \left( k \right)}},l=pk \end{aligned} ​i=1∑n​j=1∑m​σ[gcd(i,j)]=k=1∑min{n,m}​i=1∑n​j=1∑m​σ(k)[gcd(i,j)=k]=k=1∑min{n,m}​p=1∑min{⌊kn​⌋,⌊km​⌋}​μ(p)⋅⌊pkn​⌋⌊pkm​⌋=l=1∑min{m,n}​⌊ln​⌋⌊lm​⌋k∣l∑​σ(k)μ(k),l=pk​
如果没有 $a$ 的限制，则本题的做法与 $1.2$ 节相同，则其时间复杂度为 $O(q\cdot (\sqrt{m}+\sqrt{n}))$ 。

加上 $a$ 的限制之后，本题就不能像上面那样容易去做了。

我们考虑离线处理所有的查询，之后一起输出答案。

对所有的查询，按照 $a$ 从小到大进行排序。

同时，对于任意的 $x\in [1,10^5]$，我们需要预处理出所有的 $\mu (x)$ 和 $\sigma (x)$，用结构体存储自变量和函数值，之后按照函数值从小到大进行排序。

设 $f\left(x\right)={\sum }_{k|x}\sigma \left(k\right)\mu \left(k\right)\left[\sigma \left(k\right)\le a\right]$，其前缀和为 $prefix[x]$。

对于查询 $i$，$a_i\ge a_{i-1}$，在执行到该查询之前，所有满足 $\sigma (k)\le a_{i-1}$ 的 $\sigma \left(k\right)\mu \left(k\right)$ 均已经被贡献到 $prefix$ 中。

此时，我们仅需将所有满足 $a_{i-1}<\sigma (k)\le a_i$ 的 $\sigma \left(k\right)\mu \left(k\right)$ 贡献到 $prefix$ 中即可，注意要处理 $k$ 的所有小于等于 $10^5$ 的倍数。

$prefix$ 可使用树状数组进行维护。

设 $N=10^5$，则总的时间复杂度为 $O(N{\log }^{2}N+q\sqrt{N}\log N)$。

对于取模问题，直接让 int 自然溢出即可，最终答案若为负数，则直接对 $2^{31}-1 =2147483647 $ 进行与运算即可。

具体实现细节可参考 $14.3$ 节的代码。

14.3 代码

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e5+5;

bool isnotprime[maxn];
int prime[maxn], prefix[maxn];
int mu[maxn];
int cnt = 0;
struct node
{
    int id;
    int value;
} sigma[maxn];
struct data
{
    int id, n, m, a;
} query[maxn];
bool operator < (data a, data b)
{
    return a.a < b.a;
}
bool operator < (node a, node b)
{
    return a.value < b.value;
}
void mobiussieve(int range)
{
    mu[1] = 1;
    for(int i = 2; i <= range; i++)
    {
        if(!isnotprime[i])
        {
            prime[++cnt] = i;
            mu[i] = -1;
        }
        for(int j = 1; j <= cnt; j++)
        {
            if(i * prime[j] > range)    break;
            isnotprime[i * prime[j]] = true;
            if(i % prime[j] == 0)
            {
                mu[i * prime[j]] = 0;
                break;
            }
            else mu[i * prime[j]] = -mu[i];
        }
    }
}

int lowbit(int x)
{
    return x & (-x);
}

void add(int pos, int value)
{
    while(pos <= 100000)
    {
        prefix[pos] += value;
        pos += lowbit(pos);
    }
    return;
}

int querysum(int pos)
{
    int res = 0;
    while(pos >= 1)
    {
        res += prefix[pos];
        pos -= lowbit(pos);
    }
    return res;
}

int ans[maxn];
int main()
{
    int q;
    scanf("%d", &q);

    for(int i = 1; i <= q; i++)
    {
        scanf("%d%d%d", &query[i].n, &query[i].m, &query[i].a);
        query[i].id = i;
    }

    sort(query + 1, query + 1 + q);

    mobiussieve(100000);

    for(int i = 1; i <= 100000; i++)
    {
        for(int j = 1; j <= 100000 / i; j++)
        {
            sigma[i * j].value += i;
        }
    }

    sort(sigma + 1, sigma + 1 + 100000);

    int pos = 0;
    for(int i = 1; i <= q; i++)
    {
        int a = query[i].a;
        int n = query[i].n;
        int m = query[i].m;
        int id = query[i].id;
        while(pos + 1 <= 100000 && sigma[pos + 1].value <= a)
        {
            pos++;
            for(int k = 1; k * sigma[pos].id <= 100000; k++)
            {
                add(k * sigma[pos].id, sigma[pos].value * mu[k]);
            }
        }
        int answer = 0;
        int pos2 = 0, last = 0;
        int range = min(n, m);
        while(pos2 + 1 <= range)
        {
            pos2 = min(n / (n / (pos2 + 1)), m / (m / (pos2 + 1)));
            answer += (n / pos2) * (m / pos2) * (querysum(pos2) - querysum(last));
            last = pos2;
        }
        answer &= 2147483647;
        ans[id] = answer;
    }

    for(int i = 1; i <= q; i++) printf("%d\n", ans[i]);
    
    return 0;
}

15. 2019 ICPC - 南昌邀请赛 - 网络预选赛 - G tsy’s number

15.1 题面

求
$$\sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^n\frac{\phi \left(i\right)\phi \left(j^2\right)\phi \left(k^3\right)}{\phi \left(i\right)\phi \left(j\right)\phi \left(k\right)}\phi \left(\text{gcd}\left(i,j,k\right)\right)$$
对 $2^{30}$ 取模之后的值。

共有 $T(1\le T\le 10^4)$ 组查询，$1\le n\le 10^7$。

15.2 解题过程

15.2.1 Tricks

$(1)$ 对于任意的整数 $m,n$，都有
$$\phi \left(mn\right)=\phi \left(m\right)\phi \left(n\right)\cdot \frac{\text{gcd}\left(m,n\right)}{\phi \left(\text{gcd}\left(m,n\right)\right)}$$
因此，
$$\phi \left(j^2\right)={\phi }^2\left(j\right)\cdot \frac{j}{\phi \left(j\right)}=j\cdot \phi \left(j\right)$$

$$\phi \left(k^3\right)=\phi \left(k\right)\cdot \phi \left(k^2\right)\cdot \frac{k}{\phi \left(k\right)}=k\cdot \phi \left(k^2\right)=k^2\cdot \phi \left(k\right)$$

$(2)$ 幂和公式
$$\sum _{i=1}^{n}i=\frac{n\left(n+1\right)}{2}$$

$$\sum _{i=1}^n{i}^2=\frac{n\left(n+1\right)\left(2n+1\right)}{6}$$

15.2.2 推导过程

$$\begin{array}{rl}& \sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^n\frac{\phi \left(i\right)\phi \left(j^2\right)\phi \left(k^3\right)}{\phi \left(i\right)\phi \left(j\right)\phi \left(k\right)}\phi \left(\text{gcd}\left(i,j,k\right)\right)\\ =& \sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^n\frac{\phi \left(j^2\right)\phi \left(k^3\right)}{\phi \left(j\right)\phi \left(k\right)}\phi \left(\text{gcd}\left(i,j,k\right)\right)\\ =& \sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^n\frac{j\cdot \phi \left(j\right)\cdot k^2\cdot \phi \left(k\right)}{\phi \left(j\right)\phi \left(k\right)}\phi \left(\text{gcd}\left(i,j,k\right)\right)\\ =& \sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^{n}j\cdot k^2\cdot \phi \left(\text{gcd}\left(i,j,k\right)\right)\\ =& \sum _{t=1}^n\sum _{i=1}^n\sum _{j=1}^n\sum _{k=1}^{n}j\cdot k^2\cdot \phi \left(t\right)\left[\text{gcd}\left(i,j,k\right)=t\right]\\ =& \sum _{t=1}^n\sum _{i=1}^{\lfloor \frac{n}{t}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{t}\rfloor }\sum _{k=1}^{\lfloor \frac{n}{t}\rfloor }{t}^3\cdot j\cdot k^2\cdot \phi \left(t\right)\left[\text{gcd}\left(i,j,k\right)=1\right]\\ =& \sum _{t=1}^n{t}^3\phi \left(t\right)\sum _{i=1}^{\lfloor \frac{n}{t}\rfloor }\sum _{j=1}^{\lfloor \frac{n}{t}\rfloor }j\sum _{k=1}^{\lfloor \frac{n}{t}\rfloor }{k}^2\left[\text{gcd}\left(i,j,k\right)=1\right]\\ =& \sum _{t=1}^n\begin{array}{c}{t}^3\phi \left(t\right){\sum }_{i=1}^{\lfloor \frac{n}{t}\rfloor }{\sum }_{j=1}^{\lfloor \frac{n}{t}\rfloor }j{\sum }_{k=1}^{\lfloor \frac{n}{t}\rfloor }{k}^2{\sum }_{d=1}^{\lfloor \frac{n}{t}\rfloor }\mu \left(d\right)\left[d|i\right]\left[d|j\right]\left[d|k\right]\end{array}\\ =& \sum _{t=1}^n\begin{array}{c}{t}^3\phi \left(t\right){\sum }_{i=1}^{\lfloor \frac{n}{t}\rfloor }\left[d|i\right]{\sum }_{j=1}^{\lfloor \frac{n}{t}\rfloor }j\left[d|j\right]{\sum }_{k=1}^{\lfloor \frac{n}{t}\rfloor }{k}^2\left[d|k\right]{\sum }_{d=1}^{\lfloor \frac{n}{t}\rfloor }\mu \left(d\right)\end{array}\\ =& \sum _{t=1}^n\begin{array}{c}{t}^3\phi \left(t\right){\sum }_{d=1}^{\lfloor \frac{n}{t}\rfloor }\mu \left(d\right)\lfloor \frac{n}{td}\rfloor \cdot \frac{\lfloor \frac{n}{td}\rfloor \left(\lfloor \frac{n}{td}\rfloor +1\right)}{2}\cdot \frac{\lfloor \frac{n}{td}\rfloor \left(\lfloor \frac{n}{td}\rfloor +1\right)\left(2\cdot \lfloor \frac{n}{td}\rfloor +1\right)}{6}\end{array}\\ =& \sum _{l=1}^n\lfloor \frac{n}{l}\rfloor \cdot \frac{\lfloor \frac{n}{l}\rfloor \left(\lfloor \frac{n}{l}\rfloor +1\right)}{2}\cdot \frac{\lfloor \frac{n}{l}\rfloor \left(\lfloor \frac{n}{l}\rfloor +1\right)\left(2\cdot \lfloor \frac{n}{l}\rfloor +1\right)}{6}\cdot l^3\sum _{d|l}\mu \left(d\right)\phi \left(\frac{l}{d}\right),l=td\end{array}$$

设 $f\left(l\right)=l^3{\sum }_{d|l}\mu \left(d\right)\phi \left(\frac{l}{d}\right)$，显然 $f(l)$ 为积性函数，我们可以用线性筛 $O(n)$ 地求出其前缀和。

15.3 需要注意的点

为防止向下取整，平方和应该用下面的方法去求：

ll sum2(ll n)
{
    ll tt1 = n, tt2 = n + 1, tt3 = (2LL * n + 1);
    if(tt1 % 2 == 0)        tt1 /= 2;
    else if(tt2 % 2 == 0)   tt2 /= 2;
    else                    tt3 /= 2;
    if(tt1 % 3 == 0)        tt1 /= 3;
    else if(tt2 % 3 == 0)   tt2 /= 3;
    else                    tt3 /= 3;
    return tt1 * tt2 % Mod * tt3 % Mod;
}