码迷,mamicode.com
首页 > 其他好文 > 详细

旋转卡壳解poj3608Bridge Across Islands

时间:2014-09-04 20:57:40      阅读:397      评论:0      收藏:0      [点我收藏+]

标签:os   io   for   问题   代码   sp   amp   c   ef   

继续上一次的旋转卡壳的问题,这次是求两个凸包的最短距离,

其实选择卡壳就是只要找到"当前向量面积不小于下一个向量面积"即可,

满足这个条件,当前的两个ymin和ymax点就是一对对踵点了。

代码如下,欢迎参考:

#include<cmath>
#include<cstdio>
#include<iostream>
using namespace std;
typedef struct{float x,y;} Dot;
Dot operator -(Dot a,Dot b){Dot c={a.x-b.x,a.y-b.y};return c;}
float operator *(Dot a,Dot b){return a.x*b.y-b.x*a.y;}
float operator /(Dot a,Dot b){return a.x*b.x+a.y*b.y;}//点积判断这个点是否在直线外 
float dis(Dot a,Dot b){return sqrt(pow(a.x-b.x,2)+pow(a.y-b.y,2));}
float dtl(Dot a,Dot b,Dot c)//c点到线段ab的最短距离 
{
	if((b-a)/(c-a)<0) return dis(c,a);
	if((a-b)/(c-b)<0) return dis(c,b);
	return fabs((a-b)*(c-b)/dis(a,b));
}
float ltl(Dot a,Dot b,Dot c,Dot d){return min(min(dtl(a,b,c),dtl(a,b,d)),min(dtl(c,d,a),dtl(c,d,b)));}
float work(Dot * a,Dot * b,int n,int m)
{
	Dot t;
	float ans=1e9;
	int i,ymin=0,ymax=0;
	for(i=0;i<n;i++)
		if(a[i].y<a[ymin].y)
			ymin=i;
	for(i=0;i<m;i++)
		if(b[i].y>b[ymax].y)
		   ymax=i;
	a[n]=a[0];b[m]=b[0];
	for(i=0;i<n;i++)
	{
		t=a[ymin+1]-a[ymin];
		while(t*(b[ymax]-a[ymin])<t*(b[ymax+1]-a[ymin]))//找到当前向量面积不小于下一个向量面积即可 
			ymax=(ymax+1)%m;
		ans=min(ans,ltl(a[ymin],a[ymin+1],b[ymax],b[ymax+1]));
		ymin=(ymin+1)%n;
	}
	return ans;
}
int main()
{
	int i,n,m;
	Dot a[10010],b[10010];
	while(scanf("%d%d",&n,&m)&&n!=0||m!=0)
	{
		for(i=0;i<n;i++)
			scanf("%f%f",&a[i].x,&a[i].y);//单精度跑得快,而且精度是1e-6,对于此题,够用了 
		for(i=0;i<m;i++)
			scanf("%f%f",&b[i].x,&b[i].y);
		printf("%.3f\n",work(a,b,n,m));//精度在1e-3就够了 
	}
}


旋转卡壳解poj3608Bridge Across Islands

标签:os   io   for   问题   代码   sp   amp   c   ef   

原文地址:http://blog.csdn.net/stl112514/article/details/39057849

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!