[COGS 347]地震

Posted kakakakakaka

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[COGS 347]地震相关的知识,希望对你有一定的参考价值。

时间限制:4 s   内存限制:128 MB

问题描述

某国地形狭长,中部有一列山脉,由于多发地震,山脉在不断变化中。地震发生时,山脉有可能发生如下变化:局部海拔升高或降低,板块运动产生地裂而出现一段新的山脉,或板块挤压迫使一段山脉消失。

该国的科学考察队已经预测出了近期内将要发生的一次地震的全过程,他们得到的山脉变化的信息数据格式如下。

R a b c

a,b,c为三个整数,表示[a,b]这段山脉的海拔升高了c(或降低了-c)。

I a k t1 t2 ... tk

a,k,t1,t2,...,tk为整数,在第a个位置之后出现了一段新的山脉,长度为k,各处海拔依次为t1,t2,...,tk。

M a b

a,b为两个整数,表示[a,b]这段山脉消失,后面的山脉会移动过来。

查询请求格式为

Q a b

a,b为两个整数,查询[a,b]这段山脉的最高峰的海拔。

现在他们想知道这次地震中任意时刻任意一段山脉的最高峰的海拔,请你设计程序帮助他们。

输入格式

第一行,两个整数N,Q,表示地震前山脉的长度为N,地震中有Q个事件。

第二行,N个整数,表示初始时山脉各处的海拔。 接下来Q行,每行为一个山脉变化信息或查询请求,格式如上。

输出格式

对于每一个查询请求,输出一行,为该查询请求的结果。

样例输入

10 7
1 3 4 6 3 5 9 1 4 5
R 1 4 2
Q 3 7
I 1 2 2 3
M 6 9
Q 2 5
R 1 6 -4
Q 1 3

样例输出

9
6
-1

样例说明

初始时刻,山脉各处海拔为

1 3 4 6 3 5 9 1 4 5

经过 R 1 4 2,山脉各处海拔为

3 5 6 8 3 5 9 1 4 5

查询Q 3 7,结果为9

I 1 2 2 3,在1第个位置后插入了长度为2的山脉2 3,之后山脉各处海拔为

3 2 3 5 6 8 3 5 9 1 4 5

M 6 9后,山脉各处海拔为

3 2 3 5 6 1 4 5

查询Q 2 5,结果为6

经过 R 1 6 -4,山脉各处海拔为

-1 -2 -1 1 2 -3 4 5

查询Q 1 3,结果为-1

数据规模

  • 对于40%的数据,Q<=10000。
  • 对于70%的数据,1<=Q<=100000。
  • 对于100%的数据,1<=Q<=300000。

在任何时刻,山脉的长度在[1,100000]之内,山脉各处海拔在[-2^31,2^31]之内。初始山脉长度连同插入的山脉的总长度不超过1000000。在所有请求中,插入山脉(I)、压入地下(M)、查询(Q)、升降(R)、的比例约为2:3:3:4。

题解

fhq_treap,调了蛮久,最后的总结就是如果左右儿子都为空,那么不要对它push_down,因为这样会使编号为0的点被赋值。

