扩展欧几里得求乘法逆元

参考http://blog.csdn.net/zhjchengfeng5/article/details/7786595

#include <stdio.h>
#include <math.h>
__int64 ans_gcd(__int64 a,__int64 b,__int64 *x,__int64 *y)
{
    if (b==0)
    {
        *x=1;
        *y=0;
        return a;
    }
    __int64 ans=ans_gcd(b,a%b,x,y);
    __int64 temp=*x;
    *x=*y;
    *y=temp-a/b*(*y);
    return ans;
}
__int64 calculate(__int64 a,__int64 b)
{
    __int64 x,y,ans;
    __int64 gcd=ans_gcd(a,b,&x,&y);
    if (1%gcd!=0)//如果a与b不互质则没有乘法逆元
        return -1;
    x=x*(1/gcd);
    b=abs(b);
    ans=x%b;
    if (ans<=0)
        ans=ans+b;
    return ans;
}
int main()
{

    __int64 i,j,a,b,ans;
    int n;
    scanf("%d",&n);
    while(n--)
    {
        scanf("%I64d%I64d",&a,&b);
        ans=calculate(a,b);//求a关于b的乘法逆元
        if (ans==-1)
            printf("Not Exist\n");
        else
            printf("%I64d\n",ans);

    }
    return 0;
}