积性函数

原文：http://blog.csdn.net/skywalkert/article/details/50500009
以下是本人整理~

常用公式：

①$\sum_{d|n}{φ(n)}=n\to φ(n)=n-\sum_{d|n,d<n}φ(d)$
②$\sum_{i=1}^{n}[gcd(n,i)=1]*i=\frac{n*φ(n)+[n=1]}{2}$
③$[n=1]=\sum_{d|n}{μ(d)}$
④$n=\sum_{d|n}{φ(d)}$
⑤
若$S(n)=\sum_{i=1}^{n}\phi(i)$
则$S(n)=\frac{n*(n+1)}{2}-\sum_{i=2}^{n}S(\lfloor\frac{n}{i}\rfloor)$

⑥
若$S(n)=\sum_{i=1}^{n}i*\phi(i)$
则$S(n)=\frac{n*(n+1)*(2n+1)}{6}-\sum_{i=2}^{n}i*S(\lfloor\frac{n}{i}\rfloor)$

⑦
若$S(n)=\sum_{i=1}^{n}μ(i)$
则$S(n)=1-\sum_{i=2}^{n}S(\lfloor\frac{n}{i}\rfloor)$

⑧
若$S(n)=\sum_{i=1}^{n}i*μ(i)$
则$S(n)=1-\sum_{i=2}^{n}S(\lfloor\frac{n}{i}\rfloor)$

一些函数定义：
元函数$e(n)=[n=1]$，狄利克雷卷积的乘法单位元，完全积性。
恒等函数$I(n)=1$，完全积性。
单位函数$id(n)=n$，完全积性。
幂函数$id^k(n)=n^k$，完全积性。

积性函数的定义

　　1.若$f(n)$的定义域为正整数域，值域为复数，即$f:Z^+→C$，则称$f(n)$为数论函数（$I$为数论函数集合）。
　　2.若$f(n)$为数论函数，且$f(1)=1$，对于互质的正整数$p,q$有$f(p⋅q)=f(p)⋅f(q)$，则称其为积性函数。
　　3.若$f(n)$为积性函数，且对于任意正整数$p,q$都有$f(p⋅q)=f(p)⋅f(q)$，则称其为完全积性函数。

积性函数的性质与例子

　　　　$σ_k(n)=\sum_{d|n}{d^k}$，表示n的约数的k次幂和
　　　　　　特殊：约数个数$k=0$，约数和$k=1$
　　
　　　　欧拉函数：
　　　　$φ(n)=\sum_{i=1}^{n}[(n,i)=1]=\prod_{i=1}^{k}(1-\frac{1}{p_i})$
　　　　$\sum_{d|n}{φ(n)}=n\to φ(n)=n-\sum_{d|n,d<n}φ(d)$
　　　　表示不大于n且与n互质的正整数个数(n>2时φ(n)为偶数)
　　　
　　　　$\sum_{i=1}^{n}[gcd(n,i)=1]*i=\frac{n*φ(n)+[n=1]}{2}$
　　　　表示不大于n且与n互质的正整数总和

　　　　莫比乌斯函数：
　　　　$μ(n)= \begin{cases} (-1)^t　　　n=\prod_{i=1}^{t}{p_i}\\ 0　　　　　有平方因子 \end{cases}$

　　　　莫比乌斯经典公式：
　　　　$[n=1]=\sum_{d|n}{μ(d)}$　　排列组合后二项式定理转换即可证明
　　　　$n=\sum_{d|n}{φ(d)}$　　　　将$\frac{i}{n}(1≤i≤n)$化为最简分数统计个数即可证明

　　　　若$f(n)$为积性函数，对于$n=\prod_{i=1}^{t}{p_i^{k_i}}$
　　　　则：$f(n)=\prod_{i=1}^{t}{f(p_i^{k_i})}$

　　　　若$f(n)$为完全积性函数，对于$n=\prod_{i=1}^{t}{p_i^{k_i}}$
　　　　则：$f(n)=\prod_{i=1}^{t}{f(p_i)^{k_i}}$

狄利克雷卷积与莫比乌斯反演

