hdu 1576 A/B （求乘法逆元——扩展欧几里得）

A/B

Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)

Total Submission(s): 3745 Accepted Submission(s): 2869



Problem Description

要求(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<iostream>
#include<vector>
#include<map>
#include<set>
#include<string>
#include<string.h>
#include<algorithm>
#include<cstdio>
using namespace std;
#define sz(v) (int)(v.size())
#define rep(i,a,n) for(int i=a;i<n;++i)
#define per(i,a,n) for(int i=n-1;i>=a;--i)
#define mem(a,t) memset(a,t,sizeof(a))
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define N 1000005
#define LL __int64
const int mod=1e9+7;

LL exgcd(LL a,LL b,LL &x,LL &y)
{
LL d=a;
if(a==0&&b==0)
return -1;
if(b)
return exgcd(b,a%b,y,x),y-=a/b*x;
else{
x=1;
y=0;
}
return d;
}
int mod_reverse(LL a,LL mod)
{
LL x,y,d;
d=exgcd(a,mod,x,y);
if(d==-1)
return d;
return (x%mod+mod)%mod;
}

int main()
{
//freopen("in.txt","r",stdin);
int T;
int n,b,x;
scanf("%d",&T);
while(T--){
scanf("%d%d",&n,&b);
x=mod_reverse(b,9973);
printf("%d\n",n*x%9973);
}
return 0;
}