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;
}

 

  

posted @ 2014-07-28 16:13  默默无语敲代码  阅读(227)  评论(0编辑  收藏  举报