积性函数的筛法

下面的所有知识均属数论范畴。

---

欧拉筛

可以做到 $O(n)$ 的线性复杂度！

可以筛所有的积性函数（即对于所有的 $a\perp b$ 都满足 $f(ab)=f(a)\ast f(b)$ 的函数)。

举两个例子：

线性筛欧拉函数：

inline void Sieve()
{
   
    phi[1] = 1;
    for(register int i = 2; i <= MAXN; i++)	//MAXN表示值域
    {
   
        if(!NotPrime[i])	//是质数
        {
   
            phi[i] = i - 1;
            Prime[++len] = i;
        }
        for(register int j = 1; j <= len && i * Prime[j] <= MAXN; j++)
        {
   
            NotPrime[i * Prime[j]] = true;
            if(i % Prime[j] == 0)	//此时 Prime[j] 为 i 的最小质因子
			{
   
            	phi[i * Prime[j]] = phi[i] * prime[j];
				break;
            }
            phi[i * Prime[j]] = phi[i] * (Prime[j] - 1);
        }
    }
    return;
}

线性筛莫比乌斯函数：

inline void Sieve()
{
   
    miu[1] = 1;
    for(register int i = 2; i <= MAXN; i++)
    {
   
        if(!NotPrime[i])
        {
   
            miu[i] = -1;
            Prime[++len] = i;
        }
        for(register int j = 1; j <= len && i * Prime[j] <= MAXN; j++)
        {
   
            NotPrime[i * Prime[j]] = true;
            if(i % Prime[j] == 0) break;	//此时 miu 函数值为 0，可以直接跳过
            miu[i * Prime[j]] = -miu[i];
        }
    }
    return;
}

可以发现，每个数最多只会被自己的最小质因子筛到一次，所以时间复杂度是线性的 $O(n)$。

不过要注意，每个数被自己的最小质因子筛到时，对应的 i * Prime[j]仍然是需要计算函数值的！（当然 $\mu$ 函数此时就等于 $0$，可以直接跳过）

---

杜教筛

在非线性时间内求解积性函数的前缀和。

设现在我们要求的是 $f$ 的前缀和，即 $S(n)={\sum }_{i=1}^{n}f(i)$。

我们考虑把 $f$ 卷上一个 $g$ 的前缀和，即：
$$\begin{array}{rl}{S}^{\prime }(n)& =\sum _{i=1}^n(f\ast g)(i)\\ & =\sum _{i=1}^n\sum _{d|i}f(d)g(\frac{i}{d})\\ & =\sum _{d=1}^{n}g(d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }f(i)\\ & =\sum _{d=1}^{n}g(d)S(\lfloor \frac{n}{d}\rfloor )\end{array}$$
然后再推导一下，发现：
$$g(1)S(n)=\sum _{d=1}^{n}g(d)S(\lfloor \frac{n}{d}\rfloor )-\sum _{d=2}^{n}g(d)S(\lfloor \frac{n}{d}\rfloor )$$
即：
$$S(n)=\frac{{\sum }_{i=1}^n(f\ast g)(i)-{\sum }_{d=2}^{n}g(d)S(\lfloor \frac{n}{d}\rfloor )}{g(1)}$$
所以，只要我们找到合适的 $g$，就能够快速求出 $S(n)$ 了。

（这个时候建议看看下文的莫反公式）

对于 $\mu$ 的前缀和，我们发现有 $\mu \ast 1=ϵ$，即定义 $g=1$，代入公式有 $S(n)=1-{\sum }_{i=2}^{n}S(\frac{n}{i})$。

对于 $\phi$ 的前缀和，我们发现有 $\phi \ast 1=id$，即定义