BZOJ 3289: Mato的文件管理 【莫队 + 树状数组】

任意门：https://www.lydsy.com/JudgeOnline/problem.php?id=3289

3289: Mato的文件管理

Time Limit: 40 Sec  Memory Limit: 128 MB
Submit: 4917  Solved: 2041
[Submit][Status][Discuss]

Description

Mato同学从各路神犇以各种方式（你们懂的）收集了许多资料，这些资料一共有n份，每份有一个大小和一个编号

。为了防止他人偷拷，这些资料都是加密过的，只能用Mato自己写的程序才能访问。Mato每天随机选一个区间[l,r

]，他今天就看编号在此区间内的这些资料。Mato有一个习惯，他总是从文件大小从小到大看资料。他先把要看的

文件按编号顺序依次拷贝出来，再用他写的排序程序给文件大小排序。排序程序可以在1单位时间内交换2个相邻的

文件（因为加密需要，不能随机访问）。Mato想要使文件交换次数最小，你能告诉他每天需要交换多少次吗？

Input

第一行一个正整数n，表示Mato的资料份数。

第二行由空格隔开的n个正整数，第i个表示编号为i的资料的大小。

第三行一个正整数q，表示Mato会看几天资料。

之后q行每行两个正整数l、r，表示Mato这天看[l,r]区间的文件。

n,q <= 50000

Output

q行，每行一个正整数，表示Mato这天需要交换的次数。

Sample Input

4
1 4 2 3
2
1 2
2 4

Sample Output

0
2
//样例解释：第一天，Mato不需要交换
第二天，Mato可以把2号交换2次移到最后。

题意概括：

其实就是查询区间的逆序数。

解题思路：

首先根据数据范围，可以先对数据进行离散化。

用树状数组维护当前区间出现的数，那么 X 的逆序数就是在 X 之前已存在当前区间的比 X 大的数，通过树状数组这个很容易实现，然后通过莫队离线操作进行区间的转移。

AC code:

 #include <bits/stdc++.h>
 #define INF 0x3f3f3f3f
 #define LL long long
 using namespace std;
 const int MAXN = 5e4+;

 int N, M;
 struct Query
 {
     int l, r, idx;
 }Q[MAXN];
 int unit;
 int a[MAXN], b[MAXN], ans[MAXN];
 int t[MAXN];
 int lowbit(int x){return x&(-x);}
 bool cmp(struct Query a, struct Query b)
 {
     int la = a.l/unit, lb = b.l/unit;
     if(la == lb) return a.r < b.r;
     return la < lb;
 }

 void add(int x, int val)
 {
     for(int i = x; i < MAXN; i+=lowbit(i))
         t[i] += val;
 }

 int ask(int x)
 {
     int res = ;
     for(int i = x; i; i-=lowbit(i))
         res+=t[i];
     return res;
 }

 void solve()
 {
     int L = , R = , cur = ;
     for(int i = ; i <= M; i++){
         while(L < Q[i].l) {cur-=ask(a[L]-); add(a[L], -);L++;}
         while(L > Q[i].l) {cur+=ask(a[L-] - ); add(a[L-], ); L--;}
         while(R < Q[i].r) {cur+=R-L+-ask(a[R+]); add(a[R+], ); R++;}
         while(R > Q[i].r) {cur-=R-L+-ask(a[R]); add(a[R], -); R--;}
         ans[Q[i].idx] = cur;
     }

     for(int i = ; i <= M; i++){
         printf("%d\n", ans[i]);
     }
 }

 int main()
 {
     scanf("%d", &N);
     unit = sqrt(N);
     for(int i = ; i <= N; i++){ scanf("%d", &a[i]); b[i] = a[i];}
     scanf("%d", &M);
     for(int i = ; i <= M; i++){
         scanf("%d %d", &Q[i].l, &Q[i].r);
         Q[i].idx = i;
     }
     sort(b+, b+N+);                   //离散化
     int cnt = unique(b+, b++N)-b-;
     for(int i = ; i <= N; i++){
         a[i] = lower_bound(b+, b++cnt, a[i]) - b;
     }

     sort(Q+, Q++M, cmp);              //莫队分区

     solve();

     return ;
 }