[NOIP2017考前复习] Euler函数与Mobius函数

本博客中部分文章已被迁移至新博客，但本文由于内容较简单并未迁移且不在迁移计划中。

---

前排提醒：如无特殊说明，本文中所有符号与**《初等数论》（第三版，潘承洞 潘承彪著），北京大学出版社**一书中意义相同

关于$Euler$函数$\phi (x)$

$Euler$函数$\phi (x)$的定义与基本性质

本着定义是重中之重的原则，我们先来回顾一下$Euler$函数的定义：

对于任意正整数$n$，我们记小于$n$且与$n$互素的正整数的个数$\phi (m)$为$m$的Euler函数。即：$$\phi (m)=\sum _{d\le m}[(m,d)=1]$$

---

由定义立即可推出它的一些简单的性质：

1. $\phi (m)=m$的既约剩余系的个数。特别地，对于素数$p,\phi (p)=p-1$
2. 设$m=m_1{m}_2$，若$m_1$和$m_2$的素因子集合相同，那么$\phi (m)=m_2\phi (m_1)$
3. 设$m=m_1{m}_2$，若$(m_1,m_2)=1,$那么$\phi (m)=\phi (m_1)\phi (m_2),$即$Euler$函数$\phi (x)$是积性函数
4. 记$m$的标准分解为$m=\prod p_i^{{\alpha }_i},$则有$\phi (m)=m\prod (1-\frac{1}{p_i})$
5. 对于任意正整数$m$，都有$$m=\sum _{d|m}\phi (d)$$

这些性质都比较显然，在这里不作证明

$\phi (x)$的计算方法

如果只是询问一个数的函数值，最显然的思路就是O($\sqrt{n}$)求解。

inline int phi(int x)
{
	int ans=x;
	for(int i=2;i<=sqrt(x);i++)  if(x%i==0)
	{
		ans-=ans/i;
		while(x%i==0)  x/=i;
	}
	return ans;
}

但更多的时候我们是需要求很多很多的$\phi (x)$的值，这样的做法就显得不是很妙。
我们发现可以在筛素数的时候顺手处理$\phi (x)$的值，于是就有了下面这两种很常用的预处理$\phi (x)$的方法

//Eratosthenes筛法，时间复杂度O(n log log n)
void init()
{
	for(int i=1;i<=maxn;i++)  phi[i]=i;
	for(int i=2;i<=maxn;i++)  if(!vis[i])
	{
		prime[++cnt]=i;  phi[i]--;
		for(int j=i*2;j<=maxn;j+=i)  vis[j]=1,phi[j]-=phi[j]/i;
	}
}

//线性筛法，时间复杂度O(n)
void init()
{
	for(i=2;i<=maxn;i++)
	{
        if(!vis[i])  prime[++cnt]=i,phi[i]=i-1;
        for(int j=1;j<=cnt&&i*prime[j]<=maxn;j++)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j])  phi[i*prime[j]]=phi[i]*(prime[j]-1);
            else
            {
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
        }
    }
}

关于$Mobius$函数$\mu (x)$

$Mobius$函数$\mu (x)$的定义与性质

废话不多说，直接上定义：

