设\(l[i]\)为i左边第一个比i大的数的下标。\(r[i]\)为i右边第一个比i大的数的下标。

我们把\(p1,p2\)分开考虑。

当产生贡献为\(p1\)时\(i\)和\(j\)一定满足,分别为\(l[x],r[x]\)枚举每一个值为\(i\),\(j\)之间最大值可证。

党产生贡献为\(p2\)时\(i\)和\(j\)满足分别为\(l[x],[x+1,r[x]-1]\)或\([l[x]+1,x-1],r[x]\),此时\(a[x]\)为\(i\),\(j\)之间最大值,\(i\),\(j\)一个比\(a[x]\)大,一个比\(a[x]\)小。

然后就把问题转化为二维数点问题。产生贡献的点对对应坐标系中的一个点(为了避免重复计数如\((r[x],l[x])\)和\((l[x],r[x])\),可以把小的作为横坐标,大的作为纵坐标)。然后我们每一次询问就是横坐标在\([l,r]\)之间,纵坐标在\([l,r]\)之间的权值和。

然后就可以用主席树做了。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<vector>
using namespace std;
#define int long long
const int N=201000;
struct line{ int l,r,w; };
vector<line> vec[N];
int tot,root[N],lazy[N*70],ch[N*70][2],sum[N*70];
int n,m,p1,p2,a[N],stack[N],top,l[N],r[N];
void add(int l,int r,int L,int R,int w,int pre,int &now){
now=++tot;
ch[now][0]=ch[pre][0];
ch[now][1]=ch[pre][1];
lazy[now]=lazy[pre];
sum[now]=sum[pre]+(R-L+1)*w;
if(l==L&&r==R){
lazy[now]+=w;
return;
}
int mid=(l+r)>>1;
if(L>mid)add(mid+1,r,L,R,w,ch[pre][1],ch[now][1]);
else if(R<=mid)add(l,mid,L,R,w,ch[pre][0],ch[now][0]);
else{
add(l,mid,L,mid,w,ch[pre][0],ch[now][0]);
add(mid+1,r,mid+1,R,w,ch[pre][1],ch[now][1]);
}
}
int check(int l,int r,int L,int R,int pre,int now){
if(l==L&&r==R)return sum[now]-sum[pre];
int mid=(l+r)>>1;
if(L>mid)return (lazy[now]-lazy[pre])*(R-L+1)+check(mid+1,r,L,R,ch[pre][1],ch[now][1]);
else if(R<=mid)return (lazy[now]-lazy[pre])*(R-L+1)+check(l,mid,L,R,ch[pre][0],ch[now][0]);
else return (lazy[now]-lazy[pre])*(R-L+1)
+check(l,mid,L,mid,ch[pre][0],ch[now][0])
+check(mid+1,r,mid+1,R,ch[pre][1],ch[now][1]);
}
int read(){
int sum=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){sum=sum*10+ch-'0';ch=getchar();}
return sum*f;
}
signed main(){
n=read(),m=read(),p1=read(),p2=read();
for(int i=1;i<=n;i++)a[i]=read();
for(int i=1;i<=n;i++){
while(a[i]>a[stack[top]]&&top)r[stack[top--]]=i;
stack[++top]=i;
}
while(top)r[stack[top--]]=n+1;
for(int i=n;i>=1;i--){
while(a[i]>a[stack[top]]&&top)l[stack[top--]]=i;
stack[++top]=i;
}
while(top)l[stack[top--]]=0;
for(int i=1;i<=n;i++){
line x;
if(i!=n&&i+1<=r[i]-1){
x.l=i+1;x.r=r[i]-1;x.w=p2;
vec[l[i]].push_back(x);
}
if(l[i]+1<=i-1&&i!=1){
x.l=l[i]+1;x.r=i-1;x.w=p2;
vec[r[i]].push_back(x);
}
x.l=i;x.r=i;x.w=p1;
vec[l[i]].push_back(x);
vec[r[i]].push_back(x);
}
for(int i=1;i<=n;i++){
root[i]=root[i-1];
for(int j=0;j<vec[i].size();j++)
add(1,n,vec[i][j].l,vec[i][j].r,vec[i][j].w,root[i],root[i]);
}
for(int i=1;i<=m;i++){
int l=read(),r=read();
printf("%lld\n",check(1,n,l,r,root[l-1],root[r]));
}
return 0;
}

