BZOJ3744: Gty的妹子序列

3744: Gty的妹子序列

Time Limit: 20 Sec Memory Limit: 128 MB
Submit: 111 Solved: 38
[Submit][Status]

Description

我早已习惯你不在身边，

人间四月天 寂寞断了弦。

回望身后蓝天，

跟再见说再见……

某天,蒟蒻Autumn发现了从 Gty的妹子树(bzoj3720) 上掉落下来了许多妹子,他发现

她们排成了一个序列,每个妹子有一个美丽度。

Bakser神犇与他打算研究一下这个妹子序列,于是Bakser神犇问道:"你知道区间

[l,r]中妹子们美丽度的逆序对数吗?"

蒟蒻Autumn只会离线乱搞啊……但是Bakser神犇说道:"强制在线。"

请你帮助一下Autumn吧。

给定一个正整数序列a,对于每次询问,输出al...ar中的逆序对数,强制在线。

Input

第一行包括一个整数n(1<=n<=50000),表示数列a中的元素数。

第二行包括n个整数a1...an(ai>0,保证ai在int内)。

接下来一行包括一个整数m(1<=m<=50000),表示询问的个数。

接下来m行,每行包括2个整数l、r(1<=l<=r<=n),表示询问al...ar中的逆序

对数(若ai>aj且i<j,则为一个逆序对)。

l,r要分别异或上一次询问的答案(lastans),最开始时lastans=0。

保证涉及的所有数在int内。

Output

对每个询问,单独输出一行,表示al...ar中的逆序对数。

Sample Input

4

1 4 2 3

1

2 4

Sample Output

2

HINT

Source

By Autumn

题解：

出题人卡常数丧心病狂。。。

我原来思考过这个问题，还想什么时候出道题呢，结果已经被出出来了。。。

我们分块，然后预处理出i到j块的逆序对数f[i][j]，然后查询的时候不在块内的就暴力查询区间内<或>该数的个数，然后就可以n√n*logn了。。

然后出题人卡常数啊。。。

预处理的时候我偷懒用了主席树结果就T成翔啊。。。

无奈写了BIT。。。

代码：TLE

#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<string>
#define inf 1000000000
#define maxn 50000+5
#define maxm 3000000+5
#define maxk 1000+5
#define eps 1e-10
#define ll long long
#define pa pair<int,int>
#define for0(i,n) for(int i=0;i<=(n);i++)
#define for1(i,n) for(int i=1;i<=(n);i++)
#define for2(i,x,y) for(int i=(x);i<=(y);i++)
#define for3(i,x,y) for(int i=(x);i>=(y);i--)
#define mod 1000000007
#define l(x) ((x-1)*block+1)
#define r(x) ((x)==b
?n:(x)*block)
using namespace std;
inline int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
return x*f;
}
int n,m,tot,a[maxn],b[maxn],v[maxn],s[maxm],ss[maxn],ls[maxm],rs[maxm],rt[maxn],f[maxk][maxk];
inline bool cmp1(int x,int y){return a[x]<a[y];}
inline void update(int l,int r,int x,int &y,int z)
{
y=++tot;
s[y]=s[x]+1;
if(l==r)return;
ls[y]=ls[x];rs[y]=rs[x];
int mid=(l+r)>>1;
if(z<=mid)update(l,mid,ls[x],ls[y],z);else update(mid+1,r,rs[x],rs[y],z);
}
inline int query(int x,int y,int z)
{
int l=1,r=m,mid,xx=rt[x-1],yy=rt[y],res=0;
while(l!=r)
{
mid=(l+r)>>1;
if(z<=mid){xx=ls[xx];yy=ls[yy];r=mid;}
else {res+=s[ls[yy]]-s[ls[xx]];xx=rs[xx];yy=rs[yy];l=mid+1;}
}
return res;
}
inline int sum(int x){int t=0;for(;x;x-=x&(-x))t+=ss[x];return t;}
inline void add(int x,int y){for(;x<=m;x+=x&(-x))ss[x]+=y;}
int main()
{
freopen("input.txt","r",stdin);
freopen("output.txt","w",stdout);
n=read();
for1(i,n)a[i]=read(),b[i]=i;
sort(b+1,b+n+1,cmp1);
for1(i,n)
{
if(i==1||(a[b[i]]!=a[b[i-1]]))m++;
v[b[i]]=m;
}
m++;
for1(i,n)update(1,m,rt[i-1],rt[i],v[i]);
int block=sqrt(n);
for1(i,n)b[i]=(i-1)/block+1;
for1(i,b
)
{
int tmp=0;
for2(j,l(i),n)
{
add(v[j],1);
tmp+=j-l(i)+1-sum(v[j]);
f[i][b[j]]=tmp;
}
for2(j,l(i),n)add(v[j],-1);
}
//for1(i,b
)for1(j,b
)cout<<i<<' '<<j<<' '<<f[i][j]<<endl;
int cs=read(),ans=0;
while(cs--)
{
int x=read()^ans,y=read()^ans;if(x>y)swap(x,y);int bx=b[x],by=b[y];
if(by-bx<2)
{
ans=0;
for2(i,x,y-1)ans+=query(i,y,v[i]);
}
else
{
ans=f[bx+1][by-1];
for2(i,x,r(bx))ans+=query(i,r(by-1),v[i]);
for2(i,l(by),y)ans+=i-x+1-query(x,i,v[i]+1);
}
printf("%d\n",ans);
}
return 0;
}

View Code