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-8using 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;}

转http://www.cnblogs.com/nanke/archive/2011/09/12/2174032.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;}