hdu 1576 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

 

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#define LL __int64
const int maxn = 1e3 + 10;
const int mo = 9973;
using namespace std;

void exgcd(LL a, LL b, LL &d, LL &x, LL &y)
{
    if(!b) {d = a; x = 1; y = 0;}
    else{ exgcd(b, a%b, d, y, x); y -= x*(a/b); }
}

LL mo_reve(LL a, LL n)
{
    LL x, y;
    LL d;
    exgcd(a, n, d, x, y);
    if(d==1) return (x%n+n)%n;
    else return -1;
}

int main()
{
     int t, n, b;
     scanf("%d", &t);
     while(t--)
     {
         scanf("%d%d", &n, &b);
         int x = mo_reve(b, mo);
         printf("%d\n", n*x%mo);
     }
     return 0;
}