首页 > 其他 > 详细

Codeforces Round #454 (Div. 1, based on Technocup 2018 Elimination Round 4) CodeForces 906D Power Tower (欧拉降幂)

时间:2019-05-07 13:34:47      阅读:271      评论:0      收藏:0      [点我收藏+]

题目链接:http://codeforces.com/contest/906/problem/D

题目大意:给定n个整数w[1],w[2],……,w[n],和一个数m,然后有q个询问,每个询问给出一个l,r,求w[l]^w[l+1]^w[l+2]……w[r]  %m  ,即a[l]到a[r]的幂次方

解题思路:利用欧拉降幂公式

技术分享图片

第一个要求a和p互质,第2个和第3个为广义欧拉降幂,不要求a和p互质,用在这题刚好。

因为有两种情况,所以我们需要自定义一下降幂取模公式。

我们对整个区间进行递归处理,每一个数的指数是它后一个数到右端点的幂。

递归终止条件为到右端点或者p的欧拉函数值为1,再求欧拉函数值的时候我们需要进行记忆化,否则会超时

代码:

#include<iostream>
#include<cstdio>
#include<map>
using namespace std;
#define ll long long
#define MOD(a,b) a>=b?a%b+b:a
#define N 100005
map<ll,ll> mp;
int n,l,r,q;
ll mod,w[N];
ll qpow(ll a,ll b,ll p){
    ll res=1;
    while(b){
        if(b&1) res=MOD(res*a,p); //为保证指数结果正确,应该用自定义取模 
        b>>=1;
        a=MOD(a*a,p);
    }
    return res;
}
ll phi(ll x){
    if(mp[x]) return mp[x];
    ll tmp=x,res=x;
    for(int i=2;i*i<=x;i++){
        if(x%i==0){
            res=res*(i-1)/i;
            while(x%i==0) x/=i;
        }
    }
    if(x>1) res=res*(x-1)/x;
    return mp[tmp]=res;
}
ll solve(int l,int r,ll m){
    if(l==r||m==1) return MOD(w[l],m);
    else return qpow(w[l],solve(l+1,r,phi(m)),m);
}
int main() {
    scanf("%d%I64d",&n,&mod);
    for(int i=1;i<=n;i++) scanf("%I64d",&w[i]);
    scanf("%d",&q);
    while(q--){
        scanf("%d%d",&l,&r);
        printf("%I64d\n",solve(l,r,mod)%mod);
    }
    return 0;
}

 

 

bzoj 3884 上帝与集合的正确用法

题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=3884

题目大意:和上题很像,只不过所有数都是2,且次方是无穷的了,给定一个正整数p,求2^(2^(2^(2^(2^...)))) mod p的值

解题思路:方法几乎是一样的,因为每次递归幂的模数就会变成原来的欧拉函数值,所以最多经过log(p),模数就会变成1,然后后面结果都一样的了,没必要递归下去,直接结束就好了。

代码:

#include<iostream>
#include<cstdio>
#include<map>
using namespace std;
#define ll long long
#define MOD(a,b) a>=b?a%b+b:a
#define N 100005
map<ll,ll> mp;
int n,l,r,q;
ll mod;
ll qpow(ll a,ll b,ll p){
    ll res=1;
    while(b){
        if(b&1) res=MOD(res*a,p);  
        b>>=1;
        a=MOD(a*a,p);
    }
    return res;
}
ll phi(ll x){
    if(mp[x]) return mp[x];
    ll tmp=x,res=x;
    for(int i=2;i*i<=x;i++){
        if(x%i==0){
            res=res*(i-1)/i;
            while(x%i==0) x/=i;
        }
    }
    if(x>1) res=res*(x-1)/x;
    return mp[tmp]=res;
}
ll solve(ll m){
    if(m==1) return 1;
    else return qpow(2,solve(phi(m)),m);
}
int main() {
    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%lld",&mod);
        printf("%lld\n",solve(mod)%mod);
    }
    return 0;
}

 

Codeforces Round #454 (Div. 1, based on Technocup 2018 Elimination Round 4) CodeForces 906D Power Tower (欧拉降幂)

原文:https://www.cnblogs.com/zjl192628928/p/10824842.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!