SPOJ DQUERY (離線數狀陣列||線上主席樹)
轉載請註明出處,謝謝http://blog.csdn.net/ACM_cxlove?viewmode=contents by---cxlove
題目:給出一個序列,查詢區間內有多少個不同的樹
跟 著島娘,適妞一起學主席樹。。。
解法一:離線做法
將查詢區間按左端點排序
對於相同的數,先更新最左邊的位置
然後根據查詢區間,不斷更新next,保證查詢區間內
只更新過一個位置
維護字首和用樹狀陣列,時空效率都高
int s[N],a[N];
int n,q,ans[M];
int nxt[N];
map<int,int >mp;
struct Question{
int l,r,id;
bool operator<(const Question q)const{
return l<q.l;
}
}Q[M];
void add(int x,int val){
for(int i=x;i<=n;i+=lowbit(i))
s[i]+=val;
}
int query(int x){
int ret=0;
for(int i=x;i>0;i-=lowbit(i))
ret+=s[i];
return ret;
}
int main(){
while(scanf("%d",&n)!=EOF){
mp.clear();
for(int i=1;i<=n;i++){
scanf("%d",&a[i]);
if(mp.find(a[i])==mp.end()){
mp[a[i]]=i;
add(i,1);
}
}
mp.clear();
for(int i=n;i;i--){
if(mp.find(a[i])==mp.end())
nxt[i]=n+1;
else
nxt[i]=mp[a[i]];
mp[a[i]]=i;
}
scanf("%d",&q);
for(int i=0;i<q;i++){
scanf("%d%d",&Q[i].l,&Q[i].r);
Q[i].id=i;
}
sort(Q,Q+q);
int t=1;
for(int i=0;i<q;i++){
while(t<=n&&t<Q[i].l)
add(nxt[t++],1);
ans[Q[i].id]=query(Q[i].r)-query(Q[i].l-1);
}
for(int i=0;i<q;i++)
printf("%d\n",ans[i]);
}
return 0;
}
對於每一個起始位置,都建立一棵主席樹
從後往前更新,更新的時候就是在next位置先減掉,然後更新在當前位置。
主席樹又稱為可持久化線段樹,其本質上是儲存了所有的歷史資訊。
也可以理解為多棵線段樹,在更新的時候,充分共同利用歷史資訊
對於點更新來說,更新的總是某一棵子樹,另外一棵子樹,則指向原來的結點即可。
可是即使如此,主席樹還是需要大量的空間。
map<int,int>mp;
int a[N],tot,n,q;
int T[M],lson[M],rson[M],val[M];
int bulid(int l,int r){
int root=tot++;
val[root]=0;
int m=(l+r)>>1;
if(l!=r){
lson[root]=bulid(l,m);
rson[root]=bulid(m+1,r);
}
return root;
}
int update(int root,int pos,int v){
int newroot=tot++,tmp=newroot;
int l=1,r=n;
val[newroot]=val[root]+v;
while(l<r){
int m=(l+r)>>1;
//更新的時候需要充分利用歷史資訊
//更新原來的左子樹,右子樹不變
if(pos<=m){
lson[newroot]=tot++;rson[newroot]=rson[root];
newroot=lson[newroot];root=lson[root];
r=m;
}
//更新右子樹
else{
rson[newroot]=tot++;lson[newroot]=lson[root];
newroot=rson[newroot];root=rson[root];
l=m+1;
}
val[newroot]=val[root]+v;
}
return tmp;
}
int query(int root,int pos){
int ret=0;
int l=1,r=n;
while(pos<r){
int m=(l+r)>>1;
if(pos<=m){
r=m;
root=lson[root];
}
else{
ret+=val[lson[root]];
root=rson[root];
l=m+1;
}
}
return ret+val[root];
}
int main(){
while(scanf("%d",&n)!=EOF){
tot=0; //結點數
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
T[n+1]=bulid(1,n);
for(int i=n;i;i--){
int nxt;
map<int,int>::iterator it=mp.find(a[i]);
if(it==mp.end()) nxt=n+1;
else nxt=it->second;
//如果這是第一次出現,也就是最後一個位置上,則直接更新
if(nxt>n)
T[i]=update(T[i+1],i,1);
//在原來的位置上擦掉,在當前位置更新
else{
int t=update(T[i+1],nxt,-1);
T[i]=update(t,i,1);
}
mp[a[i]]=i;
}
scanf("%d",&q);
while(q--){
int l,r;
scanf("%d%d",&l,&r);
printf("%d\n",query(T[l],r));
}
}
return 0;
}
相關文章
- D-query SPOJ - DQUERY (主席樹)
- HDU 5862 Counting Intersections(樹狀陣列+掃描線+離散化)陣列
- HDU 1556 Color the ball(線段樹|樹狀陣列)陣列
- bzoj3110: [Zjoi2013]K大數查詢(主席樹+樹狀陣列)陣列
- 【主席數】可持續化線段樹
- 【Leetcode每日一題】327. 區間和的個數(線段樹/樹狀陣列)LeetCode每日一題陣列
- 樹狀陣列陣列
- bzoj1146: [CTSC2008]網路管理Network(dfs序+主席樹+樹狀陣列)陣列
- bzoj2588: Spoj 10628. Count on a tree(主席樹+LCA)
- 解析樹狀陣列陣列
- LeetCode C++ 劍指 Offer 51. 陣列中的逆序對【歸併排序/樹狀陣列/線段樹】LeetCodeC++陣列排序
- POJ 3468 【區間修改+區間查詢 樹狀陣列 | 線段樹 | 分塊】陣列
- 2024年暑假關於線段樹和樹狀陣列的小知識點陣列
- 浪潮伺服器陣列offine離線伺服器陣列
- 二維樹狀陣列陣列
- 樹狀陣列詳解陣列
- 樹狀陣列基礎陣列
- bzoj1803: Spoj1487 Query on a tree III(DFS序+主席樹)
- Leetcode 327. 區間和的個數 (字首和 + 離散化 + 樹狀陣列)LeetCode陣列
- 樹狀陣列模板題 & (樹狀陣列 1:單點修改,區間查詢)陣列
- 求區間不同數的個數【樹狀陣列求解】陣列
- 樹狀陣列快速入門陣列
- 樹狀陣列和逆序對陣列
- 【筆記/模板】樹狀陣列筆記陣列
- 數列求和【線段樹基礎】
- 浮點數(float)線上轉位元組陣列的工具陣列
- 位元組陣列轉浮點數線上轉換器陣列
- POJ-2352 Stars(樹狀陣列)陣列
- 樹狀陣列模板+習題集陣列
- 樹狀陣列(我是真小白)陣列
- 資料結構——樹狀陣列資料結構陣列
- 樹狀陣列3種基本操作陣列
- 演算法隨筆——主席樹(可持久化線段樹)演算法持久化
- 【主席樹】P3919 【模板】可持久化線段樹 1持久化
- float浮點數線上轉位元組陣列轉換工具陣列
- LeetCode 493. 翻轉對(歸併排序 || 離散化+樹狀陣列)LeetCode排序陣列
- HDU 1541 & POJ 2352 Stars (樹狀陣列)陣列
- bzoj3155: Preprefix sum(樹狀陣列)陣列
- 【luogu3368】模板 樹狀陣列 2陣列