HDU 1576 A/B

题目链接：A/B

解题思路：数论真是费脑筋，网上看了好多感觉都是互相抄，都不到为什么扩展GCD之后的系数就是解。

 ≡

我这里给出我的证明

A ≡N(MOD 9973) ① 

B*B^-1 ≡1(MOD 9973) ②

这里我们称B^-1为B的逆元而不是真的1 / B。

1.首先逆元怎么求，我们知道若gcd(a, b)=1则一定存在p, q使得ap+bq=1

B*B^-1 = 9973*K + 1  => B*B^-1 - 9973*K = 1,则如果 gcd(B, K)=1则B在K的逆元一定存在。而题设条件gcd(9973, B)=1使得B的逆元存在

2.原式变形

A = KB

KB ≡ NBB^-1  (MOD 9973)

根据百度百科上同余条目有这么一条

若ac ≡ bc (mod m) c≠0 则 a≡ b (mod m/gcd(c,m)) 其中gcd(c,m)表示c,m的最大公约数，特殊地 ,gcd(c,m)=1 则a ≡ b (mod m)。
那原式中的B就可以依此消除。
得到K≡NB^-1(MOD 9973)，所以答案就是N * B^-1 % 9973

#include<cstdio>
#include<cmath>
#define ll __int64
#define MOD 9973

ll ext_gcd(ll a, ll b, ll &x, ll &y){
    if(a == 0 && b == 0){
        return -1;
    }
    if(b == 0){
        x = 1; y = 0;
        return a;
    }
    int r = ext_gcd(b, a % b, y, x);
    y -= x * (a / b);
    return r;
}

ll mod_inverser(ll a, ll m){
    ll x, y;
    ll d = ext_gcd(a, m, x, y);
    if(d == 1)
        return (m + x % m) % m;
    else{
        return -1;
    }
}


int main(){
    int i, j, k, t;
    ll n, b;
    scanf("%d", &t);
    for(i = 0; i < t; i++){
        scanf("%I64d%I64d", &n, &b);
        b = mod_inverser(b, MOD);
        printf("%I64d\n", (b * n + MOD) % MOD);
    }
    return 0;
}