洛谷 P3919 【模板】可持久化陣列(可持久化線段樹/平衡樹)-可持久化線段樹(單點更新,單點查詢)
P3919 【模板】可持久化陣列(可持久化線段樹/平衡樹)
題目背景
UPDATE : 最後一個點時間空間已經放大
標題即題意
有了可持久化陣列,便可以實現很多衍生的可持久化功能(例如:可持久化並查集)
題目描述
如題,你需要維護這樣的一個長度為 NN 的陣列,支援如下幾種操作
-
在某個歷史版本上修改某一個位置上的值
-
訪問某個歷史版本上的某一位置的值
此外,每進行一次操作(對於操作2,即為生成一個完全一樣的版本,不作任何改動),就會生成一個新的版本。版本編號即為當前操作的編號(從1開始編號,版本0表示初始狀態陣列)
輸入輸出格式
輸入格式:
輸入的第一行包含兩個正整數 N, MN,M, 分別表示陣列的長度和操作的個數。
第二行包含NN個整數,依次為初始狀態下陣列各位的值(依次為 a_iai,1 \leq i \leq N1≤i≤N)。
接下來MM行每行包含3或4個整數,代表兩種操作之一(ii為基於的歷史版本號):
-
對於操作1,格式為v_i \ 1 \ {loc}_i \ {value}_ivi 1 loci valuei,即為在版本v_ivi的基礎上,將 a_{{loc}_i}aloci 修改為 {value}_ivaluei
-
對於操作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}^31≤N,M≤103
對於50%的資料:1 \leq N, M \leq {10}^41≤N,M≤104
對於70%的資料:1 \leq N, M \leq {10}^51≤N,M≤105
對於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}^91≤N,M≤106,1≤loci≤N,0≤vi<i,−109≤ai,valuei≤109
經測試,正常常數的可持久化陣列可以通過,請各位放心
資料略微凶殘,請注意常數不要過大
另,此題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 }
沒了,菜到變形。