【模板】欧拉函数

https://oi-wiki.org/math/euler/

费马小定理

如果 $p$ 是一个质数，对于任意小于 $p$ 的正整数 $a$ 均有 ： $a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

运用1：模意义下，$a$ 的逆元 $=a^{p-2}(\mathrm{m}\mathrm{o}\mathrm{d}p)$

$\because a^{p-1}=a\times a^{p-2}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$

$\therefore \frac{1}{a}\equiv a^{p-2}(\mathrm{m}\mathrm{o}\mathrm{d}p)$

ps: 当 $p$ 不是质数时，需要用 扩展卢卡斯

二次探测定理

若 $p$ 为素数，$a^2\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$ 小于 $p$ 的正整数解只有 $1$ 和 $p-1$

满足费马小定理和二次探测定理的数可以确定是素数

扩展：

求解 $a^2=n\mathrm{m}\mathrm{o}\mathrm{d}p$ -【模板】 二次剩余

欧拉函数

欧拉函数、欧拉定理

对正整数n，欧拉函数是小于或等于n的正整数中与n互质的数的数目

设 $N=p_1^{a_1}\ast p_2^{a_2}\ast ...\ast p_n^{a_n}$，则

$$\phi (N)=N(1-\frac{1}{p_1})(1-\frac{1}{p_2})...(1-\frac{1}{p_n})$$

若n是质数p的k次幂，$\phi (n)=p^k-p^{k-1}=(p-1)p^{k-1}$

当n为奇质数时 $\phi (2n)=\phi (n)$

当n为质数时 $\phi (n)=n-1$

若 $a$ $b$ 互质时，$\phi (ab)=\phi (a)\phi (b)$

扩展：

$n$ 的所有因子之和 $=(p_1^0+p_1^1+p_1^2...+p_1^{a_1})\ast (p_2^0+p_2^1+p_2^2...+p_2^{a_2})...(p_n^0+p_n^1+p_n^2...+p_n^{a_n})$

欧拉筛

当 i 为质数时，可以直接求得 $\phi (i)=i-1$

当 i 不为质数时，其通过它的质因子得到 $\phi$ 值，

通过任意数 i，其倍数 i*prime[j] 的$\phi$ 值，有两种情况：

1. 当i%prime[j]==0，$\phi (i\ast prime[j])=\phi (i)\ast prime[j]$
2. 当i%prime[j]!=0，$\phi (i\ast prime[j])=\phi (i)\ast (prime[j]-1)$

相关证明 - AcWing 874 筛法求欧拉函数

欧拉定理

1. 当 $(a,n)$ 互质时，可以得出 $a^{\phi (n)}=1(\mathrm{m}\mathrm{o}\mathrm{d}n)$
2. 设三角形的外接圆半径为$R$，内切圆半径为 $r$，外心与内心的距离为 $d$，则 $d^2=R^2-2Rr$
3. 如果一个联通平面图 $G$ 有个 $V$ 顶点、$E$ 条边、$F$ 个面，那么 $V-E+F=2$
4. $e^{i\theta }=\cos \theta +i\sin \theta$

欧拉降幂

$$a^b\mathrm{m}\mathrm{o}\mathrm{d}p=a^{b\mathrm{m}\mathrm{o}\mathrm{d}\phi (p)+\phi (c)}\mathrm{m}\mathrm{o}\mathrm{d}p$$

---

板子

欧拉函数

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int n, m, K;

int phi(int n) { //欧拉函数 板子
    int res = n;
    for (int i = 2; i * i <= n; ++i) {
        if (n % i == 0) {
            res = res / i * (i - 1);//防止溢出
            while (n % i == 0)n /= i;
        }
    }
    if (n > 1) res = res / n * (n - 1);//质数
    return res;
}

int main() {
    ios::sync_with_stdio(0);

    int T;
    cin>>T;
    for (int cs = 1; cs <= T; ++cs) {
        cin >> n;
        cout << phi(n) << endl;
    }
    return 0;
}

欧拉筛

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 1e6 + 10;
int n, m, K;

int phi[N], prime[N], cnt = 0;
bool vis[N];

void get_euler(int n) {//欧拉筛
    phi[1] = 1;
    for (int i = 2; i <= n; ++i) {
        if (!vis[i]) {
            prime[cnt++] = i;
            phi[i] = i - 1;
        }
        for (int j = 0; i * prime[j] <= n; ++j) {
            vis[i * prime[j]] = true;
            if (i % prime[j] == 0) {
                phi[i * prime[j]] = phi[i] * prime[j];
                break;
            }
            phi[i * prime[j]] = phi[i] * (prime[j] - 1);
        }
    }
}

int main() {
    ios::sync_with_stdio(0);
    cin >> n;
    get_euler(n);
    ll res = 0;
    for (int i = 1; i <= n; ++i) {
        res += phi[i]; //求 1~n的欧拉数之和
    }
    cout <<res << endl;
    return 0;
}

欧拉降幂

ll phi(ll n) { //欧拉函数 板子
    int res = n;
    for (ll i = 2; i * i <= n; ++i) {
        if (n % i == 0) {
            res = res / i * (i - 1);//防止溢出
            while (n % i == 0)n /= i;
        }
    }
    if (n > 1) res = res / n * (n - 1);//质数
    return res;
}

ll qpow(ll a, ll b) {//快速幂
    a %= mod;
    ll res = 1;
    while (b) {
        if (b & 1)res = res * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return res;
}


ll euler(ll a, ll b,ll p) {// a^b % p
    return qpow(a, ((b % phi(p)) + phi(p));
}