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. leetcode-14-basic-breadthFirstSearch

    BFS: breadth first search 107. Binary Tree Level Order Traversal II 解题思路: 本来我是用map<int,int>存所有 ...

  2. __vet_atags

    参考:atags--__vet_atags标签    arch/arm/include/asm/setup.h /* * linux/include/asm/setup.h * * Copyright ...

  3. 笔记-python-urllib

    笔记-python-urllib 1.      简介 PYTHON3中将urllib,urllib2整合到URLLIB中 包括以下模块 urllib.request 请求模块(核心) urllib. ...

  4. 同一条sql在mysql5.6和5.7版本遇到的问题。

    之前用的是mysql 5.6版本,执行select * from table group by colunm 是可以出结果的, 但是切换的5.7版本,这条sql就报错, Expression #1 o ...

  5. App架构经验总结

    作者:李纪钢,网名 Keegan小钢,博客地址:http://keeganlee.me.目前在广州日报新媒体有限公司,负责移动产品的研发工作. 关于:本文整理自CSDN架构主题月子活动金牌架构师微课堂 ...

  6. BNU 3692 I18n 模拟

    题意: 在一篇文章中,单词可以缩写.例如单词Internationalization可以缩写为I18n,缩写的两端是原单词的首尾字母,中间的数字是被省略的字母的个数. 现在给你一篇缩写的文章,输出展开 ...

  7. 洛谷P1085不高兴的津津

  8. Leetcode 435.无重叠区间

    无重叠区间 给定一个区间的集合,找到需要移除区间的最小数量,使剩余区间互不重叠. 注意: 可以认为区间的终点总是大于它的起点. 区间 [1,2] 和 [2,3] 的边界相互"接触" ...

  9. BZOJ2326 [HNOI2011]数学作业 【矩阵快速幂】

    题解 我们设f[i]表示前i个数模M意义下的答案 则f[i] = f[i - 1] * 100...0 + i[i是几位就有几个0] 可以写出矩阵递推式: 之后按位数分组矩乘就好了 #include& ...

  10. JS实现并集,交集和差集

    var set1 = new Set([1,2,3]);var set2 = new Set([2,3,4]); 并集let union = new Set([...set1, ...set2]); ...