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
已知gcd(B,9973)=1,所以B和9973互质,设A=Bx,n=A-(A/9973)9973=A%9973.
所以Bx1+9973y1=n,又Bx+9973y=1;用exgcd求出x,即x1=x*n=A/B。
ac代码:
#include<iostream>
#include<algorithm>
using namespace std;
int exgcd(int a,int b,int &x,int &y){
if(b==0){
x=1;
y=0;
return a;
}
int r=exgcd(b,a%b,x,y);
int t=x;
x=y;
y=t-(a/b)*y;
return r;
}
int main(){
int t;
cin>>t;
while(t--){
int n,b,a,x,y;
cin>>n>>b;
exgcd(b,9973,x,y);
cout<<(n*(x%9973+9973))%9973<<endl;
}
// exgcd(b,n,x,y);
// cout<<x<<" "<<y<<endl;
return 0;
}