然后空节点的最大值要赋值为-inf。

  1 #include<iostream>
  2 #include<cstdio>
  3 #include<cstring>
  4 #include<cstdlib>
  5 #include<cmath>
  6 #include<algorithm>
  7 #include<ctime>
  8 #include<queue>
  9 #define inf (2147483647)
 10 #define ll(x) tre[x].child[0]
 11 #define rr(x) tre[x].child[1]
 12 #define son(x,t) tre[x].child[t]
 13 using namespace std;
 14 int n,m,cnt,d[100005],root;
 15 int gi(){
 16   int ans=0,f=1;char i=getchar();
 17   while(i<0||i>9){if(i==-)f=-1;i=getchar();}
 18   while(i>=0&&i<=9){ans=ans*10+i-0;i=getchar();}
 19   return ans*f;
 20 }
 21 struct Treap{
 22   int child[2],x,mmax,lazy,size;
 23 }tre[100005];
 24 queue<int>mem;
 25 void add_lazy(int root,int c){
 26   tre[root].mmax+=c;
 27   tre[root].x+=c;
 28   tre[root].lazy+=c;
 29 }
 30 void push_up(int root){
 31   int l=ll(root),r=rr(root);
 32   tre[root].mmax=max(tre[root].x,max(tre[l].mmax,tre[r].mmax));
 33   tre[root].size=tre[l].size+tre[r].size+1;
 34 }
 35 void push_down(int root){
 36   if(!tre[root].lazy)return;
 37   if(ll(root))add_lazy(ll(root),tre[root].lazy);
 38   if(rr(root))add_lazy(rr(root),tre[root].lazy);
 39   tre[root].lazy=0;
 40 }
 41 int newnode(int x){
 42   int pos;
 43   if(!mem.empty())pos=mem.front(),mem.pop();
 44   else pos=++cnt;
 45   tre[pos].mmax=tre[pos].x=x;
 46   tre[pos].lazy=0;
 47   tre[pos].child[1]=tre[pos].child[0]=0;
 48   tre[pos].size=1;
 49   return pos;
 50 }
 51 void build(int &root,int left,int right){
 52   int mid=(left+right)>>1;
 53   root=newnode(d[mid]);
 54   if(left<mid)build(ll(root),left,mid-1);
 55   if(mid<right)build(rr(root),mid+1,right);
 56   push_up(root);
 57 }
 58 void split(int now,int k,int &x,int &y){
 59   if(!now)x=y=0;
 60   else{
 61     push_down(now);
 62     if(tre[ll(now)].size+1<=k){
 63       x=now;
 64       split(rr(now),k-tre[ll(now)].size-1,rr(now),y);
 65       push_up(x);
 66     }
 67     else{
 68       y=now;
 69       split(ll(now),k,x,ll(now));
 70       push_up(y);
 71     }
 72   }
 73 }
 74 int merge(int x,int y){
 75   if(!x||!y)return x+y;
 76   if(rand()&1){
 77     push_down(x);
 78     rr(x)=merge(rr(x),y);
 79     push_up(x);
 80     return x;
 81   }
 82   else{
 83     push_down(y);
 84     ll(y)=merge(x,ll(y));
 85     push_up(y);
 86     return y;
 87   }
 88 }
 89 void delet(int root){
 90   if(root){
 91     mem.push(root);
 92     delet(ll(root));
 93     delet(rr(root));
 94   }
 95 }
 96 int find(int root,int k){
 97   int y=tre[ll(root)].size;
 98   if(y+1==k)return root;
 99   else if(y>=k)return find(ll(root),k);
100   else return find(rr(root),k-y-1);
101 }
102 int main(){
103   freopen("equake.in","r",stdin);
104   freopen("equake.out","w",stdout);
105   srand(time(0));
106   int i,j,x,y,z,rt;
107   n=gi();m=gi();
108   for(i=1;i<=n;i++)
109     d[i]=gi();
110   build(root,1,n);
111   tre[0].mmax=-inf;
112   tre[0].x=-inf;
113   while(m--){
114     char s[2];
115     int a,b,c;
116     scanf("%s",s);a=gi();b=gi();
117     if(s[0]==R){
118       scanf("%d",&c);
119       split(root,b,x,z);
120       split(x,a-1,x,y);
121       add_lazy(y,c);
122       root=merge(merge(x,y),z);
123     }
124     if(s[0]==I){
125       split(root,a,x,y);
126       for(i=1;i<=b;i++)
127     d[i]=gi();
128       build(rt,1,b);
129       root=merge(merge(x,rt),y);
130     }
131     if(s[0]==M){
132       split(root,b,x,z);
133       split(x,a-1,x,y);
134       delet(y);
135       root=merge(x,z);
136     }
137     if(s[0]==Q){
138       split(root,b,x,z);
139       split(x,a-1,x,y);
140       printf("%d\n",tre[y].mmax);
141       root=merge(merge(x,y),z);
142     }
143   }
144   return 0;
145 }

 

以上是关于[COGS 347]地震的主要内容,如果未能解决你的问题,请参考以下文章

的海地地震数据代码

检索地震 JSON 结果时出现问题?

解题报告:hdu2191汶川地震 - 多重背包模板

数学建模基于matlab 地震勘探marmousi模型含Matlab源码 1977期

COGS2479 HZOI2016—偏序

代码随想录算法训练营第13天 | ● 239. 滑动窗口最大值 ● 347.前 K 个高频元素 ● 总结