杜教筛

杜教筛

求积性函数$f(x)$的前缀和$S(n) = \sum_{i=1}^n f(i)$

狄利克雷卷积：

$$(g*f)(x) = \sum_{d|x} g(d) f(\dfrac{x}{d})$$

$$\sum_{x = 1}^n (g*f)(x) = \sum_{x=1}^n \sum_{d|x} g(d) f(\dfrac{x}{d})$$

$$= \sum_{d = 1}^n \sum_{d|x}^n g(d) f(\dfrac{x}{d})$$

$$= \sum_{d=1}^n g(d) \sum_{x=1}^{n/d}f(x)$$

$$= \sum_{d=1}^n g(d) S(\lfloor \dfrac{n}{d} \rfloor)$$

即

$$\sum_{x=1}^n (g*f)(x) = \sum_{d=1}^n g(d) S(\lfloor \dfrac{n}{d} \rfloor)$$

我们有

$$g(1)S(n) = \sum_{i=1}^n g(i)S(\lfloor \dfrac{n}{i} \rfloor) - \sum_{i=2}^n g(i)S(\lfloor \dfrac{n}{i} \rfloor)$$

$$g(1)S(n) = \sum_{i=1}^n (g*f)(i) - \sum_{i=2}^n g(i)S(\lfloor \dfrac{n}{i} \rfloor)$$

---

假设我们要求$S(n) = \sum_{i=1}^n \mu(i)$

我们只需要找出$g(x)$使得$(g*f)(x)$的前缀和好算就可以了
我们知道

$$\sum_{d|n} \mu(d) = [n = 1]$$

那么我们可以设$g(x) = 1$
这样的话

$$(g*f)(x) = \sum_{d|x} f(d) g(\dfrac{d}{x}) = [x = 1]$$

$$\sum_{i=1}^x (g*f)(x) = 1$$

$$g(1)S(n) = 1 - \sum_{i=2}^n S(\lfloor \dfrac{n}{i} \rfloor)$$

后面的部分可以数论分块
首先先线性筛出前面一部分（$1000000$）$\mu(x)$的前缀和，后面记忆化搜索即可

---

再假设我们要求$S(n) = \sum_{i=1}^n \phi(i)$

我们知道

$$\sum_{d|n} \phi(d) = n$$

$$(g*f)(x) = \sum_{d|x} f(d)g(\dfrac{x}{d})$$

我们发现如果令$g(x) = 1$

$$(g*f)(x) = \sum_{d|x} f(d) = x$$

$$g(1)S(n) = \sum_{i=1}^n i - \sum_{i=2}^n g(i)S(\lfloor\dfrac{n}{i} \rfloor)$$

$$S(n) = \dfrac{x(x+1)}{2} - \sum_{i=2}^n S(\lfloor\dfrac{n}{i} \rfloor)$$

---

51nod1238 最小公倍数之和V3

求

$$\sum_{i=1}^n \sum_{j=1}^n lcm(i, j)$$

---

$$ans = \sum_{i=1}^n \sum_{j=1}^n \dfrac{ij}{gcd(i, j)}$$

$$= \sum_{d=1}^n \dfrac{1}{d} \sum_{i=1}^n \sum_{j=1}^n ij \ [gcd(i, j) = d]$$

$$= \sum_{d=1}^n \dfrac{1}{d} \sum_{i=1}^{n/d} \sum_{j=1}^{n/d} id \ jd \ [gcd(i, j) = 1]$$

$$= \sum_{d=1}^n d \sum_{i=1}^{n/d} \sum_{j=1}^{n/d} ij \ [gcd(i, j) = 1]$$

设

$$f(x) = \sum_{i=1}^x\sum_{j=1}^xij[gcd(i,j)=1]$$

$$= 2\Big(\sum_{i=1}^x i \sum_{j=1}^i j[gcd(i, j) = 1]\Big) - 1$$

$$= 2 \Big( \sum_{i=1}^x i \ \dfrac{i\phi(i) +[i=1]}{2} \Big) - 1$$

$$= \Big( \sum_{i=1}^x i^2 \ \phi(i) +[i=1] \Big) - 1$$

$$= \sum_{i=1}^x i^2 \ \phi(i)$$

代回原式子

$$ans = \sum_{d=1}^n d \sum_{i=1}^{n/d} i^2 \ \phi(i)$$

设

$$S(n) = \sum_{i=1}^n i^2 \ \phi(i)$$

$$S(n) = \sum_{i=1}^n (g*f)(i) - \sum_{i=2}^n g(i)S(\dfrac{n}{i})$$

推一下卷积

$$\sum_{i=1}^n (g*f)(x) = \sum_{i=1}^n \sum_{d|i} f(d)g(\dfrac{i}{d})$$

$$= \sum_{i=1}^n \sum_{d|i} d^2 \phi(d) g(\dfrac{i}{d})$$

我们发现设$g(x) = x^2$就可以把$d^2$约掉

$$= \sum_{i=1}^n \sum_{d|i} d^2 \phi(d) \dfrac{i^2}{d^2}$$

$$= \sum_{i=1}^n i^2 \sum_{d|i} \phi(d)$$

$$= \sum_{i=1}^n i^3$$

$$= \Big( \dfrac{n(n+1)}{2}\Big)^2$$

$$g(1)S(n) = \sum_{i=1}^n (g*f)(x) - \sum_{i=2}^n g(i)S(\dfrac{n}{i})$$

$$S(n) = \Big( \dfrac{n(n+1)}{2}\Big)^2 - \sum_{i=2}^n i^2 S(\dfrac{n}{i})$$

我们知道

$$\sum_{i=1}^n i^2 = \dfrac{n(n+1)(2n+1)}{6}$$

我们可以杜教筛$S(n)$

$$ans = \sum_{d=1}^n d \ S(\dfrac{n}{d})$$