bzoj 3744: Gty的妹子序列 分块+树状数组

题意

给出n个数，要求资瓷求区间逆序对，并要求强制在线。
n,m<=50000

分析

先把所有数离散化，这样就可以方便用树状数组啦。
然后分一波块。
然后我们在$O(n\sqrt nlog)$的时间内预处理出两个数组：
ans[i,j]表示第i块到第j个数之间有多少个逆序对。
sum[i,j]表示前j块里面有多少个数不大于i。
那么设现在要询问的区间为[l,r]
假如l和r在同一块内则用树状数组直接统计即可。
若不在同一块内则先求出l所在的块+1到r的逆序对数即为ans[pos[i]+1][r]，然后我们再统计l到end[l]这几个数对答案的贡献就好啦。

1A就是爽！

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 50005
#define M 305
using namespace std;

int n,ans[M][N],sum[N][M],a[N],c[N],rank,pos[N],sta[M],end[M],block;
struct data{int num,id;}b[N];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

bool cmp(data a,data b)
{
    return a.num<b.num;
}

void ins(int x,int y)
{
    while (x<=rank)
    {
        c[x]+=y;
        x+=x&(-x);
    }
}

int find(int x)
{
    int ans=0;
    while (x)
    {
        ans+=c[x];
        x-=x&(-x);
    }
    return ans;
}

int query(int l,int r)
{
    if (pos[l]==pos[r])
    {
        int tot=0;
        for (int i=l;i<=r;i++)
        {
            tot+=find(rank)-find(a[i]);
            ins(a[i],1);
        }
        for (int i=l;i<=r;i++)
            ins(a[i],-1);
        return tot;
    }
    int tot=ans[pos[l]+1][r];
    for (int i=l;i<=end[pos[l]];i++)
    {
        tot+=sum[a[i]-1][pos[r]-1]-sum[a[i]-1][pos[l]]+find(rank)-find(a[i]);
        ins(a[i],1);
    }
    for (int i=sta[pos[r]];i<=r;i++)
        tot+=find(rank)-find(a[i]);
    for (int i=l;i<=end[pos[l]];i++)
        ins(a[i],-1);
    return tot;
}

int main()
{
    n=read();
    block=sqrt(n);
    for (int i=1;i<=n;i++)
    {
        a[i]=read();
        pos[i]=(i+block-1)/block;
        if (!sta[pos[i]]) sta[pos[i]]=i;
        end[pos[i]]=i;
        b[i].num=a[i];
        b[i].id=i;
    }
    sort(b+1,b+n+1,cmp);
    rank=0;
    for (int i=1;i<=n;i++)
    {
        if (b[i].num!=b[i-1].num) rank++;
        a[b[i].id]=rank;
    }
    for (int i=1;i<=pos[n];i++)
    {
        for (int j=sta[i];j<=n;j++)
        {
            ans[i][j]=ans[i][j-1]+find(rank)-find(a[j]);
            ins(a[j],1);
        }
        for (int j=sta[i];j<=n;j++)
            ins(a[j],-1);
    }
    for (int i=1;i<=n;i++)
        sum[a[i]][pos[i]]++;
    for (int i=1;i<=rank;i++)
        for (int j=1;j<=pos[n];j++)
            sum[i][j]+=sum[i][j-1];
    for (int j=1;j<=pos[n];j++)
        for (int i=1;i<=rank;i++)
            sum[i][j]+=sum[i-1][j];
    int lastans=0;
    int m=read();
    for (int i=1;i<=m;i++)
    {
        int l=read(),r=read();
        l^=lastans;r^=lastans;
        lastans=query(l,r);
        printf("%d\n",lastans);
    }
    return 0;
}