原题:

给你一个长度为1e5的序列ai,问你它的所有子序列的最大值与最小值之差的1000次方的和是多少

即∑_{p是a的子序列}(max{p}-min{p})^1000

这题难点在于(max-min)^1000这个东西外边套了次方,不太好处理

首先需要注意一点的性质是因为只考虑最大值和最小值,所以考虑子序列实际上是考虑集合的子集

那么可以对原序列排序,然后枚举两个端点,对应区间的包含端点的所有子集就是最大值和最小值分别为端点的子集

即此区间对答案的贡献为((a[j]-a[i])*2^(j-i-1))*1000

接下来处理(max-min)^1000

一个想法是差分转差为和,然而并没有什么卵子用

另一个想法是枚举区间时不枚举端点,而是枚举长度,这样可以把2的次幂提出来

但是仍旧无法把(max-min)^1000展开

最后的做法是暴力把上式展开,或者手玩小数据,可以发现规律

(涉及到区间问题时常用的一个思路是固定一个端点,考虑与另一个端点有关的变化情况)

所有左端点为1的区间的贡献和为

∑((a[j]-a[1])*2^(j-2))*1000

=2^1000*(a[2]-a[1])^1000+2^1001*(a[3]-a[1])^1000+...+2^(1000+n-1)*(a[n]-a[1])^1000

=2^1000*((a[2]-a[1])^1000+2*(a[3]-a[1])^1000+...+2^(n-1)*(a[n]-a[1])^1000)

把2^1000提出来,令S=(a[2]-a[1])^1000+2*(a[3]-a[1])^1000+...+2^(n-1)*(a[n]-a[1])^1000

暴力展开1000次方

=(a[2]^1000+1000*a[2]^999*a[1]+...+a[1]^1000)+2*(a[3]^1000+...+a[1]^1000)+...+2^(n-1)*(a[n]^1000+...+a[1]^1000)

合并同次项

=(a[2]^1000+2*a[3]^1000+4*a[4]^1000+...+2^(n-1)*a[n]^1000)

-1000*a[1]*(a[2]^999+2*a[3]^999+4*a[4]^999+...+2^(n-1)*a[n]^999)

+C(1000,2)*a[1]^2*(...)

...

+a[1]^1000*(1+2+4+...+2^(n-1))

发现规律了木有!

随着右端点下标增加,同次项的系数每次乘2

而a[j]^0,a[j]^1,a[j]^2,...,a[j]^1000可以预处理

那么如果从左往右枚举左端点i,(a[i+1]^k+2*a[i+2]^k+...+2^p*a[n]^k)可以快速由(a[i+2]^k+...+2^(p-1)*a[n]^k))乘2+a[i+1]^k得到

C可以预处理,a[i]^k也预处理好了,那么对于每一个左端点我们都可以1000次O(1)操作得到它对答案的贡献

耗时1e5*1e3,有点紧,要善用递推性质,盲目龟速乘会挂

还需要注意的一个问题是序列中的相同元素

实际上手玩小数据或者直接证明都可以发现,相同的数按照上述方法处理仍然能得到正确的结果

这个留给读者证明233

没有代码

