HDU 6588 Function 题解

题意：给定 $n$，求 $\sum _{i=1}^n\gcd \left(\lfloor \sqrt[3]{i}\rfloor ,i\right)\mathrm{m}\mathrm{o}\mathrm{d}998244353$，其中 $n\le 1\times 10^{21}$。

首先考虑分块，根据 $\sqrt[3]{i}$ 的值可以将 $[1,n]$ 的整数分成 $\sqrt[3]{n}$ 块。记 $m=\lfloor \sqrt[3]{n}\rfloor$，原式可化为 $\left(\sum _{i=1}^m\sum _{j=i^3}^{(i+1{)}^3-1}\gcd (i,j)\right)+\sum _{i=m^3}^n\gcd (m,i)$。首先来考虑固定 $i$，整块的做法。

首先容易注意到，$\gcd (i,j)=\gcd (j,i\mathrm{m}\mathrm{o}\mathrm{d}j)$，因而对于一个固定的 $d$，$\sum _{i=d^3}^{d^3+3d^2+3d}\gcd (d,i)=\sum _{i=d^3}^{d^3+3d^2+3d}\gcd (d,i\mathrm{m}\mathrm{o}\mathrm{d}d)$。而集合 $[d^3,d^3+3d^2+3d]$ 可以拆分为 $3d^2+1$ 个完整的 $\sum _{i=1}^d\gcd (i,d)$ 和最后的一项 $\gcd (d,d^3+3d^2+3d)=d$，因而原式前半部分可化为：
$$\sum _{d=1}^m\sum _{i=d^3}^{d^3+3d^2+3d}\gcd (i,d)=\sum _{d=1}^m\left((3d+1)\sum _{i=1}^d\gcd (d,i)\right)+d$$
对于最后一部分，式子也是形如 $\sum _{i=m^3}^n\gcd (m,i)$，考虑 $\sum _{i=1}^n\gcd (i,a)$ 的求法。
$$\begin{array}{rl}\sum _{i=1}^n\gcd (i,a)=& \sum _{d|a}d\sum _{i=1}^n[\gcd (i,a)=d]\\ =& \sum _{d|a}d\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{ed|(di),e|\frac{a}{d}}\mu (e)\\ \stackrel{\underset{t=ed}{}}{=}& \sum _{d|a}d\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }\sum _{t|di,t|a}\mu \left(\frac{t}{d}\right)\\ =& \sum _{t|a}\lfloor \frac{n}{t}\rfloor \sum _{d|t}d\mu \left(\frac{t}{d}\right)\\ =& \sum _{t|a}\lfloor \frac{n}{t}\rfloor \phi (t)\end{array}$$
若 $a\ne n$，对于单个的 $a$ 可以考虑在 $\mathcal{O}(\sqrt{a})$ 的时间内枚举因子求得。对于大部分的情况——$a=n$，原式可化简为 $\sum _{t|n}\frac{n}{t}\phi (t)=\phi \ast id$，因而是一个积性函数，不妨设为 $f(n)$，考虑欧拉筛线性的求得。

欧拉筛的难点在于 $f(p^k)\to f(p^{k+1})$ 的转化。
$$\begin{array}{rl}f(p^k)=\sum _{t|p^k}\frac{p^k}{t}\phi (t)=& \sum _{i=0}^k{p}^{k-i}\phi (p^i)\\ =& p^k+\sum _{i=1}^k{p}^{k-i}{p}^{i-1}(p-1)\\ =& p+(p-1)\sum _{i=1}^k{p}^{k-1}\\ =& p+(p-1)k{p}^{k-1}\\ =& p(f(p^{k-1})+\phi (p^{k-1}))\end{array}$$

因而维护一个 $\phi$ 数组，最小质因子 $\mathcal{M}$，在该数中最小质因子 $\mathcal{T}$ 的次方数（例如，${\mathcal{T}}_{12}=2^2=4$）即可。对于 $p\nmid n$，则 $f(np)=f(n)f(p)=f(n)(2p-1)$；若 $p|n$，考虑 ${\mathcal{T}}_n=p^i$。
$$\begin{array}{rl}f(np)=& f\left(\frac{n}{p^i}\times p^{i+1}\right)\\ =& f\left(\frac{n}{p^i}\right)f(p^{i+1})\\ =& f\left(\frac{n}{p^i}\right)\times p(f(p^i)+\phi (p^i)))\\ =& p\left[f(n)+f\left(\frac{n}{p^i}\right)\phi (p^i)\right]\end{array}$$
因而同步的转移即可。

#include <bits/stdc++.h>
using namespace std;
const long long mod = 998244353;
const int N = 10000000;
char a[30];
int sum[N + 5];
int f[N + 5], prime[N + 5], phi[N + 5], minfac[N + 5], mintimes[N + 5];
int tot;
long long power(long long a, long long x)
{
    long long ans = 1;
    while(x)
    {
        if (x & 1)
            ans = ans * a % mod;
        a = a % a % mod;
        x >>= 1;
    }
    return ans;
}
long long inv(long long a)
{
    return power(a, mod - 2);
}
inline void sieve(int n)
{
    phi[1] = f[1] = 1;
    for (int i = 2; i <= n; i++)
    {
        if(!minfac[i])
        {
            minfac[i] = prime[++tot] = mintimes[i] = i;
            phi[i] = i - 1;
            f[i] = i + i - 1;
        }
        for (int j = 1; j <= tot && (long long)prime[j] * i <= n; j++)
        {
            long long num = prime[j] * i;
            if (i % prime[j])
            {
                minfac[num] = mintimes[num] = prime[j];
                phi[num] = (long long)phi[i] * (prime[j] - 1);
                f[num] = (long long)f[i] * (2 * prime[j] - 1) % mod;
            }
            else
            {
                minfac[num] = prime[j];
                mintimes[num] = (long long)prime[j] * mintimes[i];
                phi[num] = (long long)phi[i] * prime[j];
                f[num] = (f[i] + (long long)f[i / mintimes[i]] * phi[mintimes[i]] % mod) % mod * (long long)prime[j] % mod;
            }
        }
    }
    for (int i = 1; i <= n; i++)
        sum[i] = ((3ll * i + 3) * f[i] % mod + i + sum[i - 1]) % mod;
}
long long g(long long n, long long m)
{
    long long ans = 0;
    for (int i = 1; i * i <= m;i++)
        if (m % i == 0)
        {
            int j = m / i;
            ans = (ans + phi[i] * (n / i) % mod) % mod;
            if (i != j)
                ans = (ans + phi[j] * (n / j) % mod) % mod;
        }
    return ans;
}
int main()
{
    sieve(N);
    int t;
    scanf("%d", &t);
    while(t--)
    {
        __int128_t n = 0;
        scanf("%s", a);
        int len = strlen(a);
        for (int i = 0; i < len;i++)
            n = n * 10 + a[i] - 48;
        __int128_t left = 0, right = N, lim = 0;
        while (left <= right)
        {
            __int128_t mid = (left + right) >> 1;
            if (mid * mid * mid <= n)
            {
                lim = mid;
                left = mid + 1;
            }
            else
                right = mid - 1;
        }
        long long ans = (((n / lim - lim * lim) % mod * f[lim] % mod + lim) % mod + g(n % lim, lim)) % mod;
        ans = (ans + sum[lim - 1]) % mod;
        printf("%lld\n", ans);
    }
    return 0;
}