影魔
  • 这么简单的方法尽然想不到,我是真的菜
  • 对每个点,用单调栈的方式处理出他左右第一个比他大的数的位置,你可以把\(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影魔的更多相关文章

  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. 4826: [Hnoi2017]影魔

    4826: [Hnoi2017]影魔 https://lydsy.com/JudgeOnline/problem.php?id=4826 分析: 莫队+单调栈+st表. 考虑如何O(1)加入一个点,删 ...

  3. 【LG3722】[HNOI2017]影魔

    [LG3722][HNOI2017]影魔 题面 洛谷 题解 先使用单调栈求出\(i\)左边第一个比\(i\)大的位置\(lp_i\),和右边第一个比\(i\)大的位置\(rp_i\). 考虑\(i\) ...

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

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

  5. 【BZOJ4826】[Hnoi2017]影魔 单调栈+扫描线

    [BZOJ4826][Hnoi2017]影魔 Description 影魔,奈文摩尔,据说有着一个诗人的灵魂.事实上,他吞噬的诗人灵魂早已成千上万.千百年来,他收集了各式各样的灵魂,包括诗人.牧师.帝 ...

  6. [bzoj4826][Hnoi2017]影魔_单调栈_主席树

    影魔 bzoj-4826 Hnoi-2017 题目大意:给定一个$n$个数的序列$a$,求满足一下情况的点对个数: 注释:$1\le n,m\le 2\cdot 10^5$,$1\le p1,p2\l ...

  7. bzoj4826 [Hnoi2017]影魔

    Description 影魔,奈文摩尔,据说有着一个诗人的灵魂.事实上,他吞噬的诗人灵魂早已成千上万.千百年来,他收集了各式各样的灵魂,包括诗人.牧师.帝王.乞丐.奴隶.罪人,当然,还有英雄.每一个灵 ...

  8. BZOJ:4826: [Hnoi2017]影魔

    Description 影魔,奈文摩尔,据说有着一个诗人的灵魂.事实上,他吞噬的诗人灵魂早已成千上万.千百年来,他收集了各式各样的灵魂,包括诗人.牧师.帝王.乞丐.奴隶.罪人,当然,还有英雄.每一个灵 ...

  9. [AH/HNOI2017]影魔

    题目背景 影魔,奈文摩尔,据说有着一个诗人的灵魂. 事实上,他吞噬的诗人灵魂早已成千上万. 千百年来,他收集了各式各样的灵魂,包括诗人. 牧师. 帝王. 乞丐. 奴隶. 罪人,当然,还有英雄. 题目描 ...

随机推荐

  1. React-Route的属性exact

    exact是Route下的一条属性,一般而言,react路由会匹配所有匹配到的路由组价,exact能够使得路由的匹配更严格一些. exact的值为bool型,为true是表示严格匹配,为false时为 ...

  2. (转)Xpath语法格式整理

    原文 经常在工作中会使用到XPath的相关知识,但每次总会在一些关键的地方不记得或不太清楚,所以免不了每次总要查一些零碎的知识,感觉即很烦又浪费时间,所以对XPath归纳及总结一下. 在这篇文章中你将 ...

  3. 从harbor部署到在k8s中使用

    一.概述 harbor是什么呢?英文单词的意思是:港湾.港湾用来存放集装箱(货物的),而docker的由来正是借鉴了集装箱的原理,所以harbor是用于存放docker的镜像,作为镜像仓库使用.官方的 ...

  4. Luogu P1514 引水入城

    我承认我有点懒(洛谷已经发过题解了,但我发誓要坚持写博客) 这道题坑了我3天…… 首先一看就与染色问题类似,果断BFS(写DFS炸了) 先将最上面(靠近水)的一行全部扔进队列里,做一遍BFS 再对最下 ...

  5. bitcoin 源码解析 - 交易 Transaction(二) - 原理篇

    这篇文章我断断续续写了呃···· 应该快三个星期了? 所以前后的风格可能差别相当大.真是十分的怠惰啊··· 最近实在是不够努力.用python重写bitcoin的项目也卡在网络编程部分(这方面真是我的 ...

  6. Centos下安装破解confluence6.3的操作记录

    confluence是一个专业的企业知识管理与协同软件,可以用于构建企业wiki.通过它可以实现团队成员之间的协作和知识共享.现在大多数公司都会部署一套confluence,用作内部wiki.现在co ...

  7. vue项目环境搭建

    安装node.js $ npm install -g vue-cli $ vue init webpack my-project ?Project name ?Project description ...

  8. Spring方法级别的验证

    设置验证点及验证方式(1)Spring方法级别的验证有多种验证方式,比较常用的有 @NotBlank:主要是对字符串的验证,不为null且去除空白符之后长度大于0 @NotNull:主要是对对象的验证 ...

  9. linux内核期中总结

    20135132陈雨鑫 + 原创作品转载请注明出处 + <Linux内核分析>MOOC课程http://mooc.study.163.com/course/USTC-1000029000  ...

  10. JQuery监听页面滚动总结

    1.当前滚动的地方的窗口顶端到整个页面顶端的距离: var winPos = $(window).scrollTop(); 2.获取指定元素的页面位置: $(val).offset().top; 3. ...