HDU1576 A/B 扩展欧几里得

A/B

Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 10303 Accepted Submission(s): 8240

Problem 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整除，所以设A/B=x,则A=Bx。
n=A%9973，所以A=9973k+n，即Bx-9973k=n。
又gcd(B,9973)=1，由扩展欧几里得定理，有:Bx1+9973y1=1，等号两边同乘以n，有nBx1+9973y1n=n。
所以Bx-9973k=nBx1+9973y1n，x=nx1。
所以问题变为x%9973,x1用扩展欧几里得原理就可以求出来了。
注意结果要非负。

c++程序如下

#include<cstdio>
using namespace std;
const int  mod = 9973;
void exgcd(int a, int b, int &x, int &y)//扩展欧几里得算法
{
	if (b == 0)
	{
		x = 1;
		y = 0;
		return;
	}
	int x1, y1;
	exgcd(b, a%b, x1, y1);
	x = y1;
	y = x1 - (a / b)*y1;
}
int main()
{
	int t, n, b;
	scanf_s("%d", &t);
	while (t--)
	{
		scanf_s("%d%d", &n, &b);
		int x, y, z;
		exgcd(b, mod, x, y);
		x *= n;
		z = (x%mod + mod) % mod;//使得结果不是负数
		printf("%d\n", z);
	}
	return 0;
}