Solution -「COCI 2014-2015 #2」「洛谷 P6406」Norma
\(\mathcal{Description}\)
Link.
给定 \(\{a_n\}\),求:
\]
答案对 \(10^9\) 取模。
\(\mathcal{Solution}\)
挺可爱的一道题 w。
静态序列计数问题,可以考虑分治:对于 \([l,r]\)(\(l<r\)),令分割点 \(p=\lfloor\frac{l+r}2\rfloor\),对满足左端点在 \([l,p]\),右端点在 \((p,r]\) 的区间进行计算,然后递归两个区间继续求解。
对于本题,可以枚举区间左端点 \(i=p,p-1,\cdots,l\),记 \(s=\min_{u=i}^p\{a_u\}\),\(t=\max_{u=i}^p\{a_u\}\),\(j=\max_{u\in(p,r]}\{u|\min_{v=i}^j\{a_v\}=s\}\),\(k=\max_{u\in(p,r]}\{u|\max_{v=i}^j\{a_v\}=t\}\),只讨论 \(j\le k\),发现对于将要计数的区间 \([i,x]\),分三种情况:
\(x\in(p,j]\):最小值、最大值都在 \([l,p]\) 中取到,那么只需要关心 \(x\) 的位置,故此情况对答案的贡献为:
\[pq\sum_{x\in(p,j]}(x-i+1)
\]\(x\in(j,k]\):最大值在 \([l,p]\) 中取到,而最小值会在 \((p,x]\) 中取到,则要求:
\[q\sum_{x\in(j,k]}(x-i+1)\min_{u=p+1}^x\{a_u\}
\]求和内是一个距离 \(\times\) 权值的形式,尝试预处理出值的前缀和和下标 \(\times\) 值的前缀和,就能 \(\mathcal O(1)\) 计算了,这里略过,详见代码。
\(x\in(k,r]\):最小值、最大值都在 \((p,x]\) 中取到,类似地求:
\[\sum_{x\in(k,r]}(x-i+1)\min_{u=p+1}^x\{a_u\}\max_{u=p+1}^x\{a_u\}
\]仍然预处理出两种前缀和,\(\mathcal O(1)\) 计算。
综上,在分治的 \(\mathcal O(n\log n)\) 的时间内解决本题。
\(\mathcal{Code}\)
/* Clearink */
#include <cstdio>
inline int rint () {
int x = 0; char s = getchar ();
for ( ; s < '0' || '9' < s; s = getchar () );
for ( ; '0' <= s && s <= '9'; s = getchar () ) x = x * 10 + ( s ^ '0' );
return x;
}
const int MAXN = 5e5, MOD = 1e9;
int n, a[MAXN + 5], ans;
int smx[MAXN + 5], vmx[MAXN + 5];
int smn[MAXN + 5], vmn[MAXN + 5];
int sxn[MAXN + 5], vxn[MAXN + 5];
inline void chkmin ( int& a, const int b ) { b < a && ( a = b, 0 ); }
inline void chkmax ( int& a, const int b ) { a < b && ( a = b, 0 ); }
inline int mul ( const long long a, const int b ) { return a * b % MOD; }
inline int sub ( int a, const int b ) { return ( a -= b ) < 0 ? a + MOD : a; }
inline int add ( int a, const int b ) { return ( a += b ) < MOD ? a : a - MOD; }
inline void addeq ( int& a, const int b ) { ( a += b ) >= MOD && ( a -= MOD, 0 ); }
inline int sum ( const int l, const int r ) {
return r < l ? 0 : ( ( l + r ) * ( r - l + 1ll ) >> 1 ) % MOD;
}
inline void solve ( const int l, const int r ) {
if ( l == r ) return addeq ( ans, mul ( a[l], a[l] ) );
int mid = l + r >> 1;
smx[mid] = vmx[mid] = smn[mid] = vmn[mid] = sxn[mid] = vxn[mid] = 0;
for ( int i = mid + 1, mn = a[mid + 1], mx = a[mid + 1];
i <= r; ++i, chkmin ( mn, a[i] ), chkmax ( mx, a[i] ) ) {
vmx[i] = add ( vmx[i - 1], mx ),
smx[i] = add ( smx[i - 1], mul ( i - mid, mx ) );
vmn[i] = add ( vmn[i - 1], mn ),
smn[i] = add ( smn[i - 1], mul ( i - mid, mn ) );
vxn[i] = add ( vxn[i - 1], mul ( mn, mx ) ),
sxn[i] = add ( sxn[i - 1], mul ( i - mid, mul ( mn, mx ) ) );
}
for ( int i = mid, mn = a[i], mx = a[i], j = mid + 1, k = mid + 1;
i >= l; --i, chkmin ( mn, a[i] ), chkmax ( mx, a[i] ) ) {
for ( ; j <= r && mn <= a[j]; ++j );
for ( ; k <= r && a[k] <= mx; ++k );
int p = j < k ? j : k, q = j ^ k ^ p; // [mid+1,p),[p,q),[q,r].
addeq ( ans, mul ( mul ( mn, mx ), sum ( mid + 2 - i, p - i ) ) );
if ( j <= k ) { // max is constant while min will change.
addeq ( ans, mul ( mx, add (
mul ( mid - i + 1, sub ( vmn[k - 1], vmn[j - 1] ) ),
sub ( smn[k - 1], smn[j - 1] ) ) ) );
} else { // min is constant while max will change.
addeq ( ans, mul ( mn, add (
mul ( mid - i + 1, sub ( vmx[j - 1], vmx[k - 1] ) ),
sub ( smx[j - 1], smx[k - 1] ) ) ) );
}
addeq ( ans, add (
mul ( mid - i + 1, sub ( vxn[r], vxn[q - 1] ) ),
sub ( sxn[r], sxn[q - 1] ) ) );
}
solve ( l, mid ), solve ( mid + 1, r );
}
int main () {
n = rint ();
for ( int i = 1; i <= n; ++i ) a[i] = rint ();
solve ( 1, n );
printf ( "%d\n", ans );
return 0;
}
Solution -「COCI 2014-2015 #2」「洛谷 P6406」Norma的更多相关文章
- 「区间DP」「洛谷P1043」数字游戏
「洛谷P1043」数字游戏 日后再写 代码 /*#!/bin/sh dir=$GEDIT_CURRENT_DOCUMENT_DIR name=$GEDIT_CURRENT_DOCUMENT_NAME ...
- Solution -「JSOI 2019」「洛谷 P5334」节日庆典
\(\mathscr{Description}\) Link. 给定字符串 \(S\),求 \(S\) 的每个前缀的最小表示法起始下标(若有多个,取最小的). \(|S|\le3\time ...
- Solution -「洛谷 P4372」Out of Sorts P
\(\mathcal{Description}\) OurOJ & 洛谷 P4372(几乎一致) 设计一个排序算法,设现在对 \(\{a_n\}\) 中 \([l,r]\) 内的元素排 ...
- Solution -「POI 2010」「洛谷 P3511」MOS-Bridges
\(\mathcal{Description}\) Link.(洛谷上这翻译真的一言难尽呐. 给定一个 \(n\) 个点 \(m\) 条边的无向图,一条边 \((u,v,a,b)\) 表示从 ...
- Solution -「APIO 2016」「洛谷 P3643」划艇
\(\mathcal{Description}\) Link & 双倍经验. 给定 \(n\) 个区间 \([a_i,b_i)\)(注意原题是闭区间,这里只为方便后文描述),求 \(\ ...
- Solution -「POI 2014」「洛谷 P5904」HOT-Hotels 加强版
\(\mathcal{Description}\) Link. 给定一棵 \(n\) 个点的树,求无序三元组 \((u,v,w)\) 的个数,满足其中任意两点树上距离相等. \(n\le1 ...
- Solution -「JLOI 2015」「洛谷 P3262」战争调度
\(\mathcal{Description}\) Link. 给定一棵 \(n\) 层的完全二叉树,你把每个结点染成黑色或白色,满足黑色叶子个数不超过 \(m\).对于一个叶子 \(u\), ...
- Solution -「CTS 2019」「洛谷 P5404」氪金手游
\(\mathcal{Description}\) Link. 有 \(n\) 张卡牌,第 \(i\) 张的权值 \(w_i\in\{1,2,3\}\),且取值为 \(k\) 的概率正比于 \ ...
- Solution -「CEOI 2017」「洛谷 P4654」Mousetrap
\(\mathscr{Description}\) Link. 在一个含 \(n\) 个结点的树形迷宫中,迷宫管理者菈米莉丝和一只老鼠博弈.老鼠初始时在结点 \(y\),有且仅有结点 \(x\ ...
随机推荐
- jsp使用${}语法,对应的spring后台使用Model
实现的前提,一定要在jsp文件顶部加上 <%@page isELIgnored="false" %> 如上图${aa}与${bb},是直接使用Model的key值对应的 ...
- Go语言系列之自定义实现日志库
日志库logo gitee地址传送门:https://gitee.com/zhangyafeii/logo 日志库需求分析 1. 支持往不同的地方输出日志 2. 日志分级别 Debug Trace I ...
- 新增访客数量MR统计之MR数据输出到MySQL
关注公众号:分享电脑学习回复"百度云盘" 可以免费获取所有学习文档的代码(不定期更新)云盘目录说明:tools目录是安装包res 目录是每一个课件对应的代码和资源等doc 目录是一 ...
- test_6 python的列表去重
1.使用内置函数set() set() 函数创建一个无序不重复元素集,可进行关系测试,删除重复数据,还可以计算交集.差集.并集等. 2.创建一个空的列表进行比较,把不重复的元素添加到新的列表中 #co ...
- 《设计模式面试小炒》策略和工厂模式替代业务场景中复杂的ifelse
<设计模式面试小炒>策略和工厂模式替代业务场景中复杂的ifelse 我是肥哥,一名不专业的面试官! 我是囧囧,一名积极找工作的小菜鸟! 囧囧表示:小白面试最怕的就是面试官问的知识点太笼统, ...
- 【pwn】攻防世界 pwn新手区wp
[pwn]攻防世界 pwn新手区wp 前言 这几天恶补pwn的各种知识点,然后看了看攻防世界的pwn新手区没有堆题(堆才刚刚开始看),所以就花了一晚上的时间把新手区的10题给写完了. 1.get_sh ...
- csapp lab2 拆炸弹
1. 实验内容 包含一个二进制应用bomb,需要根据该应用猜测程序的运行过程.程序主体包含了六个函数phase_1到phase_6,每个函数会根据用户的输入做出反应,当输入符合要求时,会炸弹拆解成功, ...
- 【刷题-LeetCode】191 Number of 1 Bits
Number of 1 Bits Write a function that takes an unsigned integer and return the number of '1' bits i ...
- Error: xz compression not available解决办法
centos6升级php时误安装: rpm -Uvh https://mirror.webtatic.com/yum/el7/epel-release.rpm rpm -Uvh https://mir ...
- Integer缓冲区相关问题--valueOf()方法
今天在学习过程中了解到一个现象,代码如下: Integer num1 = 100; Integer num2 = 100; System.out.println(num1==num2?true:fal ...