　　　　一、
　　　　数论函数$f$和$g$狄利克雷卷积定义为$(f∗g)(n)=\sum _{d|n}{f(d)⋅g(\frac{n}{d})}$。
　　　　狄利克雷卷积满足交换律、结合律，对加法满足分配律。
　　　　存在单位元函数$e(n)=[n=1]$使得$f∗e=f=e∗f$。
　　　　若$f$和$g$为积性函数则$f∗g$也为积性函数。

　　　　二、
　　　　狄利克雷卷积的一个常用技巧是对于积性函数$f$与恒等函数$I$的卷积的处理。比如：
　　　　已知 $n=\prod_{i=1}^{t}{p_i^{k_i}}$
　　　　若　$g(n)=\sum_{d|n}{f(d)}$
　　　　则　$g(n)=\prod_{i=1}^{t}\sum_{j=0}^{k_i}{f(p_i^j)}$

　　　　三、
　　　　莫比乌斯反演也是对于$g(n)=\sum_{d|n}{f(d)}$的讨论。
　　　　☆但是不要求$f$是积性函数☆（适用于已知$g(n)$求$f(n)$的情况）
　　　　$f(n)=\sum_{d|n}{g(d)⋅μ(\frac{n}{d})}$
　　　　证明：
　　　　∵$I∗μ=e$
　　　　∴$g∗μ=f∗I∗μ=f∗e=f$
　　　　
　　　　欧拉函数和莫比乌斯函数之间的关系
　　　　∵$\sum_{d|n}{φ(d)}=id(n)$
　　　　∴$φ(n)=\sum_{d|n}μ(d)\frac{n}{d}$
　　　　∴$\frac{φ(n)}{n}=\sum_{d|n}\frac{μ(d)}{d}$

杜教筛

　　一个简单的例子：
　　　　求前n个正整数的约数之和，即$\sum_{i=1}^{n}{σ(i)}$其中$n\le 10^{12}$
　　推导：
　　　　$\sum_{i=1}^{n}{σ(i)}=\sum_{i=1}^{n}\sum_{j=1}^{n}[j|i]*j=\sum_{i=1}^{n}i*\sum_{j=1}^{n}[i|j]=\sum_{i=1}^{n}{i*\lfloor\frac{n}{i}\rfloor}$
　　当$i\le\sqrt{n}$时$\lfloor\frac{n}{i}\rfloor$只有$\sqrt{n}$个值
　　当$i\le\sqrt{n}$时
　　∵$\lfloor\frac{n}{i}\rfloor\le\sqrt{n}$
　　∴$\lfloor\frac{n}{i}\rfloor$只有$\sqrt{n}$个值
　　对于固定的$\lfloor\frac{n}{i}\rfloor$，$i$的取值是一段连续的区间，是$[\lfloor\frac{n}{\lfloor\frac{n}{i}\rfloor+1}\rfloor+1,\lfloor\frac{n}{\lfloor\frac{n}{i}\rfloor}\rfloor]$

　　练习：求前n个正整数的约数个数之和
　　$\sum_{i=1}^{n}{i*\lfloor\frac{n}{i}\rfloor}=\sum_{i=1}^{n}{\frac{\lfloor\frac{n}{i}\rfloor*(\lfloor\frac{n}{i}\rfloor+1)}{2}}$(一个常见的表达形式)

　　另一个例子：
　　求前n个正整数的欧拉函数之和，即$\sum_{i=1}^{n}φ(i)$，其中$n≤10^{11}$
　　
　　推导：
　　我们利用$\sum_{d|n}{φ(n)}=n$来化简。
　　它的另一个形式：$φ(n)=n-\sum_{d|n,d<n}φ(d)$
　　记$\phi(n)=\sum_{i=1}^{n}{φ(i)}$
　　那么：
　　

