[bzoj3744] Gty的妹子序列

　　还是分块>_<

　　强制在线求区间逆序。。如果没有空间限制的话显然是树套树？

　　分成根号n块，预处理出sum[i][j]表示从第i块的左端，到第j个点 这段区间内的逆序对个数。这部分可以直接树状数组O(n^1.5 logn)

　　那么对于每个询问，我们只要考虑左侧的零散节点对答案的贡献了。。

　　因此需要快速求出整块以及右侧零散节点插入后的树状数组。

　　再预处理出t[i][]表示前i块的数  的树状数组。。所以整块整块的树状数组就可以由前缀和得出。

　　再暴力插入右侧零散节点。。最后左侧零散节点就边插边统计了。。

　　总的时间复杂度O(n^1.5 logn)

　　块的大小就是根号n吧。。。不过具体来说比根号n小一些比较优

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn=50003;
struct zs{
    int v;unsigned short pos;
}b[maxn];int cnt;
int id[maxn],mp[maxn],l[305],r[305];
int sm[305][maxn],t[305][maxn],tt[maxn];
int i,j,k,n,m,x,y,kuai,la;
  
int ra,fh;char rx;
inline int read(){
    rx=getchar(),fh=1,ra=0;
    while((rx<'0'||rx>'9')&&rx!='-')rx=getchar();
    if(rx=='-')fh=-1,rx=getchar();
    while(rx>='0'&&rx<='9')ra*=10,ra+=rx-48,rx=getchar();return ra*fh;
}
  
bool cmp(zs a,zs b){return a.v<b.v;}
inline int query(int x,int y){
    register int j;int i,ans=0,k=id[y]-1;
    if(x==l[id[x]])return sm[id[x]][y];
    if(id[x]==id[y]){
        for(i=x;i<=y;i++){
            for(j=mp[i],ans+=i-x;j;j-=j&-j)ans-=tt[j];
            for(j=mp[i];j<=cnt;j+=j&-j)++tt[j];
        }
        for(i=x;i<=y;i++)
            for(j=mp[i];j<=cnt&&tt[j]!=0;j+=j&-j)tt[j]=0;
        return ans;
    }
    ans=sm[id[x]+1][y];
    for(i=y;i>=l[id[y]];i--)for(j=mp[i];j<=cnt;j+=j&-j)++tt[j];
    if(id[x]+1<id[y])
        for(i=r[id[x]];i>=x;i--){
            for(j=mp[i]-1;j;j-=j&-j)ans+=tt[j]+t[k][j]-t[id[x]][j];
            for(j=mp[i];j<=cnt;j+=j&-j)++tt[j];
        }
    else{
        for(i=r[id[x]];i>=x;i--){
            for(j=mp[i]-1;j;j-=j&-j)ans+=tt[j];
            for(j=mp[i];j<=cnt;j+=j&-j)++tt[j];
        }
    }
    for(i=r[id[x]];i>=x;i--)for(j=mp[i];j<=cnt&&tt[j]!=0;j+=j&-j)tt[j]=0;
    for(i=y;i>=l[id[y]];i--)for(j=mp[i];j<=cnt&&tt[j]!=0;j+=j&-j)tt[j]=0;
    return ans;
}
int main(){
    n=read();kuai=min(n,166);//kuai=1;
    for(i=1;i<=n;i++)b[i].v=read(),b[i].pos=i,id[i]=(i+kuai-1)/kuai;
      
    for(i=1;i<=id[n];i++)l[i]=(i-1)*kuai+1,r[i]=i<id[n]?kuai*i:n;
    sort(b+1,b+1+n,cmp);
    for(i=1;i<=n;mp[b[i].pos]=cnt,i++)
        if(b[i].v!=b[i-1].v||i==1)cnt++;
    for(i=1;i<=id[n];i++){
        for(j=l[i];j<=r[i];j++)++t[i][mp[j]];
        for(j=1;j<=cnt;t[i][j]+=t[i-1][j],j++)if(j+(j&-j)<=cnt)t[i][j+(j&-j)]+=t[i][j];
    }
     
    for(i=1;i<=id[n];i++){
        if(i>1)memset(tt,0,(cnt+1)<<2);
        for(j=l[i];j<=n;j++){
            for(sm[i][j]=sm[i][j-1]+j-l[i],k=mp[j];k;k-=k&-k)sm[i][j]-=tt[k];
            for(k=mp[j];k<=cnt;k+=k&-k)++tt[k];
        }
    }memset(tt,0,(cnt+1)<<2);
     
    for(m=read();m;m--)
        x=read()^la,y=read()^la,
        la=query(x,y),printf("%d\n",la);
    return 0;
}

　　一开始没删调试挂了两发= =我果然是傻逼QAQ