HDU_1576 A/B

A/B（扩展欧几里得算法）

 

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
  
  

 

分析：

扩展欧几里得的基本式是： ax1+by1=gcd(a,b)

分析此题题意：

n=A%9973=A-A/9973*9973。又因为题中已知A%B=0，所以令A/B=X,那么A=BX。

所以Bx-A/9973*9973=n。令A/9973=y，则Bx-9973y=n。

只要求出x的值，即可算出x%9973，也就是(A/B)%9973了。

题中又已知。。。。gcd(B,9973)=1,所以n=1.  那么gcd(B,9973)=Bx1+9973y1

等式两边同乘以n，得B(nx1)-9973(-ny1)=n。可知nx1就是Bx-9973y=n的解了。即x=nx1。

但是求出的x可能是负数，而题中X必须为正数。

可以做这样的转化：(x%9973+9973)%9973

代码：

#include<stdio.h>
#define mod 9973
void extend_gcd(int a,int b,int &x,int &y)  //扩展欧几里得算法的子函数 
{
   if(b==0)
   {
      x=1;
      y=0;
   }
   else 
   {
      extend_gcd(b,a%b,y,x);
      y-=x*(a/b);
   }
}
int main()
{
   int T,n,b,x,temp,y;
   scanf("%d",&T);
   while(T--)
   {
      scanf("%d%d",&n,&b);
      extend_gcd(b,mod,x,y);
      x*=n;                   //由题意分析出来的x的求法 
      temp=(x%mod+mod)%mod;      //题中x必为正，但是上面求出的可能为负，那么就把它转变为必正 
      printf("%d\n",temp);
   }
   return 0;
}