16.7.7 T1:YL杯超级篮球赛

一年一度的高一YL杯超级篮球赛开赛了。当然，所谓超级，意思是参赛人数可能多余5人。小三对这项篮球非常感兴趣，所以一场都没有落下。每个中午都准时守侯在篮球场看比赛。经过一个星期的研究，小三终于对篮球的技战术找到了一丝丝感觉了。他发现打YL杯的每个班都有一套相似的进攻战术：
　　1 ：控球后卫带球到前场，找到一个最佳攻击点 ( x , y )
　　2 ：所有除控卫以外的队员都从各自的当前位置迅速向 ( x , y ) 移动
　　3 ：控球后卫根据场上情况组织进攻
　　这个战术对于一般情况是非常奏效的，但是每个队员毕竟不像小三一样每天精力过剩，每个队员都有一个疲劳指数W，显然对于每个队员的移动需要消耗一些能量。
　　假设一个队员从位置 ( x1 , y1 ) 移动到 ( x , y )的能量消耗为 w * (ABS ( x - x1 ) +ABS ( y - y1 ) ), 这里ABS为绝对值函数。那么我们希望整个队伍一次进攻的能量消耗当然是越少越好。显然能量消耗的多少直接取决于控球后卫对于攻击点 ( x , y )的选择。
　　因为参赛人数众多，所以小三希望你能编写一个程序，即帮他找出某个时刻的最佳攻击点。

Input

　　第一行：一个整数N，表示篮球队人数
　　第二行：一共N个整数，其中的第i个数Wi表示第i个队员的疲劳指数。
　　第3~N+2行：每一行两个整数X和Y，其中的第i+2行，表示第i个队员的当前位置的横坐标和纵坐标。

Output

　　一个实数。表示所有队员集合到最佳攻击位置的能量消耗总和，答案保留两位小数。

Sample Input

1
1
0 0

Sample Output

0.00

Data Constraint

第一题N<=50000,x,y都在longint范围内。

解题思路：这道题目解起来有些难度，如果要打正解，那就是三分了。
三分如此困难，肯定用水法。
一开始我是毫无思绪的，因为这N的范围实在太大，令人无法下手，可是，我们有分治法，把最佳位置（x，y）分成两次来求，然后嘛，就用中位数来求最佳咯。另外要注意精度问题，其实，这道题目的答案的小数点全是0，因为整数减整数之后才乘以一个整数，肯定还是整数，但用小数的话更容易答案对（不知道为什么，因为本人改成实数后就对了），然后就是如果小数部分无法变成00，建议取证然后输出个为“.00”的字符串，或者乘以10除以10然后向下取整。
代码如下：（仅供参考）

#include<cmath>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#define F(i,x,y) for(int i=x;i<=y;i++)
using namespace std;
const int maxn=50005;
const double p=0.00001;
int a[maxn][3];
int w[maxn];
int k,t,n,m;
double l,r,ans,num,m1,m2;
double f(double x)
{
        double s=0;
        F(i,1,n) s+=(w[i]*abs(x-a[i][1]));
        return s;
}
double g(double y)
{
        double s=0;
        F(i,1,n) s+=(w[i]*abs(y-a[i][2]));
        return s;
}

int main()
{
    freopen("data10.in","r",stdin);
    scanf("%d",&n);
    F(i,1,n) scanf("%d",&w[i]);

    F(i,1,n) scanf("%d %d",&a[i][1],&a[i][2]);

    l=-0x7fffffff;
    r=0x7fffffff;

    while(r-l>=p)
    {
        m1=(r-l)/3+l;
        m2=((r-l)/3)*2+l;
        if(f(m1)<f(m2)) r=m2; 
        else l=m1;
    }  
    ans+=f(l);
    l=-0x7fffffff;
    r=0x7fffffff;
    while(r-l>=p)
    {
        m1=(r-l)/3+l;
        m2=((r-l)/3)*2+l;
        if(g(m1)<g(m2)) r=m2;
        else l=m1;
    }
    ans+=g(l);

    printf("%.2lf\n",floor(ans*10)/10);
    return 0;
}