[hihoCoder 1150] 基站选址(几何题，中心点)

描述

需要在一个N × M的网格中建立一个通讯基站，通讯基站仅必须建立在格点上。

网格中有A个用户，每个用户的通讯代价是用户到基站欧几里得距离的平方。

网格中还有B个通讯公司，维护基站的代价是基站到最近的一个通讯公司的路程（路程定义为曼哈顿距离）。

在网格中建立基站的总代价是用户通讯代价的总和加上维护基站的代价，最小总代价。

输入

第一行为一个整数T，表示数据组数。

每组数据第一行为四个整数：N, M, A, B。

接下来的A+B行每行两个整数x, y，代表一个坐标，前A行表示各用户的坐标，后B行表示各通讯公司的坐标。

输出

对于每组数据输出一行"Case #X: Y"，X代表数据编号（从1开始），Y代表所求最小代价。

数据范围

1 ≤ T ≤ 20

1 ≤ x ≤ N

1 ≤ y ≤ M

1 ≤ B ≤ 100

小数据

1 ≤ N, M ≤ 100

1 ≤ A ≤ 100

大数据

1 ≤ N, M ≤ 107

1 ≤ A ≤ 1000

样例输入

2
3 3 4 1
1 2
2 1
2 3
3 2
2 2
4 4 4 2
1 2
2 4
3 1
4 3
1 4
1 3

样例输出

Case #1: 4
Case #2: 13

题解：根据用户考察四个中心点，距离平方比曼哈顿距离更重要。想法大胆，可否证明？

#include<iostream>
#include<cstdio>
#include<cmath>
#include<climits>
#include<algorithm>
struct node {
    int x,y;
}p[1001],q[101];
using namespace std;
int main() {
    int t;
    scanf("%d",&t);
    for(int cnt=1;cnt<=t;cnt++) {
        int n,m,a,b;
        scanf("%d%d%d%d",&n,&m,&a,&b);
        long tmp_x=0,tmp_y=0;
        for(int i=0;i<a;i++) {
            scanf("%d%d",&p[i].x,&p[i].y);
            tmp_x+=p[i].x;
            tmp_y+=p[i].y;
        }
        for(int i=0;i<b;i++) {
            scanf("%d%d",&q[i].x,&q[i].y);
        }
        tmp_x=tmp_x/a;
        tmp_y=tmp_y/a;
        int dir[][2]={0,0,1,0,0,1,1,1};
        int center_x,center_y;
        long long res=-1;
        for(int i=0;i<4;i++) {
            int ne_x=tmp_x+dir[i][0];
            int ne_y=tmp_y+dir[i][1];
            long long dis=0;
            for(int k=0;k<a;k++) {
                dis=dis+(long long)(p[k].x-ne_x)*(p[k].x-ne_x)+(long long)(p[k].y-ne_y)*(p[k].y-ne_y);
            }
            long long dis2=-1;
            for(int k=0;k<b;k++) {
                if(dis2==-1||abs(q[k].x-ne_x)+abs(q[k].y-ne_y)<dis2) {
                    dis2=abs(q[k].x-ne_x)+abs(q[k].y-ne_y);
                }
            }
            if(res==-1||dis+dis2<res) {
                center_x=ne_x;
                center_y=ne_y;
                res=dis+dis2;
            }
        }
        printf("Case #%d: %lld\n",cnt,res);
    }
    return 0;
}