[BZOJ1901][luogu2617]Dynamic Rankings(树状数组+主席树)

题面

单点修改,区间求第k大

分析

首先,这道题卡权值线段树套treap的做法,所以只能用主席树做

对于静态的查询,root[i]对应的主席树的区间[l,r]保存的是a[1]~a[i]有多少个值落在区间[l,r]内。如果我们要修改a[i],则要修改O(n)棵主席树。那么我们像树状数组那样维护n棵主席树,不同的是每棵主席树里保存的是,a[i-lowbit(i)+1]~a[i]有多少个值落在区间[l,r]内

对于查询[ql,qr]时的做差
我们要像树状数组求和那样,把root[i],root[i-lowbit(i)],....共O(logn)棵主席树的值加起来,才能得到a[1]~a[i]有多少个数落在1~qr里面,1~(ql-1)同理,维护两个临时数组x,y存储这O(logn)棵主席树的根即可

对于单点修改,我们像树状数组修改那样,往O(logn)棵主席树中插入即可

时间复杂度\(O(n \log ^2 n)\),空间复杂度略大,为\(O(n\log^2n)\)

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<ctime>
#include<algorithm>
#define maxn 200005
#define maxlogn 27
using namespace std;
inline void qread(int& x) {
    x=0;
    char c=getchar();
    int sign=1;
    while(c<'0'||c>'9') {
        if(c=='-') sign=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9') {
        x=x*10+c-'0';
        c=getchar();
    }
    x*=sign;
}
inline void qread(char &x){
    char c=getchar();
    while(c<'A'||c>'Z'){
        c=getchar();
    }
    x=c; 
}
inline void qprint(int x) {
    if(x<0) {
        putchar('-');
        qprint(-x);
    } else if(x==0) {
        putchar('0');
        return;
    } else {
        if(x/10>0) qprint(x/10);
        putchar('0'+x%10);
    }
}

int n,m;
struct node {
    int ls;
    int rs;
    int sum;
} tree[maxn*maxlogn*maxlogn];
int root[maxn*maxlogn];
int ptr;
void push_up(int x) {
    tree[x].sum=tree[tree[x].ls].sum+tree[tree[x].rs].sum;
}
void update(int &x,int upos,int uval,int last,int l,int r) {
    x=++ptr;
    tree[x]=tree[last];
    if(l==r) {
        tree[x].sum+=uval;
        return;
    }
    int mid=(l+r)>>1;
    if(upos<=mid) update(tree[x].ls,upos,uval,tree[last].ls,l,mid);
    else update(tree[x].rs,upos,uval,tree[last].rs,mid+1,r);
    push_up(x);
}
int totx,x[maxn*maxlogn],toty,y[maxn*maxlogn];
int query(int k,int l,int r) {
    if(l==r) return l;
    int lcnt=0;
    int mid=(l+r)>>1;
    for(int i=1;i<=totx;i++){
        lcnt-=tree[tree[x[i]].ls].sum;
    }
    for(int i=1;i<=toty;i++){
        lcnt+=tree[tree[y[i]].ls].sum;
    } 
    if(k<=lcnt){
        for(int i=1;i<=totx;i++) x[i]=tree[x[i]].ls;
        for(int i=1;i<=toty;i++) y[i]=tree[y[i]].ls;
        return query(k,l,mid);
    }else{
        for(int i=1;i<=totx;i++) x[i]=tree[x[i]].rs;
        for(int i=1;i<=toty;i++) y[i]=tree[y[i]].rs;
        return query(k-lcnt,mid+1,r);
    }
}
inline int lowbit(int x){
    return x&-x;
}


struct oper {
    int type;
    int l,r,k;
    int pos,val;
    int dval;
} q[maxn];
int a[maxn];
int sz;
int b[maxn*2];

void add(int xpos,int xval,int v){
    for(int i=xpos;i<=sz;i+=lowbit(i)){
        update(root[i],xval,v,root[i],1,sz); 
    }
}
int answer(int l,int r,int k){
    totx=toty=0;
    for(int i=l-1;i;i-=lowbit(i)) x[++totx]=root[i];
    for(int i=r;i;i-=lowbit(i)) y[++toty]=root[i];
    return query(k,1,sz);
}
int main() {
    char cmd;
    qread(n);
    qread(m);
    for(int i=1; i<=n; i++) {
        qread(a[i]);
        b[++sz]=a[i];
    }
    for(int i=1; i<=m; i++) {
        qread(cmd);
        if(cmd=='Q') {
            q[i].type=0;
            qread(q[i].l);
            qread(q[i].r);
            qread(q[i].k);
        } else {
            q[i].type=1;
            qread(q[i].pos);
            qread(q[i].val);
            b[++sz]=q[i].val;
        }
    }
    sort(b+1,b+1+sz);
    sz=unique(b+1,b+1+sz)-b-1;
    for(int i=1; i<=n; i++) a[i]=lower_bound(b+1,b+1+sz,a[i])-b;
    for(int i=1; i<=m; i++) {
        if(q[i].type==1) {
            q[i].dval=lower_bound(b+1,b+1+sz,q[i].val)-b;
        }
    }
    for(int i=1; i<=n; i++) {
        add(i,a[i],1);
    }
//  T2.debug(1);
    for(int i=1; i<=m; i++) {
        if(q[i].type==0) {
            qprint(b[answer(q[i].l,q[i].r,q[i].k)]);
            putchar('\n');
        } else {
            add(q[i].pos,a[q[i].pos],-1);
            add(q[i].pos,q[i].dval,1);
            a[q[i].pos]=q[i].dval;

//          T2.debug(1);
        }
    }
}

转载于:https://www.cnblogs.com/birchtree/p/10851759.html

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
SQLAlchemy 是一个 SQL 工具包和对象关系映射(ORM)库,用于 Python 编程语言。它提供了一个高级的 SQL 工具和对象关系映射工具,允许开发者以 Python 类和对象的形式操作数据库,而无需编写大量的 SQL 语句。SQLAlchemy 建立在 DBAPI 之上,支持多种数据库后端,如 SQLite, MySQL, PostgreSQL 等。 SQLAlchemy 的核心功能: 对象关系映射(ORM): SQLAlchemy 允许开发者使用 Python 类来表示数据库表,使用类的实例表示表中的行。 开发者可以定义类之间的关系(如一对多、多对多),SQLAlchemy 会自动处理这些关系在数据库中的映射。 通过 ORM,开发者可以像操作 Python 对象一样操作数据库,这大大简化了数据库操作的复杂性。 表达式语言: SQLAlchemy 提供了一个丰富的 SQL 表达式语言,允许开发者以 Python 表达式的方式编写复杂的 SQL 查询。 表达式语言提供了对 SQL 语句的灵活控制,同时保持了代码的可读性和可维护性。 数据库引擎和连接池: SQLAlchemy 支持多种数据库后端,并且为每种后端提供了对应的数据库引擎。 它还提供了连接池管理功能,以优化数据库连接的创建、使用和释放。 会话管理: SQLAlchemy 使用会话(Session)来管理对象的持久化状态。 会话提供了一个工作单元(unit of work)和身份映射(identity map)的概念,使得对象的状态管理和查询更加高效。 事件系统: SQLAlchemy 提供了一个事件系统,允许开发者在 ORM 的各个生命周期阶段插入自定义的钩子函数。 这使得开发者可以在对象加载、修改、删除等操作时执行额外的逻辑。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值