HDU-#4355 Party All the Time（三分法+凸函数求极值）

       题目大意：给出函数 S³*W，在给出该函数的坐标以及对应的值，求在左右区间中的极小值。

      解题思路：分析函数可知，其满足凸函数的性质，因此可以利用三分法进行求解。做法是在区间[l,r]的两个三分点m1和m2，比较f(m1)和f(m2)的大小。如果前者小于后者，则解在[l,m2]中，反之在[m1,r]中。就这样进行三分，直到求出极值点，不理解可以画下图模拟下。详见code。

      题目来源：http://acm.hdu.edu.cn/showproblem.php?pid=4355

      code：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
using namespace std;

const int MAXN = 50000+10;
const double eps = 1e-12;
int T,n,k;
double s[MAXN],w[MAXN];

double f(double x){ //所要计算的函数
    double ans=0;
    for(int i=0;i<n;i++)
        ans+= fabs((s[i]-x)*(s[i]-x)*(s[i]-x))*w[i];
    return ans;
}

int main(){
    k=0;
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        for(int i=0;i<n;i++)
            scanf("%lf%lf",&s[i],&w[i]);
        double l=s[0],r=s[n-1]; //左右边界的设置，题目给出是递增的坐标
        while(l+eps<r){ //三分，eps为精度
            double m1=(l+r)/2;
            double m2=(m1+r)/2;
            if(f(m1)<f(m2)) r=m2;
            else l=m1;
        }
        printf("Case #%d: %.0f\n",++k,f(l));
    }
    return 0;
}