1. 程式人生 > >A/B 擴充套件歐幾里德

A/B 擴充套件歐幾里德

要求(A/B)%9973,但由於A很大,我們只給出n(n=A%9973)(我們給定的A必能被B整除,且gcd(B,9973) = 1)。

Input

資料的第一行是一個T,表示有T組資料。 
每組資料有兩個數n(0 <= n < 9973)和B(1 <= B <= 10^9)。

Output

對應每組資料輸出(A/B)%9973。

Sample Input

2
1000 53
87 123456789

Sample Output

7922
6060

#include<cstdio>
typedef long long ll;
const ll mod=9973;
ll exgcd(ll a,ll b,ll &x,ll &y)  //擴充套件歐幾里德 
{
    if(!b)
    {
        x=1;
        y=0;
        return a;
    }
    ll r=exgcd(b,a%b,x,y);
    ll t=x;
    x=y;
    y=t-(a/b)*y;
    return r;

ll inv(ll a,ll mod)
{
    ll x,y;
    exgcd(a,mod,x,y);
return (mod+x%mod)%mod;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        ll n,b;
        scanf("%d%d",&n,&b);
        printf("%d\n",(n*inv(b,mod)%mod)%mod);
    }
}