exgcd之求逆元

  下面我们总结一下exgcd的另一个用途---求逆元
  首先什么是逆元呢？
  a*x=1(mod m)这里，我们称 x 是 a 关于 m 的乘法逆元,其实也就是解a*x+m*y==1这个方程的解
  再来举个例子详细解释一下逆元的用法
  打个比方吧，就是给你个式子：(a*b/c)%MOD(1<a,b,c<=10^18),那么我们怎么求呢，如果是a*b*c%MOD是好求的，我们可以先化为：
(a%MOD)*(b%MOD)*(c%MOD)%MOD;但是(a*b/c)%MOD(1<a,b,c<=10^18)这个式子，如果我们先进行取模的话，(a%MOD)*(b*MOD)/(c%MOD)的话，我们不确定 a 和 b 取完模之后还能整除 c ，所以就引入了逆元。
(a/b)%MOD就可以化为：( (a%MOD)*_b )%MOD
其中_b就是b%MOD的逆元，那么逆元是怎么求的呢，
如果 b*_b % MOD == 1的话，那么就称 _b 是 b 的逆元，
说到这，其实求逆元就变得非常容易了，套用exgcd就可，那求最小乘法逆元就和最小正整数解
如出一辙了

例题：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=4712求最小逆元

#include <iostream>  
#include <cstdio>  
#include <cstring>  
#include <cmath>  


using namespace std;  
  
typedef long long LL;  
typedef double DB;  
  
LL e_gcd(LL a,LL b,LL &x,LL &y)  
{  
    if(b==0)  
    {  
        x=1;  
        y=0;  
        return a;  
    }  
    LL ans=e_gcd(b,a%b,x,y);  
    LL temp=x;  
    x=y;  
    y=temp-a/b*y;  
    return ans;  
}  
  
LL cal(LL a,LL b,LL c)  
{  
    LL x,y;  
    LL gcd=e_gcd(a,b,x,y);  
    if(c%gcd!=0) return -1;  
    x*=c/gcd;  
    b/=gcd;  
    if(b<0) b=-b;  
    LL ans=x%b;  
    if(ans<=0) ans+=b;  
    return ans;  
}  
  
int main()  
{  
    LL a,b,t;  
    scanf("%lld",&t);  
    while(t--)  
    {  
        scanf("%lld%lld",&a,&b);  
        LL ans=cal(a,b,1);  
        if(ans==-1) printf("Not Exist\n");  
        else printf("%lld\n",ans);  
    }  
    return 0;  
}