做题感悟:理解错题意了,我以为图象只有两段,所以我分了两段分别三分,两组样例只对了一个,看了别人的才知道需要对每一段三分。
解题思路:根据题意:可以将图象分成 m - 1 段 a0 ~ a1 , a1 ~ a2 , ai ~ ai+1 , an-2 ~ an-1 对每段三分即可,为什么要这样呢?
假设 : 最小值为 F( x ) , so F( x ) = n / (x * x ) + k * x + b ;这个函数在某段区间可能是单调递增的,也可能先递减后递增,这样都没关系三分都可以解决,如果单纯的用二分先递减再递增的情况就解决不了。
代码:
#include<stdio.h> #include<iostream> #include<map> #include<stack> #include<string> #include<string.h> #include<stdlib.h> #include<math.h> #include<vector> #include<queue> #include<algorithm> using namespace std ; #define LEN sizeof(struct node) #define pret(a,b) memset(a,b,sizeof(a)) #define lld __int64 const double PI = 3.1415926535898 ; const double INF = 999999999 ; const double esp = 1e-8 ;// 精确到 7 位即可 const lld md= 2810778 ; const int MX = 10005 ; struct node { double x,y ; }T[MX] ; int n ; double m,k,b,mx ; double find(double x) { return m/(x*x)+k*x+b ; } double third_search(double le,double rt) { double m1,m2,mid,midd ; while(le+esp<=rt) { mid=(rt+le)/2.0 ; midd=(mid+rt)/2.0 ; m1=find(mid) ; m2=find(midd) ; m1 > m2 ? le=mid : rt=midd ; } return find(mid) ;// 切记不要把 mid 换成 le } int main() { while(~scanf("%d%lf",&n,&m)) { for(int i=0 ;i<n ;i++) scanf("%lf%lf",&T[i].x,&T[i].y) ; double best=INF ; for(int i=1 ;i<n ;i++) // 对每一段三分 { k=(T[i].y-T[i-1].y)/(T[i].x-T[i-1].x) ; b=T[i].y-k*T[i].x ; mx=third_search(T[i-1].x,T[i].x) ; if(mx<best) best=mx ; } printf("%.3lf\n",best) ; } return 0 ; }