杜教筛模板题

链接

一般求法：

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

找出两个积性函数$h(x)$和$g(x)$，满足$g(x)=h(x)\ast f(x)$，（*表示卷积）

然后有${\sum }_{i=1}^{n}g(x)={\sum }_{i=1}^n{\sum }_{d|i}h(d)\cdot f(\frac{n}{d})$

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

记$s(x)={\sum }_{i=1}^{x}f(x)$

${\sum }_{i=1}^{n}g(x)={\sum }_{d=1}^{n}h(d)\cdot s(\lfloor \frac{n}{d}\rfloor ))$

右边提取第一项$h(1)\ast s(1)$，把剩余的移到左边等式变为：

$h(1)\ast s(n)={\sum }_{i=1}^{n}g(x)-{\sum }_{d=2}^{n}h(d)\cdot s(\lfloor \frac{n}{d}\rfloor ))$

通常${\sum }_{i=1}^{n}g(x)$和$h(1)$很容易快速得到，所以只要解决${\sum }_{d=2}^{n}h(d)\cdot s(\lfloor \frac{n}{d}\rfloor ))$这部分就可以得到$s(n)$
这部分可以使用分块快速得到

一般是先对可打表的范围预处理，然后求解的时候如果n较小，直接输出答案即可；n较大时采用递归求解的方式来解决。需要保存那些不在打表范围的结果，也就是记忆化。

模板题中：
求${\sum }_{i=1}^n\varphi (i)$
容易想到$\varphi \ast I=id$，其中$I(x)=1,id(x)=x$

求${\sum }_{i=1}^n\mu (i)$
容易想到$\mu \ast I=e$，其中$I(x)=1,e(x)=(x=1?1:0)$

后面照着过程一步一步下来很容易可以解出

常用数论函数：

函数
$d(n)={\sum }_{d\Vert n}1$	因数个数
$\sigma (n)={\sum }_{d\Vert n}d$	因数和
$e(n)=n==1?1:0$	幺元函数
$I(n)=1$	恒等函数
$\varphi (n)$	欧拉函数
$\mu (n)$	莫比乌斯函数
$id(n)=n$	单位函数
$idk(n)=n^k$	幂函数

$\mu \ast 1=ϵ$
$\varphi \ast 1=Id$
$\varphi =Id\ast \mu$
$d=1\ast 1$
$1=\mu \ast d$

参考：

#include <bits/stdc++.h>
#include <unordered_map>
using namespace std;
typedef long long ll;
const int N=5e6+5;
const int M=5e6;
int mu[N];
ll phi[N];
bool vis[N];
int prime[N],tot;

void work_pre(int maxn)
{
    phi[1]=mu[1]=1;
    for(int i=2;i<=maxn;i++)
    {
        if(!vis[i])
        {
            prime[++tot]=i;
            mu[i]=-1;phi[i]=i-1;
        }
        for(int j=1;j<=tot&&prime[j]*i<=maxn;j++)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else mu[i*prime[j]]=-mu[i],phi[i*prime[j]]=phi[i]*(prime[j]-1);
        }
    }
    for(int i=1;i<=maxn;i++){
        phi[i]+=phi[i-1];
        mu[i]+=mu[i-1];
    }
}


unordered_map<int,int>mp1;
unordered_map<int,ll>mp2;

int djs_mu(int n){
    if(n<=M){
        return mu[n];
    }
    if(mp1[n]){
        return mp1[n];
    }
    int ret=1;
    for(int l=2,r;n>=l;l=r+1){
        r=n/(n/l);
        ret-=(r-l+1)*djs_mu(n/l);
    }
    return mp1[n]=ret;
}

ll djs_phi(int n){
    if(n<=M){
        return phi[n];
    }
    if(mp2[n]){
        return mp2[n];
    }
    ll ret=1ll*n*(n+1)/2;
    for(int l=2,r;n>=l;l=r+1){
        r=n/(n/l);
        ret-=1ll*(r-l+1)*djs_phi(n/l);
    }
    return mp2[n]=ret;
}


int main(){
    work_pre(M);
    int t;
    scanf("%d",&t);
    for(int i=1;i<=t;i++){
        int n;
        scanf("%d",&n);
        printf("%lld %d\n",djs_phi(n),djs_mu(n));
    }
    return 0;
}