51nod 1256 乘法逆元

51nod 1256 乘法逆元

http://www.51nod.com/Challenge/Problem.html#!#problemId=1256

theory

https://baike.baidu.com/item/%E6%89%A9%E5%B1%95%E6%AC%A7%E5%87%A0%E9%87%8C%E5%BE%B7%E7%AE%97%E6%B3%95/1053275?fr=aladdin

对于不完全为 0 的非负整数 a，b，gcd（a，b）表示 a，b 的最大公约数，必然存在整数对 x，y ，使得 gcd（a，b）=ax+by。
求解 x，y的方法，设 a>b。
1，显然当 b=0，gcd（a，b）=a。此时 x=1，y=0；
2，a>b>0 时
设 ax1+ by1= gcd(a,b);
bx2+ (a mod b)y2= gcd(b,a mod b);
根据朴素的欧几里德原理有 gcd(a,b) = gcd(b,a mod b);
则:ax1+ by1= bx2+ (a mod b)y2;
即:ax1+ by1= bx2+ (a - [a / b] * b)y2=ay2+ bx2- [a / b] * by2;
也就是ax1+ by1 == ay2+ b(x2- [a / b] *y2);
根据恒等定理得：x1=y2; y1=x2- [a / b] *y2;
这样我们就得到了求解 x1,y1 的方法：x1，y1 的值基于 x2，y2.
上面的思想是以递归定义的，因为 gcd 不断的递归求解一定会有个时候 b=0，所以递归可以结束。

扩展欧几里德算法是用来在已知a, b求解一组x，y使得ax+by = Gcd(a, b) =d(解一定存在，根据数论中的相关定理)。扩展欧几里德常用在求解模线性方程及方程组中。

思路

K * M % N = 1等价于 K*M=N*X+1  即 K*M+N*(-X)=1 
根据扩展欧几里德算法，求出K和（-X）；
而K是应为正整数，即需要若K为负整数，则需要将之化正，即与负数取模同理，将K加上N，直至K>0为止，所得的数即为最小的乘法逆元；
若K为正整数，则直接输出即可

例子：
输入
2 3
输出
2

#include <iostream>
using namespace std;

//扩展欧几里得
int exgcd(int a,int b,int &x,int &y){
    if(b==0){
        x = 1;
        y = 0;
        return a;
    }
    int ans = exgcd(b,a%b,x,y);
    int t;
    t = x;
    x = y;
    y = t - a/b*y;
    return ans;
}

int main() {
    int m,n;
    int x,y;
    cin>>m>>n;
    exgcd(m,n,x,y);
    while(x<0) x+=n;
    cout<<x<<endl;
    return 0;
}