3289: Mato的文件管理
Time Limit: 40 Sec Memory Limit: 128 MBSubmit: 303 Solved: 122
[ Submit][ Status]
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
1 4 2 3
2
1 2
2 4
Sample Output
0
2
2
【题解】就是给你若干询问的区间,让你求出区间内的值完成排序的相邻交换需要多少次(即求区间内逆序对的数量)。求逆序对数,有 二分 和 树状数组(我使用了树状数组,此题中两个O()理论上都是一样的,但我觉得树状数组好写,常数又小)这里面可能有重复的情况,不去理他就可以了。至于到有多组询问,暴力应该是要萎的,所以我用莫队算法做。
基本步骤具体如何推,详见代码。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define maxn 51001
using namespace std;
const int S=222;
struct pp{int id,v;}a[maxn];
int n,w[maxn],f[maxn],L[maxn],R[maxn],p[maxn],num[maxn];
long long val[maxn];
inline void insert(int x,int y){for(;x;f[x]+=y,x-=x&(-x));}
inline int query(int x){int tot=0;for(;x<=n;tot+=f[x],x+=x&(-x));return tot;}
bool cmp(pp x,pp y){
if(x.v<y.v)return 1;
else if(x.v>y.v)return 0;
else if(x.id<y.id)return 1;
else return 0;
}
bool cmpx(int x,int y)
{return L[x]/S==L[y]/S?((L[x]/S)&1)^(R[x]<R[y]):L[x]<L[y];}
inline void read(int &x)
{
x=0;char ch;
while(ch=getchar(),ch==' ' || ch=='\r' || ch=='\n');
while(ch!=' ' && ch!='\r' && ch!='\n' )x=x*10+(int)ch-48,ch=getchar();
}
int main()
{
read(n);
for(int i=1;i<=n;i++)
{
read(a[i].v);a[i].id=i;
}
//离散;
sort(a+1,a+1+n,cmp);
int l=1;w[a[1].id]=1;
for(int i=2;i<=n;i++)
if(a[i].v==a[i-1].v)w[a[i].id]=l;
else w[a[i].id]=++l;
//莫队思想之分块实现;
int q;read(q);
for(int i=1;i<=q;i++)read(L[i]),read(R[i]),p[i]=i;
sort(p+1,p+1+q,cmpx);
memset(f,0,sizeof(f));
memset(num,0,sizeof(num));
int x=L[p[1]],y=R[p[1]];
long long ans=0;
for(int i=x;i<=y;i++)ans+=query(w[i]),insert(w[i],1),ans-=num[w[i]],num[w[i]]++;
val[p[1]]=ans;
for(int i=2;i<=q;i++)
{
while(x<L[p[i]])num[w[x]]--,insert(w[x],-1),ans-=query(1)-query(w[x++]);
while(x>L[p[i]])ans+=query(1)-query(w[--x]),insert(w[x],1),num[w[x]]++;
while(y<R[p[i]])ans+=query(w[++y]),insert(w[y],1),ans-=num[w[y]],num[w[y]]++;
while(y>R[p[i]])num[w[y]]--,insert(w[y],-1),ans-=query(w[y])-num[w[y--]];
val[p[i]]=ans;
}
for(int i=1;i<=q;i++)printf("%lld\n",val[i]);
return 0;
}