hdoj1576A/B问题求逆元的多种方法

hdoj1576
题意：
输入n，B求(A/B)%9973，其中n=A%9973，B与9973互质。
思路：
（A/B）%9973可以转变成（A* B^-1）%9973，其中B^-1是B的逆元。所以这题只要求出B的逆元即可。
求逆元的方法：
1.因为9973是质数，所以可以用费马小定理 (a^-1=a^p-2(mod p) ) 加上快速幂求出B的逆元。
2.因为B与9973互质，利用欧拉定理（a^φ(n) ）=1 (mod n)，其中φ(n)是1~n上与n互质的数的总数，可以利用欧拉函数
（φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)……(1-1/pn) 其中pi是n的所以质因子）求出。同样的加上快速幂就可以把B的逆元求出来了。
其实费马小定理是欧拉定理的一种特殊情况。
3.利用扩展欧几里得求出B的逆元。 扩展欧几里德算法是用来在已知a, b求解一组x，y，使它们满足贝祖等式： ax+by = gcd(a, b) =d。扩展欧几里德常用在求解模线性方程及方程组中。如果d不等于1，说明没有ax+by=1的解，但是本题中没有这种情况。
下面代码
利用方法1或2：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

int euler_phi(int n)
{
    int sum=n;
    for(int i=2;i<=n;i++)
    {
        if(n%i==0)
        {
            sum-=sum/i;
            while(n%i==0)
                n/=i;
        }
    }
    if(n>1)
        sum-=sum/n;
    return sum;
}

ll qpow(ll a,ll m)//利用欧拉定理加上快速幂计算
{
    ll b=1;
    ll p=euler_phi(m)-1;//如果利用费马小定理直接p=m-2计算
    while(p)//计算 a^φ(n)-1^
    {
        if(p%2)
            b=(a*b)%m;
        a=(a*a)%m;
        p/=2;
    }
    return b%m;
}

int main()
{
    ll n,b,t;
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld %lld",&n,&b);
        printf("%d\n",n*qpow(b,9973)%9973);
    }
    return 0;
}

利用方法3核心代码：

int ex_gcd(int a,int b,int &x,int &y)
{
    int d=a;
    if(b==0)
    {
        x=1;
        y=0;
    }
    else
    {
        d=ex_gcd(b,a%b,x,y);
        int temp=y;
        y=x-(a/b)*y;
        x=temp;
    }
    return d;
}

int mod_inverse(int a,int m)//计算逆元
{
    int x,y;
    int k=ex_gcd(a,m,x,y);//如果k不是1证明a和m不是互质，说明a的逆元不存在，但是本题中不存在这种情况
    if(k!=1)
        return -1;
    else
        return (m+x)%m;//这就是a的逆元
}