HDU 4355-三分

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

题目 给n个人的位置 和 权值w

选 数轴上一实数点， 求每个人到该点的距离dis的三次方*w，的和

让这个和最小。

分析得，为单峰函数，三分解决

#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <queue>
#include <map>
#include <set>
#include <vector>
#include <iostream>
//#include <time.h>

using namespace std;

const double pi=acos(-1.0);
double eps=0.000001;

struct node
{
double x,w;
node() {}
node( double a,double b)
{
x=a,w=b;
}
};
node aa[50005];

int n;
double F(double pos)
{
double ans=0;
for (int i=1; i<=n; i++)
{
double sb=fabs(aa[i].x-pos);
ans+=sb*sb*sb*aa[i].w;
}
return ans;
}
//double aaa[10005];
int main()
{
int t;
cin>>t;
// srand( (unsigned)time( NULL ) );
int cnt=1;
while(t--)
{
cin>>n;
double x,w;
int i;
/* for (i=1;i<=n;i++)
aaa[i]=rand()%1000*1.5;*/
//sort(aaa+1,aaa+1+n);

for (i=1; i<=n; i++)
{
scanf("%lf %lf",&aa[i].x,&aa[i].w);
// aa[i].x=aaa[i]*pow(-1.0,1.0*i);
// aa[i].w=rand()%15*1.0;
}

double l=-1e6;
double r=1e6;
for (i=0;i<100;i++)
{
double m1=l+(r-l)/3;
double m2=r-(r-l)/3;
if (F(m2)>F(m1))
r=m2;
else
l=m1;
}

printf("Case #%d: %.0lf\n",cnt++,F(l));
/*
double minn=F(1);
for (i=1;i<=n;i++)
{
minn= min(minn,F(i));
}
long long oo=(long long) (minn);
printf("ans2::%lld\n",oo);*/

}

return 0;

}