hdu4355Party All the Time(三分查找模板题)

题目链接：
http://acm.hdu.edu.cn/showproblem.php?pid=4355
思路：

看完题没证出所求值在区间内是单峰函数。但是看到这种在连续区间上求最值的，应该就要想到二分或三分，而所求值必不可能在区间上是单调的，所以三分试试（三分可求单峰函数最值）。

代码：

#include <bits/stdc++.h>

using namespace std;
const double esp = 1e-4;
const int maxn = 5e4+7;
double arr[maxn];
double pos[maxn];
int n;

double getSum(double x)
{
double ans = 0.0;
for(int i = 0;i<n;i++)
{
double p = fabs(x-pos[i]);
p = p*p*p;
ans += p*arr[i];
}
return ans;
}

double trisearch(double l,double r)
{
while(r-l>esp)
{
double ml = l+(r-l)/3.0;
double mr = r-(r-l)/3.0;
if(getSum(ml)<getSum(mr))
r = mr;
else
l = ml;
}
return getSum(l);
}
int main()
{
int t,kase = 1;
scanf("%d",&t);
while(t--)
{
double l = 1000000,r = -1000000;
scanf("%d",&n);
for(int i = 0;i<n;i++)
{
scanf("%lf%lf",&pos[i],&arr[i]);
if(pos[i]<l)
l = pos[i];
if(pos[i]>r)
r = pos[i];
}
printf("Case #%d: %.0lf\n",kase++,trisearch(l,r));
}
return 0;
}