Description
我早已习惯你不在身边，
人间四月天 寂寞断了弦。
回望身后蓝天，
跟再见说再见……
某天,蒟蒻Autumn发现了从 Gty的妹子树上掉落下来了许多妹子,他发现
她们排成了一个序列,每个妹子有一个美丽度。
Bakser神犇与他打算研究一下这个妹子序列,于是Bakser神犇问道:"你知道区间([l,r])中妹子们美丽度的逆序对数吗?"
蒟蒻Autumn只会离线乱搞啊……但是Bakser神犇说道:"强制在线。"
请你帮助一下Autumn吧。
给定一个正整数序列(a),对于每次询问,输出(a_{l}...a_{r})中的逆序对数,强制在线。
Input
第一行包括一个整数(n(1 le n le 50000)),表示数列(a)中的元素数。
第二行包括(n)个整数(a_{1}...a_{n})((a_{i}>0),保证(a_{i})在(int)内)。
接下来一行包括一个整数(m(1 le m le 50000)),表示询问的个数。
接下来m行,每行包括(2)个整数(l,r(1 le l le r le n)),表示询问(a_{l}...a_{r})中的逆序
对数(若(a_{i}>a_{j})且(a_{i}<a_{j}),则为一个逆序对)。
(l,r)要分别异或上一次询问的答案((lastans)),最开始时(lastans=0)。
保证涉及的所有数在(int)内。
Output
对每个询问,单独输出一行,表示(a_{l}...a_{r})中的逆序对数。
Sample Input
4
1 4 2 3
1
2 4
Sample Output
2
一道很典型的分块题目。
首先将序列离散化，之后分块记录这些值：(have[i][j])表示第(i)个块中小于等于(j)的元素个数；(rev[i])表示第(i)个块的逆序对个数；(ref2[i][j])表示(i)块与下标小于等于(j)的序列元素构成的逆序对数((A_{j})不在第(i)块内)；(ref3[i][j])第(i)块与前(j)块所构成逆序对数目。（预处理在(O(nsqrt{n}))的时间内完成）
询问的话，块与块之间的我们可以用(ref3)在(O(sqrt{n}))的时间内算出，块与散部可以用(ref2)在(O(sqrt{n}))的时间内算出。散与散的只能用树状数组暴力求解，在(O(sqrt{n}logn))的时间内算出结果。
综上，本题时间复杂度(O(nsqrt{n}+msqrt{n}logn))，空间复杂度(O(nsqrt{n}))。有不理解的参考代码。
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
using namespace std;
#define maxn 50010
const int len = 230;
int n,m,ans,A[maxn],bac[maxn],have[len][maxn];
int rev[maxn],tree[maxn],ref2[len][maxn],ref3[len][len],tot;
inline int begin(int ord) { return (ord-1)*len+1; }
inline int end(int ord) { return min(ord*len,n); }
inline int belong(int ord) { return (ord+len-1)/len; }
inline int size(int ord) { return end(ord)-begin(ord)+1; }
inline int lowbit(int x) { return x & -x; }
inline void ins(int x,int y) { for (;x <= n;x += lowbit(x)) tree[x] += y; }
inline int calc(int x) { int ret = 0; for (;x;x -= lowbit(x)) ret += tree[x]; return ret; }
inline void ready()
{
for (int i = 1;i <= n;++i) bac[i] = A[i];
sort(bac+1,bac+n+1); tot = unique(bac+1,bac+n+1)-bac-1;
for (int i = 1;i <= n;++i) A[i] = lower_bound(bac+1,bac+tot+1,A[i])-bac;
tot = (n + len-1)/len;
for (int i = 1;i <= tot;++i)
{
int l = begin(i),r = end(i);
for (int j = l;j <= r;++j) have[i][A[j]]++;
for (int j = 1;j <= n;++j) have[i][j] += have[i][j-1];
for (int j = l;j <= r;++j) for (int k = j+1;k <= r;++k) rev[i] += A[j] > A[k];
}
for (int i = 1;i <= tot;++i)
{
int l = begin(i),r = end(i);
for (int j = 1;j < l;++j)
ref2[i][j] += ref2[i][j-1],ref2[i][j] += have[i][A[j]-1];
for (int j = r+1;j <= n;++j)
ref2[i][j] += ref2[i][j-1],ref2[i][j] += size(i)-have[i][A[j]];
}
for (int i = 1;i <= n;++i)
{
int now = belong(i);
for (int j = 1;j < now;++j)
ref3[now][j] += size(j)-have[j][A[i]];
}
for (int i = 1;i <= tot;++i)
for (int j = 2;j < i;++j)
ref3[i][j] += ref3[i][j-1];
}
inline int work(int l,int r)
{
int tot = 0,L = belong(l),R = belong(r),p = end(L),q = begin(R);
for (int i = L+1;i < R;++i)
{
tot += rev[i];
tot += ref3[i][i-1] - ref3[i][L];
tot += ref2[i][p]-ref2[i][l-1];
tot += ref2[i][r]-ref2[i][q-1];
}
if (L != R)
{
for (int i = r;i >= q;--i) tot += calc(A[i]-1),ins(A[i],1);
for (int i = p;i >= l;--i) tot += calc(A[i]-1),ins(A[i],1);
for (int i = r;i >= q;--i) ins(A[i],-1);
for (int i = p;i >= l;--i) ins(A[i],-1);
}
else
{
for (int i = r;i >= l;--i) tot += calc(A[i]-1),ins(A[i],1);
for (int i = r;i >= l;--i) ins(A[i],-1);
}
return tot;
}
int main()
{
freopen("3744.in","r",stdin);
freopen("3744.out","w",stdout);
scanf("%d",&n);
for (int i = 1;i <= n;++i) scanf("%d",A+i);
ready();
scanf("%d",&m);
for (int i = 1;i <= m;++i)
{
int l,r; scanf("%d %d",&l,&r);
l ^= ans; r ^= ans;
printf("%d
",ans = work(l,r));
}
fclose(stdin); fclose(stdout);
return 0;
}