tjut 3609

#include <iostream>  
using namespace std;  
#define LL unsigned __int64  
#define M 205  
  
int phi[M];  
  
int Euler (int n)    //求n的欧拉值  
{  
    int i, res = n;  
    for (i = 2; i * i <= n; i++)  
    {  
        if (n % i == 0)  
        {  
            res = res - res/i;  
            while (n % i == 0)  
                n /= i;  
        }  
    }  
    if (n > 1)  
        res = res - res/n;  
    return res;  
}  
  
LL qmod (LL a, LL b, int c)    //快速幂取模  
{  
    LL res = 1;  
    for ( ; b; b >>= 1)  
    {  
        if (b & 1)  
            res = res * a % c;  
        a = a * a % c;  
    }  
    return res;  
}  
  
LL isok (LL a, LL b, int c)    //目的是判断a^b是否>=c  
{  
    LL res = 1, i;  
    for (i = 0; i < b; i++)  
    {  
        res *= a;  
        if (res >= c)  
            return res;  
    }  
    return res;  
}  
  
LL upup (LL a, int k, int num)    //按照题意递归地使用公式  
{  
    if (phi[num] == 1) return 1;    //显然符合公式条件,很容易套公式知道是1,这实际上是剪枝  
    if (k == 1) return a % phi[num];//返回上一层的a的幂  
    LL b = upup (a, k-1, num+1);    //得到a的幂b  
    LL x = isok (a, b, phi[num]);  
    if (x >= phi[num])    //使用公式的条件,满足条件可以进行降幂,返回的是上层a的幂  
        return qmod (a % phi[num], b, phi[num]) + phi[num];  
    else return x;    //不使用公式,直接返回a^b,也属于上层a的幂  
}  
  
int main()  
{  
    LL a;  
    int k;  
    phi[0] = 100000000;  
    for (k = 1; k < M; k++)  
        phi[k] = Euler (phi[k-1]);    //预处理所需欧拉值  
    while (~scanf ("%I64u%d", &a, &k))  
    {  
        printf ("%I64u\n", upup (a, k, 0) % phi[0]);  
    }  
    return 0;  
} 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值