欧拉定理与扩展欧拉定理

欧拉函数

欧拉函数的求法等请看这里，本篇只讲述推论

本文大部分来自知乎Pecco，供自己复习使用

1、若$p$是质数则$\phi (p^n)=p^{n-1}(p-1)$

在不大于$p^n$的正整数中，只有$p、2p、3p、...、p\ast p^{n-1}$总共$p^{n-1}$个所以$\phi (p^n)=p^n-p^{n-1}=p^{n-1}(p-1)$

2、若$a\mid x$，则$\phi (ax)=a\phi (x)$

因为$a$是$x$的约数，所以$a$的所有质因数都被$x$包含，根据欧拉函数的定义易得结论$\phi (ax)=a\phi (x)$

3、若$a、b$互质，则$\phi (ab)=\phi (a)\ast \phi (b)$

因为$gcd(a,b)=1$，所以$a$和$b$的所有质因数都不相同，根据欧拉函数定义易得$\phi (ax)=a\phi (x)$，且符合这种性质的数论函数被称为积性函数

欧拉函数可以单个求，也可以使用埃氏筛法、和线性筛求，上边链接里有，这里不说了

欧拉定理

若$a$与$m$互质，则$a^{\phi (m)}\equiv 1(modm)$

当$m$是质数时，欧拉定理退化成费马小定理$a^{m-1}\equiv 1(modm)$

重要推论：当$a$与$m$互质时，$a^b\equiv a^{bmod\phi (m)}(modm)$

因为$a^b=a^{\phi (m)\lfloor b/\phi (m)\rfloor +bmod\phi (m)}\equiv 1\ast a^{bmod\phi (m)}$

---

那么$a$和$m$不互质时我们怎么办呢？

这时我们引入扩展欧拉定理

若$b\ge \phi (m)$，则$a^b\equiv a^{bmod\phi (m)+\phi (m)}(modm)$

当$b<\phi (m)$时，直接用快速幂计算即可

当$b\ge \phi (m)$时，分两种情况讨论

1、$a$与$m$互质，

$a^b\equiv a^{bmod\phi (m)}(modm)\ast 1\equiv a^{bmod\phi (m)+\phi (m)}(modm)$
显然成立

2、$a$与$m$不互质

我们先把m分解质因数，即$m={p_1}^{q_1}{p_2}^{q_2}...{p_n}^{q_n}$

我们只需要证明对于每个${p_i}^{q_i}$都有
$a^b\equiv a^{bmod\phi (m)+\phi (m)}(mod{p_i}^{q_i})$

假设$m_1、m_2$互质，且$x\equiv y(mod{m}_1)，x\equiv y(mod{m}_2)，则x\equiv y(mod{m}_1{m}_2)$

因为$x-y\mid m_1$且$x-y\mid m_2$，所以两者可以合并

现在开始分类讨论${p_i}^{q_i}$

1、$gcd(a,{p_i}^{q_i})=1$

上面已经证明，根据欧拉定理必然有
$a^b=a^{\phi (m)\lfloor b/\phi (m)\rfloor +bmod\phi (m)}\equiv a^{bmod\phi (m)+\phi (m)}(mod{p_i}^{q_i})$

因为欧拉函数是积性函数，所以$\phi (m)\lfloor b/\phi (m)\rfloor$必然是$\phi ({p_i}^{q_i})$的倍数

2、$gcd(a,{p_i}^{q_i})\ne 1$

因为$p_i$是一个质数，所以$a$必然是$p_i$的倍数，设$a=n{p}_i$(以下$p_i$用$p$表示)

由欧拉定理可知

$\phi ({p_i}^{q_i})={p_i}^{q_i-1}(p_i-1)$，易知 $\phi ({p_i}^{q_i})\ge q_i$

粗略证明：
$\because p_i是一个质数$
$\therefore p_i\ge 2，p_i-1>0$

$已知q_i\ge 1$(都是质因数的幂了，这不废话吗 )

$现在问题已经转化成了证明a^{x-1}\ge x(x\ge 1)$

$动手画一下函数就知道了$

$综上所述，\phi ({p_i}^{q_i})\ge q_i$

则$b\ge \phi (m)\ge \phi ({p_i}^{q_i})\ge q_i$

所以${p_i}^{q_i}$既是$a^{bmod\phi (m)+\phi (m)}$的因数，又是$a^b$的因数

所以$a^b\equiv a^{bmod\phi (m)+\phi (m)}\equiv 0(mod{p_i}^{q_i})$

综上所述：$a^b\equiv a^{bmod\phi (m)+\phi (m)}(modm)(b\ge \phi (m))$

P5091 【模板】扩展欧拉定理

具体代码如下

#include<bits/stdc++.h>

using namespace std;

typedef long long LL;

LL qmi(LL a, LL b, LL p)
{
	LL res = 1;
	while(b)
	{
		if(b & 1) res = res * a % p;
		a = a * a % p;
		b >>= 1;
	}
	return res;
}

int main()
{
	LL a, m;
	cin >> a >> m;
	
	LL phi = m, temp = m; 
	
	for(LL i = 2; i * i <= temp; ++ i)
		if(temp % i == 0)
		{
			phi = phi / i * (i - 1);
			while(temp % i == 0) temp /= i;
		}
	
	if(temp > 1)
	{
		phi = phi / temp * (temp - 1);	
	}	
	
	LL b = 0;
	char c;
	bool flag = false;
	while(!isdigit(c = getchar()));
	for(; isdigit(c); c = getchar())
	{
		b = b * 10 + (c - '0');
		if(b >= phi)
		{
			flag = true;
			b %= phi;
		}
	}
	
	if(flag) b += phi;
	LL ans = qmi(a, b, m);
	
	cout << ans << endl;
	
	return 0;
}