hdu 4036 暴力二分 高精度 Rolling Hongshu

题意：http://acm.hdu.edu.cn/showproblem.php?pid=4036

神一般的二分。一个sweet potato想要去约会，从起点滚到终点，途经一个mountain，有m个峰顶，第一个峰顶，第m个峰顶分别为起点和终点，忽略摩擦力，在每一个斜坡上分布着bitter potatoes，每经过一个bitter potato，bitter potato就会以一定的初速度追赶sweet potato，问，sweet potato应该以最小多大的初速度滚动，保证不被bitter potato追到，完成约会；

解：

不仅要考虑到苦patato的时候速度要大，而且要考虑过山峰的时候速度要大于等于0.因为这个wa了甚久。然后全删了，本来还以为是我以前写过，于是就把以前代码贴过来改。改。但是事实是那是昨天朱老板用我的账号交的。。。不管怎样，题意以及解法是懂了的。

还有，尽量少开方。

网上还有代码是直接求的，不二分。求出到每个山峰需要多少速度，求出不被苦potato追到需要多少速度，取最大的。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <iostream>
#define eps 1.0e-8
using namespace std;

struct node{
double x,y,w,v;
bool flag;
}point[5555];

struct line{
double k,b;
};

line getl(node a,node b){
line ans;
ans.k = (b.y - a.y)/(b.x - a.x);
ans.b = a.y - a.x*ans.k;
return ans;
}

bool cmp(node a,node b){
if(a.x != b.x)
return a.x < b.x;
return a.y < b.y;
}
int main(){
int T,n,m,w;
scanf("%d",&T);
int id = 1;
double g = 20.0;
while(T--){
scanf("%d%d%d",&n,&m,&w);
for(int i = 0 ; i < n ; i ++){
scanf("%lf%lf",&point[i].x,&point[i].y);
point[i].flag = 0;
}
for(int i = n ; i < n + m ; i ++){
point[i].y = -1000000000;
scanf("%lf%lf%lf",&point[i].x,&point[i].v,&point[i].w);
point[i].flag = 1;
}
sort(point,point+n+m,cmp);
for(int i = 0 ; i < n + m ; i ++){
if(point[i].flag == 1){
int j = i;
while(j--){
if(point[j].flag == 0)
break;
}
int k = i;
while(k++){
if(point[k].flag == 0)
break;
}
line tp = getl(point[j],point[k]);
point[i].y = point[i].x*tp.k + tp.b;
}
}
double u,d,mid;
u = 999999999999LL;
d = 0.00;
while(d <= u){
mid = (u+d)/2;
bool flag = 1;
node s = point[0];
s.v = mid*mid;
for(int i = 1 ; i < n + m ; i ++){
s.v = 2*(g*(s.y - point[i].y) + s.v/2);
s.y = point[i].y;
s.x = point[i].x;
if(s.v < 0){
flag = 0;
break;
}
if(point[i].flag){
if(s.v <= point[i].v*point[i].v){
flag = 0;
break;
}
}
}
if(flag){
u = mid - eps;
}
else
d = mid + eps;
}
printf("Case %d: %.2f\n",id++,fabs((u + d) / 2.0));
}
return 0;
}

转/article/6683775.html

#include<stdio.h>
#include<math.h>
#include<string.h>
struct  PEAK
{
double x,h;
}p[1100];
struct bit
{
double x,v;
}pat[1100];
int main()
{
int G=20;
double h,weight;
int t,cases=1,i,j,n,m,w;
scanf("%d",&t);
while(t--)
{
double v0=0,temp;
scanf("%d%d%d",&n,&m,&w);
for(i=1;i<=n;i++) //求出到达任意一个峰顶需要的最小速度
{
scanf("%lf%lf",&p[i].x,&p[i].h);
if(p[i].h-p[1].h<=0)  //峰顶在起始点一下，初速度可以为零
temp=0;
else temp=sqrt(2*G*(p[i].h-p[1].h));
if(temp>v0) v0=temp;
}
double vx,v=0;
for(i=1;i<=m;i++)   //根据能量守恒，求出到达每一个bitter potato需要的最小速度
{
scanf("%lf%lf%lf",&pat[i].x,&pat[i].v,&weight);
pat[i].x+=p[1].x;     //将bitter potato的坐标改为相对于原点
for(j=1;j<n;j++)
{
if(pat[i].x>=p[j].x&&pat[i].x<=p[j+1].x)  //找出当前bitter potato所在的斜坡
{
h=(1.0*(p[j+1].h-p[j].h)/(p[j+1].x-p[j].x))*(pat[i].x-p[j].x) + p[j].h;  //bitter potatoe相对于原点的高度
break;
}
}
vx=sqrt(1.0*pat[i].v*pat[i].v+1.0*2*G*(h-p[1].h)); //求出到达当前bitter potato需要的速度
if(vx>v)
v=vx;
}
if(v>v0)
printf("Case %d: %.2lf\n",cases++,v);
else
printf("Case %d: %.2lf\n",cases++,v0);
}
return 0;
}