hdu_5286_wyh2000 and sequence(分块)

题目链接：hdu_5286_wyh2000 and sequence

题意：

给一段长度为N的序列，每次询问l-r（l和r和上一次询问的答案有关）内 不同的数的 出现次数的次方 的和。强制在线

题解：

这里贴个达哥的题解：

大体思路就是，把n个数分成sqrt（n）块，每块sqrt(n)个数，然后求出任意两块i到j（包含第i块和第j块）的ans，对于每次询问l和r，找到刚好包含l和r的一段连续的块，因为我们已经知道了这段连续块的答案，所以只要再去掉两段多出来的一部分，就好了。

为了得到任意两块间的答案G[i][j],我们可以枚举每个块的起点，然后for到最后，暴力一边，就好了，

对于多出来的一部分，我们只要知道这一部分里的数，在包含它的连续块中已经出现了几次，然后就可以很容易地更新答案了。于是我们可以先预处理出一个后缀和cnt[i][j]:表示从第i块开始，第j小的数出现了多少次，然后，对于每次询问q，我们 只要 只要 只要 （后缀和减一下）找到这一部分里的数在连续的块中出现了次数就好了，别的数就算了，都找的话，肯定超时。

之前还要预处理各种数次方的结果，还要排序，去重。

cnt[i][j]:从第i块开始到最后第j小的数出现的次数

pos[i]:在A中下标为i的数是第几小

 #include<bits/stdc++.h>
 #define F(i,a,b) for(int i=a;i<=b;++i)
 using namespace std;
 typedef long long ll;

 const int N=5e4+,P=1e9+;

 int sqr,n,m,l,r,len,t,pos[N],cnt[][N];
 ll A[N],a[N],num[N],vis[N];
 vector<ll>v[N];//第i小的数的j次方
 ll G[][];//块i到j的答案

 int T_T(int l,int r)
 {
     int L=l/sqr,R=r/sqr;
     ll ans=G[L][R];
     int LL=L*sqr,RR=min(n-,(R+)*sqr-);
     F(i,LL,l-)num[pos[i]]=cnt[L][pos[i]]-cnt[R+][pos[i]];
     F(i,r+,RR)num[pos[i]]=cnt[L][pos[i]]-cnt[R+][pos[i]];
     while(LL<l)
     {
         ans-=v[pos[LL]][num[pos[LL]]];
         ans+=v[pos[LL]][--num[pos[LL]]];
         LL++;
     }
     while(RR>r)
     {
         ans-=v[pos[RR]][num[pos[RR]]];
         ans+=v[pos[RR]][--num[pos[RR]]];
         RR--;
     }
     return (ans%P+P)%P;
 }

 int main()
 {
     scanf("%d",&t);
     while(t--)
     {
         scanf("%d%d",&n,&m);
         F(i,,n-)scanf("%lld",a+i),A[i]=a[i];
         sort(a,a+n),len=unique(a,a+n)-a;
         F(i,,len-)v[i].clear(),v[i].push_back(),vis[i]=;
         F(i,,n-)vis[pos[i]=lower_bound(a,a+len,A[i])-a]++;
         F(i,,len-)
         {
             ll p=;
             F(j,,vis[i])v[i].push_back(p=p*a[i]%P);
         }
         sqr=sqrt(n+0.5);
         memset(cnt,,sizeof(cnt));
         for(int i=;i*sqr<n;i++)
         {
             ll ans=;
             memset(num,,sizeof(num));
             for(int j=i*sqr;j<n;j++)
             {
                 cnt[i][pos[j]]++;
                 ans-=v[pos[j]][num[pos[j]]];
                 ans+=v[pos[j]][++num[pos[j]]];
                 if((j+)%sqr==||j==n-)G[i][j/sqr]=(ans%P+P)%P;
             }
         }
         int la=;
         while(m--)
         {
             int L,R;
             scanf("%d%d",&L,&R);
             l=(L^la)%n,r=(R^la)%n;
             if(l>r)l^=r,r^=l,l^=r;
             printf("%d\n",la=T_T(l,r));
         }
     }
     return ;
 }