bzoj 3203: [Sdoi2013]保护出题人 凸包
题目大意:
题解
首先我们考虑对一大波僵尸来袭的情况进行分析
假设来袭的僵尸是\(\{ a_1,a_2,a_3,...,a_n\}\)
那么距离分别为\(\{dis,dis+d,dis+d*2,...,dis+d*(n-1)\}\)
那么我们知道,在这波僵尸中,我们应该取到的最小的攻击速度为
\(max{\frac{sum_i}{dis_i}}\)
其中\(sum_n = \sum_{i=1}^{n}a_i,dis_i = dis + (i-1)*d\)
然后我们把这个式子转化到原序列中,发现式子变成了下面这个样子
假设排头的僵尸是\(p\) :
那么有\(ans_p = max_{i=1}^n\{\frac{sum_p - sum_{i-1}}{(dis_p+p*d - i*d)}\}\)
而我们的目的就是最小化\(\sum_{i=1}^nans_i\)
我们观察上面关于\(ans_p\)的式子,发现这实际上是两个点的斜率
即\((dis_p+p*d,sum_p)\)和\((i*d,sum_{i-1})\)两个点的斜率
所以我们枚举每个p,然后对于每个p都计算出所有的\((i*d,sum_{i-1})\)
然后再取斜率最大值即可...
TLE
我们发现每次针对不同的p计算的时候\((i*d,sum_{i-1})\)居然都是相同的!
所以我们根本不需要每次都重新计算
那么每次当我们枚举到一个p的时候,计算所有满足\(i\leq p\)的点中与其最大斜率即可
那怎么计算这个最大斜率呢??
我们知道每次加入的点在横纵坐标上一定都比上一次加入的点要大
说起来比较矛盾,但是我们需要维护一个不会将任何点包括的"上凸壳"
然后在这个"上凸壳"上的斜率明显是一个单峰函数,所以我们三分即可.
#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
inline void read(ll &x){
x=0;char ch;bool flag = false;
while(ch=getchar(),ch<'!');if(ch == '-') ch=getchar(),flag = true;
while(x=10*x+ch-'0',ch=getchar(),ch>'!');if(flag) x=-x;
}
inline int cat_max(const int &a,const int &b){return a>b ? a:b;}
inline int cat_min(const int &a,const int &b){return a<b ? a:b;}
const int maxn = 100010;
const double eps = 1e-9;
inline int dcmp(const double &x){
if(x < eps && x > -eps) return 0;
return x > 0 ? 1 : -1;
}
struct Point{
double x,y;
Point(const double &a = 0,const double &b = 0){x=a;y=b;}
void print(){
printf("Point : (%lf,%lf)\n",x,y);
}
};
typedef Point Vector;
inline Vector operator + (const Vector &a,const Vector &b){
return Vector(a.x+b.x,a.y+b.y);
}
inline Vector operator - (const Vector &a,const Vector &b){
return Vector(a.x-b.x,a.y-b.y);
}
inline double cross(const Vector &a,const Vector &b){
return a.x*b.y - a.y*b.x;
}
inline double slope(const Point &a,const Point &b){
return (a.y - b.y)/(a.x - b.x);
}
Point p[maxn];int m = 0;
inline void insert(const Point &x){
while(m > 1 && dcmp(cross(p[m] - p[m-1],x - p[m])) <= 0) -- m;
p[++m] = x;
}
ll a[maxn],dis[maxn],sum[maxn];
double solve(const Point &x){
int l = 1,r = m,midx,midy;
while(l <= r-3){
midx = (l+l+r)/3;midy = (l+r+r)/3;
if(slope(p[midx],x) > slope(p[midy],x)) r = midy;
else l = midx;
}
double ans = .0;
for(int i=l;i<=r;++i) ans = max(ans,slope(p[i],x));
return ans;
}
int main(){
ll n,d;read(n);read(d);
for(int i=1;i<=n;++i){
read(a[i]);read(dis[i]);
sum[i] = sum[i-1] + a[i];
}double ans = .0;
for(int i=1;i<=n;++i){
insert(Point(i*d,sum[i-1]));
ans += solve(Point(dis[i]+i*d,sum[i]));
}printf("%.0f\n",ans);
getchar();getchar();
return 0;
}
bzoj 3203: [Sdoi2013]保护出题人 凸包的更多相关文章
- BZOJ 3203 [SDOI2013]保护出题人 (凸包+三分)
洛谷传送门 题目大意:太长略 每新加入一个僵尸,容易得到方程$ans[i]=max{\frac{sum_{i}-sum_{j-1}}{s_{i}+d(i-j)}}$ 即从头开始每一段僵尸都需要在规定距 ...
- BZOJ 3203 Luogu P3299 [SDOI2013]保护出题人 (凸包、斜率优化、二分)
惊了,我怎么这么菜啊.. 题目链接: (bzoj)https://www.lydsy.com/JudgeOnline/problem.php?id=3203 (luogu)https://www.lu ...
- 【BZOJ】3203: [Sdoi2013]保护出题人(几何+三分+特殊的技巧)
http://www.lydsy.com/JudgeOnline/problem.php?id=3203 wa无数次QAQ,犯sb错....一是数组没有引用...二是输出转成了int(越界了sad). ...
- 【bzoj3203】[Sdoi2013]保护出题人 凸包+二分
题目描述 输入 第一行两个空格隔开的正整数n和d,分别表示关数和相邻僵尸间的距离.接下来n行每行两个空格隔开的正整数,第i + 1行为Ai和 Xi,分别表示相比上一关在僵尸队列排头增加血量为Ai 点的 ...
- [BZOJ3203][SDOI2013]保护出题人(凸包+三分)
https://www.cnblogs.com/Skyminer/p/6435544.html 先不要急于转化成几何模型,先把式子化到底再对应到几何图形中去. #include<cstdio&g ...
- 【BZOJ3203】[Sdoi2013]保护出题人 二分+凸包
[BZOJ3203][Sdoi2013]保护出题人 Description Input 第一行两个空格隔开的正整数n和d,分别表示关数和相邻僵尸间的距离.接下来n行每行两个空格隔开的正整数,第i + ...
- [BZOJ3203] [SDOI2013]保护出题人(二分+凸包)
[BZOJ3203] [SDOI2013]保护出题人(二分+凸包) 题面 题面较长,略 分析 对于第i关,我们算出能够打死前k个个僵尸的最小能力值,再取最大值就可以得到\(y_i\). 前j-1个僵尸 ...
- 洛谷 P3299 [SDOI2013]保护出题人 解题报告
P3299 [SDOI2013]保护出题人 题目描述 出题人铭铭认为给SDOI2012出题太可怕了,因为总要被骂,于是他又给SDOI2013出题了. 参加SDOI2012的小朋友们释放出大量的僵尸,企 ...
- [SDOI2013]保护出题人
题目 出题人铭铭认为给SDOI2012出题太可怕了,因为总要被骂,于是他又给SDOI2013出题了. 参加SDOI2012的小朋友们释放出大量的僵尸,企图攻击铭铭的家.而你作为SDOI2013的参赛者 ...
随机推荐
- Tornado 安装及简单程序示例
1.安装步骤:tar xvzf tornado-3.2.tar.gz cd tornado-3.2 python setup.py build sudo python setup.py install ...
- 利用certutil.exe 传文件
certutill.exe 在Windows 7 及其之后的所有Windows Server和Workstation版本均预装 1. Encode file: certutil -encode kk. ...
- 统计Visual Studio项目的代码行数
原文转自 https://blog.csdn.net/tyc129/article/details/74279806 使用Visual Studio 自带的在文件中查找功能中的正则表达式实现代码统计功 ...
- CentOS7手动编译安装内核4.11.7
1. 进入/usr/src/目录 cd /usr/src 2. 下载内核源码,网址:https://www.kernel.org wget https://cdn.kernel.org/pub/lin ...
- python实战===爬取所有微信好友的信息
''' 爬取所有T信好友的信息 ''' import itchat from pandas import DataFrame itchat.login() friends=itchat.get_fri ...
- mybatis-plus的学习
1.mybatisplus 提供了比较齐全的crud即增删改查,不需要在mapper.xml里写sql可以直接调用 原文链接:http://blog.csdn.net/u014519194/artic ...
- idea+spring+springmvc+mybatis+mybatis+maven
使用SSM(Spring,SpringMVC和Mybatis) 1.1.Spring Spring是一个开源框架,Spring是于2003 年兴起的一个轻量级的Java 开发框架,由Rod Johns ...
- php的设计模式------工厂模式
1.工厂模式简介 属于创建型模式.定义一个创建对象的接口,让其子类自己决定实例化哪一个工厂类,工厂模式使其创建过程延迟到子类进行主要解决的问题:接口选择的问题. 2.分类 2.1 简单工厂模式 接口: ...
- Myeclipse实用快捷键总结
alt+shift+J 为选中的类/方法添加注释 ctrl+T 显示选中类的继承树 ctrl+shift+X/Y 将选中的字符转换为大写/小写 ctrl+shift+R 打开资源 ctrl+shift ...
- linux命令(46):chgrp命令
在lunix系统里,文件或目录的权限的掌控以拥有者及所诉群组来管理.可以使用chgrp指令取变更文件与目录所属群组,这种方式采用群组名称或群组识别码都可以.Chgrp命令就是change group的 ...