Party All the Time
题目链接:HDU - 4355
题意:在x轴上有n个点,每个点有一个权值w,定义f(x)=(fabs(x-a))^3 * w(a是x轴上一点); 求出一点a使得x轴上所有点对应的f(x)和最小;
图像是下凸的双曲线,三分解决;
注意精度,不要太高,太高会T;(上次做了个精度低了会WA的,,,,一般是1e-6, 1e-8, 1e-10中选,也不一定,,,)
#include <bits/stdc++.h>
#define eps 1e-6//注意了精度高了会T;
using namespace std;
int N;
double x[50010], w[50010];
double cal(double a){
double ans=0.0;
int i=0;
for(i=0; i<N; i++){
ans+=pow(fabs(x[i]-a), 3)*w[i];
}
return ans;
}
int main(){
int T, cas=0;
scanf("%d", &T);
while(T--){
cas++;
scanf("%d", &N);
for(int i=0; i<N; i++){
scanf("%lf%lf", &x[i], &w[i]);
}
double l=x[0], r=x[N-1];
while(r-l>eps){
double mid=l+(r-l)/3;
double midd=r-(r-l)/3;
if(cal(mid)>cal(midd)) l=mid;
else r=midd;
}
double temp=cal(l);
long long ans=(int)temp;
if(temp-ans>=0.5) ans++;
printf("Case #%d: %lld\n", cas, ans);
}
return 0;
}