[Ynoi2019模拟赛]Yuno loves sqrt technology I

题目描述

给你一个长为n的排列，m次询问，每次查询一个区间的逆序对数，强制在线。

题解

MD不卡了。。TMD一点都卡不动。

强制在线的话也没啥好一点的方法，只能分块预处理了。

对于每个块，我们设lef[i]表示i到这个i这个元素所在块的块头的区间逆序对，rig[i]表示到块尾的逆序对。

在设一个cnt[i][j]表示从第i个块到第j个块的逆序对。

然后考虑如何处理询问。

整块之间的可以直接O(1)取答案，对于散块，我们需要求出散块内部的答案，散块对整块的答案，散块对散块的答案。

对于散块内的，因为散块在当前块中一定是它的前缀或后缀，所以我们直接O(1)取答案。

对于散块对整块的，我们可以在维护一个tag[i][j]表示前i个块，大于/小于j的元素有多少个，然后询问就扫描散块，前缀和统计即可。

对于散块对散块的，可以在每个块内维护元素的相对大小，然后对于两个不相交的区间，可以用桶排+扫描求出逆序对。

但还有一个特殊情况，就是lr在同一个块里，前面的方法不是很实用。

我们设当前块头尾h，那么答案可以表示为ans(h-r)-ans(h-l)-ans(h~l,l~r)，前两个部分已经求过了，后面的用桶排+扫描即可。

不过常数略大，开2s应该能过。

代码

#pragma GCC optimize(2)
#pragma GCC optimize(3)
#pragma GCC optimize ("Ofast")
#pragma GCC optimize ("Ofast","inline","unroll-loops")
#pragma GCC optimize ("no-stack-protector")
#include <cstdio>
#include <algorithm>
#include <iostream>
#include <cstring>
#include <cassert>
#include <cmath>
#define N 100002
#define M 340
using namespace std;
typedef long long ll;
int n,m,tr[N],mis[M],rnk[N],n1,ji[M],be[N],a[N],b[M],lef[N],tag1[M][N],tag2[M][N],rig[N],tong[N];
int q[N],pos[N];
ll ans,cnt[M][M];
inline void add(int x,int y){while(x<=n)tr[x]+=y,x+=x&-x;}
inline int query(int x){int ans=;while(x)ans+=tr[x],x-=x&-x;return ans;}
inline int merge(int l2,int r2,int l1,int r1){
    int ans=;
    for(int i=l1;i<=r1;++i)ji[rnk[i]]=a[i];
    q[]=;
    for(int i=;i<=n1;++i)if(ji[i]){
      q[++q[]]=ji[i];ji[i]=;
    }
    int p=;
    for(int i=l2;i<=r2;++i)ji[rnk[i]]=a[i];
    for(int i=;i<=n1;++i)if(ji[i]){
        while(p<=q[]&&(!p||q[p]<=ji[i]))p++;p--;
        ans+=p;ji[i]=;
    }
    return ans;
}
inline int rd(){
    int x=;char c=getchar();bool f=;
    while(!isdigit(c)){if(c=='-')f=;c=getchar();}
    while(isdigit(c)){x=(x<<)+(x<<)+(c^);c=getchar();}
    return f?-x:x;
}
int main(){
    n=rd();m=rd();
    n1=;
    for(int i=;i<=n;++i)a[i]=rd(),be[i]=(i-)/n1+,pos[a[i]]=i;
    for(int i=;i<=be[n];++i){
        int l=(i-)*n1+,r=min(n,i*n1);int top=;
        for(int j=l;j<=r;++j)b[++top]=a[j];
        sort(b+,b+top+);
        for(int j=;j<=top;++j)rnk[pos[b[j]]]=j;
        for(int j=l;j<=r;++j){
            lef[j]=query(n-a[j]+);
            if(j!=l)lef[j]+=lef[j-];
            add(n-a[j]+,);tong[n-a[j]+]=;
        }
        int moss=;
        for(int j=n;j>=;--j){
          moss+=tong[n-j+];
          tag2[i][j]=tag2[i-][j]+moss;
        }
        for(int j=l;j<=r;++j)add(n-a[j]+,-),tong[n-a[j]+]=;
        for(int j=r;j>=l;--j){
            rig[j]=query(a[j]);
            if(j!=r)rig[j]+=rig[j+];
            add(a[j],);tong[a[j]]=;
        }
        moss=;
        for(int j=;j<=n;++j){
          moss+=tong[j];
          tag1[i][j]=tag1[i-][j]+moss;
        }
        for(int j=l;j<=r;++j)add(a[j],-),tong[a[j]]=;
    }
    for(int i=;i<=be[n];++i){
        cnt[i][i]=lef[min(n,i*n1)];
        for(int j=i+;j<=be[n];++j){
            cnt[i][j]=cnt[i][j-];int x=min(n,n1*j);
            for(int k=(j-)*n1+;k<=x;++k){
                cnt[i][j]+=tag2[j-][a[k]]-tag2[i-][a[k]];
            }
            cnt[i][j]+=lef[x];
        }
    }
    int l,r;
    while(m--){
        l=rd();r=rd();l^=ans;r^=ans;
        if(l>r)l^=r^=l^=r;ans=;
      //  if(r>n||l<1){ans=0;puts("0");continue;}
        if(be[l]==be[r]){
              ans=lef[r];if(l!=(be[l]-)*n1+)ans-=lef[l-],ans-=merge((be[l]-)*n1+,l-,l,r);
        }
        else{
            int st=be[l]+,en=be[r]-;
           for(int i=l;i<=be[l]*n1;++i)ans=ans+tag1[en][a[i]]-tag1[st-][a[i]];
           for(int i=en*n1+;i<=r;++i)ans=ans+tag2[en][a[i]]-tag2[st-][a[i]];
           ans+=cnt[st][en]+rig[l]+lef[r];
           ans+=merge(l,(st-)*n1,en*n1+,r);
        }
         printf("%lld\n",ans);
    }
    return ;
}