题意

如果对一个数操作\(k\)次,那么这个数会变成\(c^{c^{...^{a_i}}}\),其中\(c\)有\(k\)个。

根据P4139 上帝与集合的正确用法这道题,我们可以知道一个数不断变为自己的欧拉函数,大约\(log\)次就会变成1,而任何数模\(1\)都是\(0\),于是我们可以用势能线段树解决。

因为模数不变,因此我们可以预处理所有\(\varphi(\varphi(...\varphi(p)...))\),之后在线段树上记录操作次数。

这样是三个\(log\)的,因为还要快速幂,可以对每个\(\varphi(\varphi(...\varphi(p)...))\)预处理,用光速幂解决。

注意,扩展中国剩余定理\(a_k\equiv a^{k\%\varphi(p)+\varphi(p)}\pmod{p}\)适用当且仅当\(k\geqslant \varphi(p)\),因此我们在求值时用一个\(flag\)表示是否要\(+\varphi(p)\)。

code:

#include<bits/stdc++.h>
using namespace std;
#define int long long
#define ls(p) (p<<1)
#define rs(p) (p<<1|1)
const int maxn=50010;
const int maxt=10010;
int n,m,mod,C,maxtim;
int a[maxn];
int pw1[60][maxt],pw2[60][maxt];
bool flag;
bool flag1[60][maxt],flag2[60][maxt];
vector<int>ve;
struct Seg
{
#define sum(p) (seg[p].sum)
#define cnt(p) (seg[p].cnt)
int sum,cnt;
}seg[maxn<<2];
inline int read()
{
char c=getchar();int res=0,f=1;
while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
while(c>='0'&&c<='9')res=res*10+c-'0',c=getchar();
return res*f;
}
inline int phi(int x)
{
int res=x,tmp=x;
for(int i=2;i*i<=tmp;i++)
{
if(tmp%i)continue;
res=res/i*(i-1);
while(tmp%i==0)tmp/=i;
}
if(tmp>1)res=res/tmp*(tmp-1);
return res;
}
inline void pre_work()
{
int tmp=mod;
ve.push_back(mod);
while(tmp>1)tmp=phi(tmp),ve.push_back(tmp);
ve.push_back(1);
for(unsigned int i=0;i<ve.size();i++)
{
pw1[i][0]=1;
for(int j=1;j<=10000;j++)
{
pw1[i][j]=pw1[i][j-1]*C;
if(pw1[i][j]>=ve[i])pw1[i][j]%=ve[i],flag1[i][j]=1;
flag1[i][j]|=flag1[i][j-1];
}
}
for(unsigned int i=0;i<ve.size();i++)
{
pw2[i][0]=1;flag2[i][1]=flag1[i][10000];
for(int j=1;j<=10000;j++)
{
pw2[i][j]=pw2[i][j-1]*pw1[i][10000];
if(pw2[i][j]>=ve[i])pw2[i][j]%=ve[i],flag2[i][j]=1;
flag2[i][j]|=flag2[i][j-1];
}
}
}
inline void up(int p)
{
sum(p)=(sum(ls(p))+sum(rs(p)))%mod;
cnt(p)=min(cnt(ls(p)),cnt(rs(p)));
}
void build(int p,int l,int r)
{
if(l==r){sum(p)=a[l];return;}
int mid=(l+r)>>1;
build(ls(p),l,mid);build(rs(p),mid+1,r);
up(p);
}
inline int power(int x,int id)
{
flag=0;
int res=pw1[id][x%10000]*pw2[id][x/10000];
if(res>=ve[id])res%=ve[id],flag=1;
flag|=flag1[id][x%10000]|flag2[id][x/10000];
return res;
}
int calc(int x,int dep,int k)
{
flag=0;
if(dep==k)
{
if(x>=ve[dep])flag=1,x%=ve[dep];
return x;
}
int tmp=calc(x,dep+1,k);
return power(flag?tmp+ve[dep+1]:tmp,dep);
}
void change(int p,int l,int r,int ql,int qr)
{
if(cnt(p)>=(int)ve.size()-1)return;
if(l==r)
{
cnt(p)++;
sum(p)=calc(a[l],0,cnt(p));
return;
}
int mid=(l+r)>>1;
if(ql<=mid)change(ls(p),l,mid,ql,qr);
if(qr>mid)change(rs(p),mid+1,r,ql,qr);
up(p);
}
int query(int p,int l,int r,int ql,int qr)
{
if(l>=ql&&r<=qr)return sum(p);
int mid=(l+r)>>1,res=0;
if(ql<=mid)res=(res+query(ls(p),l,mid,ql,qr))%mod;
if(qr>mid)res=(res+query(rs(p),mid+1,r,ql,qr))%mod;
return res;
}
signed main()
{
//freopen("test.in","r",stdin);
//freopen("test.out","w",stdout);
n=read();m=read();mod=read();C=read();
for(int i=1;i<=n;i++)a[i]=read();
build(1,1,n);
pre_work();
for(int i=1;i<=m;i++)
{
int op=read(),l=read(),r=read();
if(!op)change(1,1,n,l,r);
else printf("%lld\n",query(1,1,n,l,r));
}
return 0;
}

