P2179-[NOI2012]骑行川藏【导数,二分】
正题
题目链接:https://www.luogu.com.cn/problem/P2179
题目大意
给出EEE和nnn个si,ki,uis_i,k_i,u_isi,ki,ui求一个序列viv_ivi满足
∑i=1nkisi(vi−ui)2≤E\sum_{i=1}^nk_is_i(v_i-u_i)^2\leq Ei=1∑nkisi(vi−ui)2≤E
的情况下最小化
∑i=1nsivi\sum_{i=1}^n\frac{s_i}{v_i}i=1∑nvisi
1≤n≤1041\leq n\leq 10^41≤n≤104
解题思路
洛谷题解上一个十分神奇的做法看起来。(主要是看不懂拉格朗日乘数法/kk)
首先考虑对于段路的行驶时间ti=sivit_i=\frac{s_i}{v_i}ti=visi,我们可以画出消耗的能量EEE和tit_iti的函数。
对于函数f(E)=tif(E)=t_if(E)=ti不难发现的是在vi≥uiv_i\geq u_ivi≥ui的情况下EEE越小这个函数对应位置的导数越小。
也就是消耗单位能量减少的时间也就越少,性价比就越低。而我们现在要给每段路分配一个tit_iti使得消耗能量和等于EEE且tit_iti和最小的话。
根据贪心的思想有选出若干个的tit_iti满足对应位置的导数相等。
那么我们就找到了所有路的共性,考虑二分这个导数,但是我们先对这个函数f(v)=tEf(v)=\frac{t}{E}f(v)=Et求个导。
t′=−svi2,E′=2kisi(vi−ui)t'=-\frac{s}{v_i^2},E'=2k_is_i(v_i-u_i)t′=−vi2s,E′=2kisi(vi−ui)
f′(v)=t′E′=−s2kisivi2(vi−ui)f'(v)=\frac{t'}{E'}=-\frac{s}{2k_is_iv_i^2(v_i-u_i)}f′(v)=E′t′=−2kisivi2(vi−ui)s
然后我们二分出f′(vi)=xf'(v_i)=xf′(vi)=x然后再二分出对应的速度viv_ivi就好了。
code
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1e4+10;
int n;double E,s[N],k[N],v[N];
double getv(double x,int p){double l=max(v[p],0.0),r=100000;for(int i=1;i<=100;i++){double V=(l+r)/2.0;if(-2.0*k[p]*V*V*x*(V-v[p])<1.0)l=V;else r=V;}return (l+r)/2.0;
}
double check(double x){double E=0;for(int i=1;i<=n;i++){double V=getv(x,i);E+=k[i]*s[i]*(V-v[i])*(V-v[i]);}return E;
}
int main()
{scanf("%d",&n);scanf("%lf",&E);for(int i=1;i<=n;i++)scanf("%lf%lf%lf",&s[i],&k[i],&v[i]);double l=-1e5,r=0;for(int i=1;i<=100;i++){double mid=(l+r)/2.0;if(check(mid)<=E)l=mid;else r=mid;}double mid=(l+r)/2.0,ans=0;for(int i=1;i<=n;i++)ans+=s[i]/getv(mid,i);printf("%.12lf\n",ans);return 0;
}
P2179-[NOI2012]骑行川藏【导数,二分】相关推荐
- 【洛谷】P2179 [NOI2012]骑行川藏
题解 感谢小迪给我讲题啊,这题小迪写挺好的我就不写了吧 小迪的题解 代码 #include <iostream> #include <cstdio> #include < ...
- [NOI2012]骑行川藏
[NOI2012]骑行川藏 思路一:二分导数 http://www.cnblogs.com/RabbitHu/p/9019762.html 考虑"性价比"即花费单位能量缩短的时间. ...
- 2876: [Noi2012]骑行川藏
2876: [Noi2012]骑行川藏 Time Limit: 20 Sec Memory Limit: 128 MBSec Special Judge Submit: 1243 Solved: ...
- bzoj 2876: [Noi2012]骑行川藏 拉格朗日数乘
2876: [Noi2012]骑行川藏 Time Limit: 20 Sec Memory Limit: 128 MBSec Special Judge Submit: 1033 Solved: ...
- bzoj 2876: [Noi2012]骑行川藏 二分+拉格朗日乘数法
题意 蛋蛋非常热衷于挑战自我,今年暑假他准备沿川藏线骑着自行车从成都前往拉萨.川藏线的沿途有着非常美丽的风景,但在这一路上也有着很多的艰难险阻,路况变化多端,而蛋蛋的体力十分有限,因此在每天的骑行前设 ...
- bzoj2876 [Noi2012]骑行川藏 [二分+拉格朗日乘数法]
Description 蛋蛋非常热衷于挑战自我,今年暑假他准备沿川藏线骑着自行车从成都前往拉萨.川藏线的沿途有着非常美丽的风景,但在这一路上也有着很多的艰难险阻,路况变化多端,而蛋蛋的体力十分有限,因 ...
- bzoj2876: [Noi2012]骑行川藏 :拉格朗日乘数法
题目链接 bzoj2876 题目描述 Description 蛋蛋非常热衷于挑战自我,今年暑假他准备沿川藏线骑着自行车从成都前往拉萨.川藏线的沿途有着非常美丽的风景,但在这一路上也有着很多的艰难险阻, ...
- 2876: [Noi2012]骑行川藏 - BZOJ
Description 蛋蛋非常热衷于挑战自我,今年暑假他准备沿川藏线骑着自行车从成都前往拉萨.川藏线的沿途有着非常美丽的风景,但在这一路上也有着很多的艰难险阻,路况变化多端,而蛋蛋的体力十分有限,因 ...
- BZOJ2876 [Noi2012]骑行川藏
本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000 作者博客:http://www.cnblogs.com/ljh2000-jump/ ...
- [BZOJ2876] [NOI2012]骑行川藏
Description 蛋蛋非常热衷于挑战自我,今年暑假他准备沿川藏线骑着自行车从成都前往拉萨.川藏线的沿途有着非常美丽的风景,但在这一路上也有着很多的艰难险阻,路况变化多端,而蛋蛋的体力十分有限,因 ...
最新文章
- 记录一次简单、高效、无错误的linux上安装pytorch的过程
- SNF快速开发平台MVC-高级查询组件
- Postfix无法正常发送邮件故障的排除
- Android零基础入门第44节:ListView数据动态更新
- const 的学习(转载)
- Markdown中数学公式练习
- better-scroll 上拉加载,下拉刷新(解决移动端长页面卡顿)
- 抓wifi包工具linux,测试工具之在Ubuntu使用Wireshark抓无线数据包
- 怎么理解 IaaS、SaaS 和 PaaS 的区别?
- 在Simplicity Studio下创建适用于EFR32的工程项目
- 程序员要记住的几条建议
- python什么软件开发好_python怎样才能学好?python软件开发什么
- 为计算机新建一个用户,win10
- 学习记录——背包问题基础公式解释回顾
- MySQL存储过程中游标使用
- Python+OpenCV:Hough圆检测(Hough Circle Transform)
- poj 无向图最小环问题 详解,模板
- python函数传递列表_python传递列表作为函数参数
- saltstack的探索-利用脚本增加用户
- 程序设计语言的基本概念
热门文章
- qt4.7 mysql_详解Qt 4.7编译和访问Mysql驱动
- mysql数据库的安装和配置文件_MySQL 数据库安装与配置详解
- 京瓷1020怎么打印自检页_[建筑]喷墨打印机如何打印自检页 详细
- c++11 string u8_深入理解C++11:C++11新
- 服务器不知道怎么回事安卓系统很卡,为什么安卓系统很容易变卡?该怎么解决?看完长知识了...
- 360浏览器5兼容模式吗_个独模式真能将企业总体税负降低至5%吗,为何大家都热衷核定征收...
- php 常用rpc框架,php的轻量级rpc框架yar
- 机器学习——决策树学习
- 算法题目——读书知识点统计问题(POJ-3320)(尺取法)
- java使用教程——组件及事件处理——窗口(设置窗口的颜色和背景)