自己的一点（求积性函数前n项和）杜教筛总结（其实也是瞎抄的）

杜教筛可以解决哪一类的问题：

例如题目给你一个函数，让你求它的前项和 ${\sum }_{i=1}^{n}f(i)$ ，但是这个n很大（n=1e10），显然在秒数级是不能算出来的，之前我们有学过拉格朗日插值定理求前n项和，它的复杂度与这个函数最高阶k紧密相连，差不多为O(k)，但是你要想用拉格朗日插值法，你得知道这个函数f为多项式函数才行。比如让你求前n项欧拉函数和，莫比乌斯前n项和，你就用不上这个拉格朗日插值法了。

此时，为了解决这个问题，杜教筛就出来了，它就是用来求这些的。我们先大概说下他的原理：我们先在毫秒复杂度的基础上预处理前T项和的结果（T<n），然后前n项和我们通过我们构造的函数表达式，用整数分块来去求解，此时复杂度会降为O（$n^{\frac{2}{3}}$）左右。


杜教筛也就是积性函数前n项和的参考博客：
神犇1
神犇2
下面写的如果有跟网上雷同，请不要介意，我只是写来给我以后回顾总结用的。

先认识下几个函数先：


函数f和g的卷积为：$(f\ast g)(n)={\sum }_{d|n}f(d)\ast g(\frac{n}{d})$
元函数 $ϵ(n)$，满足 $ϵ(n)=[n==1]$
恒等函数 $I(n)$ ，满足恒等为1
单位函数 $id(n)$，满足 $id(n)=n$

常用公式：

1，$[n==1]={\sum }_{d|n}u(d)$
2，${\sum }_{i=1}^n[gcd(n,i)=1]\ast i=\frac{n\ast \phi (n)+[n==1]}{2}$
3，$\phi (n)={\sum }_{d|n}u(d)\frac{n}{d}$
4，$\frac{\phi (n)}{n}={\sum }_{d|n}\frac{u(d)}{d}$
5，$n={\sum }_{d|n}\phi (d)$

杜教筛的通用表达式：
已知我们要求函数f的前n项和，我们设 $S(n)={\sum }_{i=1}^{n}f(i)$，我们要找到函数g，并且函数f与函数g的卷积为 $h=f\ast g$，那么我们有如下表达式（证明见参考博客链接）：
$g(1)S(n)={\sum }_{i=1}^{n}h(i)-{\sum }_{d=2}^{n}g(d)\ast S(\lfloor \frac{n}{d}\rfloor )$




通常我们会拿恒等函数$I$给g去和f卷积。

例如1:求莫比乌斯前n项和:${\sum }_{i=1}^{n}u(i)$
$h=u\ast I={\sum }_{d|n}u(d)=ϵ$（莫比乌斯函数性质）
故最终式子就为：$S(n)=1-{\sum }_{d=2}^{n}S(\lfloor \frac{n}{d}\rfloor )$


例如2：求欧拉函数前n项和：${\sum }_{i=1}^n\phi (i)$
我们函数拿恒等函数给g去和f卷积，
$h=\phi \ast I={\sum }_{d|n}\phi (d)=id(n)=n$ (欧拉函数性质)
故最终的式子就为：$S(n)={\sum }_{i=1}^{n}i-{\sum }_{d=2}^{n}S(\lfloor \frac{n}{d}\rfloor )$


例如2：题目链接：哆啦A梦传送门
题意：已知 $n^2-3\ast n+2={\sum }_{d|n}f(d)$
求 ${\sum }_{i=1}^{n}f(i)mod(1e9+7)$ ( $n\le 1e9$)


这题是很模板的一道题，很显然我们也拿恒等函数$I(n)$去和f卷积，得：
$h=f\ast I={\sum }_{d|n}f(d)=n^2-3\ast n+2$


故最终的式子就为：
$S(n)={\sum }_{i=1}^n(i^2-3\ast n+2)-{\sum }_{d=2}^{n}S(\lfloor \frac{n}{d}\rfloor )=\frac{n\ast (n-1)\ast (n-2)}{3}-{\sum }_{d=2}^{n}S(\lfloor \frac{n}{d}\rfloor )$

我们设 $g(n)={\sum }_{d|n}f(i)=i^2-3\ast i+2$，由莫比乌斯反演可得：
$f(n)={\sum }_{d|n}u(d)g(\frac{n}{d})$
我们先预处理函数f前600000项来。
代码：

#include<bits/stdc++.h>
#include<tr1/unordered_map>
using namespace std;

typedef long long LL;

const LL mod=1e9+7;
LL n;
#define N 600000

LL sum[N+10],f[N+10];
bool vis[N+10];
int mu[N+10],tot,prime[N+10];

tr1::unordered_map<LL,LL> w;

void init()
{
    tot=0;
    mu[1]=1;
    for(int i=2;i<=N;i++)
    {
        if(!vis[i]){
            prime[++tot]=i;
            mu[i]=-1;
        }
        for(int j=1;j<=tot&&prime[j]*i<=N;j++)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0) break;
            else mu[i*prime[j]]=-mu[i];
        }
    }
    LL item,t;
    for(int d=1;d<=N;d++)///枚举n的因子d
    {
        for(int n=d;n<=N;n+=d){ ///计算f(n)
                t=(LL)(n/d);
            item=(t*t%mod-3*t%mod+2+mod)%mod;///计算g(n/d)的值
            f[n]=(f[n]+item*mu[d]%mod+mod)%mod;

        }
    }

    ///预处理
    sum[0]=0;
    for(int i=1;i<=N;i++)
        sum[i]=(sum[i-1]+f[i])%mod;
}
int num=0;
LL djs(LL x)
{
    ///在预处理的范围内，直接返回
    if(x<=600000) return sum[x];
    if(w[x]) return w[x]; ///以及算好，直接返回

    LL t=x%mod;
    LL ans=t*(t-1)%mod*(t-2)%mod*333333336;
    for(LL l=2,r;l<=x;l=r+1)///整数分块
    {
        r=x/(x/l);
        ans=(ans-(r-l+1)*djs(x/l)%mod+mod)%mod;
        num++;
    }

    return w[x]=ans;
}

int main()
{

    init();

    int ncase;
    scanf("%d",&ncase);

    while(ncase--)
    {
        scanf("%lld",&n);
        printf("%lld\n",djs(n));

    }
    return 0;
}

我的标签：做个有情怀的程序员。