【初等数论】欧拉函数

欧拉函数

欧拉函数的定义

欧拉函数 $\phi (n)$ ：对于一个正整数 $n$ ，小于 $n$ 且与 $n$ 互质的正整数的个数。

欧拉函数的积性证明

笛卡尔乘积是指在数学中，两个集合 $X$ 和 $Y$ 的笛卡尔积（Cartesian product），又称直积，表示为 $X\times Y$，第一个对象是X的成员而第二个对象是Y的所有可能有序对的其中一个成员，例如 { 1 , 2 } × { a , b } = { ( 1 , a ) , ( 1 , b ) , ( 2 , a ) , ( 2 , b ) } \{1,2\}\times \{a,b\}=\{(1,a),(1,b),(2,a),(2,b)\} {1,2}×{a,b}={(1,a),(1,b),(2,a),(2,b)}.

简化剩余系也称既约剩余系或缩系，是 $m$ 的完全剩余系中与 $m$ 互素的数构成的子集，如果模m的一个剩余类里所有数都与 $m$ 互素，就把它叫做与模 $m$ 互素的剩余类。

引理：$(a,x)=1$ ， $(b,x)=1$ $⟺$ $(ab,x)=1$ ;

证明：
设 $A$ 为模 $a$ 的既约剩余系，$B$ 为模 $b$ 的既约剩余系，$C$ 为模 $ab$ 的既约剩余系，

定义映射 $f$ : $C\to A\times B$：
$$f(z)=(zmoda,zmodb)$$

因为 $\gcd (a,b)=1$，由中国剩余定理可知，对于任意的 $(x,y)\in A\times B$，存在唯一的 $z\in C$ 使得：
$$z\equiv x(moda),z\equiv y(modb)$$

这表明 $f$ 是双射。

因为 $|C|=\phi (ab)$，$\mid A\times B\mid =\phi (a)\cdot \phi (b)$ 且 $f$ 为双射；

所以$$\phi (ab)=\phi (a)\cdot \phi (b)$$

欧拉函数的性质

性质1：积性函数，$a$ 与 $b$ 互质，则
$$\phi (ab)=\phi (a)\phi (b)$$

性质2：若 $p$ 为素数，则
$$\phi (p)=p-1$$

性质3：若 $p$ 为素数，$k$ 为正整数，则
$$\phi (p^k)=p^k-p^{k-1}=p^{k-1}(p-1)$$

简单证明一下：因为 $p$ 为素数，所以只有当一个数 $q$ 的因数包含 $p$ 时， $q$ 与 $p$ 不互质，所以在小于 $p^k$ 的数中，只有 { p , 2 p , 3 p ⋯ p k − 1 p } \{p,2p,3p\cdots p^{k-1}p\} {p,2p,3p⋯pk−1p} 与 $p$ 不互质，一共 $p^{k-1}$ 个，所以与 $p^k$ 互质的有 $p^k-p^{k-1}$ 。

性质4：质因数分解形式的公式
$$\phi (n)=n\prod _{i=1}^c(1-\frac{1}{p_i})$$

证明： $\phi (n)=\phi (p_1^{k_1})\cdot \phi (p_2^{k_2})\cdots \phi (p_c^{k_c})$

$⟹\phi (n)=(p_1^{k_1}-p_1^{k_1-1})\cdot (p_2^{k_2}-p_2^{k_2-1})\cdots (p_c^{k_c}-p_c^{k_c-1})$

$⟹\phi (n)=p_1^{k_1}(1-\frac{1}{p_1})\cdot p_2^{k_2}(1-\frac{1}{p_2})\cdots p_c^{k_c}(1-\frac{1}{p_c})$

$⟹\phi (n)=n(1-\frac{1}{p_1})\cdot (1-\frac{1}{p_2})\cdots (1-\frac{1}{p_c})$

$⟹\phi (n)=n{\prod }_{i=1}^c(1-\frac{1}{p_i})$

性质5：一个重要的恒等式
$$\sum _{d|n}\phi (d)=n$$

证明：设集合 S = { 1 , 2 , 3 , ⋯   , n } S=\{1,2,3,\cdots,n\} S={1,2,3,⋯,n},

