CodeForces Gym101550E Exponial

Toooooocold

于 2017-11-24 11:16:08 发布

阅读量446

点赞数

CC 4.0 BY-SA版权

分类专栏：数论------------------------- 数论-降幂公式

本文链接：https://blog.youkuaiyun.com/qq_32506797/article/details/78622699

数论------------------------- 同时被 2 个专栏收录

22 篇文章

订阅专栏

数论-降幂公式

2 篇文章

订阅专栏

本文介绍了一个关于模指数运算的问题，并通过欧拉定理来简化计算过程。使用了模幂运算、欧拉函数等算法，适用于解决大规模整数幂次模运算问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

$Problem:$ $F(n) = n^{F_(n-1)} \bmod m$

$Idea:$ 欧拉定理降幂公式： $a^b\bmod p = a^{\phi p+b\bmod{\phi p}}\bmod p (b \ge \phi p)$

$Code:$

#include<bits/stdc++.h>
using namespace std;

#define fi first
#define se second
#define pb push_back
#define CLR(A, X) memset(A, X, sizeof(A))
#define bitcount(X) __builtin_popcountll(X)
typedef long long LL;
typedef pair<int, int> PII;
const double eps = 1e-10;
const double PI = acos(-1.0);
const auto INF = 0x3f3f3f3f;
int dcmp(double x) { if(fabs(x) < eps) return 0; return x<0?-1:1; }

LL qpow(LL a, LL b, LL M) {
    LL ret = 1;
    while(b) {
        if(b & 1) ret = ret*a%M;
        a = a*a%M;
        b >>= 1;
    }
    return ret;
}

LL phi(LL n) {
    LL ret = n;
    for(LL i = 2; i*i <= n; i++) {
        if(n%i) continue;
        while(n%i == 0) n /= i;
        ret = ret/i*(i-1);
    }
    if(n != 1) ret = ret/n*(n-1);
    return ret;
}

LL F(LL n, LL m) {
    if(m == 1) return 0;
    if(n == 1) return 1%m;
    if(n == 2) return 2%m;
    if(n == 3) return 9%m;
    if(n == 4) return 262144%m;
    LL p = phi(m);
    return qpow(n, p+F(n-1, p), m);
}

int main() {
    LL n, m;
    scanf("%lld%lld", &n, &m);
    printf("%lld\n", F(n, m));
    return 0;
}