4869: [Shoi2017]相逢是问候

Time Limit: 40 Sec  Memory Limit: 512 MB
Submit: 1205  Solved: 409
[Submit][Status][Discuss]

Description

Informatikverbindetdichundmich.
信息将你我连结。B君希望以维护一个长度为n的数组,这个数组的下标为从1到n的正整数。一共有m个操作,可以
分为两种:0 l r表示将第l个到第r个数(al,al+1,...,ar)中的每一个数ai替换为c^ai,即c的ai次方,其中c是
输入的一个常数,也就是执行赋值ai=c^ai1 l r求第l个到第r个数的和,也就是输出:sigma(ai),l<=i<=rai因为
这个结果可能会很大,所以你只需要输出结果mod p的值即可。
 

Input

第一行有三个整数n,m,p,c,所有整数含义见问题描述。
接下来一行n个整数,表示a数组的初始值。
接下来m行,每行三个整数,其中第一个整数表示了操作的类型。
如果是0的话,表示这是一个修改操作,操作的参数为l,r。
如果是1的话,表示这是一个询问操作,操作的参数为l,r。
1 ≤ n ≤ 50000, 1 ≤ m ≤ 50000, 1 ≤ p ≤ 100000000, 0 < c <p, 0 ≤ ai < p
 

Output

对于每个询问操作,输出一行,包括一个整数表示答案mod p的值。
 

Sample Input

4 4 7 2
1 2 3 4
0 1 4
1 2 4
0 1 4
1 1 3

Sample Output

0
3

HINT

鸣谢多名网友提供正确数据,已重测!

Source

 
大家肯定都知道欧拉降幂公式了,就是x^y ≡ x^(y%p + [y>=p]*p )    mod p
这个厉害就厉害在可以在x和p不互质的情况下使用。
 
而且对于本题来说,我们可以发现一个位置最多被修改 log(P) 次就会永远成为一个定值。
所以我们暴力修改动态维护一下前缀和就好了。
 
注意,我们迭代的时候
 
