P3834 【模板】可持久化線段樹 2

纯粹的發表於2024-08-10

原題連結

題解

總體上來講,就是二分 \(x\) 查詢 插入 \(l-1\) 時有多少數小於等於 \(x\),查詢插入 \(r\) 時有多少數小於等於 \(x\)
然後減一減,看看是不是 小於等於 \(k-1\)

我認為目前沒有比ai講的更清楚的了,請點選這裡

code

#include<bits/stdc++.h>
using namespace std;
/*
#define int long long
#define double long double
#define lowbit(x) ((x)&(-x))
const int inf=1e18;
const int mod=1e9+7;

const int N=4e5;
int qpow(int a,int n)
{
    int res=1;
    while(n)
    {
        if(n&1) res=res*a%mod;
        a=a*a%mod;
        n>>=1;
    }
    return res;
}
int inv(int x)
{
    return qpow(x,mod-2);
}
int fa[2000005];
int finds(int now){return now==fa[now]?now:finds(fa[now]);}

vector<int> G[200005];

int dfn[200005],low[200005];
int cnt=0,num=0;
int in_st[200005]={0};
stack<int> st;
int belong[200005]={0};

void scc(int now,int fa)
{
    dfn[now]=++cnt;
    low[now]=dfn[now];
    in_st[now]=1;
    st.push(now);

    for(auto next:G[now])
    {
        if(next==fa) continue;

        if(!dfn[next])
        {
            scc(next,now);
            low[now]=min(low[now],low[next]);
        }
        else if(in_st[next])
        {
            low[now]=min(low[now],dfn[next]);
        }
    }

    if(low[now]==dfn[now])
    {
        int x;
        num++;
        do
        {
            x=st.top();
            st.pop();
            in_st[x]=0;
            belong[x]=num;
        }while(x!=now);
    }
}
vector<int> prime;
bool mark[200005]={0};
void shai()
{
    for(int i=2;i<=200000;i++)
    {
        if(!mark[i]) prime.push_back(i);

        for(auto it:prime)
        {
            if(it*i>200000) break;

            mark[it*i]=1;
            if(it%i==0) break;
        }
    }
}
*/
const int N=200005;
const int LOG=25;
struct persist
{
    int left,right,count;
}tree[N*LOG];
int tot=0;
int a[N];
int root[N]={0};
int b[N];

int update(int pre,int l,int r,int x)
{
    int now=++tot;

    tree[now]=tree[pre];

    tree[now].count++;

    if(l==r) return now;

    int mid=(l+r)/2;
    if(mid<x) tree[now].right=update(tree[pre].right,mid+1,r,x);
    else tree[now].left=update(tree[pre].left,l,mid,x);

    return now;
}


int query(int node,int l,int r,int x)
{
    if(l>x) return 0;
    if(r<=x) return tree[node].count;

    int mid=(l+r)/2;

    return query(tree[node].right,mid+1,r,x)+query(tree[node].left,l,mid,x);
}
void solve()
{
    int n,m;
    cin>>n>>m;

    for(int i=1;i<=n;i++)
    {
        cin>>a[i];
        b[i]=a[i];
    }

    sort(b+1,b+1+n);
    int len=unique(b+1,b+1+n)-b-1;


    for(int i=1;i<=n;i++)
    {
        int pos=lower_bound(b+1,b+1+len,a[i])-b;
        root[i]=update(root[i-1],1,len,pos);
    }


    while(m--)
    {
        int x,y,k;

        cin>>x>>y>>k;

        int l=1,r=len+1;

        while(l+1<r)
        {
            int mid=(l+r)/2;
            if(query(root[y],1,len,mid-1)-query(root[x-1],1,len,mid-1)<=k-1) l=mid;
            else r=mid;
        }

        cout<<b[l]<<'\n';
    }
}
signed main()
{
    ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    int TT=1;
    //cin>>TT;
    while(TT--) solve();
    return 0;
}


相關文章