bzoj3289: Mato的文件管理

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]区间的文件。

Output

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

Sample Input

4
1 4 2 3
2
1 2
2 4

Sample Output

0
2

HINT

Hint

n,q <= 50000

样例解释：第一天，Mato不需要交换

第二天，Mato可以把2号交换2次移到最后。

莫队+树状数组

最小的交换次数肯定是逆序对数，用树状数组维和即可

附：这题应该要离散化，不过我偷懒没写了（rank快垫底了）

code：

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cmath>
 4 #include<cstring>
 5 #include<algorithm>
 6 #define maxn 50005
 7 #define lowbit(x) ((x)&(-(x)))
 8 using namespace std;
 9 char ch;
10 int n,q,siz,tmp,a[maxn],pos[maxn],ans[maxn];
11 struct query{
12     int l,r,id;
13 }list[maxn];
14 bool cmp(query a,query b){
15     if (pos[a.l]!=pos[b.l]) return pos[a.l]<pos[b.l];
16     return pos[a.r]<pos[b.r];
17 }
18 struct bit{
19     int v[maxn*20];
20     void insert(int x,int add){for (;x<=maxn*20;x+=lowbit(x)) v[x]+=add;}
21     int query(int x){
22         int ans=0;
23         for (;x;x-=lowbit(x)) ans+=v[x];
24         return ans;
25     }
26 }T;
27 bool ok;
28 void read(int &x){
29     for (ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if (ch=='-') ok=1;
30     for (x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar());
31     if (ok) x=-x;
32 }
33 void work(){
34     int l=1,r=0;
35     for (int i=1;i<=q;i++){
36         for (;r<list[i].r;r++) T.insert(a[r+1],1),tmp+=(r+1-l+1)-T.query(a[r+1]);
37         for (;r>list[i].r;r--) T.insert(a[r],-1),tmp-=(r-l)-T.query(a[r]);
38         for (;l<list[i].l;l++) T.insert(a[l],-1),tmp-=T.query(a[l]-1);
39         for (;l>list[i].l;l--) T.insert(a[l-1],1),tmp+=T.query(a[l-1]-1);        
40         ans[list[i].id]=tmp;
41     }
42 }
43 int main(){
44     read(n),siz=(int)sqrt(n);
45     for (int i=1;i<=n;i++) read(a[i]);
46     read(q);
47     for (int i=1;i<=q;i++) read(list[i].l),read(list[i].r),list[i].id=i;
48     for (int i=1;i<=n;i++) pos[i]=i/siz;
49     sort(list+1,list+q+1,cmp);
50     work();
51     for (int i=1;i<=q;i++) printf("%d
",ans[i]);
52     return 0;
53 }