ACM_数学

最新推荐文章于 2024-08-15 22:41:27 发布

duiduia001

最新推荐文章于 2024-08-15 22:41:27 发布

阅读量83

点赞数

文章标签：数学算法

本文链接：https://blog.csdn.net/qq_48099121/article/details/120679167

版权

同余

$ax\equiv b(mod m)$ 其中 $d = g c d (a, m) ∣ b$

$x$ 有 $d$ 个模 $m$ 不同余的解

欧拉降幂

在这里插入图片描述
广义欧拉降幂：只要判断b与 $\phi(p)$ 的大小关系：
$a^{b}=a^{b\%\phi(p)}$ $b<\phi(p)$
$a^{b}=a^{b\%\phi(p)+\phi(p)}$ $b\geq \phi(p)$

模板题

在这里插入图片描述

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

LL a, b, c, mod = 1e9+7;

LL qmi(LL a, LL b, LL p) {
    LL res = 1;
    while (b) {
        if (b & 1)  res = 1ll * res * a % p;
        b >>= 1;
        a = 1ll * a * a % p;
    }
    return res;
}

LL phi(LL x) {  //返回x的欧拉函数值
    LL res = x;
    for (LL i = 2; i <= x / i; i ++) {
        if (x % i == 0) {
            res = res / i * (i - 1);
            while (x % i == 0)  x /= i;
        }
    }
    if (x > 1)  res = res / x * (x - 1);
    return res;
}

LL sol(LL a, LL b, LL p) {  //欧拉降幂
    LL t = phi(p);
    if (__gcd(a, p) > 1 && b >= t)    return qmi(a, b % t + t, p);
    else return qmi(a, b % t, p) % p;

}

int main() {
     while(~scanf("%lld%lld%lld", &a, &b, &c)) {
        LL t = sol(b, c, phi(mod));
        printf("%lld\n",sol(a, t, mod));
    }
    return 0;
}

Gym - 101550E

在这里插入图片描述

$求n^{(n-1)^{(n-2)^{...^{2^{1}}}}}\%m$ 的值

$AC代码在幂数b<模数m时，直接用快速幂可以AC，但是用欧拉降幂继续递归求解会WA。有解释说b<\phi(m)时不满足欧拉降幂的条件。$ ~~这种情况不就是 $b<\phi(p)$ 的情况吗，但用这种写法确实会WA~~ 。

#include <bits/stdc++.h>


#define MOD(a,b) a>=b?a%b+b:a

using namespace std;

typedef long long LL;
typedef long long ll;

LL a, b, c, mod = 1e9+7;

LL qmi(LL a, LL b, LL p) {
    LL res = 1;
    while (b) {
        if (b & 1)  res = 1ll * res * a % p;
        b >>= 1;
        a = 1ll * a * a % p;
    }
    return res;
}

LL phi(LL x) {  //返回x的欧拉函数值
    LL res = x;
    for (LL i = 2; i <= x / i; i ++) {
        if (x % i == 0) {
            res = res / i * (i - 1);
            while (x % i == 0)  x /= i;
        }
    }
    if (x > 1)  res = res / x * (x - 1);
    return res;
}

LL sol(LL a, LL p) {   //AC
    if (p == 1) return 0;	//模数为1
    if (a == 1) return 1;	//底数为1
    else if (a == 2)    return 2 % p;
    else if (a == 3)    return 9 % p;
    else if (a == 4)    return 262144 % p;
    
    //5^262144>1e9,超过模数的最大取值，用欧拉降幂
    LL t = phi(p);
    LL b = sol(a - 1, t);
    return qmi(a, b % t + t, p);
}

LL sol2(LL a, LL p) {   //WA
    if (p == 1) return 0;
    if (a == 1) return 1;

    LL t = phi(p);
    LL b = sol(a - 1, t);
    if (b < t)return qmi(a, b % t, p);
    return qmi(a, b % t + t, p);
}

int main() {
    LL n, m;
    cin >> n >> m;
    cout << sol(n, m);
    return 0;
}