题目描述

  小$P$最近喜欢上了单调数列,他觉得单调的数列具有非常多优美的性质。经过小$P$复杂的数学推导,他计算出了一个单调增数列的艺术价值等于该数列中所有书的总和。并且以这个为基础,小$P$还可以求出任意一个数列的艺术价值,它等于将这个数列顺次划分若干个极长单调区间(相邻两个单调区间的单调性必须不相同)后,每个单调区间中元素总和的平均值。比如对于数列$3\ 7\ 9\ 2\ 4\ 5$,它将被划分为$[3\ 7\ 9]\ [2]\ [4\ 5]$,其艺术价值为$\frac{19+2+9}{3}=10$。  由于小$P$的审美观,他还要求划分出的第一个单调区间必须为单调增区间,也就是说,对于数列$10\ 9\ 8$,它将被划分为$[10]\ [9\ 8]$,而不是$[10\ 9\ 8]$。
  现在小$P$手里有一个长度为$n$的序列$\{a_i\}$,他想问你,这个序列的所有子序列中,艺术价值最大的是哪个子序列,输出其艺术价值。
  注意:本体单调数列为严格单调,也就是说数列中的数必须严格上升或严格下降。


输入格式

输入的第一行为一个正整数$n$,表示序列的长度。
接下来的一行为$n$个正整数$a_1,a_2,...,a_n$,代表序列中的$n$个数。


输出格式

输出仅有一个实数,表示子序列的最大艺术价值。
实数四舍五入,保留三位小数。


样例

样例输入1:

4
1 2 5 4

样例输出1:

8.000

样例输入2:

6
3 1 7 2 6 5

样例输出2:

10.500


数据范围与提示

对于第一个样例,最优的子序列为$1\ 2\ 5$,其价值为$8$。如果选择$5\ 4$的话,虽然和为$9$,但单调区间数为$2$,因为第一个单调区间必须为递增区间。
对于第二个样例,最优的子序列为$3\ 7\ 6\ 5$,其价值为$\frac{3+7+6+5}{2}=10.5$。


题解

首先,来看一条性质:一定存在一个最优子序列,它的单调区间长度不超过$2$。

有了这个性质,这道题就简单多了,维护三个树状数组就好了。

下面是某大神总结的树状数组的用途:

$1$:快速简洁求解前缀最值
$2$:快速简洁求解前缀和,单点修改区间和,区间修改单点值。
$3$:并不简洁的求解区间修改区间求和(但也比线段树简洁)
$4$:快速维护树上每个点到根的路径权值和
$5$:快速维护树上每个点的子树和(与$4$不兼容)
$6$:用来轻松愉快的套主席树,$Splay$等各类动态数据结构
$7$:维护带删除\插入操作的排名,不过需要二分,两个$log$不推荐,推荐线段树一个$log$

时间复杂度:$\Theta(n\log n)$。

期望得分:$100$分。

实际得分:$100$分。


代码时刻

#include<bits/stdc++.h>
using namespace std;
int n;
int a[100001],b[100001];
long long tr1[100001],tr2[100001],tr3[100001];
double ans;
int lowbit(int x){return x&-x;}
void add1(int x,long long w){for(int i=x;i;i-=lowbit(i))tr1[i]=max(tr1[i],w);}
void add2(int x,long long w){for(int i=x;i<=b[0];i+=lowbit(i))tr2[i]=max(tr2[i],w);}
void add3(int x,long long w){for(int i=x;i;i-=lowbit(i))tr3[i]=max(tr3[i],w);}
long long ask1(int x)
{
long long res=0;
for(int i=x;i<=b[0];i+=lowbit(i))res=max(res,tr1[i]);
return res;
}
long long ask2(int x)
{
long long res=0;
for(int i=x;i;i-=lowbit(i))res=max(res,tr2[i]);
return res;
}
long long ask3(int x)
{
long long res=0;
for(int i=x;i<=b[0];i+=lowbit(i))res=max(res,tr3[i]);
return res;
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
b[i]=a[i];
}
b[0]=n;
sort(b+1,b+n+1);
b[0]=unique(b+1,b+b[0]+1)-b-1;
for(int i=1;i<=n;i++)
a[i]=lower_bound(b+1,b+b[0]+1,a[i])-b+1;
b[0]++;
for(int i=n;i;i--)
{
long long w1=ask1(a[i]+1),w2=ask2(a[i]-1),w3=ask3(a[i]+1);
ans=max(ans,(1.0*w3+b[a[i]-1])/2.0);
ans=max(ans,(1.0*b[a[i]-1]+w2)/2.0);
ans=max(ans,1.0*w1+b[a[i]-1]);
add1(a[i],b[a[i]-1]+w1);
add2(a[i],b[a[i]-1]+w2);
add3(a[i],b[a[i]-1]+max(w2,w3));
}
printf("%.3lf",ans);
return 0;
}

rp++

