Party All the Time HDU - 4355(三分)

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