寒假训练数论C

题目
要求(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。
解题思路
9973是一个质数。
由费马小定理的B^(9972)≡1(mod 9973)
所以BB^(9971)≡1(mod 9973)
B^(9971)为B的逆元
所以(A/B)%9973=（A B^(9971))%9973
代码

#include <cstdio>
#include <iostream>
using namespace std;
const int MOD=9973;
long long MI(long long a,long long b)
{
	long long base=a;
	long long ans=1;
	while (b)
	{
		if (b&1) ans=(ans*base)%MOD;
		base=(base*base)%MOD;
		b>>=1;
	}
	return ans;
}
int main()
{
	int T,n,b;
	scanf("%d",&T);
	while (T--)
	{
		cin>>n>>b;
		long long k=MI(b,MOD-2);
		long long ans=(n*k)%MOD;
		cout<<ans<<endl;
	}
}