题目链接: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的数是第几小
1 #include<bits/stdc++.h> 2 #define F(i,a,b) for(int i=a;i<=b;++i) 3 using namespace std; 4 typedef long long ll; 5 6 const int N=5e4+7,P=1e9+7; 7 8 int sqr,n,m,l,r,len,t,pos[N],cnt[240][N]; 9 ll A[N],a[N],num[N],vis[N]; 10 vector<ll>v[N];//第i小的数的j次方 11 ll G[240][240];//块i到j的答案 12 13 int T_T(int l,int r) 14 { 15 int L=l/sqr,R=r/sqr; 16 ll ans=G[L][R]; 17 int LL=L*sqr,RR=min(n-1,(R+1)*sqr-1); 18 F(i,LL,l-1)num[pos[i]]=cnt[L][pos[i]]-cnt[R+1][pos[i]]; 19 F(i,r+1,RR)num[pos[i]]=cnt[L][pos[i]]-cnt[R+1][pos[i]]; 20 while(LL<l) 21 { 22 ans-=v[pos[LL]][num[pos[LL]]]; 23 ans+=v[pos[LL]][--num[pos[LL]]]; 24 LL++; 25 } 26 while(RR>r) 27 { 28 ans-=v[pos[RR]][num[pos[RR]]]; 29 ans+=v[pos[RR]][--num[pos[RR]]]; 30 RR--; 31 } 32 return (ans%P+P)%P; 33 } 34 35 int main() 36 { 37 scanf("%d",&t); 38 while(t--) 39 { 40 scanf("%d%d",&n,&m); 41 F(i,0,n-1)scanf("%lld",a+i),A[i]=a[i]; 42 sort(a,a+n),len=unique(a,a+n)-a; 43 F(i,0,len-1)v[i].clear(),v[i].push_back(0),vis[i]=0; 44 F(i,0,n-1)vis[pos[i]=lower_bound(a,a+len,A[i])-a]++; 45 F(i,0,len-1) 46 { 47 ll p=1; 48 F(j,1,vis[i])v[i].push_back(p=p*a[i]%P); 49 } 50 sqr=sqrt(n+0.5); 51 memset(cnt,0,sizeof(cnt)); 52 for(int i=0;i*sqr<n;i++) 53 { 54 ll ans=0; 55 memset(num,0,sizeof(num)); 56 for(int j=i*sqr;j<n;j++) 57 { 58 cnt[i][pos[j]]++; 59 ans-=v[pos[j]][num[pos[j]]]; 60 ans+=v[pos[j]][++num[pos[j]]]; 61 if((j+1)%sqr==0||j==n-1)G[i][j/sqr]=(ans%P+P)%P; 62 } 63 } 64 int la=0; 65 while(m--) 66 { 67 int L,R; 68 scanf("%d%d",&L,&R); 69 l=(L^la)%n,r=(R^la)%n; 70 if(l>r)l^=r,r^=l,l^=r; 71 printf("%d ",la=T_T(l,r)); 72 } 73 } 74 return 0; 75 }