#include<bits/stdc++.h>
#define ll long long
#define maxn 50005
using namespace std;
int P[30],n,m,a[maxn];
int f[maxn],len,C,opt,TP;
int le,ri,tim[maxn],val[maxn];
set<int> s;
//set<int> ::iterator it; inline int add(int x,int y,const int ha){
x+=y;
return x>=ha?x-ha:x;
} inline int mod(int x,int y,const int ha){
ll P=x*(ll)y;
return P>=ha?P%ha+ha:P;
} inline int phi(int x){
int y=1;
for(int i=2;i*(ll)i<=x;i++) if(!(x%i)){
y*=(i-1),x/=i;
while(!(x%i)) x/=i,y*=i;
if(x==1) break;
} if(x!=1) return y*(x-1);
else return y;
} inline void update(int x,int y){
for(;x<=n;x+=x&-x) f[x]=add(f[x],y,P[0]);
} inline int query(int x){
int an=0;
for(;x;x-=x&-x) an=add(an,f[x],P[0]);
return an;
} inline int ksm(int x,int y,const int ha){
int an=1;
for(;y;y>>=1,x=mod(x,x,ha)) if(y&1) an=mod(an,x,ha);
return an;
} /*
int get(int sit,int pos){
if(sit==TP) return a[pos]>=P[sit]?a[pos]%P[sit]+P[sit]:a[pos];
else return ksm(C,get(sit+1,pos),P[sit]);
}
*/ inline void TOL(){
int now=*s.lower_bound(le);
for(;now<=ri;now=*s.upper_bound(now)){
update(now,P[0]-val[now]);
tim[now]++,TP=tim[now]; if(tim[now]>len) val[now]=1;
else val[now]=(a[now]>=P[TP]?a[now]%P[TP]+P[TP]:a[now]);
for(int i=TP-1;i>=0;i--){
val[now]=ksm(C,val[now],P[i]);
}
val[now]%=P[0]; update(now,val[now]); if(tim[now]==len+1) s.erase(now);
}
} int main(){
// freopen("9.in","r",stdin);
// freopen("9.out","w",stdout); scanf("%d%d%d%d",&n,&m,P,&C);
for(len=1;;len++){
P[len]=phi(P[len-1]);
if(P[len]==1) break;
}
s.insert(0),s.insert(n+1);
for(int i=1;i<=n;i++){
scanf("%d",a+i);
update(i,a[i]);
s.insert(i);
val[i]=a[i];
} while(m--){
scanf("%d%d%d",&opt,&le,&ri);
if(opt) printf("%d\n",add(query(ri),P[0]-query(le-1),P[0]));
else TOL();
} return 0;
}

  

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

  1. 2017 [六省联考] T5 分手是祝愿

    4872: [Shoi2017]分手是祝愿 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 458  Solved: 299[Submit][Statu ...

  2. 2017 [六省联考] T1 期末考试

    4868: [Shoi2017]期末考试 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 842  Solved: 385[Submit][Status ...

  3. 2017 [六省联考] T6 寿司餐厅

    4873: [Shoi2017]寿司餐厅 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 450  Solved: 316[Submit][Status ...

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

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

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

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

  6. 六省联考2017 Day1

    目录 2018.3.18 Test T1 BZOJ.4868.[六省联考2017]期末考试 T2 T3 BZOJ.4870.[六省联考2017]组合数问题(DP 矩阵快速幂) 总结 考试代码 T1 T ...

  7. 六省联考2017 Day2

    目录 2018.3.27 Test 总结 T1 T2 T3 BZOJ.4873.[六省联考2017]寿司餐厅(最小割ISAP 最大权闭合子图) 考试代码 T1 T2 T3 2018.3.27 Test ...

  8. [luogu] P3745 [六省联考2017]期末考试 (贪心)

    P3745 [六省联考2017]期末考试 题目描述 有 \(n\) 位同学,每位同学都参加了全部的 \(m\) 门课程的期末考试,都在焦急的等待成绩的公布. 第 \(i\) 位同学希望在第 \(t_i ...

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

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

随机推荐

  1. [译]The Python Tutorial#5. Data Structures

    [译]The Python Tutorial#Data Structures 5.1 Data Structures 本章节详细介绍之前介绍过的一些内容,并且也会介绍一些新的内容. 5.1 More ...

  2. leetcode-8-pointer

    void deleteNode(ListNode* node) { *node = *node->next; }

  3. Aizu - 1386 Starting a Scenic Railroad Service (思维乱搞)

    给你n个区间,求: 1:最多有多少区间与同一个区间相交. 2:相交部分的最大区间数目. Sample Input 1 4 1 3 1 3 3 6 3 6 Sample Output 1 2 2 Sam ...

  4. Linux系统监视工具

    转自      http://bbs.51cto.com/thread-971896-1.html # 1: top – 查看活动进程的命令TOP工具能够实时显示系统中各个进程的资源占用状况.默认情况 ...

  5. 新线程 handler

    class CalculateThread extends Thread { private Handler handler; @Override public void run() { super. ...

  6. linux随笔三

    1.ps   结果输出: PID TTY TIME CMD pts/ :: bash pts/ :: ps显示了程序的进程ID,其运行的终端和进程使用的cpu时间

  7. session的工作原理、django的超时时间设置及session过期判断

    1.session原理 cookie是保存在用户浏览器端的键值对 session是保存在服务器端的键值对 session服务端中存在的数据为: session = { 随机字符串1:{ 用户1的相关信 ...

  8. 聊聊、Nginx 初始化错误信息

    这篇文章我们继续学习 main 方法,我们先来看看 ngx_debug_init() 这个方法. 从方法名我们也知道,debug初始化.我们先看看方法位置在哪.我们来断点在这个方法上面. Functi ...

  9. Spring配置文件中使用ref local与ref bean的区别

    Spring配置文件中使用ref local与ref bean的区别.在ApplicationResources.properties文件中,使用<ref bean>与<ref lo ...

  10. poj1985&&第四次CCF软件认证第4题 求树的直径

    Cow Marathon Time Limit: 2000MS   Memory Limit: 30000K Total Submissions: 4216   Accepted: 2137 Case ...