洛谷P1257 平面上的最接近点对 数学 分治 排序

Posted 2020-09-20 third2333

来自洛谷上的题解

方法一：

先求第1个点与其余n-1个点的距离；

再求第2个点与其余n-2个点的距离；

再求第3个点与其余n-3个点的距离；

…………………………………………

再求第n-1个点与其余1个点的距离；

然后找出最小值。

如此的算法复杂度为O(n^2)，显然不能满足本题的需要。但……貌似洛谷神机可以……过……

欢迎朴素的同学挑战数据加强版新

方法二：

考虑以下分治算法：

设平面上的点都在点集S中，为了将S线性分割为大小大致相等的2个子集S1和S2，我们选取一垂直线l（方程：x=m）来作为分割直线。其中m为S中各点x坐标的中位数。由此将S分割为S1={p∈S|px≤m}和S2={p∈S|px>m}。从而使S1和S2分别位于直线l的左侧和右侧，且S=S1∪S2 。由于m是S中各点x坐标值的中位数，因此S1和S2中的点数大致相等。 递归地在S1和S2上解最接近点对问题，我们分别得到S1和S2中的最小距离δ1和δ2。现设δ=min (δ1,δ2)。

若S的最接近点对(p,q)之间的距离d(p,q)<δ则p和q必分属于S1和S2。不妨设p∈S1，q∈S2。那么p和q距直线l的距离均小于δ。因此，我们若用P1和P2分别表示直线l的左边和右边的宽为δ的2个垂直长条，则p∈S1，q∈S2

此时，P1中所有点与P2中所有点构成的点对均为最接近点对的候选者。在最坏情况下有n^2/4对这样的候选者。但是P1和P2中的点具有以下的稀疏性质，它使我们不必检查所有这n2/4对候选者。考虑P1中任意一点p,它若与P2中的点q构成最接近点对的候选者，则必有d(p,q)<δ。满足这个条件的P2中的点有多少个呢？容易看出这样的点一定落在一个δ×2δ的矩形R中，

由δ的意义可知P2中任何2个S中的点的距离都不小于δ。由此可以推出矩形R中最多只有6个S中的点。事实上，我们可以将矩形R的长为2δ的边3等分，将它的长为δ的边2等分，由此导出6个（δ/2）×（2δ/3）的矩形。

因此d(u,v)≤5δ/6<δ 。这与δ的意义相矛盾。也就是说矩形R中最多只有6个S中的点。由于这种稀疏性质，对于P1中任一点p，P2中最多只有6个点与它构成最接近点对的候选者。因此，在分治法的合并步骤中，我们最多只需要检查6×n/2=3n对候选者，而不是n^2/4对候选者。

我们只知道对于P1中每个S1中的点p最多只需要检查P2中的6个点，但是我们并不确切地知道要检查哪6个点。为了解决这个问题，我们可以将p和P2中所有S2的点投影到垂直线l上。由于能与p点一起构成最接近点对候选者的S2中点一定在矩形R中，所以它们在直线l上的投影点距p在l上投影点的距离小于δ。由上面的分析可知，这种投影点最多只有6个。因此，若将P1和P2中所有S的点按其y坐标排好序，则对P1中所有点p，对排好序的点列作一次扫描，就可以找出所有最接近点对的候选者，对P1中每一点最多只要检查P2中排好序的相继6个点。

至此，我们用分治法求出平面最接近点对。

 其实这种方法我觉得应该是 nlogn^2 的  

不过可以优化成nlogn ,但其实多只log 也没多大区别233

 

 

 1 type rec=record
 2      x,y:int64;
 3      end;
 4 var a:array [0..200001] of rec;
 5     tt,ttt,tttt:array [0..200001] of rec;
 6     n,l,m,min:int64;
 7     i,j,k:longint;
 8 procedure sort(l,r:longint);
 9 var i,j:longint;
10     x,y:rec;
11 begin
12   i:=l;
13   j:=r;
14   x:=a[(l+r) div 2];
15   repeat
16     while (a[i].x<x.x) or ((a[i].x=x.x) and (a[i].y>x.y)) do inc(i);
17     while (x.x<a[j].x) or ((x.x=a[j].x) and (x.y>a[j].y)) do dec(j);
18     if not(i>j)
19     then
20       begin
21         y:=a[i];
22         a[i]:=a[j];
23         a[j]:=y;
24         inc(i);
25         j:=j-1;
26       end;
27   until i>j;
28   if l<j then sort(l,j);
29   if i<r then sort(i,r);
30 end;
31 procedure sortt(l,r:longint);
32 var i,j:longint;
33     x,y:rec;
34 begin
35   i:=l;
36   j:=r;
37   x:=tttt[(l+r) div 2];
38   repeat
39     while (tttt[i].y<x.y) do inc(i);
40     while (x.y<tttt[j].y) do dec(j);
41     if not(i>j)
42     then
43       begin
44         y:=tttt[i];
45         tttt[i]:=tttt[j];
46         tttt[j]:=y;
47         inc(i);
48         j:=j-1;
49       end;
50   until i>j;
51   if l<j then sortt(l,j);
52   if i<r then sortt(i,r);
53 end;
54 function dis(a,b:rec):int64;
55 begin
56   exit(sqr(a.x-b.x)+sqr(a.y-b.y));
57 end;
58 function minn(a,b:int64):int64;
59 begin if a<b then exit(a) else exit(b); end;
60 function work(le,ri:longint):int64;
61 var tl,tr,t1,t2,min:int64;
62     i,j,k,lenr,lenl,p1,p2:longint;
63 begin
64   if ri-le=1 then exit(dis(a[ri],a[le]));
65   if ri-le=2 then exit(minn(minn(dis(a[ri],a[le]),dis(a[ri-1],a[le])),dis(a[le+1],a[ri])));
66   tl:=work(le,(le+ri) div 2);
67   tr:=work(((le+ri) div 2)+1,ri);
68   min:=minn(tl,tr);
69   t1:=(le+ri) div 2;
70   t2:=t1+1;
71   while (t1>=le) and (sqr(a[t1].x-a[(le+ri) div 2].x)<=min) do dec(t1);
72   inc(t1);
73   while (t2<=ri) and (sqr(a[t2].x-a[(le+ri) div 2].x)<=min) do inc(t2);
74   dec(t2);
75   if (t2=(le+ri) div 2) or (t1=((le+ri) div 2)+1) then exit(min);
76   for k:=t1 to t2 do tttt[k]:=a[k];
77   sortt(t1,t2);
78   for k:=t1 to t2-1 do
79     begin
80       for j:=k+1 to minn(k+8,t2) do
81         begin
82           min:=minn(min,dis(tttt[k],tttt[j]));
83         end;
84     end;
85   exit(min);
86 end;
87 begin
88   readln(n);
89   for i:=1 to n do readln(a[i].x,a[i].y);
90   sort(1,n);
91   writeln(sqrt(work(1,n)):0:4);
92 end.