[AH2017/HNOI2017]影魔(主席树+单调栈)的更多相关文章

  1. bzoj 4826: [Hnoi2017]影魔 [主席树 单调栈]

    4826: [Hnoi2017]影魔 题意:一个排列,点对\((i,j)\),\(p=max(i+1,j-1)\),若\(p<a_i,a_j\)贡献p1,若\(p\)在\(a_1,a_2\)之间 ...

  2. 【BZOJ3956】Count 主席树+单调栈

    [BZOJ3956]Count Description Input Output Sample Input 3 2 0 2 1 2 1 1 1 3 Sample Output 0 3 HINT M,N ...

  3. [BZOJ4826][HNOI2017]影魔(主席树)

    4826: [Hnoi2017]影魔 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 669  Solved: 384[Submit][Status][ ...

  4. 洛谷P3722 [AH2017/HNOI2017]影魔(线段树)

    题意 题目链接 Sol 题解好神仙啊qwq. 一般看到这种考虑最大值的贡献的题目不难想到单调数据结构 对于本题而言,我们可以预处理出每个位置左边第一个比他大的位置\(l_i\)以及右边第一个比他大的位 ...

  5. [AH2017/HNOI2017] 影魔 - 线段树

    #include<bits/stdc++.h> #define maxn 200010 using namespace std; int a[maxn],st[maxn][2],top,L ...

  6. Codeforces 781E Andryusha and Nervous Barriers 线段树 单调栈

    原文链接https://www.cnblogs.com/zhouzhendong/p/CF781E.html 题目传送门 - CF781E 题意 有一个矩形,宽为 w ,高为 h .一开始会有 w 个 ...

  7. 洛谷P4425 转盘 [HNOI/AHOI2018] 线段树+单调栈

    正解:线段树+单调栈 解题报告: 传送门! 1551又是一道灵巧连题意都麻油看懂的题,,,,所以先解释一下题意好了,,,, 给定一个n元环 可以从0时刻开始从任一位置出发 每次可以选择向前走一步或者在 ...

  8. 线段树+单调栈+前缀和--2019icpc南昌网络赛I

    线段树+单调栈+前缀和--2019icpc南昌网络赛I Alice has a magic array. She suggests that the value of a interval is eq ...

  9. 牛客多校第四场sequence C (线段树+单调栈)

    牛客多校第四场sequence C (线段树+单调栈) 传送门:https://ac.nowcoder.com/acm/contest/884/C 题意: 求一个$\max {1 \leq l \le ...

随机推荐

  1. struts 中数据处理的3中方式

    方式一: 获取servletapi中的对象 方式二: struts中封装的对象 方式三: 实现接口 方式一和方式二的区别 方式一需要额外引入包或者是方式二实现不了的功能,比如:获取url 因为方式二只 ...

  2. Python的流程控制

    条件判断 通过`if`,`elif`,`else`关键字来实现条件判断逻辑的实现,执行改结构中的其中一个,其结构如下: if condition1: pass elif condition2: pas ...

  3. Kz.layedit-layui.layedit富文本编辑器拓展

    项目介绍 首先欢迎使用 Kz.layedit!本项目基于layui.layedit富文本编辑器,在其之上拓展而来. 新增功能 html源码模式.插入hr水平线.段落格式.字体颜色.字体背景色.批量上传 ...

  4. NOI 2018 屠龙勇士 (拓展中国剩余定理excrt+拓展欧几里得exgcd)

    题目大意:略 真是一波三折的一道国赛题,先学了中国剩余定理,勉强看懂了模板然后写的这道题 把取出的宝剑攻击力设为T,可得Ti*x=ai(mod pi),这显然是ax=c(mod b)的形式 这部分用e ...

  5. SpringBoot 消费NSQ消息

    使用监听器,来实现实时消费nsq的消息 一.目前spring boot中支持的事件类型如下 ApplicationFailedEvent:该事件为spring boot启动失败时的操作 Applica ...

  6. 利用Python对文件进行批量重命名

    最近几天工作的内容是对40个项目进行考核,每个项目都需要一个考核评分表,已经有了项目的列表. 如果用常规的方法,需要复制40个文件,并逐个修改,不光工作量大,也容易出错,后期修改也不方便. 于是想到了 ...

  7. SPOJ 4491

    不妨先把所有要求的素数的对的个数写出来 f(2)=u(1)G(2)+u(2)*G(2*2)+u(3)*G(2*3)+.....u(k2)*G(2*k2) f(3)=u(1)G(3)+u(2)*G(2* ...

  8. 使用Modernizr检测支持CSS3

    使用Modernizr检测支持CSS3 如果支持某个属性,会增加一个class,名字就是该属性: 不支持,名字是no-某属性 还提供了一个全局Modernizr对象,使用如下: <script ...

  9. C++11新特性应用--介绍几个新增的便利算法(用于分区的几个算法)

    今天继续. C++11新增的关于Non-modifying sequence operations和Modifying sequence operations的算法已经写了.具体信息见之前的博客. 以 ...

  10. phonegap(cordova) 自己定义插件代码篇(六)----android ,iOS 微信支付工具整合

    还是那句话,在使用插件代码篇的时候,请先了解插件机制(如整合原生插件先阅读原生插件文档.非常重要.非常重要!非常重要!),如未了解,请先阅读入门篇.这里就专贴关键代码 必须先把官方sdk 依照要求一步 ...