[CSP-S模拟测试]:小P的单调数列(树状数组+DP)的更多相关文章

  1. 【CSP模拟赛】奇怪的队列(树状数组 &二分&贪心)

    题目描述 nodgd的粉丝太多了,每天都会有很多人排队要签名.  今天有n个人排队,每个人的身高都是一个整数,且互不相同.很不巧,nodgd今天去忙别的事情去了,就只好让这些粉丝们明天再来.同时nod ...

  2. 【洛谷】NOIP提高组模拟赛Day2【动态开节点/树状数组】【双头链表模拟】

    U41571 Agent2 题目背景 炎炎夏日还没有过去,Agent们没有一个想出去外面搞事情的.每当ENLIGHTENED总部组织活动时,人人都说有空,结果到了活动日,却一个接着一个咕咕咕了.只有不 ...

  3. 【bzoj4548】小奇的糖果 STL-set+树状数组

    题目描述 平面上有n个点,每个点有一种颜色.对于某一条线段,选择所有其上方或下方的点.求:在不包含所有颜色的点的前提下,选择的点数最多是多少.(本题中如果存在某颜色没有相应的点,那么选择任何线段都不算 ...

  4. 牛客小白月赛13-H(单调栈+树状数组)

    题目链接:https://ac.nowcoder.com/acm/contest/549/H 题意:给一个柱状图,包括每个矩阵的宽度和高度,求能组成的最大矩阵的面积. 思路:显然最大矩阵的高一定为n个 ...

  5. 2019.01.21 bzoj2441: [中山市选2011]小W的问题(树状数组+权值线段树)

    传送门 数据结构优化计数菜题. 题意简述:给nnn个点问有多少个www型. www型的定义: 由5个不同的点组成,满足x1<x2<x3<x4<x5,x3>x1>x2 ...

  6. 「模拟赛20180307」三元组 exclaim 枚举+树状数组

    题目描述 给定 \(n,k\) ,求有多少个三元组 \((a,b,c)\) 满足 \(1≤a≤b≤c≤n\)且\(a + b^2 ≡ c^3\ (mod\ k)\). 输入 多组数据,第一行数据组数\ ...

  7. bzoj 4826: [Hnoi2017]影魔【单调栈+树状数组+扫描线】

    参考:https://www.cnblogs.com/lcf-2000/p/6789680.html 这是一个相对码量少的做法,用到了区间修改区间查询的树状数组,详见:www.cnblogs.com/ ...

  8. 算法笔记求序列A每个元素左边比它小的数的个数(树状数组和离散化)

    #include <iostream> #include <algorithm> #include <cstring> using namespace std ; ...

  9. BZOJ4237 稻草人(分治+树状数组+单调栈)

    如果要询问的某个纵坐标为inf的点左边是否有点能与其构成所要求的矩形,只要用个单调栈就可以了.可以想到用分治来制造单调性. 按横坐标排序,每次考虑跨过分治中心的矩形.考虑右边的每个点能与左边的哪些点构 ...

随机推荐

  1. Vim文本编辑工具

    4文本编辑工具Vim Vim是vi的升级版,编辑文本时vi不会显示颜色而vim会显示颜色. 安装vim工具 #yum  install  –y  vim-enhanced Vim有三种模式:一般模式. ...

  2. centos7里创建用户和组

    1.创建组distro,其GID为2019groupadd -g 2019 distro 2.创建用户mandriva, 其ID号为1005:基本组为distro useradd mandriva - ...

  3. 2019JAVA第二次实验报告

    Java实验报告 班级 计算机科学与技术二班 学号 20188442 姓名 吴怡君 完成时间 2019/9/19 评分等级 实验二 Java简单类与对象 实验目的 掌握类的定义,熟悉属性.构造函数.方 ...

  4. 分享一个linux中测试网站是否正常的shell脚本

    #! /bin/bash #Author=Harry CheckUrl() { #<==定义函数,名字为CheckUrl timeout=5 #<==定义wget访问的超时时间,超时就退出 ...

  5. Mybatis-学习笔记(4)1对1、1对多、多对多

    1.1对1 有2种方式对内嵌Bean设值: 1>关联查询就一条语句.使用association关键字,直接将嵌套对象的映射表的字段赋值内嵌对象. <association property ...

  6. C++ 统计输入的句子有多少英文字母

    // ConsoleApplication1.cpp: 定义控制台应用程序的入口点.//#include "stdafx.h"#include <iostream>#i ...

  7. Redis集群,备份,哨兵机制

    原文:https://blog.csdn.net/zy345293721/article/details/87536144 1.集群        先来简单了解下redis中提供的集群策略, 虽然re ...

  8. while与格式化的练习

    练习 判断下列逻辑语句的结果,一定要自己先分析 1)1 > 1 or 3 < 4 or 4 > 5 and 2 > 1 and 9 > 8 or 7 < 6 2)n ...

  9. [LeetCode] 210. 课程表 II

    题目链接:https://leetcode-cn.com/problems/course-schedule-ii/ 题目描述: 现在你总共有 n 门课需要选,记为 0 到 n-1. 在选修某些课程之前 ...

  10. 第三方模块:gulp模块

    一.Gulp的使用 1. 使用npm install  gulp  下载gulp库文件 2. 在项目根目录下简历gulpfile.js文件 3. 重构项目的文件夹架构src目录放置源代码文件,dist ...