欧拉函数与积性函数

数论记号

• $a|b$ 记作a整除b(b%a==0)
• $[a==1]$ 表示当表达式成立时值为真
• $\varphi :$欧拉函数
• $\mu :$莫比乌斯函数

欧拉函数

算数基本定理

$$N=p_1^{c_1}{p}_2^{c_2}...p_m^{c_m}p\in prime$$

欧拉函数$\varphi (n)$表示的是$[1,n]$之间与n互质的数的数量

即
$$\varphi (n)=\sum _i^n[gcd(n,i)==1]$$
对于素数 P,$\varphi (P)=P-1$.

若合数$C=P1\ast P2$.

对于素数P1,p1的倍数有P1,2P1,3P1,..(C/P1)*P1共计C/P1个,同理,1~C中P2的倍数有C/P2个.

如果我们把这C/P1+C/P2数去掉,那么P1*P2的倍数被排除了两次,需要再加回来一次.

因此贡献是 $\varphi (C)=C-\frac{C}{P1}-\frac{C}{P2}+\frac{C}{P1\ast P2}=C\ast (1-\frac{1}{P1}-\frac{1}{P2}+\frac{1}{P1\ast P2})=C\ast (1-\frac{1}{P1})\ast (1-\frac{1}{P2})$.

当$C=P1\ast P2\ast P3$时,

$\varphi (C)=C\ast (1-\frac{1}{P1}-\frac{1}{P2}-\frac{1}{P3}+\frac{1}{P1\ast P2}+\frac{1}{P1\ast P3}+\frac{1}{P2\ast P3}-\frac{1}{P1\ast P2\ast P3})=C\ast (1-\frac{1}{P1})\ast (1-\frac{1}{P2})\ast (1-\frac{1}{P3})$.

推广后便是
$$\begin{gathered} C=p_1^{c_1}{p}_2^{c_2}...p_m^{c_m} \\ \varphi (C)=C\ast \prod _{质数p|C}(1-\frac{1}{p}) \end{gathered}$$
$1-\frac{1}{p}=\frac{p-1}{p}$

令 $T=P_1^{c_1-1}{P}_2^{c_2-1}...P_m^{c_m-1},F=P_1{P}_2...P_m$.
$$\begin{gathered} \varphi (C)=T\ast F\ast \prod _{质数p|C}(1-\frac{1}{p}) \\ \varphi (C)=T\ast \prod _{质数p|C}(p-1) \end{gathered}$$
对于$C2=C1\ast P$.

若C1中早就存在质数P则,新增P后对C2的影响只是$T$,而不会影响$F$，因为$F$,只是出现质数种类的集合
$$\varphi (C2)=\varphi (C1)\ast P$$
否则,就会影响到$F$,但是却不会影响到$T$,因为其统计的是质数出现次数减一的的贡献,出现一次又减去就相当于不减去.
$$\varphi (C2)=\varphi (C1)\ast (P-1)$$

const int maxn = 1e6+10;
int pri[maxn/8],tot,phi[maxn];
bool vis[maxn];
void init()
{
    phi[1]=1;
    for(int i=2;i<maxn;++i){
        if(!vis[i])pri[++tot]=i,phi[i]=i-1;
        //质数的欧拉函数是质数减一
        for(int j=1;j<=tot&&i*pri[j]<maxn;++j){
            vis[i*pri[j]]=1;
            if(i%pri[j])phi[i*pri[j]]=phi[i]*(pri[j]-1)//新增质数时*(j-1)
            else {
                phi[i*pri[j]]=phi[i]*pri[j]//不是新增直接乘
                break;
            }
        }
    }
}

例题

$$\begin{gathered} 令f(x)=\sum _i^x\varphi (i) \\ 求\sum _d^n\varphi (d)\ast f(n/d) \end{gathered}$$

分块的要点:被除的数为不变值,应计算不被除的数_的前缀和

