BZOJ4869 [Shoi2017]相逢是问候 【扩展欧拉定理 + 线段树】
题目链接
题解
这题调得我怀疑人生，，结果就是因为某些地方\(sb\)地忘了取模
前置题目：BZOJ3884
扩展欧拉定理：
\[c^a \equiv c^{a \mod \varphi(p) + [a \ge p]p} \pmod p
\]
我们发现当我们进行\(0\)操作，就相当于在\(a\)底部添加一层\(c\)
当我们进行得足够多的时候，\(\varphi(p)\)就会取到\(1\)，从而不再变化
所以每个位置操作次数其实是有限的，为\(O(logp)\)次
为何是\(O(logp)\)次呢？
考虑欧拉函数：
\[\varphi(n) = n \prod\limits_{i = 1}^{k} \frac{p_i - 1}{p_i}
\]
由于质数一定是奇数，所以\(p_i - 1\)一定为偶数，所以操作一次后的\(p\)一定为偶数
偶数有\(2\)这个质因子，式子中就会存在\(\frac{1}{2}\)这一项，所以至少减少\(\frac{1}{2}\)
所以到达\(1\)是\(O(logp)\)的
用线段树进行维护
每次修改重新计算\(O(log^2p)\)，总复杂度\(O(nlognlog^3p)\)，凭信仰可过
由于每次快速幂底都是\(c\)，我们可以预处理\(c^x\)，从而做到\(O(nlognlog^2p)\)
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<map>
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define mp(a,b) make_pair<int,int>(a,b)
#define cls(s) memset(s,0,sizeof(s))
#define cp pair<int,int>
#define LL long long int
#define ls (u << 1)
#define rs (u << 1 | 1)
using namespace std;
const int maxn = 100005,maxm = 100005,INF = 1000000000;
inline int read(){
int out = 0,flag = 1; char c = getchar();
while (c < 48 || c > 57){if (c == '-') flag = -1; c = getchar();}
while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();}
return out * flag;
}
int p[maxn],pi,isn[maxn];
void init(){
for (int i = 2; i <= 10000; i++){
if (!isn[i]) p[++pi] = i;
for (int j = 1; j <= pi && i * p[j] <= 10000; j++){
isn[i * p[j]] = true;
if (i % p[j] == 0) break;
}
}
}
int n,m,P[30],Pi,c,a[maxn],cnt[maxn],low[30];
LL C[30][maxn],CC[30][maxn];
int phi(int x){
int tmp = x,ans = x;
for (int i = 1; i <= pi && p[i] <= tmp; i++){
int v = p[i];
if (tmp % v == 0){
ans = ans / v * (v - 1);
while (tmp % v == 0) tmp /= v;
}
}
if (tmp - 1) ans = ans / tmp * (tmp - 1);
return ans;
}
LL qpow(int b,int md){
return CC[md][b / 10000] * C[md][b % 10000] % P[md];
}
int cal(int x,int t){
LL last,tmp = x;
if (tmp > P[t]) tmp = tmp % P[t] + P[t];
for(int i = t; i; i--)
{
last = tmp;
tmp = qpow(tmp,i - 1);
if (last >= low[i - 1]){
tmp += P[i - 1];
}
}
return tmp;
}
int sum[maxn << 2],val[maxn << 2];
void upd(int u){
sum[u] = (sum[ls] + sum[rs]) % P[0];
val[u] = val[ls] | val[rs];
}
void build(int u,int l,int r){
val[u] = 1;
if (l == r){sum[u] = a[l] % P[0]; return;}
int mid = l + r >> 1;
build(ls,l,mid);
build(rs,mid + 1,r);
upd(u);
}
void change(int u,int l,int r){
if (!val[u]) return;
if (l == r){
cnt[l]++;
if (cnt[l] >= Pi) val[u] = 0;
sum[u] = cal(a[l],cnt[l]) % P[0];
return;
}
int mid = l + r >> 1;
change(ls,l,mid);
change(rs,mid + 1,r);
upd(u);
}
void modify(int u,int l,int r,int L,int R){
if (!val[u]) return;
if (l >= L && r <= R){
change(u,l,r);
return;
}
int mid = l + r >> 1;
if (mid >= L) modify(ls,l,mid,L,R);
if (mid < R) modify(rs,mid + 1,r,L,R);
upd(u);
}
int query(int u,int l,int r,int L,int R){
if (l >= L && r <= R) return sum[u] % P[0];
int mid = l + r >> 1;
if (mid >= R) return query(ls,l,mid,L,R);
if (mid < L) return query(rs,mid + 1,r,L,R);
return (query(ls,l,mid,L,R) + query(rs,mid + 1,r,L,R)) % P[0];
}
int main(){
init();
n = read(); m = read(); P[0] = read(); c = read();
REP(i,n) a[i] = read();
while (P[Pi] != 1){
++Pi;
P[Pi] = phi(P[Pi - 1]);
}
P[++Pi] = 1;
for (register int t = 0; t <= Pi; t++){
C[t][0] = CC[t][0] = 1 % P[t];
low[t] = 0;
for (register int i = 1; i <= 10000; i++){
C[t][i] = C[t][i - 1] * c;
if (C[t][i] >= P[t] && !low[t]) low[t] = i;
C[t][i] %= P[t];
}
CC[t][1] = C[t][10000];
for (register int i = 2; i <= 10000; i++)
CC[t][i] = CC[t][i - 1] * C[t][10000] % P[t];
}
build(1,1,n);
int opt,l,r;
while (m--){
opt = read(); l = read(); r = read();
if (!opt) modify(1,1,n,l,r);
else printf("%d\n",query(1,1,n,l,r));
}
return 0;
}
