A/B (拓展欧几里得)

A/B

点我找原题

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

很难看出这是一道拓展欧几里得题目，突破口在n=A%9973，那么久假设A-n=y*9973 （1），又因为A能被B整除，

那么再假设A=B*x（2）,那么最后就只要求x%9973就行了，将（2）带入（1）得x*B-y*9973=n,这就能看出拓展欧几里得

的样子了，但是提目交代GCD(B,9973)=1，所以x*B-y*9973=1才是真正满足拓展欧几里得的，所以只要先求出x，

在乘上n就得到了我们所要的x，还要注意一点就是最先求出的x并不太一定是正数，所以用x+=9973直到x为正数即可。

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
    
    
    
    

     
     
#include <iostream>
using namespace std;
int x,y,q;
void extend_eulid(int a,int b)//拓展欧几里得
{
	if(b==0)
	{
		x=1;y=0;q=a;
	}
	else 
	{
		extend_eulid(b,a%b);
		int temp=x;x=y;y=temp-a/b*y;
	}
}
int main(int argc, char *argv[])
{
	int a,b,n,t;
	cin>>t;
	while(t--)
	{
		cin>>n>>b;
		extend_eulid(b,9973);
		while(x<0)
		{
			x+=9973;
		}
		x=x*n;
		cout<<x%9973<<endl;
	}
	return 0;
}