对于每一个因子 $d|n$，定义子集 S d = { k ∈ S   ∣   gcd ⁡ ( k , n ) = d } S_{d}=\{k\in S \space |\space \gcd(k,n)=d\} Sd​={k∈S ∣ gcd(k,n)=d}，

也就是说 $S_d$ 是 $S$ 中与 $n$ 的最大公约数为 $d$ 的数的集合。

根据 $S_d$ 集合定义，可知 $k=d\cdot m(1\le m\le \frac{n}{d})$，

由于 $\gcd (k,n)=d$，可推出 $\gcd (m,\frac{n}{d})=1$，且满足该式的 $m$ 有 $\phi (\frac{n}{d})$ 个，

所以 $|S_d|=\phi (\frac{n}{d})$，

由定义可得，$S$ 是所有 $S_d$ 的并集，所有的 $S_d$ 都不相交，即
$$S=\bigcup _{d|n}{S}_d，S_{d_1}\cap S_{d_2}=\varnothing (d_1\ne d_2)$$

因此
$$n=|S|=\sum _{d|n}|S_d|=\sum _{d|n}\phi (\frac{n}{d})=\sum _{d|n}\phi (d)$$

欧拉函数的求法

单个数的欧拉函数值

由于只需要单个数，可以直接借助质因数分解后的公式（上文的性质4）直接计算，时间复杂度 $O(\sqrt{n})$，代码如下：

ll phi(ll x){ //typedef long long ll;
    ll res=x;
    for(ll i=2;i*i<=x;i++)//由于大于sqrt(n)的质因数最多只有一个，所以只需要枚举到sqrt(n)，再特判即可
        if(x%i==0){ //如果是质因数
            res=res/i*(i-1);//先除后乘，防止爆long long
            while(x%i==0) x/=i;//质因数分解
        }
    if(x!=1) res=res/x*(x-1);//若x最后的值不为1，证明存在比sqrt(n)大的质因数
    return res;
}

埃氏筛打表求值

根据埃氏筛的性质，我们可以利用它，逐步筛去质因子的贡献，直接得到每个数的 $\phi (n)$，时间复杂度与埃氏筛求素数一致，均为 $O(n\log \log n)$，代码如下：

void eratos(){
	for(int i=2;i<maxx;i++){ //maxx是预处理的最大上限
		if(E[i]) continue;//如果E[i]有值，证明它不是素数
		for(int j=i;j<maxx;j+=i){ //筛去质数i的倍数j
			if(!E[j]) E[j]=j;//赋成它自身的值
			E[j]=E[j]/i*(i-1);//公式计算
		}
	}
	return ;
}

欧拉筛打表求值

埃氏筛可降低时间复杂度，欧拉筛当然也可以，我们可以用欧拉函数的积性对欧拉筛进行改造，时间复杂度与欧拉筛求素数一致，均为 $O(n)$，代码如下：

void oula(){
	phi[1]=1; //初始化phi[1]
	for(int i=2;i<=n;i++){
		if(!st[i])//st[i]储存i是否被筛过，!st[i]说明i为素数
			phi[i]=i-1,pri[++cnt]=i;//pri[i]是欧拉筛的素数集合
		for(int j=1;j<=cnt;j++){
			int x=i*pri[j];
			if(x>n) break;
			st[x]=1;
			if(i%pri[j]==0){ //判断pri[j]与i是否互质
				phi[x]=phi[i]*pri[j];//不互质时的计算方式（后面有证明）
				break;//欧拉筛保证每一个合数只被筛一次的核心
			}
			phi[x]=phi[i]*(pri[j]-1);//pri[j]与i互质，直接利用积性
		}
	}
    return ;
}

证明：$\phi (n)=\phi (p_1^{k_1})\cdot \phi (p_2^{k_2})\cdots \phi (p_i^{ki})\cdots \phi (p_c^{k_c})$

$⟹\phi (p_i\cdot n)=\phi (p_1^{k_1})\cdot \phi (p_2^{k_2})\cdots \phi (p_i^{ki+1})\cdots \phi (p_c^{k_c})$

$⟹\phi (p_i\cdot n)=(p_1^{k_1}-p_1^{k_1-1})\cdot (p_2^{k_2}-p_2^{k_2-1})\cdots (p_i^{k_{i+1}}-p_i^{k_i})\cdots (p_c^{k_c}-p_c^{k_c-1})$