#include<bits/stdc++.h>
typedef long long ll;
const int maxn = 1e6+10;
int pri[maxn/8],tot,phi[maxn];
ll sum[maxn];
bool vis[maxn];
void init()
{
    phi[1]=1;
    for(int i=2;i<maxn;++i){
        if(!vis[i])pri[++tot]=i,phi[i]=i-1;
        //质数的欧拉函数是质数减一
        for(int j=1;j<=tot&&i*pri[j]<maxn;++j){
            vis[i*pri[j]]=1;
            if(i%pri[j])phi[i*pri[j]]=phi[i]*(pri[j]-1)//新增质数时*(j-1)
            else {
                phi[i*pri[j]]=phi[i]*pri[j]//不是新增直接乘
                break;
            }
        }
    }
    for(int i=1;i<=n;i++)sum[i]=sum[i-1]+phi[i];
}

ll solve(ll n)
{
    ll ans=0;
    for(ll l=1,r;l<=n;l=r+1){
        r=n/(n/l);
        ans += (sum[r]-sum[l-1])*sum[n/l];
    }
    return ans;
}

int main(){
    int n;
    init();
    while(cin>>n)
        cout<<solve(n)<<endl;
    return 0;
}

积性函数

定义一

如果当a,b互质时,有$f(ab)=f(a)\ast f(b)$,那么称$f(x)$是积性函数

如果不当a,b互质时,仍有$f(ab)=f(a)\ast f(b)$,那么称$f(x)$是完全积性函数

定义二

若$f$是积性函数,在算数基本定理中 $n={\prod }_1^m{p}_i^{ci}$,则$f(n)={\prod }_i^{m}f(p_i^{c_i})$.

积性函数性质1

若 $F(X)={\sum }_{d|x}G(d)\ast H(x/d)[gcd(d,x/d)==1]$.若函数$ G,H $都是积性函数,则$F$也是积性函数

$\varphi (x)$函数积性证明:

只要证明$f(n)={\prod }_i^{m}f(p_i^{c_i})$.就可以证明$\varphi (x)$是积性函数

根据欧拉函数的求法,对于任意一个数C,对于C2=C1*P,若P|C1,则$\varphi (C2)=\varphi (C1)\ast P$,否则$\varphi (C2)=\varphi (C1)\ast (P-1)$.

因此,对于$C=p^{c_i}$,$\varphi (C)=(1-\frac{1}{p})\ast p^{c_i}$.

那么对于
$$\begin{gathered} C=p_1^{c_1}{p}_2^{c_2}...p_m^{c_m} \\ \varphi (C)=C\ast \prod _{质数p|C}(1-\frac{1}{p}) \end{gathered}$$
那么显然$\varphi$是积性函数

$\varphi (x)$函数性质1

$$Id=\varphi \ast I=\sum _{d|x}\varphi (d)\ast I(x/d)$$

其中$I(x)=1$,记作不变函数(并且是积性函数),并且根据积性函数性质$Id$也是积性函数,那么有
$$Id(N)=\prod _i^{m}Id(p_i^{c_i})$$
同时,对于$C=p^{c_i}$有
$$Id(C)=\sum _{d|C}\varphi (d)=1+(p-1)+(p-1)\ast p+(p-1)\ast p^2+...+(p-1)\ast p^{c_i}=p^{c_i}=C$$
因此
$$Id(N)=\prod _i^m{p}_i^{c_i}=N\to Id(N)=N$$
实际上这$Id(x)=x$就是自然函数,并且证明了${\sum }_{d|x}\varphi (d)=x$.

性质2

$\forall n>1,1n中与n互质的数的和为n\ast \varphi (n)/2$.

证明:

$gcd(n,i)==1,则gcd(n,n-i)==1$.

使用反证法,若$gcd(n,n-i)==k,则k|n,k|(n-i),那么k|i$.因此成立

由此可以发现互质的数成对出现,并且和为n,共计$\varphi (x)/2$对.所以和为$n\ast \varphi (n)/2$.证毕.