传送门

线段树好题。

维护区间取两种最值,区间加,求区间两种历史最值,区间最小值。

自己的写法调了一个晚上+一个上午+一个下午+一个晚上并没有调出来,90" role="presentation" style="position: relative;">9090分弃疗。

于是我开始学习新的写法。

也就是封装再封装。

然后发现自己的方法的确有不严谨的地方,但并不好改,于是写了一个凌晨参考了xuyixuan" role="presentation" style="position: relative;">xuyixuanxuyixuan神犇的写法调过了这题。

实际上就是我们对区间标记,区间最大值,区间最小值分别用三类结构体表示(顺便来个重载运算符)。然后就到了本蒟蒻出锅的地方:区间标记需要维护三类数的(我只维护了两类),最大值的,最小值的,既不是最大值也不是最小值的(这一项我并没有维护然后gg" role="presentation" style="position: relative;">gggg),然后分个类什么的再结合下面几道经典线段树应该就能做出来这道题了。

uoj#169. 【UR #11】元旦老人与数列

uoj#164. 【清华集训2015】V

bzoj3064: Tyvj 1518 CPU监控

bzoj4695: 最假女选手

hdu5306Gorgeous Sequence

代码:

#include<bits/stdc++.h>
#define lc (p<<1)
#define rc (p<<1|1)
#define mid (T[p].l+T[p].r>>1)
#define N 500005
#define inf 2e9
using namespace std;
inline int max(int a,int b){return a>b?a:b;}
inline int min(int a,int b){return a<b?a:b;}
inline int read(){
    int ans=0,w=1;
    char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')w=-1;ch=getchar();}
    while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
    return ans*w;
}
int n,m,a[N];
struct tag{int add,addmx,addmn;};
struct Max{int his,mx,sx;}mxtmp;
struct Min{int his,mn,sn;}mntmp;
struct Node{int l,r;Max mx;Min mn;tag Amax,Bmid,Cmin;}T[N<<2];
inline tag operator+(tag a,tag b){return (tag){a.add+b.add,max(a.addmx,b.addmx+a.add),min(a.addmn,b.addmn+a.add)};}
inline tag operator+(tag a,int b){return (tag){a.add+b,max(a.addmx,a.add+b),min(a.addmn,a.add+b)};}
inline Max operator+(Max a,Max b){return (Max){max(a.his,b.his),max(a.mx,b.mx),a.mx==b.mx?max(a.sx,b.sx):(a.mx<b.mx?max(a.mx,b.sx):max(a.sx,b.mx))};}
inline Min operator+(Min a,Min b){return (Min){min(a.his,b.his),min(a.mn,b.mn),a.mn==b.mn?min(a.sn,b.sn):(a.mn>b.mn?min(a.mn,b.sn):min(a.sn,b.mn))};}
inline void pushup(int p){T[p].mn=T[lc].mn+T[rc].mn,T[p].mx=T[lc].mx+T[rc].mx;}
inline void build(int p,int l,int r){
    T[p].l=l,T[p].r=r,T[p].Amax=T[p].Bmid=T[p].Cmin=(tag){0,0,0};
    if(l==r){T[p].mn=(Min){a[l],a[l],inf},T[p].mx=(Max){a[l],a[l],-inf};return;}
    build(lc,l,mid),build(rc,mid+1,r),pushup(p);
}
inline void pushadd(int p,int v){
    T[p].mx.mx+=v,T[p].mn.mn+=v;
    if(T[p].mx.sx!=-inf)T[p].mx.sx+=v;
    if(T[p].mn.sn!=inf)T[p].mn.sn+=v;
    T[p].mx.his=max(T[p].mx.his,T[p].mx.mx),T[p].mn.his=min(T[p].mn.his,T[p].mn.mn);
    T[p].Amax=T[p].Amax+v,T[p].Bmid=T[p].Bmid+v,T[p].Cmin=T[p].Cmin+v;
}
inline void pushAmax(int p,tag v){
    T[p].mx.his=max(T[p].mx.his,v.addmx+T[p].mx.mx);
    if(mxtmp.mx==mntmp.sn)T[p].mn.sn+=v.add;
    T[p].mx.mx+=v.add;
    T[p].Amax=T[p].Amax+v;
}
inline void pushBmid(int p,tag v){
    if(mxtmp.mx!=mntmp.sn&&mntmp.sn!=inf)T[p].mn.sn+=v.add;
    if(mntmp.mn!=mxtmp.sx&&mxtmp.sx!=-inf)T[p].mx.sx+=v.add;
    T[p].Bmid=T[p].Bmid+v;
}
inline void pushCmin(int p,tag v){
    T[p].mn.his=min(T[p].mn.his,v.addmn+T[p].mn.mn);
    if(mntmp.mn==mxtmp.sx)T[p].mx.sx+=v.add;
    T[p].mn.mn+=v.add;
    T[p].Cmin=T[p].Cmin+v;
}
inline void pushdown(int p){
    int chemn=min(T[lc].mn.mn,T[rc].mn.mn),chemx=max(T[lc].mx.mx,T[rc].mx.mx);
    mxtmp=T[lc].mx,mntmp=T[lc].mn;
    if(T[lc].mx.mx==chemx)pushAmax(lc,T[p].Amax);
    else if(T[lc].mx.mx==chemn)pushAmax(lc,T[p].Cmin);
    else pushAmax(lc,T[p].Bmid);
    pushBmid(lc,T[p].Bmid);
    if(T[lc].mn.mn==chemn)pushCmin(lc,T[p].Cmin);
    else if(T[lc].mn.mn==chemx)pushCmin(lc,T[p].Amax);
    else pushCmin(lc,T[p].Bmid);
    mxtmp=T[rc].mx,mntmp=T[rc].mn;
    if(T[rc].mx.mx==chemx)pushAmax(rc,T[p].Amax);
    else if(T[rc].mx.mx==chemn)pushAmax(rc,T[p].Cmin);
    else pushAmax(rc,T[p].Bmid);
    pushBmid(rc,T[p].Bmid);
    if(T[rc].mn.mn==chemn)pushCmin(rc,T[p].Cmin);
    else if(T[rc].mn.mn==chemx)pushCmin(rc,T[p].Amax);
    else pushCmin(rc,T[p].Bmid);
    T[p].Amax=T[p].Bmid=T[p].Cmin=(tag){0,0,0};
}
inline void update(int p,int ql,int qr,int v){
    if(ql>T[p].r||qr<T[p].l)return;
    if(ql<=T[p].l&&T[p].r<=qr)return pushadd(p,v);
    pushdown(p);
    if(qr<=mid)update(lc,ql,qr,v);
    else if(ql>mid)update(rc,ql,qr,v);
    else update(lc,ql,mid,v),update(rc,mid+1,qr,v);
    pushup(p);
}
inline void updmax(int p,int ql,int qr,int v){
    if(ql>T[p].r||qr<T[p].l||T[p].mx.mx<=v)return;
    if(ql<=T[p].l&&T[p].r<=qr&&T[p].mx.sx<v){
        mxtmp=T[p].mx,mntmp=T[p].mn;
        if(T[p].mx.mx==T[p].mn.mn)pushCmin(p,(tag){v-T[p].mx.mx,v-T[p].mx.mx,0});
        return pushAmax(p,(tag){v-T[p].mx.mx,v-T[p].mx.mx,0});
    }
    pushdown(p);
    if(qr<=mid)updmax(lc,ql,qr,v);
    else if(ql>mid)updmax(rc,ql,qr,v);
    else updmax(lc,ql,mid,v),updmax(rc,mid+1,qr,v);
    pushup(p);
}
inline void updmin(int p,int ql,int qr,int v){
    if(ql>T[p].r||qr<T[p].l||T[p].mn.mn>=v)return;
    if(ql<=T[p].l&&T[p].r<=qr&&T[p].mn.sn>v){
        mxtmp=T[p].mx,mntmp=T[p].mn;
        if(T[p].mx.mx==T[p].mn.mn)pushAmax(p,(tag){v-T[p].mn.mn,0,v-T[p].mn.mn});
        return pushCmin(p,(tag){v-T[p].mn.mn,0,v-T[p].mn.mn});
    }
    pushdown(p);
    if(qr<=mid)updmin(lc,ql,qr,v);
    else if(ql>mid)updmin(rc,ql,qr,v);
    else updmin(lc,ql,mid,v),updmin(rc,mid+1,qr,v);
    pushup(p);
}
inline int query(int p,int ql,int qr){
    if(ql>T[p].r||qr<T[p].l)return inf;
    if(ql<=T[p].l&&T[p].r<=qr)return T[p].mn.mn;
    pushdown(p);
    if(qr<=mid)return query(lc,ql,qr);
    if(ql>mid)return query(rc,ql,qr);
    return min(query(lc,ql,mid),query(rc,mid+1,qr));
}
inline int querymn(int p,int ql,int qr){
    if(ql>T[p].r||qr<T[p].l)return inf;
    if(ql<=T[p].l&&T[p].r<=qr)return T[p].mn.his;
    pushdown(p);
    if(qr<=mid)return querymn(lc,ql,qr);
    if(ql>mid)return querymn(rc,ql,qr);
    return min(querymn(lc,ql,mid),querymn(rc,mid+1,qr));
}
inline int querymx(int p,int ql,int qr){
    if(ql>T[p].r||qr<T[p].l)return -inf;
    if(ql<=T[p].l&&T[p].r<=qr)return T[p].mx.his;
    pushdown(p);
    if(qr<=mid)return querymx(lc,ql,qr);
    if(ql>mid)return querymx(rc,ql,qr);
    return max(querymx(lc,ql,mid),querymx(rc,mid+1,qr));
}
int main(){
    n=read(),m=read();
    for(int i=1;i<=n;++i)a[i]=read();
    build(1,1,n);
    while(m--){
        int op=read(),l=read(),r=read();
        if(op==1){int v=read();update(1,l,r,v);}
        if(op==2){int v=read();updmin(1,l,r,v);}
        if(op==3)printf("%d\n",query(1,l,r));
        if(op==4)printf("%d\n",querymn(1,l,r));
        if(op==5){int v=read();updmax(1,l,r,v);}
        if(op==6)printf("%d\n",querymx(1,l,r));
    }
    return 0;
}

2018.07.29~30 uoj#170. Picks loves segment tree VIII(线段树)的更多相关文章

  1. picks loves segment tree I

    picks loves segment tree I 题目背景 来源: \(\text {2018 WC Segment Tree Beats}\) 原作者: \(\text {C_SUNSHINE} ...

  2. hdu 5274 Dylans loves tree(LCA + 线段树)

    Dylans loves tree Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Othe ...

  3. UOJ#77. A+B Problem [可持久化线段树优化建边 最小割]

    UOJ#77. A+B Problem 题意:自己看 接触过线段树优化建图后思路不难想,细节要处理好 乱建图无果后想到最小割 白色和黑色只能选一个,割掉一个就行了 之前选白色必须额外割掉一个p[i], ...

  4. Codeforces 446C —— DZY Loves Fibonacci Numbers(线段树)

    题目:DZY Loves Fibonacci Numbers 题意比較简单,不解释了. 尽管官方的题解也是用线段树,但还利用了二次剩余. 可是我没有想到二次剩余,然后写了个感觉非常复杂度的线段树,还是 ...

  5. ACM/ICPC 2018亚洲区预选赛北京赛站网络赛 D 80 Days (线段树查询最小值)

    题目4 : 80 Days 时间限制:1000ms 单点时限:1000ms 内存限制:256MB 描述 80 Days is an interesting game based on Jules Ve ...

  6. 【2018.07.29】(深度优先搜索/回溯)学习DFS算法小记

    参考网站:https://blog.csdn.net/ldx19980108/article/details/76324307 这个网站里有动态图给我们体现BFS和DFS的区别:https://www ...

  7. 2018.07.30 cogs2632. [HZOI 2016] 数列操作d(线段树)

    传送门 线段树基本操作 区间加等差数列,维护区间和. 对于每个区间维护等差数列首项和公差,易证这两个东西都是可合并的,然后使用小学奥数的知识就可以切掉这题. 代码: #include<bits/ ...

  8. 2018.07.01 BZOJ3295: [Cqoi2011]动态逆序对(带修主席树)

    3295: [Cqoi2011]动态逆序对 **Time Limit: 10 Sec Memory Limit: 128 MB Description 对于序列A,它的逆序对数定义为满足i<j& ...

  9. 2018.07.01洛谷P2617 Dynamic Rankings(带修主席树)

    P2617 Dynamic Rankings 题目描述 给定一个含有n个数的序列a[1],a[2],a[3]--a[n],程序必须回答这样的询问:对于给定的i,j,k,在a[i],a[i+1],a[i ...

随机推荐

  1. 一个NPOI导出到excel文件的范例记录

    '使用NPOI写入新创建的excel文件,导出文件: Private Sub Sub_WriteXls() Dim XlsBook As XSSFWorkbook Dim XlsSheet As XS ...

  2. as3 程序域

    问题我要在应用程序中载入其他域的swf文件,并且允许它访问程序中的 ActionScript 解决办法使用flash.system.Security.allowDomain( ), flash.sys ...

  3. 键值集合List转换成datatable

    /// <summary> /// 键值集合List转换成datatable /// </summary> /// <param name="data" ...

  4. 转载:canal数据库同步redis

    ref: http://blog.csdn.net/tb3039450/article/details/53928351

  5. 拓扑排序获取所有可能序列JAVA实现

    在看算法基础这本书,看到有向无环图,其中介绍到了拓扑排序,讲到了获取拓扑序列的方法,结合自己的理解,用JAVA代码实现了获取所有可能序列,水平有限,效率什么的就没有考虑,下面贴上代码: package ...

  6. sql server字符串的类型

    sql server字符串分为不同的类型,下面就将为您介绍几个sql server字符串的常见类型,希望对您学习sql server字符串能有所帮助. nchar 此数据类型可存储1~4000个定长U ...

  7. Ansible playbooks

    Playbook是Ansible的配置,部署和编排语言. 他们可以描述您希望远程系统执行的策略,或一般IT流程中的一组步骤. 如果Ansible modules是您workshop的工具,则playb ...

  8. jps command not found已解决

    根据当前版本安装devel 包 eg: yum install java--openjdk-devel -y jdk小工具jps介绍 jps(Java Virtual Machine Process ...

  9. python使用外部PY文件的变量

    在用python和selenium编写登录等脚本时,一直都是给用户名和密码直接赋值.但是考虑到这样不便于管理,而且可能多个地方用到同一个变量,所以想把变量放在一个单独的文件中进行管理. 以登录脚本为例 ...

  10. suse安装gcc,升级到4.8.5

    前面这些是挂载iso,如果iso可以使用,就不需要下面几步. cd /etc/zypp/repos.d mkdir iso chmod -R 777 iso mount -o loop /media/ ...