• [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
    ",tre[y].mmax);
    141       root=merge(merge(x,y),z);
    142     }
    143   }
    144   return 0;
    145 }
  • 相关阅读:
    前端路由hash、history原理及简单的实践下
    nginx缓存配置及开启gzip压缩
    koa2基于stream(流)进行文件上传和下载
    理解nodejs中的stream(流)
    理解http浏览器的协商缓存和强制缓存
    Nginx中的Rewrite的重定向配置与实践
    Redis入门---字符串类型
    Redis报错 Server started, Redis version 3.2.13 Can't handle RDB format version 9 Fatal error loading the DB: Invalid argument. Exiting.
    redis 基本原理及安装
    Go语言之高级篇beego框架之Controller
  • 原文地址:https://www.cnblogs.com/huangdalaofighting/p/8279912.html
Copyright © 2020-2023  润新知