HDU1576A/B【扩展欧几里德算法模板】

http://acm.hdu.edu.cn/showproblem.php?pid=1576

思路：设(A/B)%9973 = K, 则A/B = k + 9973x  (x未知）， 因此A = kB + 9973xB,

又A%9973 = n， 所以kB%9973 = n,  故kB = n + 9973y (y未知)

故(k/n)B +(-y/n)*9973 = gcd(B,9973) = 1

扩展欧几里得 求出k/n,  再乘以个n，

记得取模，（求出的这个x可能是负数，由题这显然是不正确的。可以做这样的转化(x%9973+9973)%9973）

就是answer了

A/B

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

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>
#include<cstdio>
using namespace std;
#define mod 9973
int exGcd(int a,int b,int &x,int &y)
{
    if(b==0)
    {
        x=1,y=0;
        return a;
    }
    int ans=exGcd(b,a%b,x,y);
    int t=x;
    x=y;
    y=t-a/b*y;
    return ans;
}

int main()
{
    int t,n,b,x,y;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&b);
        exGcd(b,mod,x,y);
        while(x < 0)
            x += mod;
        cout<<(x*n)%mod<<endl;
    }
    return 0;
}