P3747 [六省联考2017]相逢是问候的更多相关文章

  1. 洛谷 P3747 [六省联考2017]相逢是问候 解题报告

    P3747 [六省联考2017]相逢是问候 题目描述 \(\text {Informatik verbindet dich und mich.}\) 信息将你我连结. \(B\) 君希望以维护一个长度 ...

  2. 洛谷P3747 [六省联考2017]相逢是问候

    传送门 题解 扩展欧拉定理. 线段树维护,已经全改到底了的节点就不管,不然暴力修改下去. //Achen #include<algorithm> #include<iostream& ...

  3. [BZOJ4869][六省联考2017]相逢是问候(线段树+扩展欧拉定理)

    4869: [Shoi2017]相逢是问候 Time Limit: 40 Sec  Memory Limit: 512 MBSubmit: 1313  Solved: 471[Submit][Stat ...

  4. bzoj千题计划271:bzoj4869: [六省联考2017]相逢是问候

    http://www.lydsy.com/JudgeOnline/problem.php?id=4869 欧拉降幂+线段树,每个数最多降log次,模数就会降为1 #include<cmath&g ...

  5. BZOJ4869 六省联考2017相逢是问候(线段树+欧拉函数)

    由扩展欧拉定理,a^(a^(a^(……^x)))%p中x作为指数的模数应该是φ(φ(φ(φ(……p)))),而p取log次φ就会变为1,也即每个位置一旦被修改一定次数后就会变为定值.线段树维护区间剩余 ...

  6. 【LuoguP3747】[六省联考2017] 相逢是问候

    题目链接 题意 给定一个长度为 n 的序列 a , 给定一个正整数 c 每次修改操作是把一段区间内的数 \(x_i\) 修改为 \(c^{x_i}\) 询问区间和模 p 的结果 Sol 修改是把一个数 ...

  7. 2017 [六省联考] T2 相逢是问候

    4869: [Shoi2017]相逢是问候 Time Limit: 40 Sec  Memory Limit: 512 MBSubmit: 1205  Solved: 409[Submit][Stat ...

  8. 【BZOJ4873】[六省联考2017]寿司餐厅(网络流)

    [BZOJ4873][六省联考2017]寿司餐厅(网络流) 题面 BZOJ 洛谷 题解 很有意思的题目 首先看到答案的计算方法,就很明显的感觉到是一个最大权闭合子图. 然后只需要考虑怎么构图就行了. ...

  9. 【BZOJ4868】[六省联考2017]期末考试(贪心)

    [BZOJ4868][六省联考2017]期末考试(贪心) 题面 BZOJ 洛谷 题解 显然最终的答案之和最后一个公布成绩的课程相关. 枚举最后一天的日期,那么维护一下前面有多少天可以向后移,后面总共需 ...

随机推荐

  1. win7再分配磁盘新加卷

    磁盘在系统刚分区的时候可以做磁盘分区最好 1.右键我的电脑,选在管理 2.在此窗口下依次展开选项,点击存储->磁盘管理,右边是我已经分好的盘不用看的 3.确认一下我的电脑的各个盘的空间,选择要压 ...

  2. android studio 软件常见问题

    xml文件没有智能提示 在网上看到说是去掉省电模式就可以了 然后我试了一下并不能解决, 最终我是这么解决的

  3. spring+activemq实战之配置监听多队列实现不同队列消息消费

    摘选:https://my.oschina.net/u/3613230/blog/1457227 摘要: 最近在项目开发中,需要用到activemq,用的时候,发现在同一个项目中point-to-po ...

  4. 一种分片更新ubi卷的方式(基于ubiupdatevol,拓展fifo支持)

    ubi卷的更新方式 对于ubi卷,不能像普通块设备一样进行随机读写.每次更新需要从头写入. 具体的,需要在打开对应的设备之后,先执行一个ioctl UBI_IOCVOLUP,同时传入要更新的数据大小. ...

  5. 6、UnityConfig实现AOP

    需求:我们需要给已经开发好的服务如这里的UserService,添加额外的执行逻辑,但是又不想破坏原有的服务,如:我们需要给UserService添加监控逻辑,监控的目的是看UserService服务 ...

  6. python网络编程socketserver模块(实现TCP客户端/服务器)

    摘录python核心编程 socketserver(python3.x版本重新命名)是标准库中的网络编程的高级模块.通过将创建网络客户端和服务器所必须的代码封装起来,简化了模板,为你提供了各种各样的类 ...

  7. Python面向对象-定制方法

    Python中的class可以定义许多定制方法,可以让我们方便的生成特定的类. 我们之前介绍了__slots__.__len__(),python中还有许多这样的特殊函数: __str__ >& ...

  8. grid行高亮显示

    grid var grid = new Ext.RSEGrid({ //数据显示面板 id: 'O005_grid', title: '图纸目录', loadMask: true, autoScrol ...

  9. 微信小程序—支付宝身份验证(支付宝小程序)

    查看应用:https://open.alipay.com/platform/keyManage.htm  这里找到您调用接口的应用 支付宝身份验证快速接入:https://docs.open.alip ...

  10. 65-如何部署 Calico 网络?

    Calico 是一个纯三层的虚拟网络方案,Calico 为每个容器分配一个 IP,每个 host 都是 router,把不同 host 的容器连接起来.与 VxLAN 不同的是,Calico 不对数据 ...