bzoj4869

http://www.lydsy.com/JudgeOnline/problem.php?id=4869

终于A了。。。参考了下dalao的代码。。。

拓展欧几里得定理，改了几次就不变了，但是用的时候要在快速幂里判是不是要用。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 100010;
int n, m, cnt;
ll p, c;
ll phi[N], table[N];
namespace seg // n^x = n^(x % phi[x] + phi[x])
{struct data {ll ans, mn;} tree[N << 2];inline ll getphi(ll x){ll ret = x, lim = x;for(ll i = 2; i * i <= lim; ++i) if(x % i == 0) {ret = ret * (i - 1) / i;while(x % i == 0) x /= i;}
//      printf("ret=%d\n", ret);if(x > 1) ret = ret * (x - 1) / x;return ret;}inline ll power(ll x, ll t, ll p, bool &flag){bool big = false;ll ret = 1; for(; t; t >>= 1) {if(t & 1) {ret = ret * x ;flag |= big | (ret >= p);        ret %= p;}x = x * x; if(x >= p) big = true, x %= p; }return ret; }ll calc(ll x, int t){if(x >= phi[t]) x = x % phi[t] + phi[t];for(int i = t - 1; i >= 0; --i) {bool flag = false;x = power(c, x, phi[i], flag);if(flag) x += phi[i];}return x % phi[0];}inline void build(int l, int r, int x){if(l == r) { tree[x].ans = table[l]; return; }int mid = (l + r) >> 1;build(l, mid, x << 1); build(mid + 1, r, x << 1 | 1);tree[x].ans = (tree[x << 1].ans + tree[x << 1 | 1].ans) % phi[0];} inline void update(int l, int r, int x, int a, int b){ //如果这次的幂和上次一样就不变了 if(tree[x].mn >= cnt) return;if(l > b || r < a) return;if(l == r){++tree[x].mn;tree[x].ans = calc(table[l], tree[x].mn);return;}int mid = (l + r) >> 1;update(l, mid, x << 1, a, b); update(mid + 1, r, x << 1 | 1, a, b);tree[x].mn = min(tree[x << 1].mn, tree[x << 1 | 1].mn);tree[x].ans = (tree[x << 1].ans + tree[x << 1 | 1].ans) % phi[0]; }inline ll query(int l, int r, int x, int a, int b){if(l > b || r < a) return 0;if(l >= a && r <= b) return tree[x].ans % phi[0];int mid = (l + r) >> 1, ret = 0;ret = (ret + query(l, mid, x << 1, a, b)) % phi[0];ret = (ret + query(mid + 1, r, x << 1 | 1, a, b)) % phi[0];return ret;    } 
} using namespace seg;
int main()
{scanf("%d%d%lld%lld", &n, &m, &p, &c);phi[0] = p;ll P = p;while(P != 1) phi[++cnt] = P = getphi(P);phi[++cnt] = 1; for(int i = 1; i <= n; ++i) scanf("%lld", &table[i]);build(1, n, 1); while(m--){int opt, l, r; scanf("%d", &opt);if(opt == 0){scanf("%d%d", &l, &r); update(1, n, 1, l, r);}if(opt == 1) {scanf("%d%d", &l, &r);printf("%lld\n", query(1, n, 1, l, r));}}return 0;
}