bzoj3289: Mato的文件管理

3289: Mato的文件管理

Time Limit: 40 Sec  Memory Limit: 128 MB
Submit: 3047  Solved: 1293
[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]区间的文件。

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次移到最后。

Source

By taorunz

嗯

我的第二道莫队

今天临时补的树状数组和用它来求逆序对

还是理解不够深刻啊

问了师兄才懂

建议可以先去写一下poj2299 Ultra-QuickSort

写了那道题之后

本道题的难点就只有那四个for了

可以在草稿纸上画一下区间

然后比较好想

看代码和注释吧

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
const int N=50005;
struct node
{
	int order,v;
}a[N];
struct edgt
{
	int l,r,id;
}q[N];
int aa[N],c[N],pos[N];int n,m;long long now=0,ans[N];
int lowbit(int i)
{
	return i&-i;
}
bool cmp(node a,node b){	return a.v<b.v;	}
bool cmp2(edgt a,edgt b)
{
	if(pos[a.l]==pos[b.l])	return a.r<b.r;
	return a.l<b.l;
}
void add(int x,int add)
{
	for(int i=x;i<=n;i+=lowbit(i))	c[i]+=add;
}
long long int qsum(int x)
{
	long long int sum=0;
	for(int i=x;i>=1;i-=lowbit(i))	sum+=c[i];
	return sum;
}
void solve()
{
	for(int i=1,l=1,r=0;i<=m;i++)
	{
		while(r<q[i].r)
		r++,add(aa[r],1),now+=r-l+1-qsum(aa[r]);	
		while(r>q[i].r)
		add(aa[r],-1),now-=r-l-qsum(aa[r]),r--;//此时区间长度只有r-l //删除就直接-1啦 
		while(l>q[i].l)
		l--,add(aa[l],1),now+=qsum(aa[l]-1);//因为丢aa[l]是<=，然后区间长度减去个数就是严格大于的，而这里要求严格小于的，所以要aa[l]-1 
		while(l<q[i].l)
		add(aa[l],-1),now-=qsum(aa[l]-1),l++;
		ans[q[i].id]=now;
	}
}
int main()
{
	int block;
	scanf("%d",&n);block=sqrt(n);
	for(int i=1;i<=n;i++)	scanf("%d",&a[i].v),a[i].order=i;
	std::sort(a+1,a+1+n,cmp);
	for(int i=1;i<=n;i++)	aa[a[i].order]=i;
	scanf("%d",&m);
	for(int i=1;i<=m;i++)	scanf("%d %d",&q[i].l,&q[i].r),q[i].id=i;
	for(int i=1;i<=n;i++)	pos[i]=(i-1)/block+1;
	std::sort(q+1,q+1+m,cmp2);
	solve();
	for(int i=1;i<=m;i++)	printf("%lld\n",ans[i]);
	return 0;
}