欧拉函数笔记

一、积性函数

1 定义

对于两个互质的正整数 $n,m$ ,如果满足：
$$f(mn)=f(m)f(n)$$
那么 $f(x)$ 这个函数就是一个积性函数。
如果对于任意两个正整数 $n,m$ ,满足：
$$f(mn)=f(m)f(n)$$
那么 $f(x)$ 就是一个完全积性函数。

2 积性函数举例

1. 常数函数 $f(n)=1$
2. 单位函数 $f(n)=[n=1]$
   证明： f ( n ) ⋅ f ( m ) = { 0   e l s e 1   n = 1 , m = 1 = f ( n ) ( m ) f(n)\cdot f(m)=\{^{1\ n=1,m=1}_{0\ else}=f(n)(m) f(n)⋅f(m)={0 else1 n=1,m=1​=f(n)(m)
3. 除数函数 ${\sigma }_k(n)={\sum }_{d\mid n}{d}^k$
4. 欧拉函数 $\phi (n)={\sum }_{j=1}^n[gcd(i,n)=1]$，又能写成$\phi (n)=n(1-\frac{1}{p_1})(1-\frac{1}{p_2})\dots (1-\frac{1}{p_n})$，用来求比 n 小且与 n 互质的正整数个数。
5. 莫比乌斯函数 $\mu (n)$
   性质：
   ∑ d ∣ n μ ( d ) = [ n = 1 ] = { 0   e l s e 1   n = 1 \sum_{d|n}\mu(d)=[n=1]=\lbrace^{1\ n=1}_{0\ else} d∣n∑​μ(d)=[n=1]={0 else1 n=1​

二、欧拉函数

以下有几个性质：(p均指质数)

1. 若 $n$ 是质数， $\phi (n)=n-1$ 。
2. 若 $n=p$，$\phi (n)=p^k-k^{k-1}$
3. 若 $n$ 是奇数， $\phi (2n)=\phi (n)$
4. ${\sum }_{d|n}\phi =n$

三、欧拉定理

若a与m互质，则：
$$a^{\varphi (m)}\equiv 1(modm)$$
这就是欧拉定理。欧拉定理、威尔逊定理、费马小定理、中国剩余定理并称为数论四大定理。

如果p是质数。那么:
$$a^{p-1}\equiv 1(modp)$$
刚好我们发现从欧拉定理变成了费马小定理。也就是说,费马小定理是欧拉定理的一个特例。

而当 $a$ 与 $p$ 互质时：
$$a^b\equiv a^{bmod\varphi (p)}(modp)$$
实际上，这是扩展欧拉定理的一部分，也叫欧拉降幂：
$$a^b\equiv \{\begin{array}{c}{a}^{b}b<\varphi (m)\\ a^{(bmod\varphi (m))+\varphi (m)}b\ge \varphi (m)\end{array}modm$$

四、线性筛

1 简介

线性筛（又名欧拉筛），是一种用于找出一定范围内素数的算法。顾名思义，它能够在 $O(n)$ 的时间内完成任务。
线性筛的核心思想是每个合数只会被它的最小质因子筛去一次。我们从小到大枚举每个数，如果这个数是质数，就将它的倍数标记为非质数。同时，每个合数只会被它的最小质因子标记，避免了重复标记。

2 练习

2.1 题目描述

给定一个范围 $[1,n]$ ，有 $q$ 词询问，每次输出第 $k$ 小的数。

2.2 代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
bool vis[100000002];
int prime[10000002],cnt;
int n,q;
void sieve(ll n){
	vis[1]=1;
	for(int i=2;i<=n;i++){
		if(!vis[i]){
			prime[++cnt]=i;
		}
		for(int j=1;j<=cnt&&i*prime[j]<=n;j++){
			vis[i*prime[j]]=1;
			if(i%prime[j]==0){
				break;
			}
		}
	}
}
int main(){
    scanf("%d%d",&n,&q);
    sieve(n);
    while(q--){
    	int k;
    	scanf("%d",&k);
		printf("%d\n",prime[k]);
	}
    return 0;
}

五、欧拉函数实践

1 求欧拉函数

1.1 题目描述

输入 $n$ ，请输出 $\phi (n)$ 的值

1.2 思路

设 $ans$ 等于 $n$ ，代表了小于 $n$ 的所有正整数。当我们发现 $n$ 可以被一个数 $i$ 整除时，我们需要从 $ans$ 中减去那些也能被 $i$ 整除的数。这些数的数量可以通过 $\frac{ans}{i}$ 来计算。

1.3 代码

#include<bits/stdc++.h>
#define ll long long
#define bug printf("---OK---")
#define pa printf("A: ")
#define pr printf("\n")
#define pi acos(-1.0)
using namespace std;
ll n;
ll phi(ll n){
    ll ans=n;
    for(int i=2;i<=sqrt(n);i++){
        if(n%i==0){
            while(n%i==0){//当n能被i整除时，不断地除以i
                n/=i;
            }
            ans-=ans/i;//从结果中减去ans/i
        }
    }
    if(n>1){//如果n大于1，说明n是一个大于sqrt(n)的质数
        ans-=ans/n;//从结果中减去ans/n
    }
    return ans;
}
int main(){
	while(scanf("%lld",&n)&&n!=0){
		printf("%lld\n",phi(n));
	}
    return 0;
}

2 Farey数列

1.1 题目描述

对于一个大于 $1$ 的整数 $n$ ，定义Farey序列 $F_n$ 为一个分数集合 { a / b , 0 < a < b ≤ n \{{a/b},0<a<b\leq n {a/b,0<a<b≤n 且 $\gcd (a,b)=1\}$
例如，前几项为
F 2 = { 1 / 2 } F_2=\{1/2\} F2​={1/2}
F 3 = { 1 / 3 , 1 / 2 , 2 / 3 } F_3=\{1/3,1/2,2/3\} F3​={1/3,1/2,2/3}
F 4 = { 1 / 4 , 1 / 3 , 1 / 2 , 2 / 3 , 3 / 4 } F_4=\{1/4,1/3,1/2,2/3,3/4\} F4​={1/4,1/3,1/2,2/3,3/4}
F 5 = { 1 / 5 , 1 / 4 , 1 / 3 , 2 / 5 , 1 / 2 , 3 / 5 , 2 / 3 , 3 / 4 , 4 / 5 } F_5=\{1/5,1/4,1/3,2/5,1/2,3/5,2/3,3/4,4/5\} F5​={1/5,1/4,1/3,2/5,1/2,3/5,2/3,3/4,4/5}
每次给定 $n$ ，求 $F_n$ 中分数的个数。

1.2 思路

分数没什么用，实际上是求欧拉函数前 $n$ 项的和。

1.3 代码

#include<bits/stdc++.h>
#define ll long long
#define bug printf("---OK---")
#define pa printf("A: ")
#define pr printf("\n")
#define pi acos(-1.0)
using namespace std;
const ll N=1e6+2;
ll t,x;
ll phi[N],prime[N],tot;
void le(int n){//线性筛的同时求欧拉函数，降低时间复杂度
    phi[1]=phi[0]=0;
    for(int i=2;i<=n;i++){
        if(!phi[i]){
        	phi[i]=i-1;
			prime[++tot]=i;
		}
        for(int j=1;j<=tot&&i*prime[j]<=n;j++){
            if(i%prime[j]==0){
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
			else{
				phi[i*prime[j]]=phi[i]*(prime[j]-1);
			}
        }
    }
}
int main(){
	le(N);//预处理
	for(int i=1;i<=N;i++){
		phi[i]+=phi[i-1];//相加
	}
	while(scanf("%lld",&x)&&x!=0){
		printf("%lld",phi[x]);pr;
	}
    return 0;
}

3 【模板】欧拉降幂

3.1 题目描述

求： $a^b(modm)$。
提示： $b\le 10^{200000000}$

3.2 思路

由于b太大了，需要用欧拉降幂的公式。使用字符串处理即可。

3.3 代码

#include<bits/stdc++.h>
#define ll long long
#define bug printf("---OK---")
#define pa printf("A: ")
#define pr printf("\n")
#define pi acos(-1.0)
#define INF 0x3f3f3f3f3f3f3f3f
using namespace std;
ll a,m;
string s;
ll qpow(ll x,ll y){
	ll ans=1;
	while(y>0){
		if(y&1==1){
			ans=(ans*x)%m;
		}
		x=(x*x)%m;
		y>>=1;
	}
	return ans;
}
ll oul(ll n){
    ll ans=n;
    for(ll i=2;i<=sqrt(n);i++){
        if(n%i==0){
            while(n%i==0){
                n/=i;
            }
            ans-=ans/i;
        }
    }
    if(n>1){
        ans-=ans/n;
    }
    return ans;
}
ll downp(ll phi){
	ll len=s.size();
	ll b=0;
	bool flag=0;
	for(ll i=0;i<len;i++){
		b=b*10+s[i]-'0';
		if(b>=phi){
			b%=phi;
			flag=1;
		}
	}
	if(flag==1){
		b+=phi;
	}
	return b;
}
int main(){
	cin>>a>>m>>s;
	cout<<qpow(a,downp(oul(m)));
	return 0;
}