hdu4355

链接：点击打开链接

题意：

小精灵不愿意走路，走路的话会使他们产生不高兴值，如果小精灵走的路是S体重是W那么他走S路程产生的不高兴值是ans=sigma（S^3）*W；现在森林里要开一个庆祝会，给出各个小精灵的坐标（注意这里是一维的）和体重，求在哪里庆祝会使小精灵的不高兴值最小，输出这个最小的不高兴值。样例中第一行是样例个数，每个样例的第一行是小精灵的个数n，接下的n行是每个小精灵的坐标和体重。输出只有一个整数，即最小的不高兴值。注意输出格式。

代码：

#include <iostream>
#include <cmath>
#include <cstdio>
#include <algorithm>
using namespace std;
double x[50005],w[50005];
long long n;
double cal(double a){
    long long i;
    double sum=0;
    for(i=0;i<n;i++)
    sum+=(fabs((x[i]-a)*(x[i]-a)*(x[i]-a))*w[i]);       //每个位置的值加和
    return sum;
}
double threesearch(double l,double r){
    double mid,midmid;
    long long i;
    for(i=1;i<=100;i++){
        mid=(r-l)/3+l;
        midmid=(r-l)/3*2+l;                             //三分
        if(cal(mid)>cal(midmid))
        l=mid;
        else
        r=midmid;
    }
    return cal(l);
}
int main(){
    long long t,i,j,k;
    double ll,rr,sum;
    scanf("%I64d",&t);
    for(k=1;k<=t;k++){
        scanf("%I64d",&n);
        ll=99999999;rr=-99999999;
        for(i=0;i<n;i++){
        scanf("%lf %lf",&x[i],&w[i]);
        ll=min(ll,x[i]);
        rr=max(rr,x[i]);
        }
//        cout<<ll<<" "<<rr<<endl;
        sum=threesearch(ll,rr);
        printf("Case #%I64d: %I64d\n",k,(long long)(sum+0.5));//这要注意精度问题
    }
    return 0;
}