IT数码 购物 网址 头条 软件 日历 阅读 图书馆
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
图片批量下载器
↓批量下载图片,美女图库↓
图片自动播放器
↓图片自动播放器↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁
 
   -> 数据结构与算法 -> CF906D:Power Tower(拓展欧拉定理、暴力) -> 正文阅读

[数据结构与算法]CF906D:Power Tower(拓展欧拉定理、暴力)

解析

大降智题。
受相逢是问候的影响,第一眼直接线段树,写完WA完才发现这个玩意没有结合律。
然后开始用 x y z = ( x y ) y z ? 1 x^{y^z}=(x^y)^{y^{z-1}} xyz=(xy)yz?1 这样的东西嗯做,结果只是在玩泥巴。

lemma:迭代求 φ \varphi φ 的次数是 log 级别的。

证明:偶数的 φ \varphi φ 至少折半,对于 n > 2 n>2 n>2 φ ( n ) \varphi(n) φ(n) 恒为偶数,因为 gcd ? ( i , n ) = g c d ( n ? i , n ) \gcd(i,n)=gcd(n-i,n) gcd(i,n)=gcd(n?i,n)

所以直接暴力往后扫模数为1就返回就可以了。

//luogu
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define ull unsigned long long
#define debug(...) fprintf(stderr,__VA_ARGS__)
#define ok debug("OK\n")
using namespace std;

const int N=1e5+100;
//const int mod=998244353;
inline ll read(){
  ll x(0),f(1);char c=getchar();
  while(!isdigit(c)) {if(c=='-')f=-1;c=getchar();}
  while(isdigit(c)) {x=(x<<1)+(x<<3)+c-'0';c=getchar();}
  return x*f;
}

#define pr pair<ll,bool>
#define mkp make_pair
inline pr ksm(ll x,ll k,int mod){
  ll res(1);bool flag=0;  
  while(k){
    if(k&1){
      res=res*x;
      if(res>=mod){
	res%=mod;flag=1;
      }
    }
    x=x*x;
    if(x>=mod&&k>1){
      x%=mod;
      flag=1;
    }
    k>>=1;
  }
  return mkp(res,flag);
}

map<int,int>mp;
int n,m;
int w[N];

inline int getphi(int x){
  if(mp.count(x)) return mp[x];
  int ori=x;
  int w=sqrt(x),res=x;
  for(int i=2;i<=w;i++){
    if(x%i==0){
      res=1ll*res*(i-1)/i;
      while(x%i==0) x/=i;
    }
  }
  if(x>1) res=1ll*res*(x-1)/x;
  return mp[ori]=res;
}

pr calc(int l,int r,int mod){
  if(l==r){return mkp(w[l]%mod,w[l]>=mod);}
  if(mod==1) return mkp(0,1);
  int phi=getphi(mod);
  pr o=calc(l+1,r,phi);
  return ksm(w[l],o.first+o.second*phi,mod);
}

signed main(){
#ifndef ONLINE_JUDGE
  freopen("a.in","r",stdin);
  freopen("a.out","w",stdout);
#endif
  n=read();int mod=read();
  for(int i=1;i<=n;i++) w[i]=read();
  m=read();
  for(int i=1;i<=m;i++){
    int l=read(),r=read();
    printf("%lld\n",calc(l,r,mod).first);
  }
  return 0;
}
/*
 */
  数据结构与算法 最新文章
【力扣106】 从中序与后续遍历序列构造二叉
leetcode 322 零钱兑换
哈希的应用:海量数据处理
动态规划|最短Hamilton路径
华为机试_HJ41 称砝码【中等】【menset】【
【C与数据结构】——寒假提高每日练习Day1
基础算法——堆排序
2023王道数据结构线性表--单链表课后习题部
LeetCode 之 反转链表的一部分
【题解】lintcode必刷50题<有效的括号序列
上一篇文章      下一篇文章      查看所有文章
加:2022-05-05 11:44:27  更:2022-05-05 11:47:47 
 
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁

360图书馆 购物 三丰科技 阅读网 日历 万年历 2024年11日历 -2024/11/26 5:19:57-

图片自动播放器
↓图片自动播放器↓
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
图片批量下载器
↓批量下载图片,美女图库↓
  网站联系: qq:121756557 email:121756557@qq.com  IT数码