ICPC20昆明M——动态开点权值线段树(主席树)

是ICPC19徐州H的子问题。ICPC19徐州H需要用带修主席树。ICPC19徐州H随便找的题解

需要发现一个结论(发现不了啊QAQ):考虑当前已知[1,s[x]]都能被表示,现在记s[x+1]为[L,R]中≤s[x]+1的所有数的和,若s[x+1]>s[x],则当前能表示的区间拓展为[1,s[x+1]],否则当前查询的答案就是s[x]+1。

证明:s[x+1]相比于上一轮新加了多个数,但只需要证明新加入1个数的情况,然后反复应用即可。设新加入的数t≤s[x]+1。则对于s[x]+i,1≤i≤t,0≤s[x]+i-t≤s[x]可以被表示,则加上t以后得知s[x]+i都可以被表示。

因此维护一个s,查询≤s+1的所有数的和t,若t=s则答案为s+1,否则令s=t继续循环。最坏情况大约是s倍增。

ans = 0;
while(true){
    LL t = ans;
    //min(ans+1,lim)防止nod[]访问越界
    ans = getsum(rt[r],rt[l-1],1,lim,1,min(ans+1,1LL*lim));
    if(ans == t) break;
}

那么我们需要一个数据结构来求a[L~R]≤x的数的和,这个数据结构就是动态开点的权值线段树。权值线段树以a[i]为下标,贡献也为a[i]。

值得注意的点:

  • 权值线段树需要的空间是n*log(max(a[i])),比nlogn大
  • modify和getsum的值域范围:参数r,应该取lim>=1e9
  • ans+1作为查询范围,可能超过1e9,需要用min(ans+1,lim)防止nod[]访问越界。可以这么截断是因为a[]的值域是1e9,如果截断成功了,ans必定没有增长。
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
#define rep(i,a,b) for(int i = (a);i <= (b);++i)
#define re_(i,a,b) for(int i = (a);i < (b);++i)
#define dwn(i,a,b) for(int i = (a);i >= (b);--i)

const int N = 1e6 + 5;
const int lim = 1e9 + 7;//主席树的区间上限

int n,m,a[N];
int rt[N],cnt = 0;

void dbg(){puts("");}
template<typename T, typename... R>void dbg(const T &f, const R &... r) {
    cout << f << " ";
    dbg(r...);
}
template<typename Type>inline void read(Type &xx){
    Type f = 1;char ch;xx = 0;
    for(ch = getchar();ch < '0' || ch > '9';ch = getchar()) if(ch == '-') f = -1;
    for(;ch >= '0' && ch <= '9';ch = getchar()) xx = xx * 10 + ch - '0';
    xx *= f;
}

struct Node{
    LL v;int l,r;
    Node(){v = l = r = 0;}
}nod[N * 40];//log(1e9)

void modify(int &ti,int tj,int l,int r,int idx,int v){
    nod[ti = ++cnt] = nod[tj];nod[ti].v += v;
    if(l == r) return;
    int mid = (l+r)>>1;
    if(idx <= mid) modify(nod[ti].l,nod[tj].l,l,mid,idx,v);
    else modify(nod[ti].r,nod[tj].r,mid+1,r,idx,v);
}

LL getsum(int ti,int tj,int l,int r,int ql,int qr){
    if(ql <= l && r <= qr) return nod[ti].v - nod[tj].v;
    int mid = (l+r)>>1;
    LL ret = 0;
    if(ql <= mid) ret += getsum(nod[ti].l,nod[tj].l,l,mid,ql,qr);
    if(qr > mid) ret += getsum(nod[ti].r,nod[tj].r,mid+1,r,ql,qr);
    return ret;
}

int main(int argc, char** argv) {
    read(n);read(m);
    rep(i,1,n) read(a[i]);
    rep(i,1,n){
        modify(rt[i],rt[i-1],1,lim,a[i],a[i]);
    }
    LL ans = 0;
    while(m--){
        int l,r;read(l);read(r);
        l = (l+ans) % n + 1;r = (r+ans) % n + 1;
        if(l > r) swap(l,r);
        ans = 0;
        while(true){
            LL t = ans;
            //min(ans+1,lim)防止nod[]访问越界
            ans = getsum(rt[r],rt[l-1],1,lim,1,min(ans+1,1LL*lim));
            if(ans == t) break;
        }
        ++ans;
        printf("%lld\n",ans);
    }
    return 0;
}
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值