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