bzoj3744 Gty的妹子序列
Posted
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj3744 Gty的妹子序列相关的知识,希望对你有一定的参考价值。
Description
给定一个正整数序列a,对于每次询问,输出al...ar中的逆序对数,强制在线。
Input
第一行包括一个整数n(1<=n<=50000),表示数列a中的元素数。
第二行包括n个整数a1...an(ai>0,保证ai在int内)。
接下来一行包括一个整数m(1<=m<=50000),表示询问的个数。
接下来m行,每行包括2个整数l、r(1<=l<=r<=n),表示询问al...ar中的逆序
对数(若ai>aj且i<j,则为一个逆序对)。
l,r要分别异或上一次询问的答案(lastans),最开始时lastans=0。
保证涉及的所有数在int内。
Output
对每个询问,单独输出一行,表示al...ar中的逆序对数。
分为sqrt(n)块,预处理每两块间和块内的逆序对数,以及每个位置的数到每块的逆序对数,取前缀和。
对询问合并区间内块内和块间答案,临时计算两端的答案。
#include<cstdio> #include<cmath> #include<algorithm> using namespace std; int n,m,l,r,la; int v[50005]; int v1[50005],v2[50005],v3[50005]; int ls[230],rs[230],bp=0; int bs[230]; int bb[230][230]; int pb[230][50001]; int sz; int cal(int l,int r){ int m=l+r>>1,ans=0; if(l<m)ans+=cal(l,m); if(m+1<r)ans+=cal(m+1,r); int p1=l,p2=m+1,p3=l; while(p1<=m&&p2<=r){ if(v2[p1]<=v2[p2])v3[p3++]=v2[p1++]; else ans+=m+1-p1,v3[p3++]=v2[p2++]; } while(p1<=m)v3[p3++]=v2[p1++]; while(p2<=r)v3[p3++]=v2[p2++]; for(int i=l;i<=r;i++)v2[i]=v3[i]; return ans; } int main(){ scanf("%d",&n); for(int i=1;i<=n;i++){ scanf("%d",v+i); v1[i]=v[i]; } sz=sqrt(n); for(int i=1;i<=n;i+=sz){ ls[bp]=i; rs[bp]=i+sz-1; if(rs[bp]>n)rs[bp]=n; for(int i=ls[bp];i<=rs[bp];i++)v2[i]=v[i]; bs[bp]=cal(ls[bp],rs[bp]); for(int i=ls[bp];i<=rs[bp];i++)v1[i]=v2[i]; bp++; } for(int i=0;i<bp;i++){ for(int j=i+1;j<bp;j++){ bb[i][j+1]+=bb[i][j]; for(int a=rs[i],p=rs[j]+1;a>=ls[i];a--){ while(p>ls[j]&&v1[p-1]>=v1[a])--p; bb[i][j+1]+=p-ls[j]; } } } for(int i=0;i<bp;i++){ for(int j=1;j<=n;j++){ pb[i+1][j]+=pb[i][j]; if(j<ls[i])pb[i+1][j]+=lower_bound(v1+ls[i],v1+rs[i]+1,v[j])-v1-ls[i]; else if(j>rs[i])pb[i+1][j]+=v1+rs[i]+1-upper_bound(v1+ls[i],v1+rs[i]+1,v[j]); } } scanf("%d",&m); while(m--){ int ans=0; scanf("%d%d",&l,&r); l^=la;r^=la; int lb=(l-1)/sz,rb=(r-1)/sz; if(lb==rb){ for(int i=l;i<=r;i++)v2[i]=v[i]; ans+=cal(l,r); }else{ for(int i=l;i<=rs[lb];i++)v2[i]=v[i]; if(l<rs[lb])ans+=cal(l,rs[lb]); for(int i=ls[rb];i<=r;i++)v2[i]=v[i]; if(ls[rb]<r)ans+=cal(ls[rb],r); for(int i=rs[lb],p=r+1;i>=l;i--){ while(p>ls[rb]&&v2[p-1]>=v2[i])--p; ans+=p-ls[rb]; } } lb++; rb--; if(lb<=rb){ for(int i=l;i<=ls[lb]-1;i++)ans+=pb[rb+1][i]-pb[lb][i]; for(int i=rs[rb]+1;i<=r;i++)ans+=pb[rb+1][i]-pb[lb][i]; for(int i=lb;i<=rb;i++)ans+=bs[i]; if(lb<rb) for(int i=lb;i<rb;i++) ans+=bb[i][rb+1]-bb[i][i+1]; } printf("%d\n",la=ans); } return 0; }
以上是关于bzoj3744 Gty的妹子序列的主要内容,如果未能解决你的问题,请参考以下文章