hdu_5193_Go to movies Ⅱ(带插入删除的逆序对,块状链表)

题目链接:hdu_5193_Go to movies Ⅱ

题意:

有n个人站成一排,每个人的身高为Hi。每次有人加入或者有人离开,就要判断有多少人站反了(i < j&&Hi>Hj) 
第一行n,m,接下来n个整数(n,m<=2e4) 
接下来m行, 
0 x y 表示有一个身高为y的人插在x后面,x=0表示插在最前面。(1≤y≤n) 
1 x 表示第x个人(从左到右)离开。

题解:

官方题解:

添加或者删除一个元素时,维护逆序对时,需要知道在它之前有多少个数比它大,在它之后有多少个数比他小。有下标和权值两个维度,可以使用两个数据结构嵌套。题目中n=20000,范围不大,外层可以使用分块维护下标,这样添加和删除元素的时候, 
也很方便,直接暴力。查找权值个数时,使用树状数组比较方便。内层通过树状数组维护权值。设每块的大小为S,那么删除或者添加元素时,维护逆序对数的复杂度是O(S+P∗logn),S是块内直接暴力更新逆序对的代价, 
​P/S∗logn在前面块找比它大和在后面块中找比它小的代价,P表示当前元素的个数。为了使这两部分复杂度尽量均摊让S=P/S∗logn,S取根号Plogn 
​。直接通过分块暴力添加和删除时,块的大小会退化,需要重构,。重构一次的复杂度是S*logn, 总的重构复杂度是,与添加和删除总的复杂度一至。因此整个问题的复杂度为O(m√nlogn).

 #include<bits/stdc++.h>
#define F(i,a,b) for(int i=a;i<=b;i++)
using namespace std; const int N=,sqr=;
int ans,tot,n,flag,op,x,y,m; inline void add(int *c,int x,int v){while(x<=n)c[x]+=v,x+=x&-x;} inline int sum(int *c,int x){int an=;while(x>)an+=c[x],x-=(x&-x);return an;} struct node
{
int dt[sqr+],sz,c[N],nxt,pre;
void init(){memset(c,,sizeof(c)),sz=,pre=,nxt=;}
}b[]; int getpos(int &x)
{
int i=;
while(b[i].sz<x&&b[i].nxt)x-=b[i].sz,i=b[i].nxt;
return i;
} int cal(int x)
{
int an=;
int id=getpos(x);
F(i,,x-)if(b[id].dt[i]>b[id].dt[x])an++;
F(i,x+,b[id].sz)if(b[id].dt[i]<b[id].dt[x])an++;
for(int i=;i!=id;i=b[i].nxt)
an+=sum(b[i].c,n)-sum(b[i].c,b[id].dt[x]);
for(int i=b[id].nxt;i!=;i=b[i].nxt)
an+=sum(b[i].c,b[id].dt[x]-);
return an;
} void del(int x){
ans-=cal(x);
int id=getpos(x);
add(b[id].c,b[id].dt[x],-);
F(i,x+,b[id].sz)b[id].dt[i-]=b[id].dt[i];
b[id].sz--;
if(b[id].sz==)//删除块
{
int pre=b[id].pre,net=b[id].nxt;
b[net].pre=pre,b[pre].nxt=net;
}
} void ins(int x,int y)
{
int xx=x;
if(flag==)
{
b[].init(),b[].dt[++b[].sz]=y;
add(b[].c,b[].dt[],),flag=;
return;
}
int id=getpos(x);//在第几个块
if(b[id].sz==sqr)//块分裂
{
b[++tot].init(),b[id].sz++;
int cnt=;
for(int i=b[id].sz;i>x;i--)
b[id].dt[i]=b[id].dt[i-];
b[id].dt[x]=y;
add(b[id].c,b[id].dt[x],);
int len=b[id].sz/;
F(i,len+,b[id].sz)
{
b[tot].dt[++cnt]=b[id].dt[i];
add(b[tot].c,b[id].dt[i],);
add(b[id].c,b[id].dt[i],-);
}
b[tot].sz=cnt,b[tot].nxt=b[id].nxt;
b[id].nxt=tot,b[id].sz=len,b[tot].pre=id;
}
else{
b[id].sz++;
for(int i=b[id].sz;i>x;i--)
b[id].dt[i]=b[id].dt[i-];
b[id].dt[x]=y;
add(b[id].c,b[id].dt[x],);
}
ans+=cal(xx);
} int main(){
while(~scanf("%d%d",&n,&m))
{
ans=,tot=,b[].sz=,flag=;
F(i,,n)scanf("%d",&y),ins(i,y);
F(i,,m)
{
scanf("%d",&op);
if(op==)scanf("%d%d",&x,&y),ins(x+,y);
else scanf("%d",&x),del(x);
printf("%d\n",ans);
}
}
return ;
}
上一篇:《Java程序设计》第一周学习记录(2)


下一篇:史上最全的Java命名规范[转]