洛谷3934:Nephren Ruq Insania——题解
https://www.luogu.org/problemnew/show/P3934
题面自己读吧(滑稽。
看到这道题就能够想到BZOJ4869:[SHOI2017]相逢是问候我们曾经用过的哲学扩展欧拉定理。
(咦什么时候这个东西都普及到noip了好方啊)
也就是说,不论询问的区间[l,r]长度有多大,实际上我们暴力算只需要logp次模数就变成了1之后的询问就不需要担心了。
区间修改可以线段树/树状数组来干。
敲完发现是70pts,很难受。
有没有发现这个公式前面有一个前提?其实,当a^b<p的时候,a的指数就不应当加phi(p)了,这也就是错的原因。
我们当然可以每次传递当前的b是否曾经对phi(p)取过模来决定是否加phi(p),经过多次修改之后终于AC了……
PS:请注意我们运算的时候很可能出现x^0%x=1这种情况,而显然我们希望的答案应当为0(虽然我并不知道为什么),所以需要特判之。
#include<cmath>
#include<queue>
#include<vector>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=5e5+;
const int M=2e7+;
inline int read(){
int X=,w=;char ch=;
while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();
return w?-X:X;
}
bool ok[N];
ll qpow(ll k,ll n,int p,int l){
ll res=;
while(n){
if(n&){
res*=k;
if(res>=p)ok[l]=,res%=p;
}
k*=k;n>>=;
if(k>=p)ok[l]=,k%=p;
}
return res;
}
int n,m;
int su[M],he[M],phi[M],tot;
ll tr[N*],lz[N*],b[N];
void Euler(int n){
phi[]=;
for(int i=;i<=n;i++){
if(!he[i])su[++tot]=i,phi[i]=i-;
for(int j=;j<=tot;j++){
if(i*su[j]>n)break;
he[i*su[j]]=;
if(i%su[j]==){
phi[i*su[j]]=phi[i]*su[j];
break;
}else phi[i*su[j]]=phi[i]*phi[su[j]];
}
}
}
inline void push(int a,int l,int r){
if(!lz[a])return;
int mid=(l+r)>>;
lz[a<<]+=lz[a];lz[a<<|]+=lz[a];
if(l==mid)tr[a<<]+=lz[a];
if(mid+==r)tr[a<<|]+=lz[a];
lz[a]=;
}
void build(int a,int l,int r){
if(l==r){
tr[a]=b[l];return;
}
int mid=(l+r)>>;
build(a<<,l,mid);build(a<<|,mid+,r);
}
void add(int a,int l,int r,int l1,int r1,int x){
if(r<l1||r1<l)return;
if(l1<=l&&r<=r1){
lz[a]+=x;
if(l==r)tr[a]+=x;
return;
}
push(a,l,r);
int mid=(l+r)>>;
add(a<<,l,mid,l1,r1,x);add(a<<|,mid+,r,l1,r1,x);
}
ll ask(int a,int l,int r,int k){
if(l==r)return tr[a];
push(a,l,r);
int mid=(l+r)>>;
if(k<=mid)return ask(a<<,l,mid,k);
else return ask(a<<|,mid+,r,k);
}
ll query(int l,int r,int p){
ok[l]=;
ll ans=ask(,,n,l);
if(ans>=p)ok[l]=;
if(ans%p==)return ;
if(p==)return ;
ans%=p;
if(l==r)return ans;
ll tmp=query(l+,r,phi[p]);
return qpow(ans,tmp+(ok[l+]?phi[p]:),p,l);
}
int main(){
Euler(2e7);
n=read(),m=read();
for(int i=;i<=n;i++)b[i]=read();
build(,,n);
for(int i=;i<=m;i++){
int op=read();
if(op==){
int l=read(),r=read(),x=read();
add(,,n,l,r,x);
}else{
int l=read(),r=read(),p=read();
printf("%lld\n",query(l,r,p));
}
}
return ;
}
当然还有另外一种(正解)的做法。
考虑到假设我们暴力枚举模拟前几次操作之后,往后的操作基本都是大于phi(p)的,而没有加phi(p)的只可能是前面几项。
故暴力处理前面几项即可,常数显然要比上面的方法要大很多。
#include<cmath>
#include<queue>
#include<vector>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int N=5e5+;
const int M=2e7+;
inline int read(){
int X=,w=;char ch=;
while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();
return w?-X:X;
}
int qpow(int k,ll n,int p){
int res=;
while(n){
if(n&)res=(ll)res*k%p;
k=(ll)k*k%p;n>>=;
}
return res;
}
int n,m,tag[N];
int su[M],he[M],phi[M],tot;
ll tr[N*],lz[N*],b[N];
void Euler(int n){
phi[]=;
for(int i=;i<=n;i++){
if(!he[i])su[++tot]=i,phi[i]=i-;
for(int j=;j<=tot;j++){
if(i*su[j]>n)break;
he[i*su[j]]=;
if(i%su[j]==){
phi[i*su[j]]=phi[i]*su[j];
break;
}else phi[i*su[j]]=phi[i]*phi[su[j]];
}
}
}
inline void push(int a,int l,int r){
if(!lz[a])return;
int mid=(l+r)>>;
lz[a<<]+=lz[a];lz[a<<|]+=lz[a];
if(l==mid)tr[a<<]+=lz[a];
if(mid+==r)tr[a<<|]+=lz[a];
lz[a]=;
}
inline void build(int a,int l,int r){
if(l==r){
tr[a]=b[l];return;
}
int mid=(l+r)>>;
build(a<<,l,mid);build(a<<|,mid+,r);
}
inline void add(int a,int l,int r,int l1,int r1,int x){
if(r<l1||r1<l)return;
if(l1<=l&&r<=r1){
lz[a]+=x;
if(l==r)tr[a]+=x;
return;
}
push(a,l,r);
int mid=(l+r)>>;
add(a<<,l,mid,l1,r1,x);add(a<<|,mid+,r,l1,r1,x);
}
inline ll ask(int a,int l,int r,int k){
if(l==r)return tr[a];
push(a,l,r);
int mid=(l+r)>>;
if(k<=mid)return ask(a<<,l,mid,k);
else return ask(a<<|,mid+,r,k);
}
inline ll qry(int l){
if(tag[l]==m)return b[l];
tag[l]=m;
return b[l]=ask(,,n,l);
}
inline ll query(int l,int r,int p){
ll ans=qry(l);
if(ans%p==)return ;
if(p==)return ;
if(l==r)return (ans%p+(ans>p)*p)%p;
int f=min(r,l+);
for(int i=l+;i<=f;i++){
if(qry(i)==){
f=i;break;
}
}
ll t=qry(f),q=;
for(int i=f-;i>=l+;i--){
q=t,t=;
ll tmp=qry(i);
while(q--){
t*=tmp;
if(t>=phi[p])return qpow(ans%p,query(l+,r,phi[p])+phi[p],p);
}
}
return qpow(ans%p,t,p);
}
int main(){
Euler(2e7);memset(tag,-,sizeof(tag));
n=read(),m=read();
for(int i=;i<=n;i++)b[i]=read();
build(,,n);
while(m--){
int op=read();
if(op==){
int l=read(),r=read(),x=read();
add(,,n,l,r,x);
}else{
int l=read(),r=read(),p=read();
printf("%lld\n",query(l,r,p));
}
}
return ;
}
+++++++++++++++++++++++++++++++++++++++++++
+本文作者:luyouqi233。 +
+欢迎访问我的博客:http://www.cnblogs.com/luyouqi233/+
+++++++++++++++++++++++++++++++++++++++++++
洛谷3934:Nephren Ruq Insania——题解的更多相关文章
- [洛谷P4118][Ynoi2016]炸脖龙I([洛谷P3934]Nephren Ruq Insania)
题目大意:有$n$个数,每个数为$s_i$,两个操作: $1\;l\;r\;x:$表示将区间$[l,r]$内的数加上$x$ $2\;l\;r\;p:$表示求$s_l^{s_{l+1}^{^{s_{l+ ...
- [洛谷3934]P3934 Nephren Ruq Insania题解
先放个奈芙莲 解法 看到这种题目就知道是欧拉降幂,然后根据某玄学证明,递归欧拉降幂从l到r不会超过\(\Theta(log_n)\),所以直接递归解决,然后区间修改直接树状数组维护一下 然后就A了 代 ...
- Luogu 3934 Nephren Ruq Insania
和Ynoi2016 炸脖龙重题了. BZOJ 5394. 首先是扩展欧拉定理: 一开始傻掉了……递归的层数和区间长度无关……也就是说我们每一次直接暴力递归求解子问题一定不会超过$logP$层,因为当模 ...
- 洛谷P1783 海滩防御 分析+题解代码
洛谷P1783 海滩防御 分析+题解代码 题目描述: WLP同学最近迷上了一款网络联机对战游戏(终于知道为毛JOHNKRAM每天刷洛谷效率那么低了),但是他却为了这个游戏很苦恼,因为他在海边的造船厂和 ...
- 洛谷P4047 [JSOI2010]部落划分题解
洛谷P4047 [JSOI2010]部落划分题解 题目描述 聪聪研究发现,荒岛野人总是过着群居的生活,但是,并不是整个荒岛上的所有野人都属于同一个部落,野人们总是拉帮结派形成属于自己的部落,不同的部落 ...
- 洛谷P1155 双栈排序题解(图论模型转换+二分图染色+栈)
洛谷P1155 双栈排序题解(图论模型转换+二分图染色+栈) 标签:题解 阅读体验:https://zybuluo.com/Junlier/note/1311990 原题地址:洛谷P1155 双栈排序 ...
- 洛谷10月月赛II题解
[咻咻咻] (https://www.luogu.org/contestnew/show/11616) 令人窒息的洛谷月赛,即将参加NOIp的我竟然只会一道题(也可以说一道也不会),最终145的我只能 ...
- [洛谷P1823]音乐会的等待 题解(单调栈)
[洛谷P1823]音乐会的等待 Description N个人正在排队进入一个音乐会.人们等得很无聊,于是他们开始转来转去,想在队伍里寻找自己的熟人.队列中任意两个人A和B,如果他们是相邻或他们之间没 ...
- BZOJ2527 & 洛谷3527:[Poi2011]Meteors——题解
+++++++++++++++++++++++++++++++++++++++++++ +本文作者:luyouqi233. + +欢迎访问我的博客:http://www.cnblogs.com/luy ...
随机推荐
- 聊聊WS-Federation
本文来自网易云社区 单点登录(Single Sign On),简称为 SSO,目前已经被大家所熟知.简单的说, 就是在多个应用系统中,用户只需要登录一次就可以访问所有相互信任的应用系统. 举例: 我们 ...
- 2019年猪年海报PSD模板-第四部分
14套精美猪年海报,免费猪年海报,下载地址:百度网盘,https://pan.baidu.com/s/1WUO4L5PHIHG5hAurv52_2A
- spring + mybatis 注解 @Transactional失效
1.问题 在使用@Transactional注解管理事务的时候会出现很多错误,比如: *** was not registered for synchronization because synchr ...
- Struts2(九.初始化用户列表时显示用户照片数目)
1.userlist.jsp //显示每个用户照片的数目(遍历每个用户) $(".picture").each(function(i,e){ $.post("${page ...
- 手动在Windows上创建kafka环境
一 主要内容请移步 参考文章 二 一点小问题 安装上面的参考文章完成配置,并且尝试传输消息后,我尝试去查看kafka的消息数据,目录是{logdir}/test-0/00000000000000.lo ...
- Python3 Tkinter-Scale
1.创建 from tkinter import * root=Tk() Scale(root).pack() root.mainloop() 2.参数 from tkinter import * r ...
- HDU 4308 Saving Princess claire_(简单BFS)
求出不使用P点时起点到终点的最短距离,求出起点到所有P点的最短距离,求出终点到所有P点的最短距离. 答案=min( 不使用P点时起点到终点的最短距离, 起点到P的最短距离+终点到P的最短距离 ) #i ...
- 【IdentityServer4文档】- 支持和咨询选项
支持和咨询选项 我们为 IdentityServer 提供多个免费和商业支持及咨询选项. 免费支持 免费支持是基于社区的,而且使用的是公共论坛 StackOverflow 有越来越多的使用 Ident ...
- android BadgeView的使用(图片上的文字提醒)
BadgeView主要是继承了TextView,所以实际上就是一个TextView,底层放了一个label,可以自定义背景图,自定义背景颜色,是否显示,显示进入的动画效果以及显示的位置等等: 这是Gi ...
- 《剑指offer》---字符串的全排列
本文算法使用python3实现 1.问题一 1.1 题目描述: 输入一个字符串,按字典序打印出该字符串中字符的所有排列.例如输入字符串abc,则打印出由字符a,b,c所能排列出来的所有字符串abc ...