求逆元（拓展欧几里得）

要求(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

这里的逆元是啥？

例如给两个数，7,20。

因为7*3%20=1，则3就是7的逆元。一个数a与它逆元的乘积对b取余等于1。

求逆元的公式就是套一个模板，假如说我们求a关于b的逆元就是解这个方程：
a*x + b*y == 1，这个方程就可以直接求了。求出的x就是a关于b的逆元。

 

#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string.h>
using namespace std;
#define mod 9973
typedef long long ll;
//a存在逆元的条件：gcd(a,mod)=1;
//拓展欧几里得求出的x是a关于b的逆元，y是b关于a的逆元 ，并且gcd(a,b)=1; 
void exgcd(ll a,ll b,ll &x,ll &y){
	if(!b){
		x=1;
		y=0;
		return ;
	}
	exgcd(b,a%b,x,y);
	ll t=x;
	x=y;
	y=t-(a/b)*y;
	return ;
}

int main(){
	ll t;
	cin>>t;
	while(t--){
		ll n,b;
		cin>>n>>b;
		ll x,y;
		exgcd(b,mod,x,y);
		x=(x%mod+mod)%mod;//求出的x可能是最小正整数或最大负整数； 
		cout<<n*x%mod<<endl;
	}
	return 0;
}