题意:如下的图,给出N个点,每个点都可以沿着线进行移动,时间花费为线的长度,问你N个点到其中一个点集中,每个点都沿着最短路移动,问最后到的点最晚的时候(注意两个最),起点和最后到的点是哪两个点。(任意一组答案)
范围:点的坐标绝对值小于10的7次,点至多为20W个。
解法:在求解曼哈顿最小距离生成树的时候,有一个求解在曼哈顿距离意义下的最近点对的过程,和这个原理相同,可以求解在这张图意义下的最远点对。
(具体的步骤不再赘述,可以参考最小曼哈顿距离生成树的论文,是用线段树进行信息维护的)复杂度为4*N*logN
代码:
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<math.h>
#include<iostream>
#include<stdlib.h>
#include<set>
#include<map>
#include<queue>
#include<vector>
#include<bitset>
#pragma comment(linker, "/STACK:1024000000,1024000000")
template <class T>
bool scanff(T &ret){ //Faster Input
char c; int sgn; T bit=0.1;
if(c=getchar(),c==EOF) return 0;
while(c!='-'&&c!='.'&&(c<'0'||c>'9')) c=getchar();
sgn=(c=='-')?-1:1;
ret=(c=='-')?0:(c-'0');
while(c=getchar(),c>='0'&&c<='9') ret=ret*10+(c-'0');
if(c==' '||c=='\n'){ ret*=sgn; return 1; }
while(c=getchar(),c>='0'&&c<='9') ret+=(c-'0')*bit,bit/=10;
ret*=sgn;
return 1;
}
#define inf 1073741823
#define llinf 4611686018427387903LL
#define PI acos(-1.0)
#define lth (th<<1)
#define rth (th<<1|1)
#define rep(i,a,b) for(int i=int(a);i<=int(b);i++)
#define drep(i,a,b) for(int i=int(a);i>=int(b);i--)
#define gson(i,root) for(int i=ptx[root];~i;i=ed[i].next)
#define tdata int testnum;scanff(testnum);for(int cas=1;cas<=testnum;cas++)
#define mem(x,val) memset(x,val,sizeof(x))
#define mkp(a,b) make_pair(a,b)
#define findx(x) lower_bound(b+1,b+1+bn,x)-b
#define pb(x) push_back(x)
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
#define NN 400400
int n,m,k,b[NN],bn,dn;
struct point{
int x,y,idx;
bool operator < (const point temp)const{
if(x!=temp.x)return x<temp.x;
else return y<temp.y;
}
}a[NN],c[NN];
double sqrt2=sqrt(2.0);
struct node{
double val;
int idx;
}t[1080000];
void built(){
for(m=1;m<bn+2;m<<=1);
rep(i,1,m<<1)t[i].val=-(inf*1.0);
}
void add(int x,int y,int idx){
int pos=findx(y-x);
double val=double(y-x)+double(x)*sqrt2;
for(int i=pos+m;i>0;i>>=1){
if(val>t[i].val){
t[i].val=val;
t[i].idx=idx;
}
}
}
node query(int x,int y){
node ans;
double val=-(inf*1.0);
int l=findx(y-x);
int r=m-1,idx=0;
for(l=l+m-1,r=r+m+1;l^r^1;l>>=1,r>>=1){
if(~l&1)if(t[l^1].val>val)val=t[l^1].val,idx=t[l^1].idx;
if(r&1) if(t[r^1].val>val)val=t[r^1].val,idx=t[r^1].idx;
}
ans.val=val;
ans.idx=idx;
return ans;
}
double caldis(int x1,int y1,int x2,int y2){
int dx=abs(x2-x1);
int dy=abs(y2-y1);
int minx=min(dx,dy);
int maxx=max(dx,dy);
return double(minx)*sqrt2+double(maxx-minx);
}
void solve(){
dn=0;
rep(i,1,n)c[i]=a[i];
double ans=-(inf*1.0);
int ax,ay;
rep(fx,1,4){
rep(i,1,n){
if(fx==2||fx==4)swap(c[i].x,c[i].y);
if(fx==3)c[i].x=-c[i].x;
}
bn=0;
rep(i,1,n)b[++bn]=c[i].y-c[i].x;
sort(b+1,b+1+bn);
bn=unique(b+1,b+1+bn)-b-1;
sort(c+1,c+1+n);
built();
add(c[n].x,c[n].y,c[n].idx);
drep(i,n-1,1){
node temp=query(c[i].x,c[i].y);
int a1=c[i].idx;
int a2=temp.idx;
if(a1&&a2){
double dis=caldis(a[a1].x,a[a1].y,a[a2].x,a[a2].y);
if(dis>ans){
ans=dis;
ax=c[i].idx;
ay=temp.idx;
}
}
add(c[i].x,c[i].y,c[i].idx);
}
}
if(n==1)printf("%d %d\n",1,1);
else if(n==2)printf("%d %d\n",1,2);
else printf("%d %d\n",ax,ay);
}
int main(){
scanff(n);
rep(i,1,n)scanff(a[i].x),scanff(a[i].y),a[i].idx=i;
solve();
return 0;
}