HDU3727 - Jewel(主席树)

题目大意

对一个序列进行以下四种操作：

1、Insert x 在序列尾部插入x

2、Query_1 s t k 查询区间[s,t]的第k小

3、Query_2 x 查询x的在序列中排名

4、Query_3 k 查询序列中的第k小

题解

第一个和第二个显然是主席树能够实现的功能，第三四个操作用vector+lower_bound就可以实现，妈蛋，写完交上去一直WA，fuck,然后把存储操作类型的每个字符串的大小从10改成15居然变RE了！！！才发现用来存储操作的数组开小了，只开到100000。。。改完就A了。。。3687MS，好慢。。。假设插入操作为n，其他操作位m,时间复杂度是nlogn+mlogn 应该没有这么慢的说啊，莫非是lower_bound调用次数太多，所以耗时很多？

代码：

#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;
#define MAXN 100005
#define lson l,m,ls[s]
#define rson m+1,r,rs[s]
typedef long long LL;
int ls[20*MAXN],rs[20*MAXN];
int cnt[20*MAXN],T[MAXN];
int a[MAXN],num[MAXN];
char op[2*MAXN][15];
int L[2*MAXN],R[2*MAXN],K[2*MAXN];
int tot;
vector<int>ivec;
void build(int l,int r,int &s)
{
s=++tot;
cnt[s]=0;
if(l==r) return;
int m=(l+r)>>1;
build(lson);
build(rson);
}
void update(int last,int p,int l,int r,int &s)
{
s=++tot;
ls[s]=ls[last],rs[s]=rs[last];
cnt[s]=cnt[last]+1;
if(l==r) return;
int m=(l+r)>>1;
if(p<=m) update(ls[last],p,lson);
else     update(rs[last],p,rson);
}
int query(int ss,int tt,int l,int r,int k)
{
if(l==r) return r;
int sum=cnt[ls[tt]]-cnt[ls[ss]];
int m=(l+r)>>1;
if(sum>=k) return query(ls[ss],ls[tt],l,m,k);
else       return query(rs[ss],rs[tt],m+1,r,k-sum);
}
int main()
{
int m,kase=0;
while(scanf("%d",&m)!=EOF)
{
memset(cnt,0,sizeof(cnt));
int len=0,n;
LL qy1=0,qy2=0,qy3=0;
ivec.clear();
for(int i=0; i<m; i++)
{
scanf("%s%d",op[i],&L[i]);
if(!strcmp(op[i],"Insert"))
{
++len;
num[len]=a[len]=L[i];
}
else if(!strcmp(op[i],"Query_1"))
scanf("%d%d",&R[i],&K[i]);
}
sort(a+1,a+len+1);
n=len;
len=unique(a+1,a+len+1)-a-1;
for(int i=1; i<=n; i++) num[i]=lower_bound(a+1,a+len+1,num[i])-a;
tot=0;
build(1,len,T[0]);
for(int i=1; i<=n; i++) update(T[i-1],num[i],1,len,T[i]);
for(int i=0; i<m; i++)
{
if(!strcmp(op[i],"Insert"))
ivec.insert(lower_bound(ivec.begin(),ivec.end(),L[i]),L[i]);
else if(!strcmp(op[i],"Query_1"))
qy1+=a[query(T[L[i]-1],T[R[i]],1,len,K[i])];
else if(!strcmp(op[i],"Query_2"))
qy2+=lower_bound(ivec.begin(),ivec.end(),L[i])-ivec.begin()+1;
else
qy3+=ivec[L[i]-1];
}
printf("Case %d:\n",++kase);
printf("%I64d\n%I64d\n%I64d\n",qy1,qy2,qy3);
}
return 0;
}