【题解】Diferenc-Diferencija [SP10622]

传送门:\(\text{Diferenc-Diferencija}\) \(\text{[SP10622]}\)

【题目描述】

序列的值被定义成其中最大的元素减去最小的元素。如序列 \((3,1,7,2)\) 的值为 \(7-1=6\), 序列 \((42,42)\) 的值为 \(42-42=0\)。

现给定一长为 \(n\) 的序列 \(a\),求出所有连续子序列的值的和。

【样例】

样例输入:
3
1
2
3
样例输出:
4 样例输入:
4
7
5
7
5
样例输出:
12 样例输入:
4
3
1
7
2
样例输出:
31

【数据范围】

\(100 \%:\) \(2 \leqslant n \leqslant 3*10^5,\) \(1 \leqslant a[i] \leqslant 10^8\)


【分析】

先将子区间的右端点固定为 \(r\),此时一共有 \(r\) 个左端点 \((l \in [1,r])\) 可与之组成连续子序列,用 \(f_1[l]\) 表示 \(max \{a[j]\}(j \in [l,r])\) ,\(f_2[l]\) 表示 \(min \{a[j]\}(j \in [l,r])\) 。于是以 \(i\) 为右端点 \(r\) 的子区间贡献和为 \(\sum_{l=1}^{r} (f_1[l]-f_2[l])\),即 \(\sum_{l=1}^{r} f_1[l] - \sum_{l=1}^{r} f_2[l]\) 。我们可以分开算 \(f_1,f_2\) 的总和。

当右端点移至 \(r+1\) 时,需要用 \(a[r+1]\) 来更新 \(f_1,f_2\),可以直接扫描 \([1,r]\),但时间不过不了关。

随着 \(l\) 的减小,\(f_1[l]=max(a[l],f_1[l+1])\),可以发现 \(f_1[l]\) 是单调不下降的,同理 \(f_2[l]\) 单调不上升。

随着 \(r\) 的增大,\(f_1[l]=max(f_1[l],a[r+1])\),可以发现 \(f_1[l]\) 仍是单调不下降的,同理 \(f_2[l]\) 单调不上升。

然后我们就会发现一个现象:每次新加进来一个数 \(a[r+1]\) 时,它会将以 \(r+1\) 结尾的一段连续的区间 \(f_1[l],f_2[l](l \in [?,r+1])\) 全部赋值为 \(a[r+1]\),而且被覆盖掉的原数对这之后的区间不再有任何贡献。

于是我们可以用两个单调栈分别维护 \(f_1,f_2\),由于下标也是单调递增,所以可以将 \(f\) 值相同的合并起来,用 \(g\) 记录 \(f\) 相同的下标个数,另设一个变量 \(S\) 表示以 \(i\) 为右端点的贡献和,当加入新的 \(a[r]\) 时就不断弹走队尾直至保持单调时结束,同时更新 \(S\),最后累加答案即可。

【Code】

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#define LL long long
#define Re register int
using namespace std;
const int N=3e5+2;
int n,t1,t2,a[N],f1[N],f2[N],g1[N],g2[N];LL S,ans;
inline void in(Re &x){
int f=0;x=0;char c=getchar();
while(c<'0'||c>'9')f|=c=='-',c=getchar();
while(c>='0'&&c<='9')x=(x<<1)+(x<<3)+(c^48),c=getchar();
x=f?-x:x;
}
int main(){
in(n);
for(Re i=1;i<=n;++i)in(a[i]);
f1[++t1]=a[1],f2[++t2]=a[1],g1[t1]=g2[t2]=1,ans=S=0;//初始化入队
for(Re i=2;i<=n;++i){
Re tmp=1;//f1[i]和f2[i]都一定会被覆盖,所以初始化为1
while(t1&&f1[t1]<=a[i])S-=(LL)f1[t1]*g1[t1],tmp+=g1[t1--];//更新最大值
f1[++t1]=a[i],g1[t1]=tmp,S+=(LL)a[i]*tmp;
tmp=1;
while(t2&&f2[t2]>=a[i])S+=(LL)f2[t2]*g2[t2],tmp+=g2[t2--];//更新最小值
f2[++t2]=a[i],g2[t2]=tmp,S-=(LL)a[i]*tmp;
ans+=S;//累加答案
}
printf("%lld\n",ans);
return 0;
}

【题解】Diferenc-Diferencija [SP10622]的更多相关文章

  1. 2016 华南师大ACM校赛 SCNUCPC 非官方题解

    我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...

  2. noip2016十连测题解

    以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...

  3. BZOJ-2561-最小生成树 题解(最小割)

    2561: 最小生成树(题解) Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1628  Solved: 786 传送门:http://www.lyd ...

  4. Codeforces Round #353 (Div. 2) ABCDE 题解 python

    Problems     # Name     A Infinite Sequence standard input/output 1 s, 256 MB    x3509 B Restoring P ...

  5. 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解

    题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...

  6. 2016ACM青岛区域赛题解

    A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Jav ...

  7. poj1399 hoj1037 Direct Visibility 题解 (宽搜)

    http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...

  8. 网络流n题 题解

    学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...

  9. CF100965C题解..

    求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...

随机推荐

  1. pytest中怎么实现参数化?

    我们在组自动化项目中,肯定会遇到需要数据驱动的地方,在unittest中我们直接采用ddt的方式,但是在pytest中我们不能用ddt来实现,所以我们采用了参数化的方式来实现 那么具体怎么完成这个参数 ...

  2. v8

    V8 - 开源,由Google开发,用C ++编写 Rhin- 由Mozilla基金会开源,完全用Java开发 SpiderMonkey 第一个JavaScript引擎,Netscape Naviga ...

  3. vs在matlab生成的dll中,load treebigger错误

    提示的错误:Warning: Variable 'model' originally saved as a TreeBagger cannot be instantiated as an object ...

  4. shiro 基本知识测试

    shiro 基本知识测试 <!--shiro核心包--> <dependency> <groupId>org.apache.shiro</groupId> ...

  5. chrom代理插件

    1.SwitchyOmega 2.Resource Override

  6. testcontainers 方便的db测试框架

    testcontainers是一个强大,简单,基于容器的db测试解决方案 目前已经支持了主流的开发语言 参考资料 https://github.com/testcontainers/testconta ...

  7. string方法介绍

    #_*_coding:utf-8_*_#作者:王佃元#日期:2019/12/9#string操作print('hello'*2) #乘法操作,输出对应次数print('helloworld'[2:]) ...

  8. Philosopher(set 线段树合并)

    直接维护乘积是肯定不可行的, 精度会爆炸, 于是我们来维护对数的和, 最后来计算最高位即可 那么转换成区间求和, 区间排序 区间排序的方式可以采用线段树维护最大递增块来解决,外层用set来维护线段树的 ...

  9. Linux性能优化实战学习笔记:第九讲

    一.中断的魅力 1.中断在生活的魅力 比如你订了一份外卖,但是不确定外卖什么时候送到,也没有别的方法了解外卖的进度,但是,配送员送外卖是不等人的,到了你这儿没人取的话,就直接走人了.所以你指能苦苦等着 ...

  10. 带有连接池的Http客户端工具类HttpClientUtil

    一.背景 业务开发中,经常会遇到通过http/https向下游服务发送请求.每次都要重复造轮子写HttpClient的逻辑,而且性能.功能参差不齐.这里分享一个高性能的.带连接池的通用Http客户端工 ...