$$\phi(n)=\sum_{i=1}^{n}{φ(i)}=\sum_{i=1}^{n}[{i-\sum_{d|i,d<i}φ(d)}]=\frac{n*(n+1)}{2}-\sum_{i=2}^{n}{\sum_{d|i,d<i}φ(d)}=\frac{n*(n+1)}{2}-\sum_{i=2}^{n}\sum_{d}[d|i]*φ(d)\\=\frac{n*(n+1)}{2}-\sum_{\frac{i}{d}=2}^{n}\sum_{d=1}^{\lfloor\frac{n}{\frac{i}{d}}\rfloor}φ(d)=\frac{n*(n+1)}{2}-\sum_{i=2}^{n}\sum_{d=1}^{\lfloor\frac{n}{i}\rfloor}φ(d)=\frac{n*(n+1)}{2}-\sum_{i=2}^{n}\phi(\lfloor\frac{n}{i}\rfloor)$$
　　
　　其中 $$\sum_{i=2}^{n}\sum_{d}[d|i]*φ(d)=\sum_{\frac{i}{d}=2}^{n}\sum_{d=1}^{\lfloor\frac{n}{\frac{i}{d}}\rfloor}φ(d)$$
　　它的思想所在就是：d的条件是要能整除i，无法方便算出。于是我们便枚举 $\frac{d}{i}$的值，这样这些值就是连续的一坨一坨的了。

　　时间复杂度证明：
　　

　　小总结：
　　如果能通过狄利克雷卷积构造一个更好计算前缀和的函数，且用于卷积的另一个函数也易计算，则可以简化计算过程。例如上题就是利用了φ∗I=id的性质，但一定注意，不是所有的这一类题都只用配个恒等函数I就可以轻松完事的，有时需要更细致的观察。
　　先根据要求的函数的一些其他表达式进行简化，将简化出来的式子或原式整体感知一下（明白它在求什么），尝试从其他方式考虑（目标就是降低复杂度，那么就要观察枚举的那些数那些比较小，那些数比较有特性之类的）

　　杜教筛通式：
　　我们要求的是$S(n)=\sum_{i-1}^{n}f(i)$，其中$f\in I$
　　假设$g\in I$且$g(1)\neq0$那么
　　

$$\sum_{i=1}^{n}(f*g)(i)=\sum_{i=1}^{n}\sum_{d|i}g(d)f(\frac{i}{d})=\sum_{i=1}^{n}g(d)\sum_{1\le i\le n,d|i}f(\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)$$
　　即
　　 $$\sum_{i=1}^{n}(f*g)(i)=\sum_{d=1}^{n}g(d)S(\lfloor\frac{n}{d}\rfloor)$$
　　∴ $$g(1)S(n)=\sum_{i=1}^{n}(f*g)(i)-\sum_{i=2}^{n}g(i)S(\lfloor\frac{n}{i}\rfloor)$$
　　如果我们能知道 $g$和$f∗g$，那么我们就能快速知道 $S(n)$，问题就能愉快地解决啦~

　　大白话口胡：
　　上面的函数$g$是我们xjb选的一个数论函数。若$\sum_{i=1}^{n}(f*g)(i)是一个定值时且g(1)\neq 0$那么:
　　

$$g(1)S(n)=\sum_{i=1}^{n}(f*g)(i)-\sum_{i=2}^{n}g(i)S(\lfloor\frac{n}{i}\rfloor)$$
　　假设 $g(1)A，\sum_{i=1}^{n}(f*g)(i)=B$
　　 $AS(n)=B-\sum_{i=2}^{n}g(i)S(\lfloor\frac{n}{i}\rfloor)$
　　（能求的前提是g(i)可以快速知道的哦~）

一坨题

　
　　1.求前n个正整数的约数之和。$\sum_{i=1}^{n}{σ_1(i)}$√见上
　　
　　2.求前n个正整数的约数个数之和。$\sum_{i=1}^{n}{σ_0(i)}$√
　　
　　3.求前n个正整数的欧拉函数之和。$\sum_{i=1}^{n}φ(i)$√见上
　　
　　4.求前n个正整数的μ值之和。$\sum_{i=1}^{n}μ(i)$√
　　
　　5.令$A(n)=\sum_{i=1}^{n}\frac{i}{(n,i)},F(n)=\sum_{i=1}^{n}A(i)$，求$F(n)$
　　
　　6.若$n=\prod_{i=1}^{t}p_i^{k_i}$则$rad(n)=\prod_{i=1}^{t}p_i$
　　求$\sum_{i=1}^{n}\sum_{d|n}rad(d)*φ(\frac{d}{rad(d)})$

