求乘法逆元模板(扩展欧几里得)

void exgcb(LL a,LL b,LL &d,LL &x,LL &y){
    if(!b){d=a;x=1;y=0;return;}
    exgcb(b,a%b,d,y,x);y-=x*(a/b);
}
LL ny(LL a,LL b){  ///求a关于b的逆元(要求a，b互质)
    LL d,x,y;
    exgcb(a,b,d,x,y);
    return d==1?(x+b)%b:-1;
}

时间： 2024-04-12 20:38:33

求乘法逆元模板(扩展欧几里得)的相关文章

51nod 125乘法逆元（扩展欧几里得）

给出2个数M和N(M < N),且M与N互质,找出一个数K满足0 < K < N且K * M % N = 1,如果有多个满足条件的,输出最小的. Input 输入2个数M, N中间用空格分隔(1 <= M < N <= 10^9) Output 输出一个数K,满足0 < K < N且K * M % N = 1,如果有多个满足条件的,输出最小的. Input示例 2 3 Output示例 2 思路: 对于正整数和,如果有,那么把这个同余方程中的最小正整数解叫做

ZOJ - 3609 —— Modular Inverse 【乘法逆，扩展欧几里得】

题目:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=4712 1. 这题数据范围太小,直接暴力即可 2. 不过其实这题也是很“直白的”求乘法逆的题目,即当b=1的特殊的模线性方程问题ax≡b mod(n),可以通过扩展欧几里得算法求解: ax≡b mod(n) => (ax) mod n = b mod n => ax=k1*n+r ... (1) b=k2n+r ... (2) (1)-(2)=>ax-

【codevs 1200】【NOIP 2012】同余方程拓展欧几里德求乘法逆元模板题

模板,,, #include<cstdio> using namespace std; void exgcd(long long a,long long b,long long &x,long long &y){ if (b==0) {x=1; y=0;} else {exgcd(b,a%b,x,y); int t=y; y=x-a/b*y; x=t;} } int main(){ long long a,b,x,y; scanf("%lld %lld\n"

hdu_1576A/B(扩展欧几里得求逆元)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1576 A/B Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 4020 Accepted Submission(s): 3091 Problem Description 要求(A/B)%9973,但由于A很大,我们只给出n(n=A%99

辗转相除法与扩展欧几里得

辗转相除法:求gcd(a,b) 扩展欧几里得:解关于x和y的方程:a*x+b*y=gcd(a,b) 1 int gcd(int a,int b){ 2 if (b==0) return a; 3 return gcd(b,a%b); 4 } 5 6 7 ------------------------------------------------ 8 9 10 int extgcd(int a,int b,int& x,int& y){ 11 int d=a; 12 if (b!=0){

暑期学习日记二—利用扩展欧几里得求逆元

最近学习了扩展欧几里得和乘法逆元的关系,在这里写一下巩固一下记忆扩展欧几里得是什么呢,在这就不详解了,可以自行百度,主要来说,对于求解ax ≡ 1(mod n)来说,当gcd(a,n)=1时,证明逆元存在,若不等于1,则证明逆元不存在. 那么当逆元存在时,我们要如何求它的逆元呢? 首先是扩展欧几里得定理,先将式子转换成 ax-ny = 1 的形式,然后我们要通过扩展欧几里得定律去获得它的最大公约数,还有它的一组解 X0,Y0 1 int exgcd(int a, int b, int &x,

hiho1530(乘法逆元)(扩展欧几里得)

#1530 : 分数取模时间限制:1000ms 单点时限:10000ms 内存限制:256MB 描述给定三个正整数 a. b 和 p,满足 b 和 p 互质.这时分数 a / b 除以 p 的余数,即 a / b MOD p 可以定义为 a × b-1 MOD p. 其中b-1 是 b 的逆元,它满足 1 ≤ b-1 < p 且 b × b-1 ≡ 1 MOD p,满足上述条件的 b-1有且仅有一个. 例如 2-1 ≡ 4 MOD 7,因为2 × 4 ≡ 1 MOD 7: 3-1 ≡ 3 M

扩展欧几里得模板（洛谷1082 同余方程NOIP 2012 提高组第二天第一题）

题目描述求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解. 输入输出格式输入格式: 输入只有一行,包含两个正整数 a, b,用一个空格隔开. 输出格式: 输出只有一行,包含一个正整数 x0,即最小正整数解.输入数据保证一定有解. 输入输出样例输入样例#1: 3 10 输出样例#1: 7 说明 [数据范围] 对于 40%的数据,2 ≤b≤ 1,000: 对于 60%的数据,2 ≤b≤ 50,000,000: 对于 100%的数据,2 ≤a, b≤ 2,000,000,000

欧几里得，扩展欧几里得（模板）

1 int gcd(int a,int b) 2 { 3 return b?gcd(b,a%b):a;//最后返回的a为最大公约数 4 } 扩展欧几里得求逆元:51nod1256 1 #include <iostream> 2 #include <stdio.h> 3 #include <string.h> 4 using namespace std; 5 int d,x,y; 6 //ax+by=1(x为a的逆元,同理y为b的逆元) 7 void Exgcd(int