主席树 专题

最新推荐文章于 2021-09-17 22:49:33 发布
最新推荐文章于 2021-09-17 22:49:33 发布
阅读量646 收藏 1
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Miracle_ma/article/details/52695369
版权

poj2104
静态区间第k大，没有修改，所以时间是O(nlogn),空间也是O(nlogn)
模版，主席树可以做减法运算，区间第k大就是root[r]−root[l−1]

模板：

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#pragma comment(linker,"/STACK:102400000,102400000")

using namespace std;
#define   MAX           100010
#define   MAXN          500005
#define   maxnode       5
#define   sigma_size    30
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   lrt           rt<<1
#define   rrt           rt<<1|1
#define   middle        int m=(r+l)>>1
#define   LL            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   pii           pair<int,int>
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   limit         10000

//const int    prime = 999983;
const int    INF   = 0x3f3f3f3f;
const LL     INFF  = 0x3f3f;
const double pi    = acos(-1.0);
const double inf   = 1e18;
const double eps   = 1e-6;
const LL     mod   = 1e9+7;
const ull    mx    = 133333331;

/*****************************************************/
inline void RI(int &x) {
      char c;
      while((c=getchar())<'0' || c>'9');
      x=c-'0';
      while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
 }
/*****************************************************/

int a[MAX];
int b[MAX];
struct Node{
    int l,r,c;
}tr[MAX*20];
int root[MAX];
int sz;

int build(int l,int r){
    int k=sz++;
    if(l==r){
        tr[k].c=0;
        return k;
    }
    int mid=(l+r)/2;
    tr[k].l=build(l,mid);
    tr[k].r=build(mid+1,r);
    return k;
}

int update(int rt,int l,int r,int pos,int d){
    int k=sz++;
    tr[k].c=tr[rt].c+1;
    if(l==r) return k;
    int mid=(l+r)/2;
    if(pos<=mid){
        tr[k].r=tr[rt].r;
        tr[k].l=update(tr[rt].l,l,mid,pos,d);
    }
    else{
        tr[k].l=tr[rt].l;
        tr[k].r=update(tr[rt].r,mid+1,r,pos,d);
    }
    return k;
}

int query(int lr,int rr,int l,int r,int k){
    if(l==r) return l;
    int mid=(l+r)/2;
    if(tr[tr[rr].l].c-tr[tr[lr].l].c>=k) return query(tr[lr].l,tr[rr].l,l,mid,k);
    else return query(tr[lr].r,tr[rr].r,mid+1,r,k-(tr[tr[rr].l].c-tr[tr[lr].l].c));
}
int main(){
    //freopen("in.txt","r",stdin);
    int n,m;
    cin>>n>>m;
    for(int i=1;i<=n;i++) scanf("%d",&a[i]),b[i]=a[i];
    sort(b+1,b+n+1);
    int tot=unique(b+1,b+n+1)-b-1;
    sz=0;
    root[0]=build(1,tot);
    for(int i=1;i<=n;i++){
        int pos=lower_bound(b+1,b+tot+1,a[i])-b;
        root[i]=update(root[i-1],1,tot,pos,1);
    }
    while(m--){
        int l,r,k;
        scanf("%d%d%d",&l,&r,&k);
        printf("%d\n",b[query(root[l-1],root[r],1,tot,k)]);
    }
    return 0;
}

HDU4417
区间内比h小的数字个数
主席树维护[1,i]的前缀线段树，放在每个高度的位置上
查询直接在root[r]上比h小的减去root[l−1]上比h小的

代码:

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#pragma comment(linker,"/STACK:102400000,102400000")

using namespace std;
#define   MAX           100010
#define   MAXN          500005
#define   maxnode       5
#define   sigma_size    30
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   lrt           rt<<1
#define   rrt           rt<<1|1
#define   middle        int m=(r+l)>>1
#define   LL            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   pii           pair<int,int>
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   limit         10000

//const int    prime = 999983;
const int    INF   = 0x3f3f3f3f;
const LL     INFF  = 0x3f3f;
const double pi    = acos(-1.0);
const double inf   = 1e18;
const double eps   = 1e-6;
const LL     mod   = 1e9+7;
const ull    mx    = 133333331;

/*****************************************************/
inline void RI(int &x) {
      char c;
      while((c=getchar())<'0' || c>'9');
      x=c-'0';
      while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
 }
