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