数论--阶乘幂&扩展欧拉定理--codeforces 906d Power Tower

最新推荐文章于 2021-04-20 23:23:28 发布

Falling～

最新推荐文章于 2021-04-20 23:23:28 发布

阅读量564

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/SM_545/article/details/81141794

版权

数论专栏收录该内容

46 篇文章 3 订阅

订阅专栏

这篇博客介绍了如何利用扩展欧拉定理解决数论问题，特别是针对给定数组的指数运算查询。通过快速幂算法和递归方法，有效地计算了数组元素之间的幂运算并取模的结果。同时，文章提供了一个名为`qk_pow`的快速幂函数实现，以及一个`f`函数来处理查询过程中的递归计算。

摘要由CSDN通过智能技术生成

给定数组w,q次查询，求w[ql]^ ( w[ql + 1] ^ ( ... w[qr - 1] ^ (w[qr]) ) ) % m

扩展欧拉定理 (a,m) != 1时，

if(b >= phi(m)) a^b % m = a ^ (b % phi(m) + phi(m) ) % m

if(b < phi(m)) a^b % m = a ^ b % m（快速幂）

1.快速幂和递归时要时刻注意指数和phi(m)的大小关系

ll qk_pow(ll a,ll b,ll p)

{

ll ans = 1;

a = a >= p ? a % p + p : a;

while (b) {

if(b & 1) ans = ans * a >= p ? ans * a % p + p : ans * a % p;

b >>= 1;

a = a * a >= p ? a * a % p + p : a * a % p;

}

return ans;

}

ll f(int ql,int qr,int pos)

{

int m = v[pos];

if(m == 1) return 1;

if(ql == qr) return w[ql] < m ? w[ql] : w[ql] % m + m;

ll x = f(ql + 1,qr,pos + 1);

ll ans = qk_pow(w[ql], x , m);

return ans;

}


#include <iostream>
#include <cstdio>
#include <vector>
#include <cmath>
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 5;
int n,m;
int w[maxn];

vector<int> v;

int phi(int m)
{
    int ans = m;
    for (int i = 2; i <= sqrt(0.5 + m); i ++) {
        if(m % i == 0){
            ans = ans / i * (i - 1);
            while(m % i == 0 ) m /= i;
        }
    }
    if(m > 1) ans = ans / m * (m - 1);
    return ans;
}

void init(int m)
{
    v.clear();
//预处理phi(m)
    int x = m;
    v.push_back(x);
    while (1) {
        x = phi(x);
        v.push_back(x);
        if(x == 1) break;
    }
}
ll qk_pow(ll a,ll b,ll p)
{
    ll ans = 1;
    a = a >= p ? a % p + p : a;
    while (b) {
        if(b & 1) ans = ans * a >= p ? ans * a % p + p : ans * a % p;
        b >>= 1;
        a = a * a >= p ? a * a % p + p : a * a % p;
    }
    return ans;
}

ll f(int ql,int qr,int pos)
{
    int m = v[pos];
    if(m == 1) return 1;
    if(ql == qr) return w[ql] < m ? w[ql] : w[ql] % m + m;
    ll x = f(ql + 1,qr,pos + 1);
    ll ans = qk_pow(w[ql], x , m);
    return ans;
}
int main()
{
   
    while (scanf("%d%d",&n,&m) != EOF) {
         init(m);
        for (int i = 1; i <= n; i ++) {
            scanf("%d",&w[i]);
        }
        int q;scanf("%d",&q);
        int ql,qr;
        ll ans;
        while (q --) {
            scanf("%d%d",&ql,&qr);
             ans = f(ql,qr,0);
            ans = (ans + m) % m;
            printf("%I64d\n",ans);
        }
    }
    return 0;
}