P3515-[POI2011]Lightning Conductor【整体二分,决策单调性】

正题

题目链接:https://www.luogu.com.cn/problem/P3507


题目大意

nnn个数字的一个序列aaa,对于每个位置iii求一个pip_ipi使得对于任意jjj满足
pi+ai−∣i−j∣≥pjp_i+a_i-\sqrt{|i-j|}\geq p_jpi+aiijpj


解题思路

化简一下发现我们是需要求出max{∣i−j∣+pj}max\{\sqrt{|i-j|}+p_j\}max{ij+pj}

分成两次去掉绝对值。
因为这个根号的性质是增长的越来越小,那么对于一个位置iii若它的maxmaxmax值位置为jjj,那么i+1i+1i+1就一定不小于jjj

利用这个单调性来优化,我们每次直接对于区间正中间midmidmid暴力求出它的答案pospospos,那么[l,mid−1][l,mid-1][l,mid1]的答案就在[L,pos][L,pos][L,pos],而[mid+1,r][mid+1,r][mid+1,r]的答案就在[pos,R][pos,R][pos,R]

然后递归下去就好了。时间复杂度O(nlog⁡n)O(n\log n)O(nlogn)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<stack>
#define ll long long
using namespace std;
const ll N=5e5+10;
ll n;double a[N],f[N],sqr[N];
stack<ll> s;
double count(ll i,ll j)
{return a[j]+sqr[abs(j-i)];}
void CDQ(ll l,ll r,ll L,ll R){if(l>r)return;ll mid=(l+r)>>1,pos=L;double tmp=count(mid,L);for(int i=L+1;i<=R&&i<=mid;i++)if(count(mid,i)>tmp)pos=i,tmp=count(mid,i);f[mid]=max(tmp,f[mid]);CDQ(l,mid-1,L,pos);CDQ(mid+1,r,pos,R);return;
}
signed main()
{scanf("%lld",&n);for(ll i=1;i<=n;i++){scanf("%lf",&a[n-i+1]);sqr[i]=sqrt((double)i);}CDQ(1,n,1,n);for(ll i=1;n-i+1>i;i++)swap(a[i],a[n-i+1]),swap(f[i],f[n-i+1]);CDQ(1,n,1,n);for(ll i=1;i<=n;i++)printf("%lld\n",(ll)ceil(f[i]-a[i]));return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/319864.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

关于二项式

组合数 (nm)\dbinom{n}{m}(mn​)&#xff1a;从 nnn 个物品中选出 mmm 个的方案数。 (nm)n!m!(n−m)!nm‾m!\dbinom{n}{m}\dfrac{n!}{m!(n-m)!}\dfrac{n^{\underline{m}}}{m!}(mn​)m!(n−m)!n!​m!nm​​ &#xff08;这个式子只依靠经典的组合意义&#xff0c;所以只在 0≤m≤…

dotnet core高吞吐Http api服务组件FastHttpApi

简介是dotNet core下基于Beetlex实现的一个高度精简化和高吞吐的HTTP API服务开源组件&#xff0c;它并没有完全实现HTTP SERVER的所有功能&#xff0c;而是只实现了在APP和WEB中提供数据服务最常用两个指令GET/SET&#xff0c;满足在应用实现JSON,PROTOBUF和MSGPACK等基于HTTP…

Poj 1284 Primitive Roots

文章目录Description题意&#xff1a;题解&#xff1a;代码&#xff1a;Poj 1284Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 6485 Accepted: 3697Description We say that integer x, 0 < x < p, is a primitive root modulo odd prime p if and onl…

P3507-[POI2010]GRA-The Minima Game【dp,博弈论】

正题 题目链接:https://www.luogu.com.cn/problem/P3507 题目大意 nnn个数&#xff0c;没人轮流取若干个并获得取走的数中最小数的权值&#xff0c;两人的目标都是自己的权值−-−对方的权值最大&#xff0c;求先手的权值−-−后手的权值。 解题思路 肯定是从大往小取&#x…

【LCT】【树状数组】Matches Are Not a Child‘s Play(luogu CF1137F)

正题 luogu CF1137F 题目大意 定义一棵树的产出序列为依次删除权值最小的叶子节点的顺序 进行q此操作&#xff1a; 1.把一个点的权值改为当前树中的最大权值1 2.查询一个点在删除序列中的位置 3.给出两个点&#xff0c;查询哪个在删除序列中的位置更前 解题思路 假设已经求出…

[XSY4220] 九万步(结论)

手动构造发现 x6x6x6 时是可行的 -101234567xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx

一起开心2020暑假训练第一周

hdu 1576 A/B oj传送 题解&#xff1a; Poj 1061 青蛙的约会 oj传送 题解&#xff1a; hdu 1525 Euclid‘s Game oj传送 题解&#xff1a; Poj 3070 Fibonacci oj传送 题解&#xff1a; HDU 2504 又见GCD oj传送 题解&#xff1a; Poj 1284 Primitive Roots oj传送 …

ASP.NET Core 中的中间件

前言由于是第一次写博客,如果您看到此文章,希望大家抱着找错误、批判的心态来看。 sky!何为中间件?在 ASP.NET Framework 中应该都知道请求管道。可参考&#xff1a;浅谈 ASP.NET 的内部机制 系列&#xff0c;个人感觉超详细。题外话&#xff1a;说到请求管道&#xff0c;就想…

【bfs】重力球(luogu 7473/NOI Online 2021 普及组 T3)

正题 luogu 7473 题目大意 给出一个正方形区域&#xff0c;中间有一些障碍 现在有两个球&#xff0c;每次操作可以使两个球同时向一个方向移动&#xff0c;直到遇到障碍或边界 现在问你让两个球到同一个位置最少要多少步 解题思路 对于每次操作&#xff0c;球只有可能停在障…

CF700E-Cool Slogans【SAM,线段树合并,dp】

正题 题目链接:https://www.luogu.com.cn/problem/CF700E 题目大意 给出一个字符串SSS&#xff0c;求一个最大的kkk使得存在kkk个字符串其中s1s_1s1​是SSS的子串&#xff0c;si1s_{i1}si1​在sis_isi​中出现了至少222次。 解题思路 首先我们需要有两个结论 si1s_{i1}si1​…

[XXSY] 构树(prufer序列,树上连通块DP)

传送门 CayleyCayleyCayley公式&#xff1a;一个完全图有nn−2n^{n-2}nn−2棵无根生成树&#xff08;可用prufer序列证明&#xff09; 扩展CayleyCayleyCayley公式&#xff1a;被确定边分为大小为a1,a2,⋯,ama_1,a_2,\cdots, a_ma1​,a2​,⋯,am​的连通块&#xff0c;则有nm−…

.NET Core中的性能测试工具BenchmarkDotnet

背景介绍之前一篇博客中&#xff0c;我们讲解.NET Core中的CSV解析库&#xff0c;在文章的最后&#xff0c;作者使用了性能基准测试工具BenchmarkDotNet测试了2个不同CSV解析库的性能&#xff0c;本篇我们来详细介绍一下BenchmarkDotNet。原文链接&#xff1a;https://dotnetco…

序列求和

来源&#xff1a;牛客网 文章目录题目描述题解&#xff1a;代码&#xff1a;时间限制&#xff1a;C/C 1秒&#xff0c;其他语言2秒 空间限制&#xff1a;C/C 131072K&#xff0c;其他语言262144K 64bit IO Format: %lld题目描述 定义S(n) 12 22 … n2&#xff0c;输出S(n) …

【区间DP】摆渡线路(2017 特长生 T4)

题目大意 给出一个园&#xff0c;圆上有100个点&#xff0c;若干条弦&#xff0c;让你选择尽量多互不相交的弦&#xff08;点可以重合&#xff09; 解题思路 可以把圆展开成链&#xff0c;然后复制一遍 设fi,jf_{i,j}fi,j​为第i个位置到第j个位置的所选弦的数量 那么可以枚…

bzoj4589-Hard Nim【FWT】

正题 题目链接:https://darkbzoj.tk/problem/4589 题目大意 求有多少个长度为nnn的数列满足它们都是不大于mmm的质数且异或和为000。 解题思路 两个初始多项式F[0]1F[0]1F[0]1&#xff0c;G[prime≤m]1G[prime\leq m]1G[prime≤m]1&#xff0c;然后答案就是FxorGnF\ xor\ G^…

some useful tricks

异或题思考方向&#xff1a;01trie树&#xff0c;分治 2. 二分图最大匹配&#xff0c;最大独立集互相转换 3. Snow 4. Code 5. Code 6. Code 7. 题目 #include<iostream> #include<cstdio> using namespace std; const int N100010; int n,p,nxt[N],no[200]…

25大技术主题向您发出最后一次约【惠】邀请

一年一度的微软技术盛会即将在上海世博中心拉开大幕金秋十月&#xff0c;来自两岸三地的百名明星讲师将携 25 大技术主题&#xff0c;齐聚上海为您奉献一场无与伦比的技术视听盛宴您&#xff0c;准备好了吗&#xff1f;25大技术主题公开&#xff0c;不负期待姗姗来迟的5系25大技…

P4980-[模板]Pólya定理

正题 题目链接:https://www.luogu.com.cn/problem/P4980 题目大意 nnn个物品图上mmm种颜色&#xff0c;求在可以旋转的情况下本质不同的涂色方案。 解题思路 既然是群论基本题就顺便写一下刚刚了解到的相关知识把&#xff08;顺便消磨一下时间 一个群(G,)(G,\times )(G,)定义…

逆元的求法

逆元&#xff1a; 对于a和p&#xff0c;若 a * inv(a) % p ≡ 1&#xff0c;则称inv(a)为a%p的逆元。其中p为质数 逆元就是在mod下&#xff0c;不能直接除以一个数&#xff0c;而要乘以他的逆元 a * inv(a) 1 (mod p) x / a可以改成 x * inv(a) % p 文章目录方法一.扩展欧几里…

[CF995F] Cowmpany Cowmpensation(树形dp,拉格朗日插值)

树形DP&#xff1a; 设f[u][i]f[u][i]f[u][i]表示给uuu的子树分配工资&#xff0c;uuu点工资为iii的方案数 f[u][i]∏v∈sonu(∑j1if[v][j])f[u][i]\prod\limits_{v\in son_u}(\sum\limits_{j1}^{i}f[v][j])f[u][i]v∈sonu​∏​(j1∑i​f[v][j]) 前缀和优化&#xff1a; 设g[u…