代码改变世界

A/B(扩展欧几里得exgcd)

2019-08-14 14:24  木木王韦  阅读(150)  评论(0)    收藏  举报

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.
所以B
x1+9973
y1=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;
}