/*****************************************************/

int a[MAX];
int b[MAX];
struct Node{
    int l,r,c;
}tr[MAX*20];
int root[MAX];
int sz;

int build(int l,int r){
    int k=sz++;
    if(l==r){
        tr[k].c=0;
        return k;
    }
    int mid=(l+r)/2;
    tr[k].l=build(l,mid);
    tr[k].r=build(mid+1,r);
    return k;
}

int update(int rt,int l,int r,int pos,int d){
    int k=sz++;
    tr[k].c=tr[rt].c+1;
    if(l==r) return k;
    int mid=(l+r)/2;
    if(pos<=mid){
        tr[k].r=tr[rt].r;
        tr[k].l=update(tr[rt].l,l,mid,pos,d);
    }
    else{
        tr[k].l=tr[rt].l;
        tr[k].r=update(tr[rt].r,mid+1,r,pos,d);
    }
    return k;
}

int query(int lr,int rr,int l,int r,int L,int R){
    if(L<=l&&r<=R) return tr[rr].c-tr[lr].c;
    int mid=(l+r)/2;
    int ans=0;
    if(L<=mid) ans+=query(tr[lr].l,tr[rr].l,l,mid,L,R);
    if(R>mid) ans+=query(tr[lr].r,tr[rr].r,mid+1,r,L,R);
    return ans;
}
int main(){
    //freopen("in.txt","r",stdin);
    int t,kase=0;
    cin>>t;
    while(t--){
        int n,m;
        kase++;
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++) scanf("%d",&a[i]),b[i]=a[i];
        sort(b+1,b+n+1);
        int tot=unique(b+1,b+n+1)-b-1;
        root[0]=build(1,tot);
        for(int i=1;i<=n;i++){
            int pos=lower_bound(b+1,b+tot+1,a[i])-b;
            root[i]=update(root[i-1],1,tot,pos,1);
        }
        printf("Case %d:\n",kase);
        while(m--){
            int l,r,h;
            scanf("%d%d%d",&l,&r,&h);
            if(l>r) swap(l,r);
            l++;r++;
            int pos=upper_bound(b+1,b+tot+1,h)-b-1;
            if(pos<1) printf("0\n");
            else printf("%d\n",query(root[l-1],root[r],1,tot,1,pos));
        }
    }
    return 0;
}

SPOJDQUERY
区间不同种类的数字个数
模板题，对于[1,i]的前缀，维护主席树时，碰到a[i]如果出现过
在pre[a[i]]处−−，i处++，没出现过直接++
然后询问[l,r],就是在root[r]上[l,r]内求和

代码：

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <sstream>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#pragma comment(linker,"/STACK:102400000,102400000")

using namespace std;
#define   MAX           30010
#define   MAXN          1000005
#define   maxnode       5
#define   sigma_size    30
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   lrt           rt<<1
#define   rrt           rt<<1|1
#define   middle        int m=(r+l)>>1
#define   LL            long long
#define   ull           unsigned long long
#define   mem(x,v)      memset(x,v,sizeof(x))
#define   lowbit(x)     (x&-x)
#define   pii           pair<int,int>
#define   bits(a)       __builtin_popcount(a)
#define   mk            make_pair
#define   limit         10000

//const int    prime = 999983;
const int    INF   = 0x3f3f3f3f;
const LL     INFF  = 0x3f3f;
const double pi    = acos(-1.0);
const double inf   = 1e18;
const double eps   = 1e-6;
const LL     mod   = 1e9+7;
const ull    mx    = 133333331;

/*****************************************************/
inline void RI(int &x) {
      char c;
      while((c=getchar())<'0' || c>'9');
      x=c-'0';
      while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0';
 }
/*****************************************************/

int a[MAX];
int b[MAX];
struct Node{
    int l,r,c;
}tr[MAX*20];
int root[MAX];
int pre[MAX];
int sz;

int build(int l,int r){
    int k=sz++;
    if(l==r){
        tr[k].c=0;
        return k;
    }
    int mid=(l+r)/2;
    tr[k].l=build(l,mid);
    tr[k].r=build(mid+1,r);
    return k;
}