对于正整数$n>1$，记其标准分解为$$n=\prod _{i=1}^r{p}_i^{{\alpha }_i}$$则有KaTeX parse error: No such environment: eqnarray at position 7: \begin{̲e̲q̲n̲a̲r̲r̲a̲y̲}̲\mu(x)= \begin{…
特别地，有$\mu (1)=1$

从以上定义，可以推出以下几点性质：

1. 若$\mu(x)\not =0$，则$x$没有比$1$大的平方数因子
2. 若$(m_1,m_2)=1$，则有$\mu(m_1m_2)=\mu(m_1)\mu(m_2)$，即$Mobius$函数$\mu(x)$是积性函数
3. 对于任意正整数$n$，都有$$\begin{eqnarray}\sum_{d|n}\mu(d)=\begin{cases}1&n=1\\0&n\not=1\end{cases}\end{eqnarray}$$
4. 对于任意正整数$n$，都有$$\mu(m)=\sum_{xmodm}e^{\frac {2\pi ix} m}$$这里$xmodm$表示对$m$的任意一组既约剩余系求和

前两个性质比较显然，下面对性质3与性质4做出简略的证明

####性质3：
$n=1$时，命题显然成立。
$n>1$时，记$n$的素因子个数为$r$，则由$Mobius$函数的定义，有$$\sum _{d|n}\mu (d)=C_r^0-C_r^1+\cdots +(-1{)}^r{C}_r^r=(1-1{)}^r=0$$
####性质4：
由复数相关知识，对于任意正整数$c$，显然有KaTeX parse error: No such environment: eqnarray at position 7: \begin{̲e̲q̲n̲a̲r̲r̲a̲y̲}̲ \sum_{i=1}^m e…
记$m$的标准分解为$$m=\prod _{i=1}^r{p}_i^{{\alpha }_i}$$设$m_j=p_j^{{\alpha }_j}(1\le j\le r)$以及$m_j{M}_j=m$，由数论的相关定理，当$x^{(j)}$分别遍历模$m_j$的缩系时，$$x=\sum _{i=1}^r{M}_i{x}^{(i)}$$遍历模$m$的缩系，于是有$$\begin{gathered} S(m)=\sum _{xmodm}{e}^{\frac{2\pi ix}{m}}=\sum _{x^{(1)}mod{m}_1}\cdots \sum _{x^{(r)}mod{m}_r}{e}^{\frac{2\pi i(M_1{x}^{(1)}+M_2{x}^{(2)}+\cdots +M_r{x}^{(r)})}{m}} \\ =S(m_1)S(m_2)\cdots S(m_r) \end{gathered}$$
容易看出KaTeX parse error: No such environment: eqnarray at position 7: \begin{̲e̲q̲n̲a̲r̲r̲a̲y̲}̲S(p_j^{\alpha_j…于是性质4得证。

$\mu (x)$的计算方法

和$Euler$函数一样，如果需要单点求值，只需要$O(\sqrt{n})$分解素因子即可。

inline int miu(int x)
{
    int ans=1;
    for(int i=1;i<=sqrt(x);i++)  if(x%i==0)
    {
        ans=-ans,x/=i;
        if(x%i==0)  return 0;
    }
    return x>1?-ans:ans;
}

至于多点求值，考虑到$\mu (x)$是积性函数，我们可以用线性筛来做到$O(n)$预处理。

void init()
{
    miu[1]=1;
    for(int i=2;i<=maxn;i++)
    {
        if(!vis[i])  prime[++cnt]=i,miu[i]=-1;
        for(int j=1;j<=cnt&&i*prime[j]<maxn;j++)
        {
            int t=i*prime[j];  vis[t]=1;
            if(i%prime[j]==0)  {miu[t]=0;  break;}
            miu[t]=-miu[i];
        }
    }
}

好啦，基础知识先介绍这么多，可以开始做题了。

典型习题

BZOJ2705 欧拉函数

题意：
给定正整数$N(N\le 2^{32})$，求值$$\sum _{i=1}^{N}gcd(i,N)$$
分析：
将$gcd(i,N)$按值分类，问题转化为求满足$(i,N)=d$的$i$的个数。
由于$(i,N)=d$等价于$(\frac{i}{d},\frac{N}{d})=1$,于是所求的个数就是$\phi (\frac{N}{d})$，问题的答案即为$$\sum _{d|N}d\phi (\frac{N}{d})$$
时间复杂度：$O(\tau (N)\sqrt{N})$
下面贴出代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
 
int N,m;
LL ans;
 
inline LL read()
{
    LL ans=0;  char ch=getchar();
    while(ch<'0'||ch>'9')  ch=getchar();
    while(ch>='0'&&ch<='9')  ans=ans*10+ch-'0',ch=getchar();
    return ans;
}
inline LL phi(LL u)
{
    LL ans=u;
    for(int i=2;i<=sqrt(u);i++)  if(u%i==0)
    {
        ans-=ans/i;
        while(u%i==0)  u/=i;
    }
    if(u>1)  ans-=ans/u;
    return ans;
}
 
int main()
{
    N=read(),m=sqrt(N);
    for(int i=1;i<=m;i++)  if(N%i==0)
    {
        ans+=i*phi(N/i);
        if(N/i!=i)  ans+=N/i*phi(i);
    }
    printf("%lld",ans);
    return 0;
}

BZOJ2818 Gcd

题意：
给定$N$，求 { 1 , 2 , ⋯   , N } 2 \{1,2,\cdots,N\}^2 {1,2,⋯,N}2中有多少对$(x,y)$满足$(x,y)$为素数
分析：
还是上一题的套路，$(i,N)=p$等价于$(\frac{i}{p},\frac{N}{p})=1$，不同的是这一次$p$应该取遍所有小于$N$的素数。答案为$$2\sum _{p\le N}\varphi (\left[\frac{N}{p}\right])-\pi (N)$$
其中$$\varphi (x)=\sum _{i=1}^x\phi (i)$$
时间复杂度：$O(n)$(线性筛)或$O(nloglogn)$(埃氏筛)

BZOJ2005 能量采集

这里用到数论中的一个结论：长为$n$宽为$m$的矩形对角线上(含端点)的整点数为$(n,m)+1$
由此可知所求答案即为$$\sum _{i=1}^m\sum _{j=1}^n(2(i,j)-1)$$乍一看貌似和上面所说的并没有太大关系，但我们可以作进一步的代数变形。
我们不妨假设总有$m\le n$
$$\begin{gathered} \sum _{i=1}^m\sum _{j=1}^n(2(i,j)-1)=2\sum _{i=1}^m\sum _{j=1}^n(i,j)-mn \\ =2\sum _{i=1}^m\sum _{j=1}^n\sum _{d|(i,j)}\phi (d)-mn=2\sum _{i=1}^m\sum _{j=1}^n\sum _{d=1}^m[d|i][d|j]\phi (d)-mn \\ =2\sum _{d=1}^m\phi (d)\left(\sum _{j=1}^n[d|j]\left(\sum _{i=1}^m[d|i]\right)\right)-mn \\ =2\sum _{d=1}^m\phi (d)\left[\frac{n}{d}\right]\left[\frac{m}{d}\right]-mn \end{gathered}$$
好了，做完了。
时间复杂度同上一题。

SPOJ-VLATTICE

说白了就是求在集合$[0,1,2,\cdots ,N{]}^3$中有多少个三元组满足$(i,j,k)=1$
因为含有数字0，所以分类讨论。
坐标轴上一共有三个点可视$(0,0,1),(0,1,0),(1,0,0)$，而在坐标平面内，点的个数是$2\varphi (N)-1$，其中$\varphi (x)$的含义同上上题。因此，只需要求$$\sum _{i=1}^N\sum _{j=1}^N\sum _{k=1}^N[(i,j,k)=1]$$
由$Mobius$函数的性质3得$$\sum _{i=1}^N\sum _{j=1}^N\sum _{k=1}^N[(i,j,k)=1]=\sum _{i=1}^N\sum _{j=1}^N\sum _{k=1}^N\sum _{d=1}^N[d|i][d|j][d|k]\mu (d)=\sum _{d=1}^N\mu (d){\left[\frac{N}{d}\right]}^3$$
综上，最终答案为$$Ans(N)=6\varphi (N)+\sum _{d=1}^N\mu (d){\left[\frac{N}{d}\right]}^3$$对于每个询问$O(\sqrt{N})$计算即可。总复杂度$O(T\sqrt{N})$

未完待续。