BZOJ 3745
题解:
分治好题
首先暴力显然rmq可以做到n^2
比较容易想到是以最值分治,这样在数据随机复杂度是nlogn,不随机还是n^2的
以最值分治只有做多与较小区间复杂度相同才是nlogn的
而这题里我们直接分治
想清楚再搞个暴力对拍还是比较好写的
代码:
#include <bits/stdc++.h>
using namespace std;
#define rint register int
#define IL inline
#define rep(i,h,t) for (rint i=h;i<=t;i++)
#define dep(i,t,h) for (rint i=t;i>=h;i--)
#define ll long long
#define mid ((h+t)>>1)
const int N=6e5;
const int mo=1e9;
int a[N],sum[N],sum2[N],sum3[N],sum4[N],sum5[N],sum6[N];
ll ans=;
void fz(int h,int t)
{
if (h==t)
{
ans=(ans+1ll*a[h]*a[h])%mo;
return;
}
int mina=a[mid],maxa=a[mid];
sum[mid+]=sum2[mid+]=sum3[mid+]=sum4[mid+]=;
dep(i,mid,h)
{
mina=min(a[i],mina);
maxa=max(a[i],maxa);
sum[i]=(sum[i+]+(1ll*mina*maxa)%mo*(mid-i+))%mo;
sum2[i]=(sum2[i+]+1ll*mina*maxa)%mo;
sum3[i]=(sum3[i+]+1ll*(mid-i+)*mina)%mo;
sum4[i]=(sum4[i+]+mina)%mo;
sum5[i]=(sum5[i+]+1ll*(mid-i+)*maxa)%mo;
sum6[i]=(sum6[i+]+maxa)%mo;
}
mina=a[mid],maxa=a[mid];
int pos1=mid,pos2=mid;
rep(i,mid,t)
{
mina=min(a[i],mina);
maxa=max(a[i],maxa);
while (a[pos1]>=mina&&pos1>=h) pos1--; pos1++;
while (a[pos2]<=maxa&&pos2>=h) pos2--; pos2++;
if (pos2<pos1)
{
ans+=(1ll*(*i-mid-pos1+)*(mid-pos1+)/)%mo*mina%mo*maxa%mo;
ans+=(sum3[pos2]-sum3[pos1]+1ll*(i-mid)*(sum4[pos2]-sum4[pos1]))%mo*maxa%mo;
ans+=(sum[h]-sum[pos2]+1ll*(i-mid)*(sum2[h]-sum2[pos2]))%mo;
ans%=mo;
} else
{
ans+=(1ll*(*i-mid-pos2+)*(mid-pos2+)/)%mo*mina%mo*maxa%mo;
ans+=(sum5[pos1]-sum5[pos2]+1ll*(i-mid)*(sum6[pos1]-sum6[pos2]))%mo*mina%mo;
ans+=(sum[h]-sum[pos1]+1ll*(i-mid)*(sum2[h]-sum2[pos1]))%mo;
ans%=mo;
}
}
if (h<=mid-) fz(h,mid-);
if (mid+<=t) fz(mid+,t);
}
int main()
{
freopen("1.in","r",stdin);
freopen("1.out","w",stdout);
ios::sync_with_stdio(false);
int n;
cin>>n;
rep(i,,n) cin>>a[i];
fz(,n);
cout<<ans<<endl;
return ;
}
BZOJ 3745的更多相关文章
- BZOJ 3745: [Coci2015]Norma(分治)
题意 给定一个正整数序列 \(a_1, a_2, \cdots, a_n\) ,求 \[ \sum_{i=1}^{n} \sum_{j=i}^{n} (j - i + 1) \min(a_i,a_{i ...
- bzoj 3745 [Coci2015]Norma——序列分治
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3745 如果分治,就能在本层仅算过 mid 的区间了. 可以从中间到左边地遍历左边,给右边两个 ...
- bzoj 3745: [Coci2015]Norma
Description Solution 考虑分治: 我们要统计跨越 \(mid\) 的区间的贡献 分最大值和最小值所在位置进行讨论: 设左边枚举到了 \(i\),左边 \([i,mid]\) 的最大 ...
- [BZOJ 3745] [COCI 2015] Norma
Description 给定一个正整数序列 \(a_1,a_2,\cdots,a_n\),求 \[ \sum_{i=1}^n\sum_{j=i}^n(j-i+1)\min(a_i,a_{i+1},\c ...
- 【刷题】BZOJ 3745 [Coci2015]Norma
Description Input 第1行,一个整数N: 第2~n+1行,每行一个整数表示序列a. Output 输出答案对10^9取模后的结果. Sample Input 4 2 4 1 4 Sam ...
- bzoj 3745: [Coci2015]Norma【分治】
参考:https://blog.csdn.net/lych_cys/article/details/51203960 真的不擅长这种-- 分治,对于一个(l,r),先递归求出(l,mid),(mid+ ...
- BZOJ 2127: happiness [最小割]
2127: happiness Time Limit: 51 Sec Memory Limit: 259 MBSubmit: 1815 Solved: 878[Submit][Status][Di ...
- BZOJ 3275: Number
3275: Number Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 874 Solved: 371[Submit][Status][Discus ...
- BZOJ 2879: [Noi2012]美食节
2879: [Noi2012]美食节 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 1834 Solved: 969[Submit][Status] ...
随机推荐
- vc++基础班[28]---动态数组及动态链表的讲解
C++中也有相应的动态数组.动态链表.映射表的模板类,就是STL中的:vector.list.map 他们属于C++标准中的一部分,对于程序的移植性来说也是不错的,但是在MFC编程中使用 CArray ...
- makefile实例
#.PHONY:cleanall cleanobj cleandiff #cleanall:cleandiff cleanobj # rm program #cleanobj: # rm obj.c ...
- hibernate框架学习之对象状态
lHibernate对象共有三种状态 •瞬时状态:瞬时对象 •持久化状态:持久化对象 •托管状态:托管对象 l瞬时对象(Transient Object),简称TO l瞬时对象指的是应用程序创建出来的 ...
- new-delete-malloc-free关系总结
new-delete-malloc-free关系总结 写在前面的话 这个系列的笔记总结是根据网上的两篇基础拓展而来的 C++经典面试题(最全,面中率最高) C++面试集锦( 面试被问到的问题 ) 面试 ...
- 深入学习NAT工作原理
深入学习NAT工作原理 我们单位中的电脑很多,组成了一个局域网,网络中只有一个电脑和外网(Internet)相连,当然有一个外网地址,但仅仅一个.我以前一直不明白,我们局域网的电脑均能上网,几台同时上 ...
- socket原理-基于UDP的套接字-基于TCP实现远程命令-粘包
- 对mysql数据库中字段为空的处理
数据库中字段为空的有两种:一种为null,另一种为空字符串.null代表数值未知,空字符串是有值得,只是为空.有时间我们想把数据库中的数据以excel形式导出时 如果碰到字段为空的,为空的字段会被后面 ...
- nginx代理跨域(mac)
首先找到nginx.conf文件,修改并添加如下配置 html 文件 <!DOCTYPE html> <html lang="en"> <head&g ...
- mysql 视图 触发器 事物 存储过程 函数 流程控制
1.视图 *** 视图是有一条sql语句的查询结果构成的虚拟表 其不是物理存在的 使用方式与普通表相同 视图的作用1.简化sql语句的编写 2.限制可以查看的数据 可以使用权限来完成 权限某一个库 的 ...
- STL的基本操作指令
list :Lists将元素按顺序储存在链表中. 与 向量(vectors)相比, 它允许快速的插入和删除,但是随机访问却比较慢. assign() 给list赋值 back() 返回最后一个元素 b ...