[HNOI2017]影魔
题意:
给定 \(n\) 个数的排列,\(m\) 次询问,每次询问询问一个区间内所有子区间的贡献。
每个区间如果两个端点分别是最大值和次大值,我们就算 \(P1\) 的贡献。
如果两个端点一个是最大值,一个不是次大值,我们就算 \(P2\) 的贡献。
\(\text{Solution:}\)
将询问离线,处理出以 \(i\) 结尾的询问的答案。
考虑怎样的点 \(j(j<i)\) 满足它能与 \(i\) 组成区间 \(i\) ,\(j\) 分别是最大值与次大值或非次大值。
钦定 \(a[i]\) 比 \(a[j]\) 大,那么当 \(a[j]\) 为左边单调减的波峰的时候与 \(a[i]\) 有贡献 \(p1\) ,当 \(a[j]\) 为左边两个单调减的波峰之间的时候与 \(a[i]\) 有贡献 \(p2\) 。
维护单调减的波峰可以用单调栈。
由于钦定了 \(a[i] > a[j]\) 所以要倒过来再搞一遍。
像这种区间上不好直接维护的东西,可以考虑离线,讨论时要讨论清楚哪种点有贡献。
\(\text{Source:}\)
#include<vector>
#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <assert.h>
#include <algorithm>
using namespace std;
#define LL long long
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define GO debug("GO\n")
inline int rint() {
register int x = 0, f = 1; register char c;
while (!isdigit(c = getchar())) if (c == '-') f = -1;
while (x = (x << 1) + (x << 3) + (c ^ 48), isdigit(c = getchar()));
return x * f;
}
template<typename T> inline void chkmin(T &a, T b) { a > b ? a = b : 0; }
template<typename T> inline void chkmax(T &a, T b) { a < b ? a = b : 0; }
const int N = 2e5 + 10;
int n,m,p1,p2;
int a[N];
vector<pair<int,int> >Q,QR[N];
LL ans[N];
namespace SGT {
#define ls (x<<1)
#define rs (x<<1|1)
LL sum[N * 4], tag[N * 4];
void init() {
memset(sum, 0, sizeof sum);
memset(tag, 0, sizeof tag);
}
void pu(int x) {
sum[x] = sum[ls] + sum[rs];
}
void pd(int x, int l, int r) {
if (tag[x]) {
int mid = (l + r) >> 1;
sum[ls] += tag[x] * (mid - l + 1);
tag[ls] += tag[x];
sum[rs] += tag[x] * (r - mid);
tag[rs] += tag[x];
tag[x] = 0;
}
}
void Add(int L, int R, int val, int x = 1, int l = 1, int r = n) {
if (L > R) return;
if (L <= l and r <= R) {
sum[x] += val * (r - l + 1);
tag[x] += val;
return;
}
pd(x, l, r); int mid = (l + r) >> 1;
if (L <= mid) Add(L, R, val, ls, l, mid);
if (R > mid) Add(L, R, val, rs, mid + 1, r);
pu(x);
}
LL query(int L, int R, int x = 1, int l = 1, int r = n) {
if (L > R) return 0;
if (L <= l and r <= R) {
return sum[x];
}
pd(x, l, r); int mid = l + r >> 1;
LL ans = 0;
if (L <= mid) ans += query(L, R, ls, l, mid);
if (mid < R) ans += query(L, R, rs, mid + 1, r);
return ans;
}
}using SGT::Add; using SGT::query;
void work() {
static int stk[N], top;
SGT::init();
top = 0;
for (int i = 1; i <= n; ++ i) {
int last = i;
while (top and a[stk[top]] <= a[i]) {
if (stk[top] + 1 < last)
Add(stk[top] + 1, last - 1, p2);
Add(stk[top], stk[top], p1);
last = stk[top--];
}
if (!top) Add(1, last - 1, p2);
else Add(stk[top] + 1, last - 1, p2);
stk[++top] = i;
for (int j = 0; j < QR[i].size(); ++ j)
ans[QR[i][j].second] += query(QR[i][j].first, i);
}
}
int main(){
#ifndef ONLINE_JUDGE
freopen("xhc.in", "r", stdin);
freopen("xhc.out", "w", stdout);
#endif
n = rint(), m = rint(), p1 = rint(), p2 = rint();
for (int i = 1; i <= n; ++ i) a[i] = rint();
for (int i = 1; i <= m; ++ i){
int l = rint(), r = rint();
Q.push_back(make_pair(l, r));
}
for (int i = 0; i < m; ++ i)
QR[Q[i].second].push_back(make_pair(Q[i].first, i + 1));
work();
for (int i = 1; i <= n; ++ i)
QR[i].clear();
for (int i = 0; i < m; ++ i)
QR[n - Q[i].first + 1].push_back(make_pair(n - Q[i].second + 1, i + 1));
reverse(a + 1, a + 1 + n);
work();
for (int i = 1; i <= m; ++ i) printf("%lld\n", ans[i]);
}
[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]影魔
题目背景 影魔,奈文摩尔,据说有着一个诗人的灵魂. 事实上,他吞噬的诗人灵魂早已成千上万. 千百年来,他收集了各式各样的灵魂,包括诗人. 牧师. 帝王. 乞丐. 奴隶. 罪人,当然,还有英雄. 题目描 ...
- HNOI2017影魔
影魔 这么简单的方法尽然想不到,我是真的菜 对每个点,用单调栈的方式处理出他左右第一个比他大的数的位置,你可以把\(0\)和\(n+1\)设成\(inf\). 显然对于每对\(lef[i]\)和\(r ...
随机推荐
- 使用PLSQL客户端登录ORACLE时报ORA-12502和ORA-12545错误的解决方案
当使用plsqldev客户端登录oracle数据库时,如果对应的tnsnames.ora中是直接使用IP地址来连接,并且未在系统的hosts文件中绑定主机名时,极易出现ORA-12502及ORA-12 ...
- iOS:绘图(18-01-25更)
目录 1.UIBezierPath(贝塞尔曲线) 1).在重写 drawRect: 方法里使用 2).在普通方法里使用,需要画布. 3). 切圆角.指定位置圆角.任意形状. 4).彩色的动画加载圆圈. ...
- idea配置SpringBoot热部署之自动Build
一.pom.xml文件导入所需依赖文件 SpringBoot热部署插件 <dependency> <groupId>org.springframework.boot</g ...
- 【NXP开发板应用—智能插排】4. PWM驱动
[前言] 首先感谢深圳市米尔科技有限公司举办的这次活动并予以本人参加这次活动的机会,以往接触过嵌入式,但那都是皮毛,最多刷个系统之类的,可以说对于嵌入式系统开发这件事情是相当非常陌生的,这次活动为我提 ...
- GitLab 基本操作
登录 在浏览其中输入http://192.168.3.11:8888 如图1登录界面. 图1 注:第一次新增用户,会发送修改密码链接到用户的邮箱中,用户会收到如图2邮件. 图2 2. 修改密码 点 ...
- python学习笔记:第21天 常用内置模块之collections和time
目录 一.collections模块 二.时间模块 也可以在我的个人博客上阅读 一.collections模块 1. Counter Counter是⼀个计数器,主要⽤统计字符的数量,之前如果我们要统 ...
- 闰年相关的问题v3.0——计算有多少闰年
# include<stdio.h>int main(){ int a,b,i; int sum = 0; printf("Input your birth year:" ...
- 初识Trie_对Trie的一些认识
Trie,之前觉得一个听起来很diao的数据结构,学了发现其实是一个挺简单的数据结构. 在计算机科学中,trie,又称前缀树或字典树,是一种有序树,用于保存关联数组,其中的键通常是字符串.与二叉查找树 ...
- springBoot 自定义redisTemplate
package com.atirm.mybatismutiplesource.config.RedisConfig; import com.atirm.mybatismutiplesource.ent ...
- itop4412开发板添加开机启动程序
1. 先编写代码,以helloworld.c为例子 #include<stdio.h> #include<unistd.h> //这个文件是什么 main() { ; ) { ...