【ecfinal2019热身赛】B题的更多相关文章

  1. 暑假集训——cf热身赛部分题有感加其题解

    刚刚开始集训,集训队队长暂时还没有拉专题,而是拉了部分codeforces上过题人数在2000左右的题组成了一场热身赛(其实就是一场练习),花了一天时间终于把它刷完了,其中很多题让我学到了很多骚操作, ...

  2. hdu 4706 Children's Day 2013年ICPC热身赛A题 模拟

    题意:按字母顺序排列成n型,简单的模拟题. 当字母排到z时从a开始重新排起. 代码: /* * Author: illuz <iilluzen[at]gmail.com> * Blog: ...

  3. hdu 4707 Pet 2013年ICPC热身赛A题 dfs水题

    题意:linji的仓鼠丢了,他要找回仓鼠,他在房间0放了一块奶酪,按照抓鼠手册所说,这块奶酪可以吸引距离它D的仓鼠,但是仓鼠还是没有出现,现在给出一张关系图,表示各个房间的关系,相邻房间距离为1,而且 ...

  4. hdu 4708 Rotation Lock Puzzle 2013年ICPC热身赛A题 旋转矩阵

    题意:给出一个n*n的矩阵,旋转每一圈数字,求出对角线可能的最大值,以及转到最大时的最小距离. 只要分析每一层就可以了,本来想用地址传递二维数组,发现行不通,改了一下就行了. 这里有个坑,比如: 1 ...

  5. hdu 4715 Difference Between Primes 2013年ICPC热身赛A题 素数水题

    题意:给出一个偶数(不论正负),求出两个素数a,b,能够满足 a-b=x,素数在1e6以内. 只要用筛选法打出素数表,枚举查询下就行了. 我用set储存素数,然后遍历set里面的元素,查询+x后是否还 ...

  6. ACM之最短路径做题笔记与记录

    在这里纪念一下从4月开始一直因为事情而荒废了的最短路,多亏了jbb的帮助,我才完成了FZU热身赛一题简单的一个用模拟链表存边以及最短路径的学习,目前(6.5)已经学会使用了最简单的djstral与sp ...

  7. 2017河工大校赛补题CGH and 赛后小结

    网页设计课上实在无聊,便开始补题,发现比赛时候僵着的东西突然相通了不少 首先,"追妹"这题,两个队友讨论半天,分好多种情况最后放弃(可是我连题目都没看啊),今天看了之后试试是不是直 ...

  8. PTA刷题笔记

    PTA刷题记录 仓库地址: https://github.com/Haorical/Code/tree/master/PTA/GPLT 两周之内刷完GPLT L2和L3的题,持续更新,包括AK代码,坑 ...

  9. 2014年acm亚洲区域赛·鞍山站

    今天北京赛站的比赛也结束了···看了一天的直播之后意识到鞍山站的比赛都过去了一个多月了···这一个月比较萎靡···整天都在睡觉写报告画工图中度过··· 鞍山比哈尔滨还是暖和很多的···就是山上有奇怪的 ...

随机推荐

  1. 浅出讲解:php的socket通信

    原文地址:https://www.cnblogs.com/aipiaoborensheng/p/6708963.html 对TCP/IP.UDP.Socket编程这些词你不会很陌生吧?随着网络技术的发 ...

  2. SMOS数据产品介绍与下载方法

    1. SMOS数据介绍 The Soil Moisture and Ocean Salinity (SMOS) 卫星是欧空局发射的一颗以探测地球土壤水含量以及海表盐度为目标的卫星,卫星所搭载的唯一载荷 ...

  3. go项目部署到linxu

    环境: 在mac上编译, 编译后上传到linux, 然后运行代码 go项目打包 一.直接部署到linux 1. 在mac上, 进入到项目目录, 执行以下命令, 进行编译: CGO_ENABLED=0 ...

  4. Celery—分布式的异步任务处理系统

    Celery 1.什么是Clelery Celery是一个简单.灵活且可靠的,处理大量消息的分布式系统 专注于实时处理的异步任务队列 同时也支持任务调度 Celery架构 Celery的架构由三部分组 ...

  5. protobuf 序列化 使用

    protocol buff是一种协议,是谷歌推出的一种序列化协议 .Java序列化协议也是一种协议 两者的目的是,将对象序列化成字节数组,或者说是二进制数据 导包 <dependency> ...

  6. springboot项目在IDEA根据不同的开发人员读取不同的配置文件

    IDEA启动项目打开项目的配置文件,修改Program argument为--spring.profiles.active=developerName启动项目,即可 命令行方式启动项目 java -j ...

  7. AtCoder Grand Contest 040 A - ><

    传送门 对于某个位置,只要知道这个位置往左最多的连续 $\text{<}$ 的数量 $x$ 和往右最多的连续 $\text{>}$ 的数量 $y$ 那么这个位置最小可能的数即为 $max( ...

  8. HDU4641 K-string(后缀自动机+线段树合并)

    先考虑没有动态加字符怎么做.计算每个节点的贡献,当|right|>=k时将len-lenfa计入即可. 动态加字符后,这个东西难以用LCT维护.于是考虑离线.建完SAM后,容易发现每个节点在时间 ...

  9. ELK 日志收集系统

    传统系统日志收集的问题 在传统项目中,如果在生产环境中,有多台不同的服务器集群,如果生产环境需要通过日志定位项目的Bug的话,需要在每台节点上使用传统的命令方式查询,这样效率非常底下. 通常,日志被分 ...

  10. javascript的装载和执行

    前言 为什么要采用js来create一个script标签,设置src然后append到head中,而不是直接使用script标签,这样不是更简单点吗? javascript的装载和执行 首先,我想说一 ...