[BZOJ 3170][Tjoi 2013]松鼠聚会

题解

我们可以发现, 我们走过一个对角线的距离也记为 $1$ , 这样的话我们可以直接补偿掉任意两个点间横坐标差与纵坐标差之间的差, 也就是说两点间的距离是这样的:

\[max(|x_0-x_1|,|y_0-y_1|)\]

实际上这是切比雪夫距离的定义, 但是单纯求所有其他点到某个点的切比雪夫距离和比较困难, 我们可以考虑按照一定规则转换点的坐标来使点与点间的曼哈顿距离等于原来的点之间的切比雪夫距离, 具体转换方式为 $(x,y) \rightarrow (x+y,x-y)$ , 转化前两点的切比雪夫距离为转化后两点间的曼哈顿距离的一半. 正确性很容易证明.

转化为曼哈顿距离后, 我们可以对所有横坐标排序并求前缀和, 然后对所有纵坐标排序求前缀和, 用前缀和减去目标的最终位置再求绝对值即可得到结果.

参考代码

GitHub

#include <cstdio>
#include <cstring>
#include <climits>
#include <cstdlib>
#include <iostream>
#include <algorithm>

const int MAXN=1e5+10;

struct Point{
long long x;
long long y;
long long posx;
long long posy;
};

int n;
long long sx[MAXN];
long long sy[MAXN];

Point P[MAXN];

void Initialize();
bool cmpx(const Point&,const Point&);
bool cmpy(const Point&,const Point&);

int main(){
Initialize();
std::sort(P+1,P+n+1,cmpx);
for(int i=1;i<=n;i++){
P[i].posx=i;
sx[i]=sx[i-1]+P[i].x;
}
std::sort(P+1,P+n+1,cmpy);
for(int i=1;i<=n;i++){
P[i].posy=i;
sy[i]=sy[i-1]+P[i].y;
}
long long ans=LLONG_MAX;
for(int i=1;i<=n;i++){
long long sum=0;
sum+=P[i].x*(P[i].posx)-sx[P[i].posx];
sum+=sx
-sx[P[i].posx]-(n-P[i].posx)*P[i].x;
sum+=P[i].y*(P[i].posy)-sy[P[i].posy];
sum+=sy
-sy[P[i].posy]-(n-P[i].posy)*P[i].y;
// printf("%lld\n",sum);
ans=std::min(ans,sum);
}
printf("%lld\n",ans/2);
return 0;
}

void Initialize(){
int x,y;
scanf("%d",&n);
for(int i=1;i<=n;i++){
scanf("%d%d",&x,&y);
P[i].x=x+y;
P[i].y=x-y;
}
}

inline bool cmpx(const Point& x,const Point& y){
return x.x<y.x;
}

inline bool cmpy(const Point& x,const Point& y){
return x.y<y.y;
}

Backup