bzoj:4105: [Thu Summer Camp 2015]平方运算

Description

 

Input

第一行有三个整数N，M，p，分别代表序列的长度、平方操作与询问操作的总次数以及在平方操作中所要模的数。

 

接下来一行N个数代表一开始的序列{X1,X2,...,XN}。

 

接下来M行，每行三个整数op,l,r。其中op代表本次操作的类型。若op=0，代表这是一次平方操作，平方的区间为[l,r]；如果op=1，代表这是一次询问操作，询问的区间为[l,r]。

 

Output

对于每次的询问操作，输出一行代表这段区间内数的总和。注意：答案没有对任何数取模。

 

Sample Input

3 3 11
1 2 3
1 1 3
0 1 3
1 1 3

Sample Output

6
14

HINT

对于100%的数据，∀i,Xi∈[0,p),l,r∈[1,n]

N，M，p的范围如下：

 

编号  N  M  p

1  1000  1000  233

2  1000  1000  2332

3  100000  100000  5

4  100000  100000  8192

5  100000  100000  23

6  100000  100000  45

7  100000  100000  37

8  55000  55000  4185

9  55000  55000  5850

10  55000  55000  2975

11  55000  55000  2542

12  55000  55000  2015

13  60000  60000  2003

14  65000  65000  2010

15  70000  70000  4593

16  75000  75000  4562

17  80000  80000  1034

18  85000  85000  5831

19  90000  90000  9905

20  100000  100000  9977

 

膜了一大把题解

线段树是肯定的……

平方是会出现循环节的（听说会很短

预处理出所有环，和指向环的链，由于没有其他修改操作，这里面的数字肯定是越修改越往环跑，进了环就处理出跑k次后答案是多少，不在环上就暴力改……

论权限号的重要性

/**************************************************************
    Problem: 4105
    User: JSZX11556
    Language: C++
    Result: Accepted
    Time:25216 ms
    Memory:275236 kb
****************************************************************/

#include<cstdio>
#include<algorithm>
#define lp (p<<1)
#define rp ((p<<1)|1)
using namespace std;

int read_p,read_ca;
inline int read(){
    read_p=;read_ca=getchar();
    while(read_ca<''||read_ca>'') read_ca=getchar();
    while(read_ca>=''&&read_ca<='') read_p=read_p*+read_ca-,read_ca=getchar();
    return read_p;
}
int n,m,p,a[],ne[],i;
bool v[],f[];
struct na{
    int l,r,w,le,c,b[],pos;
    bool v;
}t[];
inline int gcd(int x,int y){return y==?x:gcd(y,x%y);}
inline int lcm(int x,int y){return x*y/gcd(x,y);}
inline void updata(int p){
    if (t[p].l==t[p].r) return;
    t[p].w=t[lp].w+t[rp].w;
    t[p].v=t[lp].v&t[rp].v;
    if (t[p].v){
        t[p].pos=;
        t[p].le=lcm(t[lp].le,t[rp].le);
        for (i=;i<t[p].le;i++) t[p].b[i]=t[lp].b[(i+t[lp].pos)%t[lp].le]+t[rp].b[(i+t[rp].pos)%t[rp].le];
    }
}
inline void build(int p,int l,int r){
    t[p].l=l;t[p].r=r;
    if (l==r){
        t[p].w=a[l];t[p].v=f[t[p].w];
        if (t[p].v) for (t[p].b[]=t[p].w,t[p].le=,i=ne[t[p].w];i!=t[p].w;i=ne[i]) t[p].b[t[p].le++]=i;
        return;
    }
    int mid=l+r>>;
    build(lp,l,mid);build(rp,mid+,r);
    updata(p);
}
inline void hb(int p,int c){
    if (!t[p].v){
        for (i=;i<=c;i++) t[p].w=ne[t[p].w];t[p].v=f[t[p].w];
        if (t[p].v) for (t[p].b[]=t[p].w,t[p].le=,i=ne[t[p].w];i!=t[p].w;i=ne[i]) t[p].b[t[p].le++]=i;
        return;
    }
    c%=t[p].le;
    t[p].c+=c;if (t[p].c>=t[p].le) t[p].c-=t[p].le;
    t[p].pos+=c;if (t[p].pos>=t[p].le) t[p].pos-=t[p].le;
    t[p].w=t[p].b[t[p].pos];
}
inline void pd(int p){
    if (t[p].c){
        if (t[p].l!=t[p].r) hb(lp,t[p].c),hb(rp,t[p].c);
        t[p].c=;
        updata(p);
    }
}
inline void ch(int p,int l,int r){
    if (t[p].l==l&&t[p].r==r&&t[p].v) return hb(p,);
    pd(p);
    if (t[p].l==t[p].r){
        t[p].w=ne[t[p].w];t[p].v=f[t[p].w];
        if (t[p].v) for (t[p].b[]=t[p].w,t[p].le=,i=ne[t[p].w];i!=t[p].w;i=ne[i]) t[p].b[t[p].le++]=i;
        return;
    }
    int mid=t[p].l+t[p].r>>;
    if (r<=mid) ch(lp,l,r);else
    if (l>mid) ch(rp,l,r);else
    ch(lp,l,mid),ch(rp,mid+,r);
    updata(p);
}
inline int MMH(int p,int l,int r){
    pd(p);
    if (t[p].l==l&&t[p].r==r) return t[p].w;
    int mid=t[p].l+t[p].r>>;
    if (r<=mid) return MMH(lp,l,r);else
    if (l>mid) return MMH(rp,l,r);else
    return MMH(lp,l,mid)+MMH(rp,mid+,r);
}
int main(){
    register int i,j,T,l,r;
    n=read();m=read();
    p=read();
    for (i=;i<=n;i++) a[i]=read();
    for (i=;i<p;i++) ne[i]=i*i%p,f[i]=;
    for (i=;i<p;i++)
    if (!v[i]){
        for (j=i;!v[j];j=ne[j]) v[j]=;
        for (T=i;T!=j;T=ne[T]) f[T]=;
    }
    build(,,n);
    for (i=;i<=m;i++){
        T=read();l=read();r=read();
        if (T)printf("%d\n",MMH(,l,r));else ch(,l,r);
    }
}