A/B 擴充套件歐幾里德
阿新 • • 發佈:2019-01-02
要求(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);
}
}