HNOI2017影魔
影魔
- 这么简单的方法尽然想不到,我是真的菜
- 对每个点,用单调栈的方式处理出他左右第一个比他大的数的位置,你可以把\(0\)和\(n+1\)设成\(inf\)。
- 显然对于每对\(lef[i]\)和\(rig[i]\)都会做出\(p1\)的贡献
- 每个\(lef[i]\)会对\(i+1\)到\(rig[i]-1\)做出\(p2\)贡献
- 同理,每个\(rig[i]\)都会给\(lef[i]+1\)到\(i-1\)做出\(p2\)贡献
- 用结构体存下来,按顺序用线段树将贡献加入即可
- 统计贡献,对于每个询问\(l-r\)在扫到\(l-1\)时将这段区间减去,在扫到\(r\)时再将这段区间贡献加上即可
// luogu-judger-enable-o2
#include<bits/stdc++.h>
using namespace std;
typedef int sign;
typedef long long ll;
#define For(i,a,b) for(register sign i=(sign)a;i<=(sign)b;++i)
#define Fordown(i,a,b) for(register sign i=(sign)a;i>=(sign)b;--i)
const int N=2e5+5;
bool cmax(sign &a,sign b){return (a<b)?a=b,1:0;}
bool cmin(sign &a,sign b){return (a>b)?a=b,1:0;}
template<typename T>inline T read()
{
T f=1,ans=0;
char ch=getchar();
while(!isdigit(ch)&&ch!='-')ch=getchar();
if(ch=='-')f=-1,ch=getchar();
while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch-'0'),ch=getchar();
return ans*f;
}
template<typename T>inline void write(T x,char y)
{
if(x==0)
{
putchar('0');putchar(y);
return;
}
if(x<0)
{
putchar('-');
x=-x;
}
static char wr[20];
int top=0;
for(;x;x/=10)wr[++top]=x%10+'0';
while(top)putchar(wr[top--]);
putchar(y);
}
void file()
{
#ifndef ONLINE_JUDGE
freopen("3722.in","r",stdin);
freopen("3722.out","w",stdout);
#endif
}
int n,m,p1,p2;
int a[N];
struct Q
{
int l,r,x,tag,id;
bool operator < (const Q &s)const {return x<s.x;}
}opt[N<<2];
int cnt;
ll ans[N];
void input()
{
int l,r;
n=read<int>();m=read<int>();p1=read<int>();p2=read<int>();
For(i,1,n)a[i]=read<int>();
For(i,1,m)
{
l=read<int>();r=read<int>();
ans[i]+=1ll*(r-l)*1ll*p1;
opt[++cnt]=(Q){l,r,l-1,-1,i};
opt[++cnt]=(Q){l,r,r,1,i};
}
}
const int inf=0x3f3f3f3f;
namespace Tree
{
#define mid ((l+r)>>1)
#define lson h<<1,l,mid
#define rson h<<1|1,mid+1,r
ll sum[N<<2],lazy[N<<2];
void push_up(int h)
{
sum[h]=sum[h<<1]+sum[h<<1|1];
}
void push_down(int h,int l,int r)
{
if(!lazy[h])return;
int ls=h<<1,rs=ls|1;
lazy[ls]+=lazy[h];lazy[rs]+=lazy[h];
sum[ls]+=lazy[h]*1ll*(mid-l+1);
sum[rs]+=lazy[h]*1ll*(r-mid);
lazy[h]=0;
}
void update(int h,int l,int r,int s,int t,ll v)
{
if(s<=l&&r<=t)
{
lazy[h]+=v;
sum[h]+=1ll*v*1ll*(r-l+1);
}
else
{
push_down(h,l,r);
if(s<=mid)update(lson,s,t,v);
if(mid<t)update(rson,s,t,v);
push_up(h);
}
}
ll query(int h,int l,int r,int s,int t)
{
if(s<=l&&r<=t)return sum[h];
push_down(h,l,r);
ll res=0;
if(s<=mid)res=query(lson,s,t);
if(mid<t)res+=query(rson,s,t);
push_up(h);
return res;
}
}
#define rg register
int lef[N],rig[N];
struct node
{
int l,r,x;
ll v;
bool operator < (const node &s)const {return x<s.x;}
}e[N<<2];
int sz;
void init()
{
int j;
sort(opt+1,opt+cnt+1);
a[0]=a[n+1]=inf;
For(i,1,n)
{
for(j=i-1;j>=0;j=lef[j])if(a[j]>a[i])break;
lef[i]=j;
}
Fordown(i,n,1)
{
for(j=i+1;j<=n+1;j=rig[j])if(a[j]>a[i])break;
rig[i]=j;
}
//For(i,1,n)cout<<lef[i]<<' '<<rig[i]<<endl;
For(i,1,n)
{
if(lef[i]+1<i&&rig[i]<=n)e[++sz]=(node){lef[i]+1,i-1,rig[i],p2};
if(1<=lef[i]&&rig[i]>i+1)e[++sz]=(node){i+1,rig[i]-1,lef[i],p2};
if(1<=lef[i]&&rig[i]<=n)e[++sz]=(node){lef[i],lef[i],rig[i],p1};
}
sort(e+1,e+sz+1);
}
void work()
{
int pos1=1,pos2=1;
while(pos2<=cnt&&!opt[pos2].x)pos2++;
//cerr<<pos2<<' '<<cnt<<endl;
For(i,1,n)
{
for(;pos1<=sz&&e[pos1].x==i;++pos1)
Tree::update(1,0,n,e[pos1].l,e[pos1].r,e[pos1].v);
for(;pos2<=cnt&&opt[pos2].x==i;++pos2)
ans[opt[pos2].id]+=Tree::query(1,0,n,opt[pos2].l,opt[pos2].r)*opt[pos2].tag;
}
For(i,1,m)write(ans[i],'\n');
}
int main()
{
file();
input();
init();
work();
return 0;
}
HNOI2017影魔的更多相关文章
- bzoj 4826: [Hnoi2017]影魔 [主席树 单调栈]
4826: [Hnoi2017]影魔 题意:一个排列,点对\((i,j)\),\(p=max(i+1,j-1)\),若\(p<a_i,a_j\)贡献p1,若\(p\)在\(a_1,a_2\)之间 ...
- 4826: [Hnoi2017]影魔
4826: [Hnoi2017]影魔 https://lydsy.com/JudgeOnline/problem.php?id=4826 分析: 莫队+单调栈+st表. 考虑如何O(1)加入一个点,删 ...
- 【LG3722】[HNOI2017]影魔
[LG3722][HNOI2017]影魔 题面 洛谷 题解 先使用单调栈求出\(i\)左边第一个比\(i\)大的位置\(lp_i\),和右边第一个比\(i\)大的位置\(rp_i\). 考虑\(i\) ...
- [BZOJ4826][HNOI2017]影魔(主席树)
4826: [Hnoi2017]影魔 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 669 Solved: 384[Submit][Status][ ...
- 【BZOJ4826】[Hnoi2017]影魔 单调栈+扫描线
[BZOJ4826][Hnoi2017]影魔 Description 影魔,奈文摩尔,据说有着一个诗人的灵魂.事实上,他吞噬的诗人灵魂早已成千上万.千百年来,他收集了各式各样的灵魂,包括诗人.牧师.帝 ...
- [bzoj4826][Hnoi2017]影魔_单调栈_主席树
影魔 bzoj-4826 Hnoi-2017 题目大意:给定一个$n$个数的序列$a$,求满足一下情况的点对个数: 注释:$1\le n,m\le 2\cdot 10^5$,$1\le p1,p2\l ...
- bzoj4826 [Hnoi2017]影魔
Description 影魔,奈文摩尔,据说有着一个诗人的灵魂.事实上,他吞噬的诗人灵魂早已成千上万.千百年来,他收集了各式各样的灵魂,包括诗人.牧师.帝王.乞丐.奴隶.罪人,当然,还有英雄.每一个灵 ...
- BZOJ:4826: [Hnoi2017]影魔
Description 影魔,奈文摩尔,据说有着一个诗人的灵魂.事实上,他吞噬的诗人灵魂早已成千上万.千百年来,他收集了各式各样的灵魂,包括诗人.牧师.帝王.乞丐.奴隶.罪人,当然,还有英雄.每一个灵 ...
- [AH/HNOI2017]影魔
题目背景 影魔,奈文摩尔,据说有着一个诗人的灵魂. 事实上,他吞噬的诗人灵魂早已成千上万. 千百年来,他收集了各式各样的灵魂,包括诗人. 牧师. 帝王. 乞丐. 奴隶. 罪人,当然,还有英雄. 题目描 ...
随机推荐
- 认识ASP.NET Windows身份认证
本文摘自:细说ASP.NET Windows身份认证 Forms身份认证虽然使用广泛,不过,如果是在 Windows Active Directory 的环境中使用ASP.NET, 那么使用Windo ...
- angularjs呼叫Web API
今早有分享一篇<创建Web API并使用>http://www.cnblogs.com/insus/p/7771428.html 接下来,我再分享一篇,怎样在angularjs去呼叫Web ...
- 表单验证(AngularJs)
这次,学习angularjs的表单的验证,angularjs提供上表几种状态验证: 状态 描述 $invalid 未通过验证 $valid 经过验证 $pristine 未修改过 $dirty 修改 ...
- Codeforces Edu Round 63(Rated for Div. 2)
感觉现在Edu场比以前的难多了…… A: 温暖人心 /* basic header */ #include <iostream> #include <cstdio> #incl ...
- FreeCAD源码初步了解
FreeCAD简介 FreeCAD是基于OpenCASCADE的开源CAD/CAE软件,完全开源(GPL的LGPL许可证),官方源码地址,详情可参考维基百科,百度百科等等. 如果要编译FreeCAD, ...
- 分布式监控系统Zabbix-图形集中展示插件Graphtree安装笔记
Zabbix想要集中展示图像,唯一的选择是screen,后来zatree解决了screen的问题,但性能不够好.Graphtree 由OneOaaS开发并开源出来,用来解决Zabbix的图形展示问题, ...
- Mesos+Zookeeper+Marathon的Docker管理平台部署记录(1)
随着"互联网+"时代的业务增长.变化速度及大规模计算的需求,廉价的.高可扩展的分布式x86集群已成为标准解决方案,如Google已经在几千万台服务器上部署分布式系统.Docker及 ...
- Websocket(一)——原理及基本属性和方法
初次接触 WebSocket 的人,都会问同样的问题:我们已经有了 HTTP 协议,为什么还需要另一个协议?它能带来什么好处? 答案很简单,因为 HTTP 协议有一个缺陷:通信只能由客户端发起. 举例 ...
- 在Mac终端显示 Git 当前所在分支
1.进入你的home目录 cd ~ 2.编辑.bashrc文件 vi .bashrc 3.将下面的代码加入到文件的最后处 function git_branch { branch="`git ...
- 如何用chrome查看post get及返回的数据
chrome浏览器按下F12打开开发者工具 点击Network,找到过滤器 筛选XHR,Method那一列会显示POST GET: