-
Notifications
You must be signed in to change notification settings - Fork 110
/
Copy pathPersistent Segment Tree 2.cpp
55 lines (44 loc) · 1.12 KB
/
Persistent Segment Tree 2.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
const int MAXN = (1 << 20);
struct node
{
int sum;
node *l, *r;
node() { l = nullptr; r = nullptr; sum = 0; }
node(int x) { sum = x; l = nullptr; r = nullptr; }
};
typedef node* pnode;
pnode merge(pnode l, pnode r)
{
pnode ret = new node(0);
ret->sum = l->sum + r->sum;
ret->l = l;
ret->r = r;
return ret;
}
pnode init(int l, int r)
{
if(l == r) { return (new node(0));}
int mid = (l + r) >> 1;
return merge(init(l, mid), init(mid + 1, r));
}
pnode update(int pos, int val, int l, int r, pnode nd)
{
if(pos < l || pos > r) return nd;
if(l == r) { return (new node(val)); }
int mid = (l + r) >> 1;
return merge(update(pos, val, l, mid, nd->l), update(pos, val, mid + 1, r, nd->r));
}
int query(int qL, int qR, int l, int r, pnode nd)
{
if(qL <= l && r <= qR) return nd->sum;
if(qL > r || qR < l) return 0;
int mid = (l + r) >> 1;
return query(qL, qR, l, mid, nd->l) + query(qL, qR, mid + 1, r, nd->r);
}
int get_kth(int k, int l, int r, pnode nd)
{
if(l == r) return l;
int mid = (l + r) >> 1;
if(nd->l->sum < k) return get_kth(k - nd->l->sum, mid + 1, r, nd->r);
else return get_kth(k, l, mid, nd->l);
}