[BZOJ3745][COCI2015]Norma[分治]
题意
分析
考虑分治,记当前分治区间为 \(l,r\) 。
枚举左端点,然后发现右端点无非三种情况:
- 极大极小值都在左边;
- 有一个在左边;
- 极大极小值都在右边;
考虑递推 \(l\) 的同时递推最靠右的满足最大最小值在左边的位置 \(p_1,p_2\).
根据不同的情况计数即可,注意计算以 \(\rm mid\) 作为右端点的情况。
总时间复杂度为\((nlogn)\)。
重点:分治算法降低复杂度的原因:根据极大极小值的不同划分数据以及?
代码
#include<bits/stdc++.h>
using namespace std;
#define go(u) for(int i=head[u],v=e[i].to;i;i=e[i].last,v=e[i].to)
#define rep(i,a,b) for(int i=a;i<=b;++i)
#define pb push_back
#define For for(int j=1;j<=6;++j)
typedef long long LL;
inline int gi(){
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
while(isdigit(ch)){x=(x<<3)+(x<<1)+ch-48;ch=getchar();}
return x*f;
}
template<typename T>inline bool Max(T &a,T b){return a<b?a=b,1:0;}
template<typename T>inline bool Min(T &a,T b){return b<a?a=b,1:0;}
const int N=5e5 + 7;
const LL mod=1e9;
int n;
LL a[N],s[7][N],ans;
void add(LL &a,LL b){a+=b;if(a>=mod) a-=mod;}
LL s1(int x){
return 1ll*x*(x+1)/2%mod;
}
void fz(int l,int r){
if(l==r){ add(ans,a[l]*a[l]%mod);return ;}
if(l>r) return;
int mid=l+r>>1;LL mx=a[mid],mi=a[mid];
For s[j][mid-1]=0;
for(int i=mid;i<=r;++i){
Min(mi,a[i]);Max(mx,a[i]);
For s[j][i]=s[j][i-1];
add(s[1][i],mi*mx%mod);
add(s[2][i],(i-mid)*mi%mod*mx%mod);
add(s[3][i],mi);
add(s[4][i],mx);
add(s[5][i],(i-mid)*mi%mod);
add(s[6][i],(i-mid)*mx%mod);
}
mx=mi=a[mid];
int p1=mid,p2=mid;
for(int i=mid;i>=l;--i){
Min(mi,a[i]);Max(mx,a[i]);
for(;p1+1<=r&&a[p1+1]>=mi;++p1);
for(;p2+1<=r&&a[p2+1]<=mx;++p2);
int k1=min(p1,p2),k2=max(p1,p2);
add(ans,(s1(k1-i+1)-s1(mid-i)+mod)*mi%mod*mx%mod);
if(k1==p1){
add(ans,(s[5][k2]-s[5][k1]+mod)*mx%mod);
add(ans,(s[3][k2]-s[3][k1]+mod)*(mid-i+1)%mod*mx%mod);
}else{
add(ans,(s[6][k2]-s[6][k1]+mod)*mi%mod);
add(ans,(s[4][k2]-s[4][k1]+mod)*(mid-i+1)%mod*mi%mod);
}
add(ans,(s[2][r]-s[2][k2]+mod)%mod);
add(ans,(s[1][r]-s[1][k2]+mod)*(mid-i+1)%mod);
}
fz(l,mid-1);fz(mid+1,r);
}
int main(){
n=gi();
rep(i,1,n) a[i]=gi();
fz(1,n);
printf("%lld\n",ans);
return 0;
}
[BZOJ3745][COCI2015]Norma[分治]的更多相关文章
- bzoj3745: [Coci2015]Norma 分治,单调队列
链接 bzoj 思路 首先\(\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{n}\sum\limits_{k=i}^{j}max(a_k)\)可以用单调队列求解.参见 ...
- 【BZOJ3745】Norma [分治]
Norma Time Limit: 20 Sec Memory Limit: 64 MB[Submit][Status][Discuss] Description Input 第1行,一个整数N: ...
- bzoj3745: [Coci2015]Norma
Description Input 第1行,一个整数N: 第2~n+1行,每行一个整数表示序列a. Output 输出答案对10^9取模后的结果. 预处理每个位置的数作为最小/大值向左延伸的最大距离, ...
- 【BZOJ3745】[Coci2015]Norma cdq分治
[BZOJ3745][Coci2015]Norma Description Input 第1行,一个整数N: 第2~n+1行,每行一个整数表示序列a. Output 输出答案对10^9取模后的结果. ...
- 【BZOJ3745】Norma(CDQ分治)
[BZOJ3745]Norma(CDQ分治) 题面 BZOJ 洛谷 题解 这种问题直接做不好做,显然需要一定的优化.考虑\(CDQ\)分治. 现在唯一需要考虑的就是跨越当前中间节点的所有区间如何计算答 ...
- 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 ...
- BZOJ3745:[COCI2015]Norma
浅谈离线分治算法:https://www.cnblogs.com/AKMer/p/10415556.html 题目传送门:https://lydsy.com/JudgeOnline/problem.p ...
- bzoj 3745 [Coci2015]Norma——序列分治
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3745 如果分治,就能在本层仅算过 mid 的区间了. 可以从中间到左边地遍历左边,给右边两个 ...
- bzoj 3745: [Coci2015]Norma【分治】
参考:https://blog.csdn.net/lych_cys/article/details/51203960 真的不擅长这种-- 分治,对于一个(l,r),先递归求出(l,mid),(mid+ ...
随机推荐
- 使用最新版SDWebImage
使用最新版SDWebImage 1. 下载源码: 2. 测试能否编译成功: 3. 用Xcode6新建一个工程,然后将文件夹拖入到工程当中: 4. 查看其主要的源码,发现之前使用版本的方法都被弃用了: ...
- Custom Settings.ini 和 bootstrap.ini 配置
[Settings]Priority=DefaultProperties=MyCustomProperty [Default] ;SkipWizard=YES 如果跳过部署向导,则即使 SkipCap ...
- HTML基础标签的综合应用案例(颜色、斜体、加粗、下划线、a标签、无序列表、有序列表)
什么是HTML l HTML(HyperText Mark-up Language)即超文本标记语言或超文本标签语言. l 何为超文本:“超文本”可以实现页面内可以包含图片.链接,甚至音乐.程序等. ...
- 高性能网站架构缓存——redis集群
相信你已经对redis有一定的了解,并能够安装上,进行简单的使用了,但是在咱们的实际应用中,使用redis肯定不会使用单机版,不光是redis不能使用单机版,其他的也不会使用,所以今天我们来说一下re ...
- Memorize and recite an important historical speech
Memorize and recite an important historical speech memorize['memәraiz]v.[亦作memorise] 记住, 记忆 historic ...
- docker 部署django项目(nginx + uwsgi +mysql)
最近在学习用docker部署Django项目,经过百折不挠的鼓捣,终于将项目部署成功,爬过好多坑,也发现很多技能需要提高.特此写下随笔与小伙伴们分享,希望能对大家有所启发. docker的理论我就不赘 ...
- 谈谈 C++ 中的右值引用
转自:https://liam0205.me/2016/12/11/rvalue-reference-in-Cpp/ 最近在改 XGBoost 的代码.XGBoost 在代码中使用了很多来自 C++1 ...
- zookeeper for windows
1,安装java环境,安装完java之后,添加用户变量.并在命令行中测试 java -version . 2,下载zookeeper http://mirrors.shu.edu.cn/apache/ ...
- MySQL主从.md
MySQL Replication 概述 Mysql内建的复制功能是构建大型,高性能应用程序的基础.将Mysql的数据分布到多个系统上去,这种分布的机制,是通过将Mysql的某一台主机的数据复制到其它 ...
- TensorFlow函数(八)tf.control_dependencies()
tf.control_dependencies(control_inputs) 此函数指定某些操作执行的依赖关系 返回一个控制依赖的上下文管理器,使用 with 关键字可以让在这个上下文环境中的操作都 ...