SCAU培训 C 平面最佳坐标点

C  平面最佳坐标点
Time Limit:500MS  Memory Limit:65535K
题型: 编程题   语言: 无限制
描述 在平面坐标上，给n个坐标点，现定义两个点（x1,y1）和（x2,y2）的“距离”为|x1-x2|3+|x1-x2|*|y1-y2|+|y1-y2|3 现要求编程找一个点，使得n个点到该点的距离累加和最小，精度去到小数点3位
输入格式
第一行一个数n(n<=100)，
此后n行，每行两个浮点数为坐标点x y(0<=x,y<=100)
输出格式
输出最佳位置坐标

输入样例
3
0 0
100 0
100 100
输出样例
58.648 41.352
/
测试2：
测试数据用例 No.1
标准输入数据：
3
1 1
2 2
9 9
标准输出答案：
1|4.557 4.557
Provider admin

#include <iostream>
#include <cmath>
#include<algorithm>
using namespace std;

const int N = 1001;//第二次居然改成1000看傻我

int n;
long double x[N],y[N];

long double dis(long double x1,long double y1)
{
    long double sum = 0.0;
    for(int i =0;i < n;i ++ ){
            long double a=x1-x[i]>0?(x1-x[i]):(x[i]-x1);
            long double b=y1-y[i]>0?(y1-y[i]):(y[i]-y1);
            sum += a*a*a+b*b*b+a*b;
    }
    return sum;
}

long double yyy;
long double check(long double x1) // 三分y，x和y都是凹函数，结构一样
{
    //long double l= y[0],r = y[n-1];
     long double l=0.0,r=1001.0;
    while(r - l >1e-8)
    {
        long double mid = (l + r) / 2.0;
        long double mmid = (mid + r) / 2.0;
        if(dis(x1, mid) > dis(x1, mmid)) l = mid;
        else r = mmid;
        //cout<<"||||yr="<<r<<endl;// 这里是r = mmid;
    }
    yyy=r;
    return dis(x1,l);
}
int main()
{
    std::ios::sync_with_stdio(false);
    scanf("%d",&n);
    for(int i = 0;i < n;i ++ ) {
            scanf("%LF%LF",&x[i],&y[i]);
    }
//    sort(x,x+n);
//    sort(y,y+n);
//    long double l=x[0] ,r=x[n-1] ;
    long double l=0.0,r=1001.0;
    while(r - l >1e-8)  // 三分x
    {
        long double mid = (l + r) / 2.0;
        long double mmid = (mid + r) / 2.0;
        if(check(mid) > check(mmid)) l = mid;
        else r = mmid; // 这里是r = mmid;
       // cout<<"||||xr="<<r<<endl;
    }
    printf("%.3LF %.3LF",l,yyy);
    return 0;
}