P3648-[APIO2014]序列分割【斜率优化】

正题

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


题目大意

nnn个数字的序列,分割kkk次,每次的权值是左右两块数字的乘积。求最大权值和分割方案。


解题思路

显然分割顺序不会影响结果,一个分割方式的答案是每一块与其他块的乘积之和。

考虑dpdpdpfi,jf_{i,j}fi,j表示第iii次分割,到第jjj个时的方案数,有转移
fi,j=max{fi−1,k+(sj−sk)∗sk}f_{i,j}=max\{f_{i-1,k}+(s_j-s_k)*s_k\}fi,j=max{fi1,k+(sjsk)sk}
考虑两个方案k′<kk'<kk<k的优劣性
fi−1,k+sj∗sk−sk2>fi−1,k′+sj∗sk′−sk′2f_{i-1,k}+s_j*s_k-s_k^2>f_{i-1,k'}+s_j*s_{k'}-s_{k'}^2fi1,k+sjsksk2>fi1,k+sjsksk2
fi−1,k+sj∗sk−fi−1,k′+sk′2>sj∗(sk′−sk)f_{i-1,k}+s_j*s_k-f_{i-1,k'}+s_{k'}^2>s_j*(s_{k'}-s_{k})fi1,k+sjskfi1,k+sk2>sj(sksk)
fi−1,k+sj∗sk−fi−1,k′+sk′2sk′−sk<sj\frac{f_{i-1,k}+s_j*s_k-f_{i-1,k'}+s_{k'}^2}{s_{k'}-s_{k}}<s_jskskfi1,k+sjskfi1,k+sk2<sj
因为sjs_jsj递增,所以单调队列维护一个上突壳即可

然后记录一下前驱输出方案
时间复杂度O(nk)O(nk)O(nk)


codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=110000;
ll n,k,q[N],fa[210][N],s[N],f[2][N];
double slope(ll z,ll x,ll y){if(s[y]==s[x])return -1e18;return (f[z][x]-f[z][y]-s[x]*s[x]+s[y]*s[y])/(double)(s[y]-s[x]);
}
int main()
{scanf("%lld%lld",&n,&k);for(ll i=1;i<=n;i++)scanf("%lld",&s[i]),s[i]+=s[i-1];memset(f,0,sizeof(f));for(ll i=1;i<=k;i++){ll head=1,tail=0;memset(f[i&1],0,sizeof(f[i&1]));for(ll j=1;j<=n;j++){while(head<tail&&slope(~i&1,q[head],q[head+1])<=s[j])head++;f[i&1][j]=0;if(head<=tail){ll x=q[head];fa[i][j]=x;f[i&1][j]=f[~i&1][x]+(s[j]-s[x])*s[x];}while(head<tail&&slope(~i&1,q[tail],q[tail-1])>=slope(~i&1,j,q[tail]))tail--;q[++tail]=j;}}printf("%lld\n",f[k&1][n]);for(ll i=k,w=fa[k][n];i;w=fa[--i][w])printf("%lld ",w);
}

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

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

相关文章

codeforces E. Jamie and Tree LCA+dfs序+线段树

题解&#xff1a; 写起来还稍微有点麻烦。 dfs序线段树可以维护子树的整体修改和查询。 因此&#xff0c;这道题我们要往子树上靠。 我们首先从1号点进行dfs遍历&#xff0c;顺便求出点的dfs序和深度&#xff0c;然后我们采用倍增的思想&#xff0c;可以预处理出每个点的祖先…

async/await 的基本实现和 .NET Core 2.1 中相关性能提升

前言这篇文章的开头&#xff0c;笔者想多说两句&#xff0c;不过也是为了以后再也不多嘴这样的话。在日常工作中&#xff0c;笔者接触得最多的开发工作仍然是在 .NET Core 平台上&#xff0c;当然因为团队领导的开放性和团队风格的多样性&#xff08;这和 CTO 以及主管的个人能…

【DP】砝码称重 (ssl 1072)

砝码称重 ssl 1072 题目大意&#xff1a; 有6种砝码&#xff0c;每种的个数分别为a1,a2,a3,a4,a5,a6,请问可以平出多少种重量(0除外) Description 设有1g、2g、3g、5g、10g、20g的砝码各若干枚&#xff08;其总重<1000&#xff09;&#xff0c; 要求&#xff1a; 输入…

P5025-[SNOI2017]炸弹【tarjan,线段树优化建图】

正题 题目链接:https://www.luogu.com.cn/problem/P5025 题目大意 .nnn个炸弹&#xff0c;每个在xxx位置处&#xff0c;范围为rrr。定义fif_ifi​表示第iii个炸弹爆炸能连锁的炸弹数&#xff0c;要求输出∑i1nfi∗i\sum_{i1}^nf_i*ii1∑n​fi​∗i 解题思路 将每个炸弹可以炸…

使用Swashbuckle构建RESTful风格文档

本次和大家分享的是Swagger to WebApi的nuget包Swashbuckle&#xff1b;因为项目需要统一api文档的风格&#xff0c;并要支持多种开发语言&#xff08;C#&#xff0c;java&#xff0c;python&#xff09;&#xff0c;所以首先想到的是swagger来构建api文档&#xff0c;本章讲解…

【dfs】【bfs】【链表】 求连通分量 (ssl 1759)

求连通分量 ssl 1759 题目大意 由n个点组成的无向图&#xff0c;求连通在一起的点数最大是多少 原题 求一个图的连通分量 Input n 顶点数(<100) 边 Output 连通分量 Sample Input 8 6 3 1 2 2 5 5 4 4 1 8 7 0 0 Sample Output 4 方法一&#xff08;dfs …

P2472-[SCOI2007]蜥蜴【网络流】

正题 题目链接:https://www.luogu.com.cn/problem/P2472 题目大意 n∗mn*mn∗m个格子&#xff0c;每个格子的石柱高度不同&#xff0c;蜥蜴可以跳到距离不超过ddd的石柱处&#xff0c;并且先前所站的石柱高度减一&#xff0c;为0则不能站&#xff0c;然后求有多少只蜥蜴不可以…

发布 Rafy .NET Standard 版本 Nuget 包

去年年中&#xff0c;Rafy 框架的源码就已经支持了 Net Standard 2.0 版本。其开源代码也已经上传到 Github 中&#xff1a;https://github.com/zgynhqf/rafy/tree/NetStandard2.0 。但是这都只是在源码层面支持 NS2.0&#xff0c;并没有发布其正式的 Nuget 包。要使用这个版本…

一些来自STL的好东西

STL: 队列: 表达式作用#include&#xff1c;queue&#xff1e;定义queue&#xff1c;int&#xff1e;x定义一个int类型的队列&#xff0c;名为xx.push(y)从队列x的对尾插入yx.pop()队列x的队头出列gx.front()g等于x的队头hx.back()h等于x的队尾x.size()队列x的长度x.empty()判…

codeforces D.MADMAX 动态规划、记忆化搜索

题意 给出一个DAG&#xff0c;每条边上有权重(权重是小写字母的ASCII码)&#xff0c;现在两位同学A和B分别位于某两点上(可以相同)&#xff0c;其中A和B轮流走&#xff0c;但是每人所走的边权不能变小&#xff0c;走到不能走为止就输。 A先走&#xff0c;询问最后谁会赢。 题解…

GDOI2020游记

DAY0DAY0DAY0 早上上课一堆作业没写完(啊这 然后中午去竞赛室之前还没找到老师请假&#xff0c;还专门又跑回去拿了准考证(还白嫖了马老师红包 然后下午就去了&#xff0c;发现没带充电线 之后发现准考证没带&#xff0c;吓死了。还好酒店有可以打印的地方&#xff0c;还好搞定…

你关心才值得分享 | K8S网络安全之访问控制技术实践

(请允许我插播下广告&#xff0c;便于其它伙伴了解趣码 Cloud Coder)还是那句话&#xff0c;你关心才值得分享~最近的一起分享就在5.10本周四晚&#xff0c;精彩千万不要错过&#xff01;Hi&#xff0c;你是不是也曾觉得K8S&#xff08; Kubernetes &#xff09;网络安全话题范…

【dfs】【链表】连通图 (ssl 1758)

连通图 ssl 1758 题目大意 有一个由n个点组成的无向图&#xff0c;检测他是否联通 原题 判断一个图是否为一个边通图 Input n 顶点 (n<100) 边 Output 1 表示连通 0 表示不边通 Sample Input 5 1 2 2 3 5 4 0 0 Sample Output 0 解题方法 用dfs链表从1开始…

洛谷P1120小木棒 爆搜+剪枝

题解 暴搜的思路容易想到&#xff0c;但是剪枝细节有很多&#xff0c;数据很强。 搜索思路&#xff1a; a. 用dfs(left_num,left_len,bound)表示当前还需要拼left_num根木棒&#xff0c;当前正在拼的木棒还剩left_len长度&#xff0c;搜索是从大往小搜索&#xff0c;并且当前搜…

P3338-[ZJOI2014]力【FFT】

正题 题目链接:https://www.luogu.com.cn/problem/P3338 题目大意 Fj∑i1j−1qi∗qj(i−j)2−∑ij1nqi∗qj(i−j)2F_j\sum_{i1}^{j-1}\frac{q_i*q_j}{(i-j)^2}-\sum_{ij1}^n\frac{q_i*q_j}{(i-j)^2}Fj​i1∑j−1​(i−j)2qi​∗qj​​−ij1∑n​(i−j)2qi​∗qj​​ EjFjqjE_j…

从Xamarin.Essentials谈Xamarin库的封装

编者语&#xff1a;Xamarin在国内的推广还需要努力&#xff0c;其实这真的是移动端开发的一大福音&#xff0c;毕竟用一份代码的时间可以生成iOS/Android/Windows/Linux/macOS/Tizen多个平台&#xff0c;而且是原生的性能。Xamarin在Build 2018发布的新功能有Xamarin.Essential…

【最短路】【图论】【Floyed】牛的旅行(ssl 1119/luogu 1522)

牛的旅行 ssl 1119 luogu 1522 题目大意 有两堆点&#xff0c;每一堆点之中的任何两个点都一定有相连的路线&#xff0c;连接两堆点中的各一个点&#xff0c;使最远的两个点的距离最短 原题 农民John的农场里有很多牧区。有的路径连接一些特定的牧区。一片所有连通的牧区称…

codefoces 939E Maximize!好题

题解 若存在一个子集s满足答案的话&#xff0c;则该子集一定包含集合S的最大值。 反证法证明&#xff1a; 假设s集合中最大的元素为x&#xff0c;S集合中最大的元素为X。则如果把x换成X&#xff0c;最大值增加了X-x&#xff0c;而平均值增量一定不大于X-x。 这样的话&#xff0…

P2050-[NOI2012]美食节【费用流,动态连边】

正题 题目链接:https://www.luogu.com.cn/problem/P2050 题目大意 nnn个菜品mmm个厨师&#xff0c;第iii种菜需要pip_ipi​份&#xff0c;第iii个人做第jjj道菜需要时间ti,jt_{i,j}ti,j​&#xff0c;求最少等待时间和。 解题思路 这题和之前修车很像&#xff0c;数据变大了。…

用ASP.NET Core 2.0 建立规范的 REST API -- 预备知识

什么是RESTREST 是 Representational State Transfer 的缩写. 它是一种架构的风格, 这种风格基于一套预定义的规则, 这些规则描述了网络资源是如何定义和寻址的.一个实现了REST这些规则的服务就叫做RESTful的服务.最早是由Roy Fielding提出的.RPC 风格/getUsers/getUser?id1/c…