bzoj3744: Gty的妹子序列

最新推荐文章于 2020-02-14 17:23:10 发布

zhouyuyang233

最新推荐文章于 2020-02-14 17:23:10 发布

阅读量368

点赞数

分类专栏：数据结构

本文链接：https://blog.csdn.net/zhouyuyang233/article/details/69371918

版权

数据结构专栏收录该内容

32 篇文章 0 订阅

订阅专栏

传送门
如果不强制在线的话就是莫队水题。
可是~~~
果断分块走起。
预处理出每一个快到其他节点的逆序对个数
用数状数组实现。
前面的就用主席树。
时间复杂度O（nsqrt(n)logn）

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
#define fi first
#define se second
#define M 50005
using namespace std;
struct node{
    node *ls,*rs;
    int num;
    void* operator new(size_t size,node *_,node *__,int ___);
}*t[M],mem[1001001],*C=mem;
int n,m,tot,sz,ans,x,y;
int a[M],l[M],r[M],be[M],cnt[250][M];
pair<int,int> b[M];
int c[M],tim[M],T;
void* node::operator new(size_t size,node *_,node *__,int ___){
    C->ls=_;
    C->rs=__;
    C->num=___;
    return C++;
}
void change(int x){
    for (;x;x-=x&-x){
        if (tim[x]!=T) tim[x]=T,c[x]=0;
        c[x]++;
    }
}
int ask(int x){
    int s=0;
    for (;x<=tot;x+=x&-x)
        if (tim[x]==T) s+=c[x];
    return s;
}
node* build(node *p,int x,int y,int v){
    int mid=(x+y)/2;
    if (x==y) return new(0x0,0x0,p->num+1)node;
    if (v<=mid) return new(build(p->ls,x,mid,v),p->rs,p->num+1)node;
    return new(p->ls,build(p->rs,mid+1,y,v),p->num+1)node;
}
int ask(node *p1,node *p2,int x,int y,int l,int r){
    int mid=(x+y)/2;
    if (p1->num==p2->num) return 0;
    if (x==l&&y==r) return p2->num-p1->num;
    if (r<=mid) return ask(p1->ls,p2->ls,x,mid,l,r);
    if (l>mid) return ask(p1->rs,p2->rs,mid+1,y,l,r);
    return ask(p1->ls,p2->ls,x,mid,l,mid)+ask(p1->rs,p2->rs,mid+1,y,mid+1,r);
}
int query(int x,int y){
    int s=0;
    if (be[x]==be[y]){
        T++;
        for (int i=x;i<=y;i++)
            s+=ask(a[i]+1),change(a[i]);
        return s;
    }
    s=cnt[be[x]+1][y];
    for (int i=x;i<=r[be[x]];i++)
        s+=ask(t[i],t[y],0,tot+1,0,a[i]-1);
    return s;
}
int main(){
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
        scanf("%d",&b[i].fi),b[i].se=i;
    sort(b+1,b+n+1);
    for (int i=1;i<=n;i++){
        if (i==1||b[i].fi!=b[i-1].fi) ++tot;
        a[b[i].se]=tot;
    }
    t[0]=new(0x0,0x0,0)node;
    t[0]->ls=t[0]->rs=t[0];
    for (int i=1;i<=n;i++)
        t[i]=build(t[i-1],0,tot+1,a[i]);
    sz=(int)(sqrt(n)+1e-7);
    for (int i=1;i<=n;i++)
        be[i]=(i-1)/sz+1;
    for (int i=1;(i-1)*sz<=n;i++)
        l[i]=(i-1)*sz+1,r[i]=min(sz*i,n);
    for (int i=1;(i-1)*sz<=n;i++){
        ++T;
        for (int j=(i-1)*sz+1;j<=n;j++){
            cnt[i][j]=cnt[i][j-1]+ask(a[j]+1);
            change(a[j]);
        }
    }
    scanf("%d",&m);
    while (m--){
        scanf("%d%d",&x,&y);
        x^=ans; y^=ans;
        printf("%d\n",ans=query(x,y));
    }
}