int update(int rt,int l,int r,int pos,int d){
    int k=sz++;
    tr[k].c=tr[rt].c+d;
    if(l==r) return k;
    int mid=(l+r)/2;
    if(pos<=mid){
        tr[k].r=tr[rt].r;
        tr[k].l=update(tr[rt].l,l,mid,pos,d);
    }
    else{
        tr[k].l=tr[rt].l;
        tr[k].r=update(tr[rt].r,mid+1,r,pos,d);
    }
    return k;
}

int query(int rt,int l,int r,int L,int R){
    if(L<=l&&r<=R) return tr[rt].c;
    int mid=(l+r)/2;
    int ans=0;
    if(L<=mid) ans+=query(tr[rt].l,l,mid,L,R);
    if(R>mid) ans+=query(tr[rt].r,mid+1,r,L,R);
    return ans;
}
int main(){
    //freopen("in.txt","r",stdin);
    int n;
    while(cin>>n){
        for(int i=1;i<=n;i++) RI(a[i]),b[i]=a[i];
        sort(b+1,b+n+1);
        mem(pre,-1);
        int tot=unique(b+1,b+n+1)-b-1;
        sz=0;
        root[0]=build(1,n);
        for(int i=1;i<=n;i++){
            int pos=lower_bound(b+1,b+tot+1,a[i])-b;
            if(pre[pos]==-1) root[i]=update(root[i-1],1,n,i,1);
            else {
                int tmp=update(root[i-1],1,n,pre[pos],-1);
                root[i]=update(tmp,1,n,i,1);
            }
            pre[pos]=i;
        }
        int m;
        cin>>m;
        while(m--){
            int l,r;
            RI(l);RI(r);
            printf("%d\n",query(root[r],1,n,l,r));
        }
    }
    return 0;
}
Miracle_ma
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
主席树.ppt
07-08
主席树.ppt
主席树专题
Hugokung的博客
10-04 121
题目描述 最近实验室正在为其管理的超级计算机编制一套任务管理系统，而你被安排完成其中的查询部分。超级计算机中的任务用三元组(Si,Ei,Pi)描述，(Si,Ei,Pi)表示任务从第Si秒开始，在第Ei秒后结束（第Si秒和Ei秒任务也在运行），其优先级为Pi。同一时间可能有多个任务同时执行，它们的优先级可能相同，也可能不同。调度系统会经常向查询系统询问，第Xi秒正在运行的任务中，优先级最小的Ki个...
主席树 常见题目类型总结(未完待续)
HeartFireY的博客
09-17 421
???? | Powered By HeartFireY 文章目录一、 主席树维护区间第K大二、主席树维护线段树区间修改(标记永久化)三、求区间内不同的数字个数/求区间大于KKK的数字有多少四、求区间小于等于KKK的数字个数(二分查询)五、求区间Mex六、求区间内出现次数大于>=k次的最前数七、主席树+树上路径 一、 主席树维护区间第K大 ???? | Require：主席树 主席树最基础的应用，对于给定的序列建立权值数组，对于每个元素建立一棵主席树，维护[1,n][1,n][1,n]的数字个数(
主席树入门题
某司机的黑车
06-01 407
主席树 全称是可持久化权值线段树（以前一直分不清可持久化线段树和主席树的区别） 但实际上写法跟可持久化trie没什么区别，维护权值就OK了 那就不讲了 jzoj1011 题目 Description JZ拥有一个很大的野生动物园。这个动物园坐落在一个狭长的山谷内，这个区域从南到北被划分成N个区域，每个区域都饲养着一头狮子。这些狮子从北到南编号为1,2,3,…,N。每头狮子都有...
主席树_主席树模板_
09-30
供信息学奥林匹克竞赛选手使用 主席树模板
权值线段树和主席树入门
09-25
权值线段树和主席树入门PPT，权值线段树，顾名思义就是记录权值的线段树，普通的线段树直接以坐标为l，r建树，而权值线段树是以大小来建树，树上寸的信息是该权值的数量，而通过建树时二分从小到大的性质，可以用这...
区间第k大 主席树代码
11-03
普通的主席树 区间第k大 代码 包含了简单的注释，和线段树很相似
工会主席述职报告专题8篇.docx
02-20
工会主席述职报告专题8篇.docx
主席树（可持久化线段树）入门专题
写东西的地方
08-16 5964
1.poj 2104 查询区间第k小。 主席树其实相当于建立了n棵线段树，第i棵线段树是根据区间【1，i】按值建立的。对于每一棵线段树我们记录它对应的区间每个数出现的次数，所以首先要对所有的数离散化。 先考虑最简单的情况，只查询【1，n】的第k小，对于【1，n】我们按值建立一棵线段树，对于a[i]我们在位置a[i]上加1。查询第k小那么先看左子区间出现了多少个数cnt，假设左区间出现的数cnt
（ 数据结构专题 ）【 主席树 】
石旭先森的博客
12-02 138
（ 数据结构专题 ）【 主席树 】 推荐阅读：https://blog.csdn.net/bestFy/article/details/78650360 推荐阅读：https://blog.csdn.net/qq_39565901/article/details/81782739 推荐阅读：https://blog.csdn.net/chenxiaoran666/article/detai...
最详细的讲解，让你一次学会主席树
热门推荐
creatorx的博客
07-19 3万+
好久以前就想学习主席树这个黑科技，一直觉得很难，然后平时上课也没有什么好的时间，所以一直搁置到现在，最近遇到了一个比较简单，比较经典的问题，求区间第k小，比如poj2104,没有更新操作，只有查询操作，因为最近一直在学习分块思想，既然没有更新操作，我觉得可以用分块搞一下，但是一直tle，我的大致思想是把原序列分为为若干块，然后对每一块进行块内排序，每次 查询操作就是二分区间第k小元素设为x,jud...
【例题】【主席树】
ANG
12-11 455
NKOJ2664 静态区间第K小数 Time Limit:40000MS Memory Limit:234567K Total Submit:56 Accepted:39 Case Time Limit:3000MSDescription 给定一个由N个数组成的序列{A1,A2,…,AN} 每次提问(l,r,k)表示提问序列中{Al,..,Ar}中第k小的数的值。 Input 第
主席树题目汇总
weixin_30780221的博客
11-16 90
bzoj 3932 [CQOI2015]任务查询系统 　√　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 洛谷 P2633 Count on a tree bzoj 1901Zju2112 Dynamic Rankings 动态第k大数 bzoj 3295 [Cqoi2011]动态逆序对 \sum| 转载于:ht...
主席树总结（题目合集）
常胜将军的博客
03-04 1655
1、HDU 2665 题意：无修改区间第k大 思路：主席树（离线算法）网上都有各种详细的解释了，就不多说了。。。主席树的核心思想包括前缀和、二分查找、空间重复利用、转化（区间表示在此范围内的数的个数，即权值线段树）。时间和空间复杂度为nlogn。 代码： #include<iostream> #include<algorithm> using namespace...
主席树学习--入门
weixin_30273501的博客
03-15 77
本篇文章，讲讲主席树入门，以及区间第K大，或者区间第K小。学习主席树前提是要学习线段树，主席树也就是多棵线段树，每次更新都需要记录之前的线段树，如果每次新建一棵线段树。比较浪费空间时间。因为每次更新一个点只会修改一条树链。其他的树节点信息不改变，所以可以公用一些节点。这就是可持久线段树思想。 大致思路我就不再重复了，整理了一些学习资料。算法讲堂-主席树，这个视频是UESTC某位大佬讲...
可持久化专题（一）——浅谈主席树：可持久化线段树
CXR的博客
08-08 6109
前言 不得不说，可持久化数据结构真是太难了！ 由于数据结构这东西真的太玄学了，学这个主席树我真的学了很久。 简介 主席树为什么叫主席树？据说因为它是一个名字缩写为HJTHJTHJT的神犇发明的，与当时主席的名字缩写一样…… 主席树实质上就是一棵可持久化线段树，它的具体实现可以看下面。 让我们从值域线段树开始说起 要学主席树，我们就要先学值域线段树。 值域线段树的区间存的并不...
正如标题，这是一道板子题。 给定一个长为n 的序列a 。 有m 次询问，求区间[l,r] 中出现次数第k 小的数，若有多个则输出最小的。 出现次数第k 小的数指恰有k-1 个数的出现次数严格小于该数。
最新发布
06-06
好的，这是一道比较典型的分治题目，可以使用分治+主席树来解决。具体的思路如下： 1.记录每个数出现的次数，并将所有数存入一个数组v中。 2.建立主席树，对于每个节点记录该区间内不同数字的数量以及数字出现的总...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值