欧拉定理 & 扩展欧拉定理

这篇具有很好参考价值的文章主要介绍了欧拉定理 & 扩展欧拉定理。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

观前提醒:「文章仅供学习和参考,如有问题请在评论区提出」

目录
  • 前置
    • 剩余类(同余类)
    • 完全剩余系(完系)
    • 简化剩余系(缩系)
    • 欧拉函数
  • 欧拉定理
  • 扩展欧拉定理
  • 参考资料

前置


剩余类(同余类)


给定一个正整数 \(n\) ,把所有的整数根据\(n\) 的余数 \(r\in [0, n - 1]\) 分为 \(n\) 类,每一类就可以被表示为 \(C_{r} = nx + r\) 。那么这类数所构成的集合就称为\(n\) 的剩余类


完全剩余系(完系)


给定一个正整数 \(n\) ,有 \(n\) 个不同的模 \(n\) 的剩余类(因为余数 \(r\in [0, n - 1]\) )。

从这 \(n\) 个不同的剩余类中各取出一个元素,总共 \(n\) 个数,将这些数构成一个新的集合,则称这个集合为\(n\) 的完全剩余系

例如:\(n = 5\) 时,\(\{ 0, 1, 2, 3, 4 \}\)\(\{ 5, 1, -3, 8, 9 \}\) 都是一个模 \(5\) 的完全剩余系。

因为他们都有 \(5\) 个不同的模 \(5\) 的剩余类 \(r \in [0, 4]\)


简化剩余系(缩系)


给定一个正整数 \(n\) ,有 \(\varphi(n)\) 个不同的模 \(n\) 的余数 \(r\)\(n\) 互质的剩余类。

从这 \(\varphi(n)\) 个剩余类中各取出一个元素,总共 \(\varphi(n)\) 个数,将这些数构成一个新的集合,则称这个集合为模 \(n\) 的简化剩余系。

\(\varphi(n)\) 为欧拉函数,\(1 \sim n\) 中与 \(n\) 互质的数的个数。

例如:\(n = 5\) 时,\(\{ 1, 2, 3, 4\}\) 是一个模 \(5\) 的简化剩余系。\(n = 10\) 时,\(1, 3, 7, 9\) 是一个模 \(10\) 的简化剩余系。

显然,\(n\) 的简化剩余系中所有的数都与 \(n\) 互质


欧拉函数


\(1 \sim n\)\(n\) 互质的数的个数称为欧拉函数,记作 \(\varphi(n)\)

\[\begin{align} n = p_{1}^{a_{1}} p_{2}^{a_{2}} p_{3}^{a_{3}} ... p_{k}^{a_{k}}\\ \varphi(n) = n \times {\textstyle \prod_{i = 1}^{k}} \frac{p_{i} - 1}{p_{i}} \end{align} \]

欧拉定理


定义:若 \(\gcd(a, n) = 1\) ,则 $a^{\varphi(n)} \equiv 1 \pmod{n} $ 。

证明

\(\{ r_{1}, r_{2},...r_{\varphi(n)}\}\) 是一个模 \(n\) 的简化剩余系。

那么 \(r_{i}\) 就是和 \(n\) 互质,又因为 \(\gcd(a, n) = 1\) ,所以 \(a\)\(n\) 也是互质的。

那么 \(\{ ar_{1}, ar_{2},...ar_{\varphi(n)}\}\) 也是一个模 \(n\) 的简化剩余系。所以

\[\begin{align} {\textstyle \prod_{i = 1}^{\varphi(n)}}r_{i} &\equiv {\textstyle \prod_{i = 1}^{\varphi(n)}} ar_{i} \pmod{n} \\ {\textstyle \prod_{i = 1}^{\varphi(n)}}r_{i} &\equiv a ^{\varphi(n)}{\textstyle \prod_{i = 1}^{\varphi(n)}} r_{i} \pmod{n} \\ a ^{\varphi(n)} &\equiv 1 \pmod{n} \\ \end{align} \]

扩展欧拉定理


\[a^{b} = \begin{cases} a ^ {b}&, b < \varphi(m), \mod (m)\\ a ^ {b \mod \varphi(m) + \varphi(m)}&, b\ge \varphi(m), \mod(m) \end{cases} \]

P5091 【模板】扩展欧拉定理 - 洛谷

给定三个正整数 \(a, b, m\) ,求 \(a ^ b \mod m\)

\(1 \le a \le 10^9, 1 \le m \le 10^9, 1 \le b \le 10^{20000000}\)

可以根据扩展欧拉定理,当 \(b < \varphi(m)\) 时,用快速幂求解,当 \(b \ge \varphi(m)\) 时,通过定理进行降幂,然后再用快速幂求解。

实现代码

#include <bits/stdc++.h>
using namespace std;

typedef long long LL;

// 快速幂, a ^ k % p
int qmi(int a, int k, int p) {
    int res = 1;
    while (k) {
        if (k & 1) res = (LL)res * a % p;
        a = (LL)a * a % p;
        k >>= 1;
    }
    return res;
}

// 求 n 的欧拉函数
int get_phi(int n) {
    int res = n;
    for (int i = 2; i <= n / i; i++) {
        if (n % i == 0) {
            while (n % i == 0) n /= i;
            res = res / i * (i - 1);
        }
    }
    if (n > 1) res = res / n * (n - 1);
    return res;
}

// 对 b 进行降幂
int de_pow(string s, int phi) {
    int res = 0;
    bool flag = false;
    for (int i = 0; s[i]; i++) {
        res = res * 10 + s[i] - '0';
        if (res >= phi) flag = true, res %= phi;
    }
    if (flag) res += phi;
    return res;
}

