[BZOJ4826][HNOI2017]影魔(主席树)
4826: [Hnoi2017]影魔
Time Limit: 20 Sec Memory Limit: 512 MB
Submit: 669 Solved: 384
[Submit][Status][Discuss]Description
影魔,奈文摩尔,据说有着一个诗人的灵魂。事实上,他吞噬的诗人灵魂早已成千上万。千百年来,他收集了各式各样的灵魂,包括诗人、牧师、帝王、乞丐、奴隶、罪人,当然,还有英雄。每一个灵魂,都有着自己的战斗力,而影魔,靠这些战斗力提升自己的攻击。奈文摩尔有 n 个灵魂,他们在影魔宽广的体内可以排成一排,从左至右标号 1 到 n。第 i个灵魂的战斗力为 k[i],灵魂们以点对的形式为影魔提供攻击力,对于灵魂对 i,j(i<j)来说,若不存在 k[s](i<s<j)大于 k[i]或者 k[j],则会为影魔提供 p1 的攻击力(可理解为:当 j=i+1 时,因为不存在满足 i<s<j 的 s,从而 k[s]不存在,这时提供 p1 的攻击力;当 j>i+1 时,若max{k[s]|i<s<j}<=min{k[i],k[j]} , 则 提 供 p1 的 攻击 力 ); 另 一 种 情 况 , 令 c 为k[i+1],k[i+2],k[i+3]......k[j-1]的最大值,若 c 满足:k[i]<c<k[j],或者 k[j]<c<k[i],则会为影魔提供 p2 的攻击力,当这样的 c 不存在时,自然不会提供这 p2 的攻击力;其他情况的点对,均不会为影魔提供攻击力。影魔的挚友噬魂鬼在一天造访影魔体内时被这些灵魂吸引住了,他想知道,对于任意一段区间[a,b],1<=a<b<=n,位于这些区间中的灵魂对会为影魔提供多少攻击力,即考虑 所有满足a<=i<j<=b 的灵魂对 i,j 提供的攻击力之和。顺带一提,灵魂的战斗力组成一个 1 到 n 的排列:k[1],k[2],...,k[n]。Input
第一行 n,m,p1,p2第二行 n 个数:k[1],k[2],...,k[n]接下来 m 行,每行两个数 a,b,表示询问区间[a,b]中的灵魂对会为影魔提供多少攻击力。1 <= n,m <= 200000;1 <= p1,p2 <= 1000Output
共输出 m 行,每行一个答案,依次对应 m 个询问。Sample Input
10 5 2 3
7 9 5 1 3 10 6 8 2 4
1 7
1 9
1 3
5 9
1 5Sample Output
30
39
4
13
16HINT
Source
比较巧妙的地方是将询问映射到二维坐标系里,虽然比较套路但考场上不一定能想到。剩下的就是主席树裸题了。
代码用时:40min。
- #include<cstdio>
- #include<vector>
- #include<algorithm>
- #define rep(i,l,r) for (int i=l; i<=r; i++)
- typedef long long ll;
- using namespace std;
- const int N=;
- int n,m,p1,p2,l,r,cnt,top,rt[N],a[N],L[N],R[N],stk[N],v[N*],ls[N*],rs[N*];
- ll s[N*];
- struct P{ int x,y; ll p; };
- vector<P>V[N];
- void mdf(int y,int &x,int L,int R,int l,int r,int k){
- v[x=++cnt]=v[y]; s[x]=s[y]; ls[x]=ls[y]; rs[x]=rs[y];
- s[x]+=1ll*(r-l+)*k;
- if (l==L && r==R){ v[x]+=k; return; }
- int mid=(L+R)>>;
- if (r<=mid) mdf(ls[y],ls[x],L,mid,l,r,k);
- else if (l>mid) mdf(rs[y],rs[x],mid+,R,l,r,k);
- else mdf(ls[y],ls[x],L,mid,l,mid,k),mdf(rs[y],rs[x],mid+,R,mid+,r,k);
- }
- ll que(int x,int L,int R,int l,int r){
- if (!x) return ;
- if (l==L && r==R) return s[x];
- int mid=(L+R)>>; ll sum=1ll*(r-l+)*v[x];
- if (r<=mid) return sum+que(ls[x],L,mid,l,r);
- else if (l>mid) return sum+que(rs[x],mid+,R,l,r);
- else return sum+que(ls[x],L,mid,l,mid)+que(rs[x],mid+,R,mid+,r);
- }
- int main(){
- scanf("%d%d%d%d",&n,&m,&p1,&p2);
- rep(i,,n) scanf("%d",&a[i]);
- rep(i,,n){
- while (top && a[stk[top]]<a[i]) top--;
- L[i]=stk[top]; stk[++top]=i;
- }
- rep(i,,top) stk[i]=; stk[top=]=n+;
- for (int i=n; i; i--){
- while (top && a[stk[top]]<a[i]) top--;
- R[i]=stk[top]; stk[++top]=i;
- }
- rep(i,,n){
- if (i<n) V[i].push_back((P){i+,i+,p1});
- if (L[i] && R[i]<=n) V[L[i]].push_back((P){R[i],R[i],p1});
- if (L[i] && i<=R[i]-) V[L[i]].push_back((P){i+,R[i]-,p2});
- if (R[i]<=n && i>=L[i]+) V[R[i]].push_back((P){L[i]+,i-,p2});
- }
- rt[]=;
- rep(i,,n){
- rt[i]=rt[i-];
- for (vector<P>::iterator it=V[i].begin(); it!=V[i].end(); it++)
- mdf(rt[i],rt[i],,n,it->x,it->y,it->p);
- }
- rep(i,,m) scanf("%d%d",&l,&r),printf("%lld\n",que(rt[r],,n,l,r)-que(rt[l-],,n,l,r));
- return ;
- }
[BZOJ4826][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\)之间 ...
- [bzoj4826][Hnoi2017]影魔_单调栈_主席树
影魔 bzoj-4826 Hnoi-2017 题目大意:给定一个$n$个数的序列$a$,求满足一下情况的点对个数: 注释:$1\le n,m\le 2\cdot 10^5$,$1\le p1,p2\l ...
- [BZOJ4826] [HNOI2017] 影魔 单调栈 主席树
题面 因为是一个排列,所以不会有重复的.如果有重复就没法做了.一开始没有仔细看题目想了半天. 发现,如果是第一种情况,那么边界\(l\)和\(r\)就应该分别是整个区间的最大值和次大值. 然后,对于那 ...
- bzoj4826 [Hnoi2017]影魔
Description 影魔,奈文摩尔,据说有着一个诗人的灵魂.事实上,他吞噬的诗人灵魂早已成千上万.千百年来,他收集了各式各样的灵魂,包括诗人.牧师.帝王.乞丐.奴隶.罪人,当然,还有英雄.每一个灵 ...
- bzoj千题计划196:bzoj4826: [Hnoi2017]影魔
http://www.lydsy.com/JudgeOnline/problem.php?id=4826 吐槽一下bzoj这道题的排版是真丑... 我还是粘洛谷的题面吧... 提供p1的攻击力:i,j ...
- BZOJ4826 [Hnoi2017]影魔 【线段树 + 单调栈】
题目链接 BZOJ4826 题解 蒟蒻智力水平捉急orz 我们会发现相邻的\(i\)和\(j\)贡献一定是\(p1\),可以很快算出来[然而我一开始忘了考虑调了半天] 我们现在只考虑不相邻的 我们只需 ...
- [BZOJ4826][HNOI2017]影魔 可持久化线段树
链接 题意:给你 \(1\) 到 \(n\) 的排列 \(k_1,k_2,\dots,k_n\) ,对 \(i,j (i<j)\)来说,若不存在 \(k_s (i<s<j)\) 大于 ...
- 洛谷P3722 [AH2017/HNOI2017]影魔(线段树)
题意 题目链接 Sol 题解好神仙啊qwq. 一般看到这种考虑最大值的贡献的题目不难想到单调数据结构 对于本题而言,我们可以预处理出每个位置左边第一个比他大的位置\(l_i\)以及右边第一个比他大的位 ...
- [AH2017/HNOI2017] 影魔 - 线段树
#include<bits/stdc++.h> #define maxn 200010 using namespace std; int a[maxn],st[maxn][2],top,L ...
随机推荐
- 【洛谷 P3338】 [ZJOI2014]力(FFT)
题目链接 \[\Huge{E_i=\sum_{j=1}^{i-1}\frac{q_j}{(i-j)^2}-\sum_{j=i+1}^{n}\frac{q_j}{(i-j)^2}}\] 设\(A[i]= ...
- MongoDB 3.4.2 配置 CentOS 6.5 远程连接
1.新建用户 db.createUser({user: 'test', pwd: 'myPassword', roles: [{role: 'readWrite', db: 'test_db'}]}) ...
- Mysql储存过程6: in / out / inout
in 为向函数传送进去的值 out 为函数向外返回的值 intout 传送进去的值, 并且还返回这个值 )) begin then select 'true'; else select 'false' ...
- 爬行百度标题&URL案例
思路: 先将需要获取的匹配出,然后可以用"永真"(即while True:)来遍历使得URL可以一直自增变化(百度点击下一页URL的pn参数就增加10)每增加10就爬行一遍URL然 ...
- 每天一条linux命令(1):ls命令
ls命令是linux下最常用的命令.ls命令就是list的缩写缺省下ls用来打印出当前目录的清单如果ls指定其他目录那么就会显示指定目录里的文件及文件夹清单. 通过ls 命令不仅可以查看linu ...
- python之微信公众号开发(基本配置和校验)
前言 最近有微信公众号开发的业务,以前没有用python做过微信公众号开发,记录一下自己的学习和开发历程,共勉! 公众号类型 订阅号 普通订阅号 认证订阅号 服务号 普通服务号 认证服务号 服务方式 ...
- django框架<二>
django框架: Models 1.基本创建 Django提供了一个抽象层("Model")的构建和管理Web应用程序的数据. Django使用一种新的方式,即:关系对象映射 ...
- Serv-U设置允许用户更改密码【转】
最近,公司上了一套Serv-U10.5.0.6的ftp软件,应该是目前最新的版本了.上的第一天就遇到了一个问题,有领导发话了,他需要自己更改密码.找了N久才找到,分享一下. 点击管理界面的用户. 进入 ...
- 理解一条语句:SELECT difference(sum("value")) FROM "mq_enqueue" WHERE "channel" =~ /ActiveMQ_TEST/ AND $timeFilter GROUP BY time($interval)
最近使用grafana在查询InfluxDB中,用到了这一条语句 SELECT difference(sum("value")) FROM "mq_enqueue&quo ...
- Builder设计模式--改善构造器多个参数时可显著改善可读性
作为一名程序开发者,设计模式其实一直有在接触,只是没有专门的去学过,所以可能对设计模式没有一个系统的理解.在一次项目中,需要使用到第三方服务商提供的功能,为了尽快的熟悉其功能代码,在官网下了demo来 ...