洛谷-4213 【模板】杜教筛（Sum）

题目描述
给定一个正整数N$(N\le 2^{31}-1)$
求
$an{s}_1={\sum }_{i=1}^n\phi (i)$
$an{s}_2={\sum }_{i=1}^n\mu (i)$
输入格式
一共T+1行 第1行为数据组数T(T<=10) 第2~T+1行每行一个非负整数N，代表一组询问
输出格式
一共T行，每行两个用空格分隔的数ans1,ans2

输入输出样例
输入 #1 复制
6
1
2
8
13
30
2333

输出 #1 复制
1 1
2 0
22 -2
58 -3
278 -3
1655470 2

以下是转的大佬的总结

积性函数
积性函数：对于任意互质的整数 a,ba,b 有 f(ab)=f(a)f(b) 则称 f(x)f(x) 的数论函数。
完全积性函数：对于任意整数 a,ba,b 有 f(ab)=f(a)f(b) 的数论函数。
常见的积性函数：$\phi ,\mu ,\sigma ,d$
常见的完全积性函数：$ϵ,I,id$
这里特殊解释一下 $ϵ,I,id$ 分别是什么意思： $ϵ(n)=[n=1],I(n)=1,id(n)=n$
狄利克雷卷积
设 f, g 是两个数论函数，它们的狄利克雷卷积卷积是：$(f\ast g)(n)=\sum _{d|n}f(d)g(\frac{n}{d})$
性质：满足交换律，结合律
单位元：$ϵ（即f\ast ϵ=f）$
结合狄利克雷卷积得到的几个性质：
$\mu \ast I=ϵ$
$\phi \ast I=id$
$\mu \ast id=\phi$
莫比乌斯反演
若
$g(n)=\sum _{d|n}f(d)$
则
$f(n)=\sum _{d|n}\mu (d)g(\frac{n}{d})$
证明：这里需要用到前面提到的性质：$\mu \ast I=ϵ$
给出的条件等价于 $g=f\ast I$
所以 $g\ast \mu =f\ast I\ast \mu =f\ast ϵ=f$ 即 $g\ast \mu =f$即 结论。
杜教筛
设现在要求积性函数 ff 的前缀和， 设 $\sum _{i=1}^{n}f(i)=S(n)$
再找一个积性函数 g ，则考虑它们的狄利克雷卷积的前缀和
$\sum _{i=1}^n(f\ast g)(i)$
$\begin{array}{rl}& =\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}$
其中得到第一行是根据狄利克雷卷积的定义。
得到第二行则是先枚举 d 提出 g 。
得到第三行则是把 $\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }f(i)$
​ f(i) 替换为 $S(\lfloor \frac{n}{d}\rfloor )$
接着考虑 g(1)S(n) 等于什么。
可以发现，他就等于
$\sum _{i=1}^{n}g(i)S(\lfloor \frac{n}{i}\rfloor )-\sum _{i=2}^{n}g(i)S(\lfloor \frac{n}{i}\rfloor )$
（可以理解成从1开始的前缀和减去从2开始的前缀和就是第一项）
前面这个式子$\sum _{i=1}^{n}g(i)S(\lfloor \frac{n}{i}\rfloor )$
根据刚才的推导，他就等于$\sum _{i=1}^n(f\ast g)(i)$
所以得到杜教筛的核心式子：
$g(1)S(n)=\sum _{i=1}^n(f\ast g)(i)-\sum _{i=2}^{n}g(i)S(\lfloor \frac{n}{i}\rfloor )$
得到这个式子之后有什么用呢？
现在如果可以找到一个合适的积性函数 g ，使得可以快速算出 $\sum _{i=1}^n(f\ast g)(i)$
$(f\ast g)(i)$ 和 g 的前缀和，便可以用数论分块递归地求解。
$(1)\mu$ 的前缀和
考虑到莫比乌斯函数的性质 $\mu \ast I=ϵ$ ，自然想到取 $f=\mu ,g=I,f\ast g=ϵ$ 。
其中 $I$ 的前缀和和 $ϵ$ 的前缀和都弱到爆了。。
$(2)$考虑到 $\phi$ 的性质 $\phi \ast I=id$，取 $f=\phi ,g=I,f\ast g=id$
$f\ast g$ 即 id 的前缀和为 $\frac{n\ast (n+1)}{2}$
$(3)$ （综合）$\sum _{i=1}^n\phi (i)\cdot i$
令 $f=\phi \cdot id,g=id$, 考虑迪利克雷卷积的形式得到 $(f\ast g)(n)=\sum _{d|n}(\phi (d)\cdot d)\cdot (\frac{n}{d})=n\sum _{d|n}\phi (d)=n^2$
即 $(f\ast g)(i)=i^2$
这样就可以快速求得 $(f\ast g)(i)$ 的前缀和 $\frac{n(n+1)(2n+1)}{6}$

​

#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <tr1/unordered_map>
#define re register
using namespace std;
using namespace std::tr1;
inline int read() {
    int X=0,w=1; char c=getchar();
    while (c<'0'||c>'9') { if (c=='-') w=-1; c=getchar(); }
    while (c>='0'&&c<='9') X=X*10+c-'0',c=getchar();
    return X*w;
}
typedef long long LL;
typedef unsigned long long ULL;
const int MAXN=5000000;
int v[MAXN+10];
int p[MAXN+10];
LL mu[MAXN+10],phi[MAXN+10];
inline void init() {
    v[1]=mu[1]=phi[1]=1;
    re int cnt=0;
    for (re int i=2;i<=MAXN;++i) {
        if (!v[i]) p[++cnt]=i,mu[i]=-1,phi[i]=i-1;
        for (re int j=1;j<=cnt&&i*p[j]<=MAXN;++j) {
            v[i*p[j]]=1;
            if (i%p[j]) mu[i*p[j]]=-mu[i],phi[i*p[j]]=phi[i]*phi[p[j]];
            else { mu[i*p[j]]=0,phi[i*p[j]]=phi[i]*p[j]; break; }
        }
    }
    for (re int i=1;i<=MAXN;++i) mu[i]+=mu[i-1],phi[i]+=phi[i-1];
}
unordered_map<int,LL> ansmu,ansphi;
inline LL S_phi(int n) {
    if (n<=MAXN) return phi[n];
    if (ansphi[n]) return ansphi[n];
    LL ans=0;
    for (re int l=2,r;r<2147483647&&l<=n;l=r+1) //特判
        r=n/(n/l),ans+=(r-l+1)*S_phi(n/l);
    return ansphi[n]=(ULL)n*(n+1ll)/2ll-ans; //转ULL避免溢出
}
inline LL S_mu(int n){
    if (n<=MAXN) return mu[n];
    if (ansmu[n]) return ansmu[n];
    LL ans=0;
    for (re int l=2,r;r<2147483647&&l<=n;l=r+1) //同上
        r=n/(n/l),ans+=(r-l+1)*S_mu(n/l);
    return ansmu[n]=1ll-ans;
}
int main() {
    init();
    int T=read();
    while (T--) {
        int n=read();
        printf("%lld %lld\n",S_phi(n),S_mu(n));
    }
    return 0;
}