HDU1576 A/B （扩展欧几里得）

 

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

 

 

 

思路：先设（A/B）%9973 = k;       可得A/B = k + 9973 * x;=>A = kB + 9973*x*B;可得到kB % 9973 = n;   则kB = n + 9973 * y;则 kB + 9973 *y = n，y属于整数，可正可负。

（k / n） * B + 9973*（y / n） = 1；又gcd（B,9973）= 1，所以求出k = x *n %9973; 

 

 

#include <stdio.h>

void exgcd(long long a, long long b, long long &x, long long &y){
	if(b == 0){
		x = 1;		y = 0;
		return ;
	}
	
	exgcd(b,a%b,x,y);
	
	long long t = x;
	x = y;
	y = t - a/b * x;
	return ;
}

int main(){
	long long t,n,b;
	
	scanf("%lld",&t);
	
	while(t --){
		scanf("%lld%lld",&n,&b);
		
		long long x,y;
		
		exgcd(b,9973,x,y);
		
		x = (x % 9973 + 9973) % 9973;
		
		printf("%lld\n",(x * n) % 9973);
	}
	return 0;
}