hdu_4355 Party All the Time(简单数学题)

http://acm.hdu.edu.cn/showproblem.php?pid=4355

题意：
直线上n个点Xi,每个点有一个权值Wi,求存在某点S使得(abs(s-xi))^3*wi的和最小，求这个最小值。
思路：
根据题意列出函数后发现没有单调性，然后想是否可以三分，但无法证明。。。想了好久才发现，枚举n个点分成的n-1个区间，函数可以去掉绝对值。（极值点通过求导求得）分别求最小值。所有区间最小值的最小值即为答案。（PS:解题报告上说这是签到题。。。）

我的代码：

#include<stdio.h>
#include<math.h>
#define M 50001
#define eps 1e-6
double a,b,c,d,s,sumA[M],sumB[M],sumC[M],sumD[M];
double x[M],w[M],ans,delta,t,x1,x2;
int T,i,j,n,k,js=1;
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);ans=1e18;
        sumA[0]=sumB[0]=sumC[0]=sumD[0]=0;
        for(i=0;i<n;i++)
        {
            scanf("%lf%lf",&x[i],&w[i]);
                sumA[i+1]=w[i]+sumA[i];
                sumB[i+1]=w[i]*x[i]+sumB[i];
                sumC[i+1]=w[i]*x[i]*x[i]+sumC[i];
                sumD[i+1]=w[i]*x[i]*x[i]*x[i]+sumD[i];
        }
        for(k=1;k<n;k++)
        {
            a=sumA[k]-(sumA[n]-sumA[k]);
            b=3.0*(sumB[n]-sumB[k]-sumB[k]);
            c=3.0*(sumC[k]-sumC[n]+sumC[k]);
            d=sumD[n]-sumD[k]-sumD[k];
            delta=4.0*b*b-12.0*a*c;//a*s^3+b*s^2+c*s+d求导
            if(delta>0) 
            {
                x1=(-2.0*b+sqrt(delta))/(6.0*a);
                x2=(-2.0*b-sqrt(delta))/(6.0*a);
                if(x1>x[k-1]+eps&&x1<x[k]-eps) t=a*x1*x1*x1+b*x1*x1+c*x1+d,ans=t<ans?t:ans;
                if(x2>x[k-1]+eps&&x2<x[k]-eps) t=a*x2*x2*x2+b*x2*x2+c*x2+d,ans=t<ans?t:ans;
            }
            t=a*x[k-1]*x[k-1]*x[k-1]+b*x[k-1]*x[k-1]+c*x[k-1]+d,ans=t<ans?t:ans;
            t=a*x[k]*x[k]*x[k]+b*x[k]*x[k]+c*x[k]+d,ans=t<ans?t:ans;
        }
        printf("Case #%d: ",js++);
        if(n==1) printf("0\n");
        else printf("%.0lf\n",ans);//这里WA了很多次，最后猛然意识到，%.0lf本来就是四舍五入的。。。 
    }
    return 0;
}