bzoj 3289 Mato的文件管理 树状数组+莫队

Mato的文件管理

Time Limit: 40 Sec  Memory Limit: 128 MB
Submit: 4325  Solved: 1757
[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次移到最后。

HINT

 

题解：树状数组维护逆序对，然后莫队降低复杂度，要换绝对要换

 #include<cstring>
 #include<cstdio>
 #include<algorithm>
 #include<iostream>
 #include<cmath>
 #include<queue>

 #define N 50007
 using namespace std;
 inline int read()
 {
     int x=,f=;char ch=getchar();
     while(!isdigit(ch)){if(ch=='-')f=-;ch=getchar();}
     while(isdigit(ch)){x=(x<<)+(x<<)+ch-'';ch=getchar();}
     return x*f;
 }

 unsigned int now,ans[N];
 int n,m;
 int a[N],disc[N],belong[N],t[N];
 struct query
 {
     int l,r,id;
 }q[N];

 bool operator<(query a,query b)
 {
     if(belong[a.l]==belong[b.l])return a.r<b.r;
     return belong[a.l]<belong[b.l];
 }
 void add(int x,int val)
 {
     for(int i=x;i<=n;i+=i&-i)
         t[i]+=val;
 }
 unsigned int query(int x)
 {
     unsigned int sum=;
     for(int i=x;i;i-=i&-i)
         sum+=t[i];
     return sum;
 }
 void solve()
 {
     sort(q+,q+m+);
     int l=,r=;
     for(int i=;i<=m;i++)
     {
         while(l<q[i].l)
             add(a[l],-),now-=query(a[l]-),l++;
         while(r>q[i].r)
             add(a[r],-),now-=r-l-query(a[r]),r--;
         while(l>q[i].l)
             l--,add(a[l],),now+=query(a[l]-);
         while(r<q[i].r)
             r++,add(a[r],),now+=r-l+-query(a[r]);
         ans[q[i].id]=now;
     }
 }
 int main()
 {
     n=read();int t=sqrt(n);
     for(int i=;i<=n;i++)disc[i]=a[i]=read();
     sort(disc+,disc+n+);
     for(int i=;i<=n;i++)a[i]=lower_bound(disc+,disc+n+,a[i])-disc;
     m=read();
     for(int i=;i<=m;i++)
         q[i].l=read(),q[i].r=read(),q[i].id=i;
     for(int i=;i<=n;i++)
         belong[i]=(i-)/t+;
     solve();
     for(int i=;i<=m;i++)
         printf("%d\n",ans[i]);
 }

By taorunz