　　7.求$\sum_{i=1}^{n}\sum_{j=1}^{n}lcm(i,j)$√

　　8.求$\sum_{i=1}^{n}\sum_{j=1}^{n}gcd(i,j)$√

　　9.求$S(n)=\sum_{i=1}^{n}i*\phi(i)$√

　　10.求$\sum_{i=1}^{n}\sum_{j=1}^{i}\frac{lcm(i,j)}{i}$√

　　11.求$\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)$√

一坨解法

　　2.求前n个正整数的约数个数之和。$\sum_{i=1}^{n}{σ_0(i)}$
　　　　

$$\sum_{i=1}^{n}{σ_0(i)}=\sum_{i=1}^{n}\sum_{d=1}^{n}[d|n]=\sum_{d=1}^{n}\lfloor\frac{n}{d}\rfloor$$
　　
　　4.求前n个正整数的μ值之和。 $S(n)=\sum_{i=1}^{n}μ(i)$
　　　　 $$1=\sum_{i=1}^{n}[i=1]=\sum_{i=1}^{n}\sum_{d|i}μ(d)=\sum_{d=1}^{n}\sum_{i=1}^{\lfloor\frac{n}{i}\rfloor}μ(i)=\sum_{i=1}^{n}S(\lfloor\frac{n}{i}\rfloor)$$
　　　　∴ $S(n)=1-\sum_{i=2}^{n}S(\lfloor\frac{n}{i}\rfloor)$
　　7.求 $\sum_{i=1}^{n}\sum_{j=1}^{n}lcm(i,j)$
　　 $$\sum_{i=1}^{n}\sum_{j=1}^{n}lcm(i,j)=\sum_{i=1}^{n}\sum_{j=1}^{n}\frac{ij}{gcd(i,j)}\\=\sum_{d=1}^{n}d\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}i\sum_{j=1}^{\lfloor\frac{n}{d}\rfloor}j*[gcd(i,j)==1]$$
　　∵ $$\sum_{i=1}^{n}i[gcd(i,n)==1]=\frac{n*\phi(n)}{2}$$
　　∴ $$原式=\sum_{d=1}^{n}d\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}i^2*\phi(i)$$
　　令 $$S(n)=\sum_{i=1}^{n}i^2*\phi(i)$$
　　 $$g(n)=n^2*\phi(n)$$
　　∵ $$n=\sum_{d|n}\phi(d)$$
　　∴ $$(\frac{n*(n+1)}{2})^2=\sum_{i=1}^{n}\sum_{d|i}\phi(d)*i^2=\sum_{i=1}^{n}\sum_{d|i}g(d)*(\frac{i}{d})^2=\sum_{i=1}^{n}\sum_{d=1}^{\lfloor\frac{n}{i}\rfloor}g(d)*i^2$$
　　∴ $$S(n)=(\frac{n*(n+1)}{2})^2-\sum_{i=2}^{n}i^2*S(\lfloor\frac{n}{i}\rfloor)$$

　　8.求$\sum_{i=1}^{n}\sum_{j=1}^{n}gcd(i,j)$
　　

$$\sum_{i=1}^{n}\sum_{j=1}^{n}gcd(i,j)=\sum_{i=1}^{n}\sum_{j=1}^{n}\sum_{d|gcd(i,j)}\phi(d)=\sum_{i=1}^{n}\sum_{j=1}^{n}\sum_{d|i,d|j}\phi(d)\\=\sum_{d=1}^{n}\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{n}{d}\rfloor}\phi(d)=\sum_{d=1}^{n}\phi(d)*\lfloor\frac{n}{d}\rfloor*\lfloor\frac{n}{d}\rfloor$$
　　转换为问题3

　　9.求$S(n)=\sum_{i=1}^{n}i*\phi(i)$
　　

