POJ 3608 旋转卡壳求解凸包间的最短距离

题意给你两个凸包,求解两个凸包间的最短距离。


#include <iostream>
#include <vector>
#include <algorithm>
#include <cmath>
#include <cstdio>
using namespace std;

#define MAX_N 10000+10
#define INF 0x3F3F3F3F
#define EPS 1E-10

struct Point
{
	double x, y;
	Point() {}
	Point(double x, double y) : x(x), y(y) {}
	Point operator + (const Point& p){ return Point(x + p.x, y + p.y); }
	Point operator - (const Point& p){ return Point(x - p.x, y - p.y); }
	Point operator * (const double& d){ return Point(x * d, y * d); }
	bool operator < (const Point& a) const
	{
		if (x != a.x) return x < a.x;
		else return y < a.y;
	}
	double dot(const Point& p) { return x * p.x + y * p.y; }
	double det(const Point& p) { return x * p.y - y * p.x; }
};

Point P[MAX_N], Q[MAX_N];


// 字典序比较
bool cmp_x(Point a, Point b)
{
	if (a.x != b.x) return a.x < b.x;
	return a.y < b.y;
}


// 向量AB 与 AC 的叉积 如果叉积大于0,那么C在向量AB的逆时针方向,叉积小于0则在AB的顺时针方向。如果叉积等于0,则ABC共线。
inline double cross(Point A, Point B, Point C)
{
	return (B - A).det(C - A);
}

// 向量AB 与 AC 的点积 如果点积的结果为0,那么这两个向量互相垂直
inline double multi(Point A, Point B, Point C)
{
	return (B - A).dot(C - A);
}

// 两点距离
inline double dist(Point A, Point B)
{
	return sqrt((A.x - B.x) * (A.x - B.x) + (A.y - B.y) * (A.y - B.y));
}

// 逆时针排序
inline void anticlockwise_sort(Point* p, int N)
{
	for (int i = 0; i < N - 2; ++i)
	{
		double tmp = cross(p[i], p[i + 1], p[i + 2]);
		if (tmp > EPS) return;
		else if (tmp < -EPS)
		{
			reverse(p, p + N);
			return;
		}
	}
}

//计算C点到线段AB的最短距离
inline double point_to_line(Point A, Point B, Point C)
{
	if (dist(A, B) < EPS) return dist(B, C);
	if (multi(A, B, C) < -EPS) return dist(A, C);
	if (multi(B, A, C) < -EPS) return dist(B, C);
	return fabs(cross(A, B, C) / dist(A, B));
}

//求一条线段的两端点到另外一条线段的距离,反过来一样,共4种情况
inline double line_to_line(Point A, Point B, Point C, Point D)
{
	return min(min(point_to_line(A, B, C), point_to_line(A, B, D)), min(point_to_line(C, D, A), point_to_line(C, D, B)));
}


//求解两凸包之间的最短距离
inline double min_dis_2(Point* P, Point* Q, int n, int m)
{
	int yminP = 0, ymaxQ = 0;
	for (int i = 0; i < n; ++i) if (P[i].y < P[yminP].y) yminP = i;	// P上y坐标最小的顶点
	for (int i = 0; i < m; ++i) if (Q[i].y > Q[ymaxQ].y) ymaxQ = i; // Q上y坐标最大的顶点
	P[n] = P[0];	// 为了方便,避免求余
	Q[m] = Q[0];
	double arg, ans = INF;
	for (int i = 0; i < n; ++i)
	{
		while (arg = cross(P[yminP + 1], Q[ymaxQ + 1], P[yminP]) - cross(P[yminP + 1], Q[ymaxQ], P[yminP]) > EPS) ymaxQ = (ymaxQ + 1) % m;
		ans = min(ans, line_to_line(P[yminP], P[yminP + 1], Q[ymaxQ], Q[ymaxQ + 1]));
		yminP = (yminP + 1) % n;
	}
	return ans;
}



int main()
{
	int N, M;
	while (scanf("%d%d", &N, &M))
	{
	    if(N==0&&M==0) break;
		for (int i = 0; i < N; ++i)
		{
			scanf("%lf%lf", &P[i].x, &P[i].y);
		}
		for (int i = 0; i < M; ++i)
		{
			scanf("%lf%lf", &Q[i].x, &Q[i].y);
		}
		anticlockwise_sort(P, N);
		anticlockwise_sort(Q, M);
		printf("%.5f\n", min_dis_2(P, Q, N, M));
	}
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值