首页 > 其他 > 详细

主席树 浅显学习

时间:2019-08-09 23:02:23      阅读:130      评论:0      收藏:0      [点我收藏+]

之前学习了不带修改的主席树,今天想学习一下待修改的主席树。在这个之前重新回顾了一下主席树,对主席树又有了新的认识。

主席树其实就是记录了每一个历史版本的树,所以这个空间要开的很大,一般开40倍就够了。

首先要建一棵空树,然后对我们要求的数进行离散化,最后就是每次更新新建一颗树。

学习地址

void build(int &id,int l,int r)
{
    id = ++tot;
    sum[id] = 0;
    int mid = (l + r) >> 1;
    if (l == r) return;
    build(lc[id], l, mid);
    build(rc[id], mid + 1, r);
}

 

每次新建树的过程中sum++ 这个很好理解,因为我们放了一个点进去,所以新树继承原来的信息之后要+1

这个就和权值线段树有点像了。这个是update操作

int update(int id,int l,int r,int val)
{
	int rt = ++tot;
	sum[rt] = sum[id] + 1;
	lc[rt] = lc[id], rc[rt] = rc[id];
	if (l == r) return rt;
	int mid = (l + r) >> 1;
	if (val <= mid) lc[rt] = update(lc[rt], l, mid, val);//如果我们要放进去的这个数比这个mid要小,就往小的这边更新
	else rc[rt] = update(rc[rt], mid + 1, r, val);//反之就是往大的这边更新。
	return rt;
}

  

我觉得比较难写的是query

int query(int u,int v,int l,int r,int k)//u 是第一个区间,v是第二个区间
{
	int mid = (l + r) >> 1;
	int x = sum[lc[u]] - sum[lc[v]];//sum 相减 表示这个区间里面的数有多少
	if (l == r) return l;//如果l==r到达叶子节点 就说明找到了
	if (k <= x) return query(lc[u], lc[v], l, mid, k);//如果第k大要比左区间的数量小,则说明在左区间
	return query(rc[u], rc[v], mid + 1, r, k - x);//如果在右区间,那么应该要减去左区间的个数
}//这个sum其实存的是数量,我们求第k大其实就是找这一段区间的第k个数是什么,如果在右区间,那么应该是右区间的 第 k-sum[]大

  

这个题目链接K-th Number

 

整体代码

技术分享图片
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <queue>
#include <vector>
#include <string>
#include <algorithm>
#include <iostream>
#include <map>
#define inf 0x3f3f3f3f
#define inf64 0x3f3f3f3f3f3f3f3f
using namespace std;
typedef long long ll;
const int maxn = 1e5 + 10;
int root[maxn * 40], lc[maxn * 40], rc[maxn * 40], sum[maxn * 40];
int tot = 0;

void build(int &id,int l,int r)
{
    id = ++tot;
    sum[id] = 0;
    int mid = (l + r) >> 1;
    if (l == r) return;
    build(lc[id], l, mid);
    build(rc[id], mid + 1, r);
}

int update(int id,int l,int r,int val)
{
    int rt = ++tot;
    sum[rt] = sum[id] + 1;
    lc[rt] = lc[id], rc[rt] = rc[id];
    if (l == r) return rt;
    int mid = (l + r) >> 1;
    if (val <= mid) lc[rt] = update(lc[rt], l, mid, val);//如果我们要放进去的这个数比这个mid要小,就往小的这边更新
    else rc[rt] = update(rc[rt], mid + 1, r, val);//反之就是往大的这边更新。
    return rt;
}

int query(int u,int v,int l,int r,int k)//u 是第一个区间,v是第二个区间
{
    int mid = (l + r) >> 1;
    int x = sum[lc[u]] - sum[lc[v]];//sum 相减 表示这个区间里面的数有多少
    if (l == r) return l;//如果l==r到达叶子节点 就说明找到了
    if (k <= x) return query(lc[u], lc[v], l, mid, k);//如果第k大要比左区间的数量小,则说明在左区间
    return query(rc[u], rc[v], mid + 1, r, k - x);//如果在右区间,那么应该要减去左区间的个数
}//这个sum其实存的是数量,我们求第k大其实就是找这一段区间的第k个数是什么,如果在右区间,那么应该是右区间的 第 k-sum[]大

int a[maxn], b[maxn];
int main()
{
    int n, m;
    scanf("%d%d", &n, &m);
    for (int i = 1; i <= n; i++) scanf("%d", &a[i]), b[i] = a[i];
    sort(b + 1, b + 1 + n);
    int len = unique(b + 1, b + 1 + n) - b - 1;
    build(root[0], 1, len);
    for (int i = 1; i <= n; i++) {
        a[i] = lower_bound(b + 1, b + 1 + len, a[i]) - b;
        root[i] = update(root[i - 1], 1, len, a[i]);
    }
    while (m--) {
        int l, r, k;
        scanf("%d%d%d", &l, &r, &k);
        int ans = query(root[r], root[l - 1], 1, len, k);
        printf("%d\n", b[ans]);
    }
    return 0;
}
View Code

 

主席树 浅显学习

原文:https://www.cnblogs.com/EchoZQN/p/11329684.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!