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

这道题好劲啊 乍一看就知道是线段树 然而不知道怎么做
于是又看题解 愉快地引用结论 （数据里的p比较特殊？？？

首先平方之后取模 这种东西显然是有环的 然后（根据结论…..）其实最大环只有60。。
所以我们只需要预处理出每个数是否在环上 以及他的next是什么
接着每次平方操作就相当于在环上走一步 （当然你要算出环上每个位置的和）
不在环上的呢？ 暴力走。。同样的 他到环的距离不会很长 （有人说<11）

故时间复杂度为 O(60∗nlogn)
注意 我说的环不只是对单个数自己平方形成的 而是对一个数列一起平方
实在不懂的话就看看代码里的upd和pushdown吧 不难懂。。

#include<bits/stdc++.h>using namespace std;const int N=100010,M=10000;inline int read(){    char ch=getchar(); int x=0,f=1;    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();}    return x*f;}inline int lcm(int a,int b){return a/__gcd(a,b)*b;}int p,s[N<<2],w[N<<2],a[N],len[N<<2],cir[N<<2][62];bool inc[M],in[N<<2]; int nx[M],lz[N<<2];void init(int x,int v){    s[x]=v; in[x]=inc[v];    if(in[x]){        len[x]=1,cir[x][w[x]=0]=v;        for(int u=nx[v];u!=v;u=nx[u]) cir[x][len[x]++]=u;    }}void upd(int x){    int lc=x<<1,rc=lc|1;    s[x]=s[lc]+s[rc],in[x]=in[lc]&in[rc];    if(in[x]){        len[x]=lcm(len[lc],len[rc]);        int tx=w[lc],ty=w[rc];        for(int i=0;i<len[x];++i){            cir[x][i]=cir[lc][tx++]+cir[rc][ty++];            if(tx==len[lc])tx=0; if(ty==len[rc])ty=0;        }        w[x]=0;    }}void bulid(int x,int l,int r){    if(l==r){        init(x,a[l]); return;    }    int lc=x<<1,rc=lc|1,mid=(l+r)>>1;    bulid(lc,l,mid),bulid(rc,mid+1,r);    upd(x);}void pushdown(int x){    int c=lz[x],lc=x<<1,rc=lc|1; lz[x]=0;    lz[lc]+=c,lz[rc]+=c; w[lc]+=c,w[rc]+=c;    if(w[lc]>=len[lc])w[lc]%=len[lc]; s[lc]=cir[lc][w[lc]];    if(w[rc]>=len[rc])w[rc]%=len[rc]; s[rc]=cir[rc][w[rc]];}void change(int x,int l,int r,int ql,int qr){    if(l==ql&&r==qr&&in[x]){        ++lz[x],++w[x]; if(w[x]==len[x])w[x]=0;        s[x]=cir[x][w[x]]; return;    }    if(l==r){init(x,nx[s[x]]); return;}    if(lz[x])pushdown(x); int mid=(l+r)>>1,lc=x<<1,rc=lc|1;    if(qr<=mid)change(lc,l,mid,ql,qr);    else if(ql>mid)change(rc,mid+1,r,ql,qr);    else change(lc,l,mid,ql,mid),change(rc,mid+1,r,mid+1,qr);    upd(x);}int query(int x,int l,int r,int ql,int qr){    if(l==ql&&r==qr)return s[x];    if(lz[x])pushdown(x); int mid=(l+r)>>1,lc=x<<1,rc=lc|1;    if(qr<=mid)return query(lc,l,mid,ql,qr);    else if(ql>mid)return query(rc,mid+1,r,ql,qr);    return query(lc,l,mid,ql,mid)+query(rc,mid+1,r,mid+1,qr);}bool v[M];int main(){    int n=read(),m=read(),i; p=read();    for(i=1;i<=n;++i)a[i]=read();    for(i=0;i<p;++i)nx[i]=(i*i)%p,inc[i]=1;    for(i=0;i<p;++i) if(!v[i]){        int y=i; while(!v[y])v[y]=1,y=nx[y];        for(int t=i;t!=y;t=nx[t])inc[t]=0;    }    bulid(1,1,n);    while(m--){        int u=read(),l=read(),r=read();        if(!u) change(1,1,n,l,r);        else printf("%d\n",query(1,1,n,l,r));    }    return 0;}