int main() {
    int a, m;
    string b;
    cin >> a >> m >> b;

    int phi = get_phi(m);
    int B = de_pow(b, phi);
    cout << qmi(a, B, m) << "\n";

    return 0;
}

参考资料


522 剩余系 欧拉定理 扩展欧拉定理_董晓算法文章来源地址https://www.toymoban.com/news/detail-646401.html

到了这里,关于欧拉定理 & 扩展欧拉定理的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 数论——中国剩余定理、扩展中国剩余定理 学习笔记

    中国剩余定理(Chinese Remainder Theorem,CRT) 求解如下形式的一元线性同余方程组(其中 (m) 两两互质): $left{begin{matrix}x equiv a_1 pmod {m_1} \\\\x equiv a_2 pmod {m_2} \\\\ dots \\\\x equiv a_k pmod {m_k}end{matrix}right.$ 计算所有模数的积 (M = prod m_i) ; 对于第 (i) 个方程: 计算: (M_i

    2024年02月08日
    浏览(40)
  • OpenCV学习笔记之Overload报错的处理(仅供参考)

    今天在练习一个文档识别的小项目时,运行后一直提示报错,可是我还不知道问题出在哪里,源代码如下: 报错内容如下: 可是错误并不在报错所提示的那一行,而是上面的findContours()函数,新版本中这个函数要有两个返回值(以前好像是有三个),猛然反应过来以后,加上

    2024年02月12日
    浏览(36)
  • 【JS逆向一】逆向某站的 加密参数算法--仅供学习参考

    逆向日期:2024.02.06 使用工具:Node.js 文章全程已做去敏处理!!!  【需要做的可联系我】 可使用AES进行解密处理(直接解密即可):在线AES加解密工具 1、打开某某网站(请使用文章开头的AES在线工具解密):T9mpG48zIdYFB40hkjsQS4+N5rr4x4mSPHlx5EoT/+s= 2、点击右上角的登录按钮,账号

    2024年02月22日
    浏览(55)
  • 深度学习:从入门到精通课后习题解答本答案仅供参考

    第一章: 1、通过本章的学习,你认为深度学习崛起的原因有哪些? 答:(1) 计算能力的发展。深度学习的起源并不晚,但是在发展初期遭遇瓶颈的最主要原因是:当时的计算资源无法支持我们实现深度学习如此庞大复杂的计算。直到我们开始使用GPU进行计算后,深度学习才终

    2024年02月07日
    浏览(53)
  • 「学习笔记」(扩展)中国剩余定理

    有物不知其数,三三数之剩二,五五数之剩三,七七数之剩二。问物几何? 该问题出自《孙子算经》,具体问题的解答口诀由明朝数学家程大位在《算法统宗》中给出: 三人同行七十希,五树梅花廿一支,七子团圆正半月,除百零五便得知。 (2 times 70 + 3 times 21 + 2 times

    2024年02月06日
    浏览(35)
  • Unity单机手游逆向破解思路(仅供学习参考,禁止用于非法行为)

    一、安卓逆向常用工具 针对安卓单机游戏逆向,尤其是逆向使用Unity引擎开发的安卓游戏,只需了解下面的工具即可。 (1)Android Killer        Android Killer是安卓通用逆向工具,其可以对apk进行反向编译,得到smail代码,用户可以更改smail代码后,对apk重新打包,以实现破解

    2024年01月19日
    浏览(37)
  • 欧拉定理公式(包括欧拉降幂)

    a b ≡ { a b   m o d   φ ( p ) , gcd ⁡ ( a , p ) = 1 a b , gcd ⁡ ( a , p ) ≠ 1 , b φ ( p ) (   m o d     p ) a b   m o d   φ ( p ) + φ ( p ) , gcd ⁡ ( a , p ) ≠ 1 , b ≥ φ ( p ) a^{b}equivleft{begin{array}{ll} a^{b bmod varphi(p)}, operatorname{gcd}(a, p)=1 \\\\ a^{b}, operatorname{gcd}(a, p) neq 1, bvarphi(p) quad(bmo

    2024年02月07日
    浏览(47)
  • 数论——欧几里得算法、裴蜀定理、扩展欧几里得算法 学习笔记

    最大公约数 最大公约数即为 Greatest Common Divisor,常缩写为 gcd。 一组整数的公约数,是指同时是这组数中每一个数的约数的数。 (pm 1) 是任意一组整数的公约数; 一组整数的最大公约数,是指所有公约数里面最大的一个。 特殊的,我们定义 (gcd(a, 0) = a) 。 最小公倍数 最

    2024年02月08日
    浏览(45)
  • Office 2019 激活-探索(仅供参考)

    1. 新建txt文件 2. 在新建txt文件中输入 @echooff (cd /d \\\"%~dp0\\\")(NET FILE||(powershell start-process -FilePath \\\'%0\\\' -verb runas)(exit /B)) NUL 21 title Office 2019 Activator r/Piracy echo Converting... mode 40,25 (if exist \\\"%ProgramFiles%Microsoft OfficeOffice16ospp.vbs\\\" cd /d \\\"%ProgramFiles%Microsoft OfficeOffice16\\\")(if exist \\\"%ProgramFiles(x

    2023年04月08日
    浏览(62)
  • restful接口设计规范[仅供参考]

    应该尽量将API部署在专用域名之下。 如果确定API很简单,不会有进一步扩展,可以考虑放在主域名下。 应该将API的版本号放入URL。 另一种做法是,将版本号放在HTTP头信息中,但不如放入URL方便和直观。Github就采用了这种做法。 因为不同的版本,可以理解成同一种资源的不

    2024年02月15日
    浏览(46)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包