[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 ...
随机推荐
- OSD仿真_MFC程序01
Windows系统具有强大的绘图功能,可以用来模拟OSD显示.接下来将设计一个简单的模拟显示终端的程序,用于后续显示功能和菜单系统的开发.说明一下,对于Windows下的MFC编程我不怎么了解,只知道 ...
- java连接linux的三种方式(附执行命令)
# 本地调用使用JDK自带的RunTime类和Process类实现 public static void main(String[] args){ Process proc = RunTime.get ...
- linux 内存操作相关命令
清理内存 echo 1 > /proc/sys/vm/drop_caches 查看内存使用情况 free –m 查看内存条数命令: dmidecode |grep -A16 "Memo ...
- JS-类型相关
typeof检测类型typeof 返回的数据类型种类:number(js不分整形,浮点等等 所有的数字都是number类型).string.boolean.undefined.object.funct ...
- Spring的入门学习笔记 (注解)
Spring的bean管理(注解) 注解 1.代码里面特殊标记(ep:@Test),使用注解完成一些相关功能 2.注解写法 @注解名称(属性名称=属性值) 3.可以用在类,方法,属性上都可以 4.Sp ...
- DBUtils 学习使用
DBUtils 学习使用 commons-dbutils简介 commons-dbutils是Apache组织提供的一个开源JDBC工具类库,它是对JDBC的简单封装,学习成本极低,并且使用dbuti ...
- 一分钟完成pip安装
很多实用Python的小伙伴都需要使用pip安装相应的包,对于初学者而已,检查遇到pip安装不成功的情况,如以下典型错误: Traceback (most recent call last): Fil ...
- zkfc的znode不存在的问题
cd /soft/hadoop/logs/hadoop-centos-zkfc-s101.log发现: 2018-09-29 12:42:03,616 FATAL org.apache.hadoop. ...
- go基础语法-循环语句
1.基础定义 for语句的条件不需要括号(同if语句) ,golang里的循环只有for,没有while sum := 0 for i=0;i<100;i++ { sum += i } 2.条件 ...
- Go语言入门(一)特性、安装、环境搭建、第一个程序等
Go语言是谷歌推出的一种全新的编程语言,可以在不损失应用程序性能的情况下降低代码的复杂性.谷歌首席软件工程师罗布派克(Rob Pike)说:我们之所以开发Go,是因为过去10多年间软件开发的难度令人沮 ...