一道有趣的数学题分享

最新推荐文章于 2025-06-06 21:48:36 发布

herunze666

最新推荐文章于 2025-06-06 21:48:36 发布

阅读量314

点赞数 1

文章标签：算法

本文链接：https://blog.csdn.net/herunze666/article/details/148049355

版权

离散对数(洛谷P12540)

题目描述：

给定正整数 a, c, p，保证 p 是素数，求 b 使得： $a^{b}\equiv b^{c}(mod p)$

我们称整数 A,B,C 有 A≡B(modC)，当且仅当存在整数 k 使得 $a-b=c*k$ 。

解题思路：

这道题若a,p互质，则根据费马小定理得到 $a^{p - 1}\equiv 1(modp)$ ，根据题中的式子就可以得出 $a^b \equiv a^{b mod (p-1)} modp.$ 所以就有 $a^{b} mod p = a^{(k*p +b mod p) mod (p-1)} mod p$ 将这个式子再度展开就可以得到 $a^{(b mod p) mod (p-1) + k*p mod (p-1)} mod p$ 。但 $kp mod (p-1)$ 等于 $k mod(p-1)*p mod(p-1).$ 因为 $p \equiv 1 mod(p-1)$ ，所以 $p = 1 mod(p-1)$ . 所以 $k p mod (p - 1) =k * 1 mod(p-1) =k mod(p-1)$ . 所以 $b= kp+s \rightarrow b mod(p-1) = (k +s) mod(p-1).$

但是我们又希望 $a^{b} mod p = a^{ (k +bmodp) mod (p-1) } mod p.$ 同时 $b^c modp = (kp +bmodp)^c modp$ 。由于 $kp \equiv 0 modp$ ，所以条件变成了 $a^{ (k +bmodp) mod (p-1) } \equiv (bmodp)^c modp.$ 要满足这个等式我们就需要选择k使得满足上述条件，我们可以令 $bmodp=amodp$ ，此时，右边的 $(bmodp)^c modp$ 等于 $a^c mod p$ ,而对于左边来说我们希望能让这一部分等于 $a^c mod p$ 所以，可以令 $(k +bmodp) mod (p-1) \equiv c mod(p-1)$ . 因为如果这样的话，左边就是 $a^c mod p$ ，而右边同样也是 $a^c mod p$ 此时等式成立。

所以，要满足 $k + bmodp \equiv c mod(p-1).$ 因此， $k \equiv (c -bmodp) mod(p-1)$ 所以，k可以等于 $(c -bmodp) mod(p-1) + t*(p-1)$ ，其中t是任意整数。那么此时我们就可以推出结论 $b = [(c -bmodp) mod(p-1)] *p +bmodp + t*(p-1)*p$

解决代码

#include<bits/stdc++.h>
#define endl "\n"

using namespace std;
typedef long long LL;

int main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);

    LL a,b,c,p;
    cin >> a >> c >> p ;

    LL n = a % p;
    LL m = c % (p - 1);
    LL k = (m - n) % (p - 1);
    if(k < 0)
        k += (p - 1);
    if(p == 1)
        cout << 1;
    else
    {
        cout <<k * p  + n;
    }

    return 0;

}