这是本蒟蒻的第一篇博客,在此手动庆祝一下~~
其实今天本来想写一道树链剖分的题但是一直没过(蒟蒻的无奈)……所以只好写前几天莫队的题解啦
这道题分为两个知识点,树状数组和莫队
1。莫队 就是最基本的套路啦,分块+离线排序
2。树状数组 我做这道题本来是奔着莫队的,没想最开始还被树状数组坑了
一定要注意分类讨论两点:
1)pl->a[i].l 和pr->a[i].r 的变化会使ans加还是减
2)
pl->a[i].l ,所影响的是增减的数和它后面的数所组成的逆序对,所以查询树状数组最后return x后面的数
pr->a[i].r ,所影响的是增减的数和它前面的数所组成的逆序对,所以查询树状数组最后return x前面的数
代码见下:
/**************************************************************
Problem: 3289
User: DD_D
Language: C++
Result: Accepted
Time:5544 ms
Memory:2400 kb
****************************************************************/
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int N=50000+5;
int siz[N],val[N],pos[N];
int tree[N],rt[N];
int n,q,sum,ans;
struct data{
int l,r,id;
}a[N];
int cmp(data a,data b)
{
if(pos[a.l]==pos[b.l]) return a.r<b.r;
return a.l<b.l;
}
int cmp_id(data a,data b)
{
return a.id<b.id;
}
int lowbit(int x)
{
return x&(-x);
}
void add(int x,int d)
{
for(int i=x;i<=n;i+=lowbit(i))
tree[i]+=d;
sum+=d;
}
int query1(int x)//统计树状数组后面
{
int tot=0;
for(int i=x;i>0;i-=lowbit(i))
tot+=tree[i];
return sum-tot;
}
int query2(int x)//统计树状数组前面
{
int tot=0;
for(int i=x-1;i>0;i-=lowbit(i))
tot+=tree[i];
return tot;
}
void init()
{
scanf("%d",&n);
int block=(int)sqrt(n);
for(int i=1;i<=n;i++){
pos[i]=(i-1)/block+1;
scanf("%d",&siz[i]);
val[i]=siz[i];
}
sort(val+1,val+1+n);
int tot=unique(val+1,val+1+n)-val;
for(int i=1;i<=n;i++)
siz[i]=lower_bound(val+1,val+tot,siz[i])-val;
scanf("%d",&q);
for(int i=1;i<=q;i++){
scanf("%d%d",&a[i].l,&a[i].r);
a[i].id=i;
}
}
void solve()
{
for(int i=a[1].l;i<=a[1].r;i++)
add(siz[i],1),ans+=query1(siz[i]);
rt[a[1].id]=ans;
int pl=a[1].l,pr=a[1].r;
for(int i=2;i<=q;i++){
for(;pl<a[i].l;pl++)//多了 统计前面的
add(siz[pl],-1),ans-=query2(siz[pl]);
for(;pl>a[i].l;pl--)//已有的数都是在它的后面,所以统计树状数组前面
add(siz[pl-1],1),ans+=query2(siz[pl-1]);
for(;pr<a[i].r;pr++)//少了 统计后面的
add(siz[pr+1],1),ans+=query1(siz[pr+1]);
for(;pr>a[i].r;pr--)//多了
add(siz[pr],-1),ans-=query1(siz[pr]);
if(a[i].r==a[i].l) {
rt[a[i].id]=0;
continue;
}
rt[a[i].id]=ans;
}
}
int main()
{
init();
sort(a+1,a+1+q,cmp);
solve();
sort(a+1,a+1+q,cmp_id);
for(int i=1;i<=q;i++)
printf("%d\n",rt[i]);
return 0;
}