指数循环节（降幂）

最新推荐文章于 2022-04-06 00:18:24 发布

斌同学斌同學

最新推荐文章于 2022-04-06 00:18:24 发布

阅读量823

点赞数 1

分类专栏：数论文章标签： math 数论数学

本文链接：https://blog.csdn.net/GreatJames/article/details/77162961

版权

数论专栏收录该内容

34 篇文章 0 订阅

订阅专栏

指数循环节

在有些题目中我们需要对指数进行降幂处理才能计算。比如计算

其中和

这里由于很大，所以需要进行降幂。那么实际上有如下降幂公式

给定，和的值，求的值，其中，

#include <iostream>
#include <string.h>
#include <stdio.h>

using namespace std;
const int N = 1000005;
typedef long long LL;

char str[N];

int phi(int n)
{
    int rea = n;
    for (int i = 2; i * i <= n; i++)
    {
        if (n % i == 0)
        {
            rea = rea - rea / i;
            while (n % i == 0)
            {
                n /= i;
            }
        }
    }
    if (n > 1)
    {
        rea = rea - rea / n;
    }
    return rea;
}

LL multi(LL a, LL b, LL m)
{
    LL ans = 0;
    a %= m;
    while (b)
    {
        if (b & 1)
        {
            ans = (ans + a) % m;
            b--;
        }
        b >>= 1;
        a = (a + a) % m;
    }
    return ans;
}

LL quick_mod(LL a, LL b, LL m)
{
    LL ans = 1;
    a %= m;
    while (b)
    {
        if (b & 1)
        {
            ans = multi(ans, a, m);
            b--;
        }
        b >>= 1;
        a = multi(a, a, m);
    }
    return ans;
}

void Solve(LL a, char str[], LL c)
{
    LL len = strlen(str);
    LL ans = 0;
    LL p = phi(c);
    if (len <= 15)
    {
        for (int i = 0; i < len; i++)
        {
            ans = ans * 10 + str[i] - '0';
        }
    }
    else
    {
        for (int i = 0; i < len; i++)
        {
            ans = ans * 10 + str[i] - '0';
            ans %= p;
        }
        ans += p;
    }
    printf("%I64d\n", quick_mod(a, ans, c));
}

int main()
{
    LL a, c;
    while (~scanf("%I64d%s%I64d", &a, str, &c))
    {
        Solve(a, str, c);
    }
    return 0;
}