杜教筛

杜教筛

杜教筛可以在 $O(n^{\frac{2}{3}})$ 的复杂度内求解积性函数的前缀和

设数论函数 $f(x)$，求

$$\sum _{i=1}^{n}f(i)$$

令 $S(n)=\sum _{i=1}^{n}f(i)$，找一个数论函数 $g(x)$ 与 $f(x)$ 做卷积

$$\sum _{i=1}^{n}f\ast g=\sum _{i=1}^{n}g\ast f$$

展开卷积

$$\sum _{i=1}^n\sum _{d\mid i}g(d)f(\frac{i}{d})$$

交换求和次序

$$\sum _{d=1}^{n}g(d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }f(i)$$

后半部分其实就是 $S(\lfloor \frac{n}{d}\rfloor )$，换一下变量名得

$$\sum _{i=1}^{n}g(i)S(\lfloor \frac{n}{i}\rfloor )$$

那么就有递推公式

$$g(1)S(n)=\sum _{i=1}^{n}f\ast g-\sum _{i=2}^{n}g(i)S(\lfloor \frac{n}{i}\rfloor )$$

假设 $\sum _{i=1}^{n}f\ast g$ 可以快速求出，后半部分则可以用数论分块求解

莫比乌斯函数前缀和

由狄利克雷卷积

$$\epsilon =\mu \ast I$$

设 $h=\epsilon ,f=\mu ,g=I$ 套用杜教筛

$$S(n)=\sum _{i=1}^n\epsilon (i)-\sum _{i=2}^{n}I(i)S(\lfloor \frac{n}{i}\rfloor )$$

化简可得

$$S(n)=1-\sum _{i=2}^{n}S(\lfloor \frac{n}{i}\rfloor )$$

欧拉函数前缀和

法一:

先对欧拉函数用莫比乌斯反演

$$\sum _{i=1}^n\phi (i)$$

展开欧拉函数

$$\sum _{i=1}^n\sum _{j=1}^n[\gcd (i,j)=1]$$

莫比乌斯反演

$$\sum _{i=1}^n\sum _{j=1}^n\sum _{d\mid \gcd (i,j)}\mu (d)$$

交换求和次序

$$\sum _{d=1}^n\mu (d)\lfloor \frac{n}{d}{\rfloor }^2$$

那么就可以用前面的莫比乌斯函数前缀和来求解了

注意要减去 $i=1,j=1$ 的情况，并且 $/2$

法二：

由狄利克雷卷积

$$Id=\phi \ast I$$

设 $h=Id,f=\phi ,g=I$ 套用杜教筛

$$S(n)=\sum _{i=1}^{n}Id(i)-\sum _{i=2}^{n}I(i)S(\lfloor \frac{n}{i}\rfloor )$$

$$S(n)=\sum _{i=1}^{n}i-\sum _{i=2}^{n}S(\lfloor \frac{n}{i}\rfloor )$$

等差数列求和

$$S(n)=\frac{n(n+1)}{2}-\sum _{i=2}^{n}S(\lfloor \frac{n}{i}\rfloor )$$

代码(两个函数一起求)：

题目：洛谷P4213

#include <bits/stdc++.h>
#define int long long
using namespace std;
const int N = 2e6 + 5;
int T, n, mobius[N], primes[N], cnt, sum[N];
bool st[N];
unordered_map<int, int> mp;
void get_mobius(int n) {
    mobius[1] = 1;
    for (int i = 2; i <= n; i ++) {
        if (!st[i]) {
            primes[cnt ++] = i;
            mobius[i] = -1;
        }
        for (int j = 0; primes[j] * i <= n; j ++) {
            int t = primes[j] * i;
            st[t] = 1;
            if (i % primes[j] == 0) {
                mobius[t] = 0;
                break;
            }
            mobius[t] = -mobius[i];
        }
    }
    for (int i = 1; i <= n; i ++) sum[i] = sum[i - 1] + mobius[i];
}
int sum_mobius(int n) {
    if (n < N) return sum[n];
    if (mp[n]) return mp[n];
    int res = 1;
    for (int l = 2, r; l <= n; l = r + 1) {
        r = n / (n / l);
        res -= sum_mobius(n / l) * (r - l + 1);
    }
    return mp[n] = res;
}
int sum_phi(int n) {
    int res = 0;
    for (int l = 1, r; l <= n; l = r + 1) {
        r = n / (n / l);
        res += (sum_mobius(r) - sum_mobius(l - 1)) * (n / l) * (n / l);
    }
    return (res - 1) / 2 + 1;
}
signed main() {
    get_mobius(N - 1);
    cin >> T;
    while (T --) {
        cin >> n;
        cout << sum_phi(n) << " " << sum_mobius(n) << endl;
    }
}