HDOJ - 3714 三分搜索

     昨天编程之美的C题搞不出来..看了下别人的方法...三分搜索水之...以前没写过...

    三分从二分上发展来的...当这个函数不是单峰时,也能找到最值~~好牛逼....

     每次不仅取中点mid=(l+r)/2..还取令一个次中点mmid=(mid+r)/2...如果mid上的值比mmid优...r=mmid...否则 l=mid...

    这题和昨天的C题可以说一模一样...出原题不太好吧....


Program:

#include<iostream>
#include<cmath>
#include<stack>
#include<queue>
#include<set>
#include<algorithm>
#include<stdio.h>
#include<string.h>
#define ll long long
#define oo 1000000007
using namespace std; 
struct node
{
       double a,b,c;
}f[10005];
int n;
double getdata(double x)
{
       double m=f[1].a*x*x+f[1].b*x+f[1].c;
       for (int i=2;i<=n;i++)  m=max(m,f[i].a*x*x+f[i].b*x+f[i].c); 
       return m;
}
int main()
{         
       int T,i; 
       double l,r,ml,mr,dl,dr;
       scanf("%d",&T);
       while (T--)
       {
              scanf("%d",&n);
              for (i=1;i<=n;i++) scanf("%lf%lf%lf",&f[i].a,&f[i].b,&f[i].c);
              l=0.0;
              r=1000.0; 
              while (r-l>1e-10)
              {
                    ml=(l+r)/2;
                    mr=(ml+r)/2;
                    dl=getdata(ml);
                    dr=getdata(mr);
                    if (dl>dr) l=ml;
                         else r=mr;
              }
              printf("%.4lf\n",getdata(r));
       }
       return 0;
}


你可能感兴趣的:(HDOJ - 3714 三分搜索)