【BZOJ 3339】Rmq Problem:线段树优化DP_bzoj dp 线段树优化-优快云博客

本文介绍了一种数据结构——主席树，并通过一个具体问题展示了如何使用主席树进行插入操作及查询。文章详细解释了主席树的基本概念、构建过程、插入算法和查询方法。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目描述

在这里插入图片描述

分析

用这道题来复习一下主席树
我们在主席树中插入每一个数的位置，没有被插入的数默认位置为0
然后我们只需要在第 $r$ 棵主席树中找到第一个位置小于 $l$ 的数即可

代码

#pragma GCC optimize(3)
#include <bits/stdc++.h>
#define debug(x) cout<<#x<<":"<<x<<endl;
#define dl(x) printf("%lld\n",x);
#define di(x) printf("%d\n",x);
#define _CRT_SECURE_NO_WARNINGS
#define pb push_back
#define mp make_pair
#define all(x) (x).begin(),(x).end()
#define fi first
#define se second
#define SZ(x) ((int)(x).size())
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;
typedef vector<int> VI;
const int INF = 0x3f3f3f3f;
const int N = 200010;
const ll mod = 1000000007;
const double eps = 1e-9;
const double PI = acos(-1);
template<typename T>inline void read(T &a) {
    char c = getchar(); T x = 0, f = 1; while (!isdigit(c)) {if (c == '-')f = -1; c = getchar();}
    while (isdigit(c)) {x = (x << 1) + (x << 3) + c - '0'; c = getchar();} a = f * x;
}
int gcd(int a, int b) {return (b > 0) ? gcd(b, a % b) : a;}
struct Node{
    int l,r;
    int cnt;
}tr [N * 40];
int root[N];
int idx,n,m;

int build(int l,int r){
    int p = ++idx;
    if(l == r) return p;
    int mid = l + r >> 1;
    tr[p].l = build(l,mid),tr[p].r = build(mid + 1,r);
    return p;
}
int insert(int p,int l,int r,int x,int s){
    int q = ++idx;
    tr[q] = tr[p];
    if(l == r){
        tr[q].cnt = s;
        return q;
    }
    int mid = l + r >> 1;
    if(x <= mid) tr[q].l = insert(tr[p].l,l,mid,x,s);
    else tr[q].r = insert(tr[p].r,mid + 1,r,x,s);
    tr[q].cnt = min(tr[tr[q].l].cnt,tr[tr[q].r].cnt);
    return q;
}

int query(int p,int l,int r,int s){
    if(l == r) return l;
    int mid = l + r >> 1;
    if(tr[tr[p].l].cnt < s) return query(tr[p].l,l,mid,s);
    else return query(tr[p].r,mid + 1,r,s);
}

int main() {
    read(n),read(m);
    root[0] = build(1,N - 1);
    for(int i = 1;i <= n;i++){
        int x;
        read(x);
        x++;
        root[i] = insert(root[i - 1],1,N - 1,x,i);
    }
    while(m--){
        int l,r;
        read(l),read(r);
        di(query(root[r],1,N - 1,l) - 1);
    }
    return 0;
}

/**
*　　┏┓　　　┏┓+ +
*　┏┛┻━━━┛┻┓ + +
*　┃　　　　　　　┃
*　┃　　　━　　　┃ ++ + + +
*  ████━████+
*  ◥██◤　◥██◤ +
*　┃　　　┻　　　┃
*　┃　　　　　　　┃ + +
*　┗━┓　　　┏━┛
*　　　┃　　　┃ + + + +Code is far away from 　
*　　　┃　　　┃ + bug with the animal protecting
*　　　┃　 　 ┗━━━┓ 神兽保佑,代码无bug　
*　　　┃ 　　　　　　 ┣┓
*　　  ┃ 　　　　　 　┏┛
*　    ┗┓┓┏━┳┓┏┛ + + + +
*　　　　┃┫┫　┃┫┫
*　　　　┗┻┛　┗┻┛+ + + +
*/