「SDOI2016」征途
征途
Pine开始了从S地到T地的征途。
从S地到T地的路可以划分成\(n\)段,相邻两段路的分界点设有休息站。
Pine计划用\(m\)天到达T地。除第\(m\)天外,每一天晚上Pine都必须在休息站过夜。所以,一段路必须在同一天中走完。
Pine希望每一天走的路长度尽可能相近,所以他希望每一天走的路的长度的方差尽可能小。
帮助Pine求出最小方差是多少。
设方差是\(v\),可以证明,\(v\times m^2\)是一个整数。为了避免精度误差,输出结果时输出\(v\times m^2\)。
对于 \(100\%\) 的数据,\(1 \le n \le 3000\)
题解
先推导一波方差。
=\frac 1m\left(m\overline{v}^2-2\overline{v}\sum_{i=1}^mv_i+\sum_{i=1}^m v_i^2\right)\\
=-\overline{v}^2+\frac 1m\sum_{i=1}^m v_i^2
\]
答案为
\]
所以要最小化的就是\(\sum_{i=1}^mv_i^2\)。
斜率优化
设\(f(i,j)\)表示前\(i\)天走了前\(j\)段路程的最小代价。
\]
化成一次函数式
\]
点坐标很明显,单调队列维护下凸包即可。
时间复杂度\(O(mn)\)。话说我怎么几次都把\(j\)写成\(i\),调了好久。
#include<bits/stdc++.h>
#define co const
#define il inline
template<class T> T read(){
T x=0,w=1;char c=getchar();
for(;!isdigit(c);c=getchar())if(c=='-') w=-w;
for(;isdigit(c);c=getchar()) x=x*10+c-'0';
return x*w;
}
template<class T>il T read(T&x){
return x=read<T>();
}
using namespace std;
#define int long long
co int N=3003;
int n,m,s[N];
int f[N][N],q[N];
int up(int i,int k1,int k2){
return f[i][k2]+s[k2]*s[k2]-f[i][k1]-s[k1]*s[k1];
}
int down(int k1,int k2){
return s[k2]-s[k1];
}
signed main(){
read(n),read(m);
for(int i=1;i<=n;++i) s[i]=s[i-1]+read<int>();
for(int i=1;i<=n;++i) f[1][i]=s[i]*s[i];
for(int i=2;i<=m;++i){
int l=0,r=0;
for(int j=1;j<=n;++j){
while(l<r&&2LL*s[j]*down(q[l],q[l+1])>=up(i-1,q[l],q[l+1])) ++l;
int k=q[l];
f[i][j]=f[i-1][k]+(s[j]-s[k])*(s[j]-s[k]);
while(l<r&&up(i-1,q[r],j)*down(q[r-1],q[r])<=up(i-1,q[r-1],q[r])*down(q[r],j)) --r;
q[++r]=j; // edit 1:i->j
}
}
printf("%lld\n",m*f[m][n]-s[n]*s[n]);
return 0;
}
导数零点二分
这道题更像一个wqs二分的模板题,我现在把它叫做导数零点二分。
为什么呢?推荐FlashHu的博客。
考虑\(最优解-段数\)的函数图像,给每一段追加的权值相当于左右移动导函数的零点。我们要做的就是让零点恰好是我们要求的段数。
时间复杂度\(O(n \log s_n)\)。
#include<bits/stdc++.h>
#define co const
#define il inline
template<class T> T read(){
T x=0,w=1;char c=getchar();
for(;!isdigit(c);c=getchar())if(c=='-') w=-w;
for(;isdigit(c);c=getchar()) x=x*10+c-'0';
return x*w;
}
template<class T>il T read(T&x){
return x=read<T>();
}
using namespace std;
#define int long long
co int N=3003;
int n,m,s[N];
int f[N],g[N];
int q[N],l,r;
int up(int a,int b){
return f[b]+s[b]*s[b]-f[a]-s[a]*s[a];
}
int down(int a,int b){
return s[b]-s[a];
}
void solve(int w){
q[l=r=0]=0;
for(int i=1;i<=n;++i){
while(l<r&&2*s[i]*down(q[l],q[l+1])>up(q[l],q[l+1])) ++l; // edit 1:>=
f[i]=f[q[l]]+w+(s[i]-s[q[l]])*(s[i]-s[q[l]]),g[i]=g[q[l]]+1;
while(l<r&&up(q[r],i)*down(q[r-1],q[r])<up(q[r-1],q[r])*down(q[r],i)) --r;
q[++r]=i;
}
}
signed main(){
read(n),read(m);
for(int i=1;i<=n;++i) s[i]=s[i-1]+read<int>();
int l=0,r=s[n]*s[n]/m,ans;
while(l<r){
int mid=(l+r)>>1;
solve(mid);
if(g[n]>m) l=mid+1;
else ans=m*(f[n]-mid*m)-s[n]*s[n],r=mid;
}
printf("%lld\n",ans);
return 0;
}
「SDOI2016」征途的更多相关文章
- 「SDOI2016」征途 题解
「SDOI2016」征途 先浅浅复制一个方差 显然dp,可以搞一个 \(dp[i][j]\)为前i段路程j天到达的最小方差 开始暴力转移 \(dp[i][j]=min(dp[k][j-1]+?)(j- ...
- 【LOJ】#2035. 「SDOI2016」征途
题解 有人管它叫带权二分,有人管它叫dp凸优化,有人管它叫wqs二分-- 延伸出来还有zgl分治,xjp¥!%#!@#¥!# 当我没说 我们拆个式子,很容易发现所求的就是 \(m\sum_{i = 1 ...
- loj2035 「SDOI2016」征途
学了斜率优化这题就能一气呵成地做出来啦qwqqwq #include <iostream> #include <cstdio> using namespace std; typ ...
- 「SDOI2016」储能表(数位dp)
「SDOI2016」储能表(数位dp) 神仙数位 \(dp\) 系列 可能我做题做得少 \(QAQ\) \(f[i][0/1][0/1][0/1]\) 表示第 \(i\) 位 \(n\) 是否到达上界 ...
- [LOJ 2070] 「SDOI2016」平凡的骰子
[LOJ 2070] 「SDOI2016」平凡的骰子 [题目链接] 链接 [题解] 原题求的是球面面积 可以理解为首先求多面体重心,然后算球面多边形的面积 求重心需要将多面体进行四面体剖分,从而计算出 ...
- liberOJ #2033. 「SDOI2016」生成魔咒 后缀数组
#2033. 「SDOI2016」生成魔咒 题目描述 魔咒串由许多魔咒字符组成,魔咒字符可以用数字表示.例如可以将魔咒字符 1 11.2 22 拼凑起来形成一个魔咒串 [1,2] [1, 2] ...
- 「SDOI2016」数字配对
「SDOI2016」数字配对 题目大意 传送门 题解 \(a_i\) 是 \(a_j\) 的倍数,且 \(\frac{a_i}{a_j}\) 是一个质数,则将 \(a_i,a_j\) 质因数分解后,其 ...
- 【LOJ】#2070. 「SDOI2016」平凡的骰子
题解 用了一堆迷之复杂的结论结果迷之好写的计算几何???? 好吧,要写立体几何了 如果有名词不懂自己搜吧 首先我们求重心,我们可以求带权重心,也就是x坐标的话是所有分割的小四面体的x坐标 * 四面体体 ...
- 【LOJ】#2069. 「SDOI2016」齿轮
题解 我一开始还努力想这道题是不是有坑,被SDOI折磨到我觉得不能有那么水的题在-- 就是带权并查集维护一下两点间距离,如果新加一条边两个点在同一集合,看看已有的路径和新加的路径是否相等 乘积可以在模 ...
随机推荐
- Springboot项目中异常拦截设计与处理
背景: 项目运行过程中会出现各种各样的问题,常见的有以下几种情况: 业务流程分析疏漏,对业务流程的反向操作.边界分析设计不充分 调用外部服务.调用外部系统出现的超时.错误.返回值与预期不符 外部资源连 ...
- Mahalanobia Distance(马氏距离)的解释
马氏距离有多重定义: 1)可以表示 某一个样本与DataSet的距离. 2)可以表示两个DataSet之间的距离. 1) The Mahalanobis distance of an observat ...
- 《PHP - 信号/基本操作/配置》
一:PHP 信号 - SIGINT / SIGTERM / SIGQUIT - 退出FPM,在master收到退出信号后将向所有的worker进程发送退出信号,然后master退出. - SIGUSR ...
- el-input 标签中密码的显示和隐藏
效果展示: 密码隐藏: 密码显示: 代码展示: 一:<el-input>标签代码 <el-form-item label="密码" prop="pass ...
- java笔记——java加载与执行
java的加载与执行过程: 编译阶段: 首先,我们先创建一个java源文件xxx.java 2. 经过javac.exe编译(如果源文件没有语法错误),生成一个或多个.class文件,注意的是此时的字 ...
- CSS中position和float的使用
近期会更新一系列博客,对基础知识再度做个巩固和梳理. 一.position定位 (一):position的属性 1.absolute:生成绝对定位的元素,相对于最近一级定位不是static的父元素来进 ...
- AttributeError: 'sys.flags' object has no attribute 'utf8_mode'
AttributeError: 'sys.flags' object has no attribute 'utf8_mode' pycharm工程的py版本是3.6,结果即使使用py3.7编译后的py ...
- UOJ46 清华集训2014玄学(线段树)
注意到操作有结合律,容易想到用一个矩形表示第i次操作对第j个位置的数的影响.那么修改是单行内的区间修改,而查询是单列内的区间查询.这样二维线段树上以列为外层行为内层直接打标记就可以维护.然后就喜闻乐见 ...
- 使用jQuery开发messager消息框插件
1.插件使用 首先引入jquery库,然后引入dialog.js.dialog.css.messager.js.messager.css,如下: <script type="text/ ...
- 纯 CSS 画 iphone
好几天没有更新了,直接上效果吧,哈哈!(我想这个应该大部分都会!哈哈哈!) 代码如下: html: <div class="container"> <div cl ...