hdu 4311 给n个二维点,选择一个点,使得其他点到这个点的曼哈顿距离最小。
hdu 4312 给n个二维点,选择一个点,使得其他点到这个点的切比雪夫距离最小。
很显然4311 是很简单的,直接一个前缀处理一下就可以了(分别考虑x,y)。
但是第 4312 就。。。
但是曼哈顿距离和切比雪夫距离是有一定的关系的!!!
参考博客: 戳我you
代码:
/// hdu 4311
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N =100005;
struct node
{
int f;
ll x,y;
ll Ans;
}a[N];
ll pre[N];
ll post[N];
int n;
bool cmp1(node a,node b)
{
return a.x<b.x;
}
bool cmp2(node a,node b)
{
return a.y<b.y;
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
memset(pre,0,sizeof(pre));
for(int i=1;i<=n;i++) scanf("%lld %lld",&a[i].x, &a[i].y);
for(int i=1;i<=n;i++) a[i].Ans=0;
sort(a+1,a+n+1,cmp1);
for(int i=1;i<=n;i++) a[i].f=0;
for(int i=1;i<=n;i++){
pre[i]=pre[i-1]+a[i].x;
}
/*
for(int i=1;i<=n;i++){
cout<<pre[i]<<" ";
}
cout<<endl;
*/
ll cnt;
for(int i=1;i<=n;i++){
cnt=i;
a[i].Ans+=(cnt*a[i].x-pre[i]);
}
for(int i=1;i<=n;i++){
if(a[i].x!=a[i+1].x){
cnt=n-i;
ll tmp=pre[n]-pre[i]-cnt*a[i].x;
for(int j=i;j>=1;j--){
if(a[j].f) break;
a[j].f=1;
a[j].Ans+=tmp;
}
}
}
//for(int i=1;i<=n;i++) cout<<"ans "<<a[i].Ans<<endl;
sort(a+1,a+n+1,cmp2);
for(int i=1;i<=n;i++){
a[i].f=0;
}
memset(pre,0,sizeof(pre));
for(int i=1;i<=n;i++){
pre[i]=pre[i-1]+a[i].y;
}
for(int i=1;i<=n;i++){
cnt=i;
a[i].Ans+=(cnt*a[i].y-pre[i]);
}
for(int i=1;i<=n;i++){
if(a[i].y!=a[i+1].y){
cnt=n-i;
ll tmp=pre[n]-pre[i]-cnt*a[i].y;
for(int j=i;j>=1;j--){
if(a[j].f) break;
a[j].f=1;
a[j].Ans+=tmp;
}
}
}
ll Ans=1e18+5;
for(int i=1;i<=n;i++){
Ans=min(Ans,a[i].Ans);
}
printf("%lld\n",Ans);
}
return 0;
}
/*
4
6
-4 -1
-1 -2
2 -4
0 2
0 3
5 -2
*/
///hdu 4312
/// 曼哈顿距离中的(x,y) 对应切比雪夫距离中的((x+y)/2,(x-y)/2)
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N =100005;
struct node
{
int f;
ll x,y;
ll Ans;
}a[N];
ll pre[N];
ll post[N];
int n;
bool cmp1(node a,node b)
{
return a.x<b.x;
}
bool cmp2(node a,node b)
{
return a.y<b.y;
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
memset(pre,0,sizeof(pre));
for(int i=1;i<=n;i++) scanf("%lld %lld",&a[i].x, &a[i].y);
for(int i=1;i<=n;i++){
ll tx,ty;
tx=a[i].x; ty=a[i].y;
a[i].x=tx+ty; a[i].y=tx-ty;
}
for(int i=1;i<=n;i++) a[i].Ans=0;
sort(a+1,a+n+1,cmp1);
for(int i=1;i<=n;i++) a[i].f=0;
for(int i=1;i<=n;i++){
pre[i]=pre[i-1]+a[i].x;
}
/*
for(int i=1;i<=n;i++){
cout<<pre[i]<<" ";
}
cout<<endl;
*/
ll cnt;
for(int i=1;i<=n;i++){
cnt=i;
a[i].Ans+=(cnt*a[i].x-pre[i]);
}
for(int i=1;i<=n;i++){
if(a[i].x!=a[i+1].x){
cnt=n-i;
ll tmp=pre[n]-pre[i]-cnt*a[i].x;
for(int j=i;j>=1;j--){
if(a[j].f) break;
a[j].f=1;
a[j].Ans+=tmp;
}
}
}
//for(int i=1;i<=n;i++) cout<<"ans "<<a[i].Ans<<endl;
sort(a+1,a+n+1,cmp2);
for(int i=1;i<=n;i++){
a[i].f=0;
}
memset(pre,0,sizeof(pre));
for(int i=1;i<=n;i++){
pre[i]=pre[i-1]+a[i].y;
}
for(int i=1;i<=n;i++){
cnt=i;
a[i].Ans+=(cnt*a[i].y-pre[i]);
}
for(int i=1;i<=n;i++){
if(a[i].y!=a[i+1].y){
cnt=n-i;
ll tmp=pre[n]-pre[i]-cnt*a[i].y;
for(int j=i;j>=1;j--){
if(a[j].f) break;
a[j].f=1;
a[j].Ans+=tmp;
}
}
}
ll Ans=1e18+5;
for(int i=1;i<=n;i++){
Ans=min(Ans,a[i].Ans);
}
printf("%lld\n",Ans/2);
}
return 0;
}
/*
*/