$⟹\phi (p_i\cdot n)=p_1^{k_1}(1-\frac{1}{p_1})\cdot p_2^{k_2}(1-\frac{1}{p_2})\cdots p_i^{k_{i+1}}(1-\frac{1}{p_i})\cdots p_c^{k_c}(1-\frac{1}{p_c})$

$⟹\phi (p_i\cdot n)=p_i\cdot n(1-\frac{1}{p_1})\cdot (1-\frac{1}{p_2})\cdots (1-\frac{1}{p_i})\cdots (1-\frac{1}{p_c})$

即： $$\phi (p_i\cdot n)=p_i\cdot n\prod _{i=1}^c(1-\frac{1}{p_i})$$

练习题

洛谷-P2023 Longge 的问题

P2303 [SDOI2012] Longge 的问题

题目背景
Longge 的数学成绩非常好，并且他非常乐于挑战高难度的数学问题。
题目描述
现在问题来了：给定一个整数 $n$，你需要求出 $\sum _{i=1}^n\gcd (i,n)$，其中 $\gcd (i,n)$ 表示 $i$ 和 $n$ 的最大公因数。
输入格式
输入只有一行一个整数，表示 $n$。
输出格式
输出一行一个整数表示答案。
对于 $60\%$ 的数据，保证 $n\le 2^{16}$。
对于 $100\%$ 的数据，保证 $1\le n<2^{32}$。

Solution： 由题意可知 $\gcd (i,n)$ 一定为 $n$ 的一个因数，我们不妨将问题转化为，求对于每一个 $d|n$，从1到 $n$ 有多少个数满足 $\gcd (i,n)=d$ ，只需将每一个 $d$ 乘上相应的次数求和即可；

可能满足上面等式的数一定有 i ∈ { d , 2 d , 3 d , ⋯   , d ⋅ n d } i\in\{d,2d,3d,\cdots,d\cdot \frac{n}{d}\} i∈{d,2d,3d,⋯,d⋅dn​}，我们令 $i=k\cdot d(1\le k\le \frac{n}{d})$，

又因为 $n=d\cdot \frac{n}{d}$，所以只有当 $k$ 与 $\frac{n}{d}$ 互质的时候，$n$ 与 $i$ 不存在更大的公因数，问题又转化成为求所有的 $d\cdot \phi (\frac{n}{d})$ 之和

$$\sum _{i=1}^n\gcd (i,n)=\sum _{d|n}d\cdot \sum _{i=1}^n[\gcd (i,n)=d]=\sum _{d|n}d\cdot \phi (\frac{n}{d})$$

由于 $2^{32}$ 以内数的因数很少，所以可以直接求单个数的欧拉函数值，代码如下：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
namespace zyx_{

ll n,ans;
ll phi(ll x){ //求取单个数的欧拉函数值
    ll res=x;
    for(ll i=2;i*i<=x;i++)
        if(x%i==0){
            res=res/i*(i-1);
            while(x%i==0) x/=i;
        }
    if(x!=1) res=res/x*(x-1); //别忘记大于sqrt(x)的质因数
    return res;
}
int my_main(){
	cin>>n;
	for(ll i=1;i*i<=n;i++)
		if(n%i==0){
			ans+=i*phi(n/i);
			if(i*i!=n) ans+=(n/i)*phi(i);
		}
	cout<<ans;
	return 0;
}

}
int main(){
	zyx_::my_main();
	return 0;
}

洛谷-P2158 仪仗队

P2158 [SDOI2008] 仪仗队

题目描述
作为体育委员，C 君负责这次运动会仪仗队的训练。仪仗队是由学生组成的 $N\times N$ 的方阵，为了保证队伍在行进中整齐划一，C 君会跟在仪仗队的左后方，根据其视线所及的学生人数来判断队伍是否整齐（如下图）。
现在，C 君希望你告诉他队伍整齐时能看到的学生人数。

输入格式
一行，一个正整数 $N$。
输出格式
输出一行一个数，即 C 君应看到的学生人数。
说明/提示
对于 $100\%$ 的数据，$1\le N\le 40000$。

Solution： 我可以建立平面直角坐标系，不难发现，当某同学的坐标 $(x,y)$ 不互质时，其与 C 君 $(0,0)$ 之间的连线上，必定存在更小的正整数对 $(x^{,},y^{,})$，挡住了 C 君的视线，使得 $(x,y)$ 的同学不能被看到，换句话说，该问题转化为求所有满足 $x,y$ 互质的点对的个数，

