Uva--10692--Huge Mods【数论】

最新推荐文章于 2017-11-12 19:16:12 发布

altair21

最新推荐文章于 2017-11-12 19:16:12 发布

阅读量1k

点赞数

分类专栏： Uva代码文章标签：数论指数循环节

本文链接：https://blog.csdn.net/zzzz40/article/details/17126471

版权

Uva代码专栏收录该内容

11 篇文章

订阅专栏

博客探讨了指数循环节的概念，介绍了公式A^x % m = A^(x%phi(m)+phi(m)) % m，并解释了如何利用欧拉函数phi(m)进行递归计算。在递归过程中，模值m随着层次变化，对应于不同层级的m的欧拉函数值。

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

2^3^4^5 mod 10 = ?

这个东西叫指数循环节，有个公式：

A^x % m = A^(x%phi(m)+phi(m)) % m (x >= phi(m))

其中phi(m)是m的欧拉函数值，然后用递归层层推下来

递归的时候mod的值会变化，因为公式里等号右边x mod的是外面m的欧拉函数值，再往深层递归m的欧拉函数值又变成了上一层的m，此时等号右边x mod的是m的欧拉值的欧拉值，所以mod 的m是变化的

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<vector>
#include<set>
#include<queue>
#include<stack>
#include<stdlib.h>
#include<ctime>
#define PI acos(-1)
using namespace std;
int a[10];
int mm;
int eular(int n)
{
    int ret=1,i;
    for(i=2; i*i<=n; i++)
        if(n%i==0)
        {
            n/=i,ret*=i-1;
            while(n%i==0)
                n/=i,ret*=i;
        }
    if(n>1)
        ret*=n-1;
    return ret;
}
int mode(int a,int n,int m)
{
    int t = a;
    int ans = 1;
    while(n)
    {
        if(n & 1)
        {
            ans = ans * t % m;
        }
        n >>= 1;
        t =  t * t % m;
    }
    return ans;
}
int solve(int c,int n,int m)
{
    if(c==n-1) return a[c]%m;
    int temp=solve(c+1,n,eular(m));
    int ans=mode(a[c],temp+eular(m),m);
    return ans;
}
int main()
{
    int ans;
    int m,n,i,k=1;
    char c;
    while(cin>>m)
    {
        c=getchar();
        if(c=='#')  break;
        mm=eular(m);
        cin>>n;
        for(i=0; i<n; i++)
            scanf("%d",&a[i]);
        ans=solve(0,n,m);
        printf("Case #%d: %d\n",k++,ans);
    }
    return 0;
}