HDU 4311 Meeting point-1 && HDU 4312 Meeting point-2

这俩个题  题意：：给出N(<1e5)个点求找到一个点作为聚会的地方，使每个点到达这里的距离最小。4311是 曼哈顿距离 4312是 切比雪夫距离；

曼哈顿距离 ：大家都知道 对于二维坐标系a(xa,yb),b(xb,yb)的曼哈顿距离是abs（xa-xb）+abs（ya-yb）；

　　看的出来这个距离和x，y 都有关，但是X,Y并不相互影响，所以可以分开计算这样，分开计算的好处如下：

　　如果 只给一个维度的坐标系 ，我们是不可以再什么养的复杂度的时间内处理出来呢？ 大难还是很好想的先排序一下，会发现就需要O（n）的时间就可以算出来我们想要的答案

　　至于俩个维度的合并运算只要O（N）记录就可以了；

切比雪夫距离：

　　我们知道对于距离一个点的曼哈顿距离为R 的点可以形成一个正方形如下图：

　　也知道对于距离一个点的切比雪夫距离为R 的点也可以形成一个正方形如下图：

　　对比俩个图片我们知道：距离一个点的 切比雪夫距离可以通过一定的变换转换成 曼哈顿距离（因为曼哈顿距离怎么算我们上面已经讲过了）！！！；

　　比较麻烦？：但其实很简单：我们知道上面的 切比雪夫的正方形 旋转45度 在变换一定的倍数会 和 曼哈顿的正方形重合！所以我饿们可以认为 切比雪夫距离是  曼哈顿距离旋转在放大的结果：反过来只要对坐标进行同样的旋转和放大就可以了；

代码如下：：：

#include <cstdio>
#include <map>
#include <vector>
#include <set>
#include <iostream>
using namespace std;
typedef long long LL;
struct info
{
    LL x,y;
    int cnt;
    info(){}
    info(int x,int y):x(x),y(y){}
};
bool cmpx(info  a,info b)
{
    if(a.x==b.x)return a.y<b.y;
    return a.x<b.x;
}
bool cmpy(info a,info b)
{
    if(a.y==b.y) return a.x<b.x;
    return a.y<b.y;
}
info ko1[100005];
info ko2[100005];
LL dp1[100005];
LL dp2[100005];
int n;
void make(int x)
{
    int tmpx,tmpy;
    tmpx=ko1[x].x;
    tmpy=ko1[x].y;
//这么写是按远点点进行逆时针旋转的
//你也可以按别的点进行旋转  不过们的旋转点必须是一个
//列
//tmpx=ko1[x].x-2;
//tmpy=ko1[x].y+12;
    ko1[x].x=tmpx-tmpy;
    ko1[x].y=tmpy+tmpx;
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
        {
           // cin>>ko1[i].x>>ko1[i].y;
            scanf("%I64d%I64d",&ko1[i].x,&ko1[i].y);
            make(i);//有这句是12 没这句是11
        }
        sort(ko1+1,ko1+1+n,cmpx);
        for(int i=1;i<=n;i++)
        {
            ko1[i].cnt=i;
            ko2[i]=ko1[i];
        }
        sort(ko2+1,ko2+1+n,cmpy);
        for(int i=1;i<=n;i++)
        {
            dp1[i]=dp1[i-1]+abs(ko1[i].x-ko1[1].x);
            dp2[i]=dp2[i-1]+abs(ko2[i].y-ko2[1].y);
        }
        LL ans=1LL<<62;
        for(int i=1;i<=n;i++)
        {
            LL tmp2=dp2[n]-dp2[i]-(ko2[i].y-ko2[1].y)*(n-i)+(ko2[i].y-ko2[1].y)*i-dp2[i];
            int cnt=ko2[i].cnt;
            LL tmp1=dp1[n]-dp1[cnt]-(ko1[cnt].x-ko1[1].x)*(n-cnt)+(ko1[cnt].x-ko1[1].x)*cnt-dp1[cnt];
            ans=min(ans,(tmp1+tmp2)/2);//  这个12 要除以2；11不需要
        }
        cout<<ans<<endl;
    }
    return 0;
}

HDU 4311 Meeting point-1 && HDU 4312 Meeting point-2,布布扣,bubuko.com