我们不妨逐层累加点的个数，观察发现，对角线两侧的点对称分布，所以只需将每列的 $\phi (i)$ 乘以2求和即可，当然，不要忘记了对角线上还有一个同学！

#include<bits/stdc++.h>
using namespace std;
const int N=4e4+10;
namespace zyx_{

int n,sum=1,cnt,pri[N],phi[N];
bool st[N];
void oula(){ //欧拉筛打表欧拉函数
    phi[1]=1;
	for(int i=2;i<=n;i++){
		if(!st[i]) pri[++cnt]=i,phi[i]=i-1;
		for(int j=1;j<=cnt;j++){
			int x=pri[j]*i;
			if(x>n) break;
			st[x]=true;
			if(i%pri[j]==0){
                phi[x]=phi[i]*pri[j];
                break;
            }
			phi[x]=phi[i]*phi[pri[j]];
		}
	}
}
int my_main(){
	cin>>n;
	if(n==1){ //注意需要特判
		cout<<0;
		return 0;
	}
	oula();
	for(int i=2;i<=n;i++){
		sum+=phi[i-1]*2;
	}
	cout<<sum;
	return 0;
}

}
int main(){
	zyx_::my_main();
	return 0;
}

洛谷-P2568 GCD

P2568 GCD

题目描述
给定正整数 $n$，求 $1\le x,y\le n$ 且 $\gcd (x,y)$ 为素数的数对 $(x,y)$ 有多少对。
输入格式
只有一行一个整数，代表 $n$。
输出格式
一行一个整数表示答案。
说明/提示
样例输入输出 1 解释
对于样例，满足条件的 $(x,y)$ 为 $(2,2)$，$(2,4)$，$(3,3)$，$(4,2)$。
数据规模与约定
对于 $100\%$ 的数据，保证 $1\le n\le 10^7$。
来源：bzoj2818

Solution： 我们可以设 $n$ 以内的素数集合为 $P$ ,则 $\gcd (x,y)=d(d\in P)$，根据最大公约数的定义，可以得到 $x,y$ 均为 $d$ 的倍数，即 x , y ∈ { d , 2 d , 3 d , ⋯   , d ⋅ n d } x,y\in\{d,2d,3d,\cdots,d\cdot\frac{n}{d}\} x,y∈{d,2d,3d,⋯,d⋅dn​}，不妨设 $x=n\cdot d，y=m\cdot d(1\le x,y\le \frac{n}{d})$，因为 $d$ 为素数，所以 $\gcd (n,m)=1$，

该问题可以转化为，对于每一个 $n$ 以内的素数 $d$ ，$\frac{n}{d}$ 范围内有多少种 $(n,m)$ 的组合，使得 $\gcd (n,m)=1$，求它们的加和，

由于 $n,m$ 地位对等，我们可以固定住 $n$，来看 $m$，根据欧拉函数的定义，所有组合的个数即为1到 $\frac{n}{d}$ 的所有 $\phi (i)$，只需在欧拉筛中加入前缀和即可，对于 $m$ 也同理，但是由于 $(1,1)$ 的存在，最后结果应乘以2再减去1；

a n s = ∑ d ∈ P { 2 [ ∑ i = 1 n d φ ( i ) ] − 1 } ans=\sum_{d\in P}\{2[\sum_{i=1}^{\frac{n}{d}}\varphi(i)]-1\} ans=d∈P∑​{2[i=1∑dn​​φ(i)]−1}

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e7+10;
namespace zyx_{

ll read(){ //快读快写
	ll x=0;int f=1;char c=getchar();
	while(c<'0'||c>'9'){if(c=='-') f=-1;c=getchar();}
	while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+(c^48);c=getchar();}
	return x*f;
}
void print(ll x){
	if(x<0) putchar('-'),x=-x;
	if(x>9) print(x/10);
	putchar(x%10+'0');
	return ;
}
ll n,ans,cnt,pri[N];
ll phi[N],sumphi[N];
bool st[N];
void oula(){ //欧拉筛打表欧拉函数
    sumphi[1]=phi[1]=1;
	for(int i=2;i<=n;i++){
		if(!st[i]) phi[i]=i-1,pri[++cnt]=i;
        sumphi[i]=sumphi[i-1]+phi[i]; //直接计算前缀和
		for(int j=1;j<=cnt;j++){
			int x=i*pri[j];
			if(x>n) break;
			st[x]=1;
			if(i%pri[j]==0){
				phi[x]=phi[i]*pri[j];
				break;
			}
			phi[x]=phi[i]*(pri[j]-1);
		}
	}
    return ;
}
int my_main(){
	n=read();
    oula();
    for(int i=1;i<=cnt;i++)
        ans+=sumphi[n/pri[i]]*2-1;
	print(ans);
	return 0;
}

}
int main(){
	zyx_::my_main();
	return 0;
}

