码农知识堂 - 1000bd
  •   Python
  •   PHP
  •   JS/TS
  •   JAVA
  •   C/C++
  •   C#
  •   GO
  •   Kotlin
  •   Swift
  • 扩展欧几里得


    扩展欧几里得,常用于计算ax+by=gcd(a,b)的一组可行解

    设ax_{1}+by_{1}=gcd(a,b),

    bx_{2}+(a%b)y_{2}=gcd(b,a%b)

    由欧几里得可知gcd(a,b)=gcd(b,a%b)

    因此ax_{1}+by_{1}=bx_{2}+(a%b)y_{2}

    又因为a mod b = a-(\left \lfloor \frac{a}{b} \right \rfloor*b)

    所以ax_{1}+by_{1}=bx_{2}+(a-\left \lfloor \frac{a}{b} \right \rfloor*b)y_{2}

    可得ax_{1}+by_{2}=ay_{2}+b(x_{2}-\left \lfloor \frac{a}{b} \right \rfloor y_{2})

    因为a=a,b=b,所以x_{1}=y_{2},y_{1}=x_{2}-\left \lfloor \frac{a}{b} \right \rfloor y_{2}

    将x_{2},y_{2}不断代入递归求解直至gcd为0递归x=1,y=0回去求解

    [NOIP2012 提高组] 同余方程 - 洛谷

    求解方程最小正整数解,先把题目的方程分解来看ax=1+by,然后因为x和y的符号都是由我们决定的,所以可以得到ax+by=1,题面指出输入数据一定有解,因此从扩展欧几里得可知,ax+by=m中如果有解,那么一定存在m%gcd(a,b)=0,因此可得出题目中a和b一定互素,该题不能直接费马小定理求\frac{1}{a}在mod b下的逆元,因为并没有保证b是质数

    AC代码:

    1. #include
    2. #define rep(i,a,n) for(int i=a;i
    3. using namespace std;
    4. using LL = long long;
    5. int main() {
    6. ios::sync_with_stdio(false);
    7. cin.tie(nullptr);
    8. int a, b, x, y;
    9. cin >> a >> b;
    10. function<int(int, int, int&, int&)> exgcd = [&](int a, int b, int &x, int &y) {
    11. if (!b) {
    12. x = 1;
    13. y = 0;
    14. return a;
    15. }
    16. int d = exgcd(b, a % b, x, y);
    17. int t = x;
    18. x = y;
    19. y = t - (a / b) * y;
    20. return d;
    21. };
    22. int g = exgcd(a, b, x, y);
    23. cout << (x + b) % b << '\n';
    24. return 0;
    25. }

  • 相关阅读:
    DateUtils 和 CalendarUtil 类来执行日期处理操作使用说明书
    Linux时间同步练习
    Python内置函数hex()详解
    使用 Electron 来替代本地调试线上代理的场景
    ArrayList为什么线程不安全以及三种解决办法【详细】
    Java新手小白,Java基础知识
    0-1背包
    云同行 月更明
    Jenkins部署前后端不分离springboot项目
    ProEasy机器人:TCP无协议通讯(socket通讯)时打印log日志
  • 原文地址:https://blog.csdn.net/eyuhaobanga/article/details/126304464
  • 最新文章
  • 攻防演习之三天拿下官网站群
    数据安全治理学习——前期安全规划和安全管理体系建设
    企业安全 | 企业内一次钓鱼演练准备过程
    内网渗透测试 | Kerberos协议及其部分攻击手法
    0day的产生 | 不懂代码的"代码审计"
    安装scrcpy-client模块av模块异常,环境问题解决方案
    leetcode hot100【LeetCode 279. 完全平方数】java实现
    OpenWrt下安装Mosquitto
    AnatoMask论文汇总
    【AI日记】24.11.01 LangChain、openai api和github copilot
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1
正则表达式工具 cron表达式工具 密码生成工具

京公网安备 11010502049817号