HDU-1576 扩展欧几里得 or 逆元

最新推荐文章于 2021-02-24 00:30:01 发布

jesHrz

最新推荐文章于 2021-02-24 00:30:01 发布

阅读量163

点赞数

分类专栏： ACM 数学

本文链接：https://blog.csdn.net/weixin_43135318/article/details/94588863

版权

ACM 同时被 2 个专栏收录

37 篇文章 0 订阅

订阅专栏

数学

8 篇文章 0 订阅

订阅专栏

题目链接

HDU-1576

题意

一句话，给定n和B，求 $\frac{A}{B}\ mod\ 9973$ ，其中n=A%9973，数据保证B能整除A

思路

发现9973是个质数，可以直接用费马定理求下B的逆元再乘A就完事了。
然后发现可以用扩欧做。记 $\frac{A}{B}\ mod \ 9973$
则有
$\frac{A}{B}=k_1*9973+ans$ $A=k_29973+n$
然后整理一下可以得到
$n + k_29973=k_1B*9973+ans*B$
再整理一下可以得到
$k_1B-k_2)9973+ans*B=n$
然后我们就得到了一个二元一次不定方程,其中B和9973是系数，ans和9973前面的系数是未知量，套上扩欧的板子就好了。

代码

#include <algorithm>
#include <cstdio>
#include <cstring>
#define inf 1ll << 62
#define cases(yjq) for (int cas = 1; cas <= (int)yjq; ++cas)
using namespace std;
const long long mod = 9973;
inline long long Abs(long long x) { return x < 0 ? -x : x; }

long long exgcd(long long a, long long b, long long& x, long long& y) {
    if (!b) {
        x = 1;
        y = 0;
        return a;
    }
    long long t = exgcd(b, a % b, y, x);
    y -= a / b * x;
    return t;
}

bool solve(long long a, long long b, long long c, long long& x, long long& y) {
    long long g = exgcd(a, b, x, y);
    if (c % g) return false;
    x *= c / g;
    y *= c / g;
    long long tmp = Abs(b / g);
    x = (x % tmp + tmp) % tmp;
    y = (c - a * x) / b;
    return true;
}
long long n, b;
int main() {
    int t;
    scanf("%d", &t);
    cases(t) {
        scanf("%d%d", &n, &b);
        long long x, y;
        solve(b, mod, n, x, y);
        printf("%lld\n", x % mod);
    }
    return 0;
}