给定两个非连接(比如不相交)的凸多边形 P 和 Q, 目标是找到拥有最小距离的点对 (p,q) (p 属于 P 且 q 属于 Q)。
事实上, 多边形非连接十分重要, 因为我们所说的多边形包含其内部。 如果多边形相交, 那么最小距离就变得没有意义了。 然而, 这个问题的另一个版本, 凸多边形顶点对间最小距离对于相交和非相交的情况都有解存在。
回到我们的主问题: 直观的, 确定最小距离的点不可能包含在多边形的内部。 与最大距离问题相似, 我们有如下结论:
两个凸多边形 P 和 Q 之间的最小距离由多边形间的对踵点对确立。 存在凸多边形间的三种多边形间的对踵点对, 因此就有三种可能存在的最小距离模式: “顶点-顶点”的情况 “顶点-边”的情况 “边-边”的情况
换句话说, 确定最小距离的点对不一定必须是顶点。 下面的三个图例表明了以上结论: 给定结果, 一个基于旋转卡壳的算法自然而然的产生了: 考虑如下的算法, 算法的输入是两个分别有 m 和 n 个顺时针给定顶点的凸多边形 P 和 Q。 计算 P 上 y 坐标值最小的顶点(称为 yminP ) 和 Q 上 y 坐标值最大的顶点(称为 ymaxQ)。 为多边形在 yminP 和 ymaxQ 处构造两条切线 LP 和 LQ 使得他们对应的多边形位于他们的右侧。 此时 LP 和 LQ 拥有不同的方向, 并且 yminP 和 ymaxQ 成为了多边形间的一个对踵点对。 计算距离(yminP,ymaxQ) 并且将其维护为当前最小值。 顺时针同时旋转平行线直到其中一个与其所在的多边形的边重合。 如果只有一条线与边重合, 那么只需要计算“顶点-边”对踵点对和“顶点-顶点”对踵点对距离。 都将他们与当前最小值比较, 如果小于当前最小值则进行替换更新。 如果两条切线都与边重合, 那么情况就更加复杂了。 如果边“交叠”, 也就是可以构造一条与两条边都相交的公垂线(但不是在顶点处相交), 那么就计算“边-边”距离。 否则计算三个新的“顶点-顶点”对踵点对距离。 所有的这些距离都与当前最小值进行比较, 若小于当前最小值则更新替换。 重复执行步骤4和步骤5, 直到新的点对为(yminP,ymaxQ)。 输出最大距离。 旋转卡壳模式保证了所有的对踵点对(和所有可能的子情况)都被考虑到。 此外, 整个算法拥有现行的时间复杂度, 因为(除了初始化), 只有与顶点数同数量级的操作步数需要执行。
最小距离和最大距离的问题表明了旋转卡壳模型可以用在不同的条件下(与先前的直径和宽度问题比较)。 这个模型可以应用于两个多边形的问题中。 “最小盒子”问题(最小面积外接矩形)通过同一多边形上两个正交切线集合展示了另一种条件下旋转卡壳的应用。
转自http://blog.csdn.net/ACMaker/archive/2008/10/29/3178696.aspx
POJ 3608 #include<stdio.h> #include<algorithm> #include<math.h> #include<stdlib.h> #define eps 1e-8 using namespace std; struct pt { double x,y; }pp[2][30000],poly[2][30000]; double cross(pt a,pt b,pt c) { return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x); } int dblcmp(double x) { if(fabs(x)<eps) return 0; else return x<0?-1:1; } double dist(pt a,pt b) { return (b.x-a.x)*(b.x-a.x)+(b.y-a.y)*(b.y-a.y); } bool cmp1(pt a,pt b) { int x=dblcmp(cross(pp[0][0],a,b)); if(x==0) return dist(pp[0][0],a)<dist(pp[0][0],b); return x>0; } bool cmp2(pt a,pt b) { int x=dblcmp(cross(pp[1][0],a,b)); if(x==0) return dist(pp[1][0],a)<dist(pp[1][0],b); return x>0; } int p,q,p0,q0,n,m; int next1[30000],next2[30000]; double dot(pt a,pt b,pt c) { return (b.x-a.x)*(c.x-a.x)+(b.y-a.y)*(c.y-a.y); } double dis(pt a,pt b,pt c) { if(dot(a,b,c)>dist(a,b)+eps||dot(a,b,c)<-eps) return min(dist(c,b),dist(c,a)); else return fabs(cross(a,b,c)*cross(a,b,c)/dist(a,b)); } double polldis(int p,int q) { return min(min(dis(poly[0][p],poly[0][next1[p]],poly[1][q]),dis(poly[0][p],poly[0][next1[p]],poly[1][next2[q]])),min(dis(poly[1][q],poly[1][next2[q]],poly[0][p]),dis(poly[1][q],poly[1][next2[q]],poly[0][next1[p]]))); } int TB(int n,int s) { int swap=0; int i; for(i=0;i<n;i++) { if(pp[s][i].y<pp[s][swap].y||(pp[s][i].y==pp[s][swap].y&&pp[s][i].x<pp[s][swap].x)) swap=i; } poly[s][0]=pp[s][0]; pp[s][0]=pp[s][swap]; pp[s][swap]=poly[s][0]; if(s==0) sort(&pp[0][0],&pp[0][n],cmp1); else sort(&pp[1][0],&pp[1][n],cmp2); poly[s][0]=pp[s][0]; poly[s][1]=pp[s][1]; int top=1; for(i=2;i<n;i++) { while(top>0&&dblcmp(cross(poly[s][top-1],poly[s][top],pp[s][i]))<=0) top--; poly[s][++top]=pp[s][i]; } int tt=1; for(i=1;i<=top;i++) { if(!(poly[s][i].x==poly[s][i-1].x&&poly[s][i].y==poly[s][i-1].y)) poly[s][tt++]=poly[s][i]; } return tt; } double RC() { double mi=1e15; p=p0=0; q=0; while(cross(poly[0][p],poly[0][next1[p]],poly[1][next2[q]])>cross(poly[0][p],poly[0][next1[p]],poly[1][q])) q=next2[q]; q0=q; while(1) { p=next1[p]; if(dis(poly[0][p],poly[0][next1[p]],poly[1][q])<mi) mi=dis(poly[0][p],poly[0][next1[p]],poly[1][q]); while(cross(poly[0][p],poly[0][next1[p]],poly[1][next2[q]])>cross(poly[0][p],poly[0][next1[p]],poly[1][q])) { q=next2[q]; if(dis(poly[0][p],poly[0][next1[p]],poly[1][q])<mi) mi=dis(poly[0][p],poly[0][next1[p]],poly[1][q]); if(p==p0&&q==q0) { return mi; } } if(fabs(cross(poly[0][p],poly[0][next1[p]],poly[1][next2[q]])-cross(poly[0][p],poly[0][next1[p]],poly[1][q]))<eps) { mi=min(mi,polldis(p,q)); } if(p==p0) break; } return mi; } double RC2() { double mi=1e15; q=q0=0; p=0; while(cross(poly[1][q],poly[1][next2[q]],poly[0][next1[p]])>cross(poly[1][q],poly[1][next2[q]],poly[0][p])) p=next1[p]; p0=p; while(1) { q=next2[q]; if(dis(poly[1][q],poly[1][next2[q]],poly[0][p])<mi) mi=dis(poly[1][q],poly[1][next2[q]],poly[0][p]); while(cross(poly[1][q],poly[1][next2[q]],poly[0][next1[p]])>cross(poly[1][q],poly[1][next2[q]],poly[0][p])) { p=next1[p]; if(dis(poly[1][q],poly[1][next2[q]],poly[0][p])<mi) mi=dis(poly[1][q],poly[1][next2[q]],poly[0][p]); if(p==p0&&q==q0) { return mi; } } if(fabs(cross(poly[1][q],poly[1][next2[q]],poly[0][next1[p]])-cross(poly[1][q],poly[1][next2[q]],poly[0][p]))<eps) { mi=min(mi,polldis(p,q)); } if(q==q0) break; } return mi; } int main() { int i; while(scanf("%d%d",&n,&m)&&(n||m)) { for(i=0;i<n;i++) { scanf("%lf %lf",&pp[0][i].x,&pp[0][i].y); } for(i=0;i<m;i++) { scanf("%lf %lf",&pp[1][i].x,&pp[1][i].y); } int t1=n,t2=m; n=TB(t1,0); m=TB(t2,1); for(i=0;i<n;i++) next1[i]=(i+1)%n; for(i=0;i<m;i++) next2[i]=(i+1)%m; printf("%.5f\n",sqrt(min(RC(),RC2()))); } }
|