洛谷-P3601 签到题

P3601 签到题

题目背景
这是一道签到题！
建议做题之前仔细阅读数据范围！
题目描述
我们定义一个函数：$\mathrm{qiandao}(x)$ 为小于等于 $x$ 的数中，与 $x$ 不互质的数的个数。
这题作为签到题，给出 $l$ 和 $r$，求出：
$$\sum _{i=l}^r\mathrm{qiandao}(i)\mathrm{mod}666623333$$
输入格式
一行两个整数，$l$、$r$。
输出格式
一行一个整数表示答案。
输入输出样例 #1
输入 #1
233 2333
输出 #1
1056499
输入输出样例 #2
输入 #2
2333333333 2333666666
输出 #2
153096296
说明/提示
对于 $30\%$ 的数据，$l,r\le 10^3$。
对于 $60\%$ 的数据，$l,r\le 10^7$。
对于 $100\%$ 的数据，$1\le l\le r\le 10^{12}$，$r-l\le 10^6$。

Solution： 虽然 $1\le l\le r\le 10^{12}$ 的数据范围过大，但由于 $r-l\le 10^6$ ，我们可以考虑，将区间缩减到 $[1,10^6]$ ，毕竟大于 $\sqrt{n}$ 的质因数最多只有一个，可以先用欧拉筛找出 $10^6$ 以内的所有素数，接着统计各个质数对区间 $[l,.r]$ 的贡献即可，具体详解请看代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=666623333;//模数
const int N=1e6+10;
namespace zyx_{

ll read(){ //快速读取
	ll x=0;int f=1;char c=getchar();
	while(c<'0'||c>'9'){if(c=='-') f=-1;c=getchar();}
	while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+(c^48);c=getchar();}
	return x*f;
}
void print(ll x){ //快速输出
	if(x<0) print('-'),x=-x;
	if(x>9) print(x/10);
	putchar(x%10+'0');
	return ;
}
ll l,r,pri[N],phi[N],val[N];
bool st[N];
ll cnt,ans;
void oula(){ //欧拉筛筛素数
	for(int i=2;i<=N;i++){
		if(!st[i]) pri[++cnt]=i;
		for(int j=1;j<=cnt&&pri[j]*i<=N;j++){
			st[pri[j]*i]=true;
			if(i%pri[j]==0) break;
		}
	}
	return ;
}
int my_main(){
	l=read(); r=read(); oula();
    for(ll i=l;i<=r;i++){
    	phi[i-l]=val[i-l]=i; //初始化欧拉函数值与自身数值
    }
    for(ll i=1;i<=cnt&&pri[i]*pri[i]<=r;i++){
    	ll x=l/pri[i]*pri[i]+pri[i];//找到区间中质数pri[i]的最小的倍数
    	if(l%pri[i]==0) x=l;//特判
    	for(ll j=x;j<=r;j+=pri[i]){
    		phi[j-l]=phi[j-l]/pri[i]*(pri[i]-1);//欧拉函数公式
    		while(val[j-l]%pri[i]==0) val[j-l]/=pri[i];//质因数分解
    	}
    }
    for(ll i=l;i<=r;i++){
    	if(val[i-l]!=1) phi[i-l]=phi[i-l]/val[i-l]*(val[i-l]-1);//特殊处理存在大于sqrt(n)的质因数时的情况
    	(ans+=(i-phi[i-l])%mod)%=mod;//一定不要忘记了取模
    }
    print(ans);
	return 0;
}

}
int main(){
	zyx_::my_main();
	return 0;
}