$$\frac{n(n+1)(2n+1)}{6}=\sum_{i=1}^{n}i^2=\sum_{i=1}^{n}\sum_{d|i}d*\phi(d)*\frac{i}{d}\\=\sum_{i=1}^{n}i\sum_{d=1}^{\lfloor\frac{n}{i}\rfloor}d*\phi(d)=\sum_{i=1}^{n}i*S(\lfloor\frac{n}{i}\rfloor)=S(n)+\sum_{i=2}^{n}i*S(\lfloor\frac{n}{i}\rfloor)$$

　　10.求$\sum_{i=1}^{n}\sum_{j=1}^{i}\frac{lcm(i,j)}{i}$
　　

$$\sum_{i=1}^{n}\sum_{j=1}^{i}\frac{lcm(i,j)}{i}=\sum_{i=1}^{n}\sum_{j=1}^{i}\frac{j}{gcd(i,j)}=\sum_{d=1}^{n}\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{i}j*[gcd(i,j)=1]\\=\sum_{d=1}^{n}\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\frac{i*\phi(i)}{2}=\sum_{i=1}^{n}\frac{i*\phi(i)}{2}*\lfloor\frac{n}{i}\rfloor$$

　　11.求$\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)$
　　不妨设$n\le m$
　　

$$\sum_{i=1}^{n}\sum_{j=1}^{m}lcm(i,j)=\sum_{d=1}^{n}d\sum_{i=1}^{\lfloor\frac{n}{d}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{d}\rfloor}ij*[gcd(i,j)=1]$$
　　令 $F(x,y)=\sum_{i=1}^{x}\sum_{j=1}^{y}ij*[gcd(i,j)=1]$
　　∴ $$原式=\sum_{d=1}^{n}d*F(\lfloor\frac{n}{d}\rfloor，\lfloor\frac{m}{d}\rfloor)$$
　　分析 $F(x,y)$
　　 $$F(x,y)=\sum_{i=1}^{x}\sum_{j=1}^{y}ij*[gcd(i,j)=1]=\sum_{i=1}^{x}\sum_{j=1}^{y}ij\sum_{d|gcd(i,j)}μ(d)\\=\sum_{d=1}^{x}μ(d)*d^2\sum_{i=1}^{\lfloor\frac{x}{d}\rfloor}i\sum_{j=1}^{\lfloor\frac{y}{d}\rfloor}j\\=\sum_{d=1}^{x}μ(d)*d^2*\frac{\lfloor\frac{x}{d}\rfloor*(\lfloor\frac{x}{d}\rfloor+1)*\lfloor\frac{y}{d}\rfloor*(\lfloor\frac{y}{d}\rfloor+1)}{4}$$
　　代入原式：
　　 $$原式=\sum_{d=1}^{n}d\sum_{dd=1}^{\lfloor\frac{n}{d}\rfloor}μ(dd)*dd^2*\frac{\lfloor\frac{n}{d*dd}\rfloor*(\lfloor\frac{n}{d*dd}\rfloor+1)*\lfloor\frac{m}{d*dd}\rfloor*(\lfloor\frac{m}{d*dd}\rfloor+1)}{4}$$
　　预处理以后分块套分块。也可以不预处理，之间分块套分块套杜教筛

//问题4----51nod1244
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>

using namespace std;
typedef long long ll;
const int N=100100;
int p[N],v[N],mu[N],tot;
ll l,r,ans,sum[N];
map<ll,ll> mp;

ll gao(ll n){
    if (n<N) return sum[n];
    if (mp.find(n)!=mp.end()) return mp[n];
    ll res=1,i,nx;
    for (i=2;i<=n;i=nx+1){
        nx=n/(n/i);
        res=res-(nx-i+1)*gao(n/i);
    }
    return mp[n]=res;
}

void pre(){
    for (int i=2;i<N;i++){
        if (!v[i]) p[++tot]=i,mu[i]=-1;
        for (int j=1;j<=tot && i*p[j]<N;j++){
            v[i*p[j]]=1;
            if (i%p[j]==0) break;
            mu[i*p[j]]=-mu[i];
        }
    }
    mu[1]=1;
    for (int i=1;i<N;i++)
        sum[i]+=sum[i-1]+mu[i];
}

int main(){
    freopen("d.in","r",stdin);
    scanf("%lld%lld",&l,&r);pre();
    printf("%lld\n",gao(r)-gao(l-1));
}