洛谷 P3919 【模板】可持久化数组(可持久化线段树/平衡树)-可持久化线段树(单点更新,单点查询)

P3919 【模板】可持久化数组(可持久化线段树/平衡树)

题目背景

UPDATE : 最后一个点时间空间已经放大

标题即题意

有了可持久化数组,便可以实现很多衍生的可持久化功能(例如:可持久化并查集)

题目描述

如题,你需要维护这样的一个长度为 NN 的数组,支持如下几种操作

  1. 在某个历史版本上修改某一个位置上的值

  2. 访问某个历史版本上的某一位置的值

此外,每进行一次操作(对于操作2,即为生成一个完全一样的版本,不作任何改动),就会生成一个新的版本。版本编号即为当前操作的编号(从1开始编号,版本0表示初始状态数组)

输入输出格式

输入格式:

输入的第一行包含两个正整数 N, MN,M, 分别表示数组的长度和操作的个数。

第二行包含NN个整数,依次为初始状态下数组各位的值(依次为 a_iai1 \leq i \leq N1iN)。

接下来MM行每行包含3或4个整数,代表两种操作之一(ii为基于的历史版本号):

  1. 对于操作1,格式为v_i \ 1 \ {loc}_i \ {value}_ivi 1 loci valuei,即为在版本v_ivi的基础上,将 a_{{loc}_i}aloci 修改为 {value}_ivaluei

  2. 对于操作2,格式为v_i \ 2 \ {loc}_ivi 2 loci,即访问版本v_ivi中的 a_{{loc}_i}aloci的值

输出格式:

输出包含若干行,依次为每个操作2的结果。

输入输出样例

输入样例#1:  复制
5 10
59 46 14 87 41
0 2 1
0 1 1 14
0 1 1 57
0 1 1 88
4 2 4
0 2 5
0 2 4
4 2 1
2 2 2
1 1 5 91
输出样例#1:  复制
59
87
41
87
88
46

说明

数据规模:

对于30%的数据:1 \leq N, M \leq {10}^31N,M103

对于50%的数据:1 \leq N, M \leq {10}^41N,M104

对于70%的数据:1 \leq N, M \leq {10}^51N,M105

对于100%的数据:1 \leq N, M \leq {10}^6, 1 \leq {loc}_i \leq N, 0 \leq v_i < i, -{10}^9 \leq a_i, {value}_i \leq {10}^91N,M106,1lociN,0vi<i,109ai,valuei109

经测试,正常常数的可持久化数组可以通过,请各位放心

数据略微凶残,请注意常数不要过大

另,此题I/O量较大,如果实在TLE请注意I/O优化

询问生成的版本是指你访问的那个版本的复制

样例说明:

一共11个版本,编号从0-10,依次为:

* 0 : 59 46 14 87 41

* 1 : 59 46 14 87 41

* 2 : 14 46 14 87 41

* 3 : 57 46 14 87 41

* 4 : 88 46 14 87 41

* 5 : 88 46 14 87 41

* 6 : 59 46 14 87 41

* 7 : 59 46 14 87 41

* 8 : 88 46 14 87 41

* 9 : 14 46 14 87 41

* 10 : 59 46 14 87 91

题意看题目。

写的时候,数组范围开的2e5,最后两组数据过不了,一个WA,一个RE,最后发现数组开小了,开成2e6就可以了。

代码:

 1 //洛谷 P3919 【模板】可持久化数组(可持久化线段树/平衡树)
 2 //可持久化线段树-单点修改,单点查询
 3 #include<iostream>
 4 #include<cstdio>
 5 #include<cstring>
 6 #include<algorithm>
 7 #include<bitset>
 8 #include<cassert>
 9 #include<cctype>
10 #include<cmath>
11 #include<cstdlib>
12 #include<ctime>
13 #include<deque>
14 #include<iomanip>
15 #include<list>
16 #include<map>
17 #include<queue>
18 #include<set>
19 #include<stack>
20 #include<vector>
21 using namespace std;
22 typedef long long ll;
23 typedef pair<int,int> pii;
24 
25 const double PI=acos(-1.0);
26 const double eps=1e-6;
27 const ll mod=1e9+7;
28 const int inf=0x3f3f3f3f;
29 const int maxn=2e6+10;
30 const int maxm=100+10;
31 #define ios ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
32 #define lson l,m
33 #define rson m+1,r
34 
35 int val[maxn*40],ls[maxn*40],rs[maxn*40];
36 int n,m,sz=0;
37 
38 void build(int &rt,int l,int r)
39 {
40     rt=++sz;
41     if(l==r){
42         scanf("%d",&val[rt]);
43         return ;
44     }
45 
46     int m=(l+r)>>1;
47     build(ls[rt],lson);
48     build(rs[rt],rson);
49 }
50 
51 void update(int pre,int &rt,int l,int r,int p,int c)
52 {
53     rt=++sz;
54     ls[rt]=ls[pre];rs[rt]=rs[pre];
55     if(l==r){
56         val[rt]=c;
57         return ;
58     }
59 
60     int m=(l+r)>>1;
61     if(p<=m) update(ls[pre],ls[rt],lson,p,c);
62     else     update(rs[pre],rs[rt],rson,p,c);
63 }
64 
65 void query(int rt,int l,int r,int p)
66 {
67     if(l==r){
68         printf("%d\n",val[rt]);
69         return ;
70     }
71     int m=(l+r)>>1;
72     if(p<=m) query(ls[rt],lson,p);
73     else     query(rs[rt],rson,p);
74 }
75 
76 int rt[maxn];
77 
78 int main()
79 {
80     int n,m;
81     scanf("%d%d",&n,&m);
82     sz=0;
83     build(rt[0],1,n);
84     for(int i=1; i<=m; i++)
85     {
86         int t,op,p,v;
87         scanf("%d%d",&t,&op);
88         if(op==1){
89             scanf("%d%d",&p,&v);
90             update(rt[t],rt[i],1,n,p,v);
91         }
92         else{
93             scanf("%d",&p);
94             query(rt[t],1,n,p);
95             rt[i]=rt[t];
96         }
97     }
98 }

没了,菜到变形。

猜你喜欢

转载自www.cnblogs.com/ZERO-/p/9822956.html