1008 - 主席树（可持久化线段树） - Poj 2104

传送门

 

题意

求区间第k大

 

分析

就是一个很主席树的板子

发现一个写的很良心的博客，但访问量出奇的少？？？不行不行，秉着应该有福同享的宗旨，我打算分享出来

大家一起学习（我真可爱O(∩_∩)O哈哈~）

快戳快戳

 

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>//poj可能有毒吧，居然不支持bits库
#define N 100009
#define in read()
using namespace std;
int n,m,a[N],b[N],rt[4*N];
inline int read(){}
int cnt=0,mm,tot=0;
struct node{int sum,lc,rc;}tr[20*N];
//注意数据范围，每次修改都会增加log（n）个新节点，有n次修改 
void insert(int pre,int &now,int l,int r,int val){
	tr[now=++tot]=tr[pre];
	tr[now].sum++;
	if(l==r) return;
	int mid=l+r>>1;
	if(val<=mid) insert(tr[pre].lc,tr[now].lc,l,mid,val);
	else insert(tr[pre].rc,tr[now].rc,mid+1,r,val);
}
int query(int x,int y,int l,int r,int k){
	if(l==r) return l;
	int mid=l+r>>1;
	int lsum=tr[tr[y].lc].sum-tr[tr[x].lc].sum;
	if(lsum>=k) return query(tr[x].lc,tr[y].lc,l,mid,k);
	else return query(tr[x].rc,tr[y].rc,mid+1,r,k-lsum);
}
int main(){
	n=in;m=in;
	int i,j,k;
	for(i=1;i<=n;++i) a[i]=b[++cnt]=in;
	sort(b+1,b+cnt+1);//离散化
	mm=unique(b+1,b+cnt+1)-b-1;
	for(i=1;i<=n;++i)
		a[i]=lower_bound(b+1,b+mm+1,a[i])-b;
	for(i=1;i<=n;++i)//每加入一个点就重新建一颗权值线段树
		insert(rt[i-1],rt[i],1,mm,a[i]);
	int l,r;
	for(i=1;i<=m;++i)
	{
		l=in;r=in;k=in;
		printf("%d\n",b[query(rt[l-1],rt[r],1,mm,k)]);
//由于之前是把a[i]对应到树上的，现在再映射回来，得到原数
	}	
	return 0;
}

听说还有可持久化trie树，嗯……最近应该会学