uva10951(欧拉降幂)

最新推荐文章于 2021-08-12 14:17:41 发布

原创最新推荐文章于 2021-08-12 14:17:41 发布 · 484 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#欧拉降幂

uva个人训练专栏收录该内容

3 篇文章

订阅专栏

本文介绍了一种利用欧拉降幂公式进行快速幂运算的方法，并通过C++代码实现了该算法。具体而言，文章首先解释了欧拉降幂公式的原理，即如何将指数运算转换为更简单的形式来加速计算过程。然后，给出了一个具体的C++实现示例，包括初始化欧拉函数表、快速幂函数以及读取输入等关键部分。

经典的欧拉降幂问题。
欧拉降幂公式：

n^x mod m = n^{φ(m)+(x mod φ(m))}

这个公式当且仅当x>φ(m)时成立。

//欧拉降幂公式
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
#define MOD(x,m) (x>m)?m+x%m:x
int n, m, a[20];
int phi[10000 + 10], kase;
void phi_table(int n)
{
    for (int i = 2; i <= n; i++)if (!phi[i])
        for (int j = i; j <= n; j += i) {
            if (!phi[j])phi[j] = j;
            phi[j] = phi[j] / i * (i - 1);
        }
}
LL fpow(LL a, LL m, LL mod)
{
    LL ret = 1;
    while (m) {
        if(m&1)ret = MOD(ret * a, mod);
        a = MOD(a * a, mod);
        m >>= 1;
    }
    return ret;
}
bool read()
{
    char ch;
    while (scanf("%c", &ch) && ch == '\n');
    if (ch == '#')return false;
    cin.putback(ch);
    scanf("%d%d", &m, &n);
    for (int i = 0; i < n; i++)scanf("%d", &a[i]);
    return true;
}
LL solve(int x, int mod)
{
    if (x + 1 == n || mod == 1)return MOD(a[x], mod);
    LL tmp = solve(x + 1, phi[mod]);
    return fpow(a[x], tmp, mod);
}
int main()
{
    phi_table(10000);
    while (read()) {
        printf("Case #%d: %lld\n",++kase, solve(0, m) % m);
    }
}