两条线段距离

线段间最小距离

时间限制(普通/Java):1000MS/3000MS          运行内存限制:65536KByte

描述

作为程序员的沈同学深深地知道大学所学课程的重要性，没有不重要的课程，当我们不知道如何应用所学知识的时候，我们就会觉得该知识无用。学习完高数后，沈同学发现现实当中的好多问题都可以用高数来求解，高数并不像很多同学想象中的那么没用。一天，他在hdu OJ上A题，发现了一种新算法叫做旋转卡壳法，这种算法对于他来说还是有点难，但是有一个问题他觉得自己已经有能力解决，就是求解两线段中两点的最小距离。

输入

输入多组数据，每组二行，分别是线段两点的x1,y1,x2,y2。（0<=x1,y1,x2,y2<10000)

输出

对于每组测试数据，输出一行,输出结果保留小数点后2位。

样例输入

0 0 0 3
1 1 1 2
0 0 0 1
1 1 1 2
0 0 1 1
0 1 1 0

样例输出

提示

1.00
1.00
0.00

题目来源

caisz

#include<iostream>
#include<cmath>
using namespace std;

struct N
{
	double x,y;
}node[4];

double Max(double a,double b)
{
	return a>b?a:b;
}

double Min(double a,double b)
{
	return a<b?a:b;
}

int jude(N a,N b,N c,N d)
{
	if(Min(a.x,b.x)>Max(c.x,d.x))
	{
		if(Min(a.y,b.y)>Max(c.y,d.y)||Max(a.y,b.y)<Min(c.y,d.y))
		{
			return 0;
		}
	}
	if(Max(a.x,b.x)<Min(c.x,d.x))
	{
		if(Min(a.y,b.y)>Max(c.y,d.y)||Max(a.y,b.y)<Min(c.y,d.y))
		{
			return 0;
		}
	}
	return 1;
}

double cj(N a,N b,N c)
{
	return (a.x-b.x)*(c.y-b.y)-(c.x-b.x)*(a.y-b.y);
}

double dis(N a,N b)
{
	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

double pan(N a,N b,N c)
{
	return (a.x-b.x)*(c.x-b.x)+(a.y-b.y)*(c.y-b.y);
}

int main()
{
	double MIN;
	while(scanf("%lf%lf%lf%lf%lf%lf%lf%lf",&node[0].x,&node[0].y,&node[1].x,&node[1].y,&node[2].x,&node[2].y,&node[3].x,&node[3].y)!=EOF)
	{
		MIN=9999999;
		if(jude(node[0],node[1],node[2],node[3]))
		{
			if(cj(node[2],node[0],node[1])*cj(node[1],node[0],node[3])>=0&&cj(node[1],node[2],node[3])*cj(node[3],node[2],node[0])>=0)
			{
				printf("0.00\n");
				continue;
			}
		}
		MIN=Min(dis(node[0],node[2]),MIN);
		MIN=Min(dis(node[0],node[3]),MIN);
		MIN=Min(dis(node[1],node[2]),MIN);
		MIN=Min(dis(node[1],node[3]),MIN);
		if(pan(node[0],node[2],node[3])>0&&pan(node[0],node[3],node[2])>0)//node[0]点
		{
			MIN=Min(abs(cj(node[2],node[0],node[3]))/dis(node[2],node[3]),MIN);
		}
		if(pan(node[1],node[2],node[3])>0&&pan(node[1],node[3],node[2])>0)
		{
			MIN=Min(abs(cj(node[2],node[1],node[3]))/dis(node[2],node[3]),MIN);
		}
		if(pan(node[2],node[0],node[1])>0&&pan(node[2],node[1],node[0])>0)
		{
			MIN=Min(abs(cj(node[0],node[2],node[1]))/dis(node[0],node[1]),MIN);
		}
		if(pan(node[3],node[0],node[1])>0&&pan(node[3],node[1],node[0])>0)
		{
			MIN=Min(abs(cj(node[0],node[3],node[1]))/dis(node[0],node[1]),MIN);
		}
		printf("%.2lf\n",MIN);
	}
	return 0;
}