【解题报告】 Ural 1348 简单的计算几何

最新推荐文章于 2021-05-26 21:08:30 发布

x314542916

最新推荐文章于 2021-05-26 21:08:30 发布

阅读量578

点赞数 1

CC 4.0 BY-SA版权

分类专栏：【计算几何】文章标签： vector c

本文链接：https://blog.youkuaiyun.com/x314542916/article/details/7824740

【计算几何】专栏收录该内容

11 篇文章

订阅专栏

这篇博客详细介绍了如何解决Ural 1348问题，涉及计算几何的运用。通过计算直线距离和向量积，确定羊为了吃到水果需要伸长绳子的最小和最大长度。代码中包含判断double型值大小的重点部分。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目连接： Ural 1348
题目大意：有一只羊被绳子栓到了,有一块水果地，视为一个线段，我们输入线段的起点坐标和终点坐标以及绳子的长度。让你计算如果羊要吃到水果和要吃完水果所需要伸长绳子的最小和最大长度。

// 计算几何
// 重点就是判断double型值的大小 
//
#include <stdio.h>
#include <math.h>
#define MIN  0.0000001
#define MINF -0.0000001

double max(double a, double b, double c){
	double max;
	if (a - b > MIN)	max = a;
	else	max = b;
	if (max - c > MIN)	return max;
	else	return c;
}

double min(double a, double b, double c){
	double min;
	if (b - a > MIN)	min = a;
	else	min = b;
	if (c - min > MIN)	return min;
	else	return c;
}

double Abs(double x){ // 计算绝对值
	if (x < MINF ) return -x;
	return x;
}

double kkk(double x1,double y1,double x2,double y2){ // 斜率
	return ( (y1 - y2) / (x1 - x2) );
}

double length(double x1,double y1,double x2,double y2,double x0,double y0){
	// (x0,y0)点到直线(两点式)的距离
	if (x1 - x2 < MIN && x1 - x2 > MINF) return Abs(x1-x0); // 如果斜率不存在
	double k = kkk(x1,y1,x2,y2);
	return (Abs(k*x0 - y0 + y2 - k*x2)) / (sqrt(k*k + 1));
}

double lengthline(double x1,double y1,double x2,double y2){
	return  sqrt( (x1-x2)*(x1-x2) + (y1-y2)*(y1-y2) );
}

double vector(double x1,double y1,double x2,double y2){ // 传进来的值为两个向量 注意方向
	// 计算向量积
	return (x1*x2 + y1*y2);
}

int main()
{
	//freopen("in.txt","r",stdin);
	double x1,x2,x0,y1,y2,y0,len;
	while(scanf("%lf%lf%lf%lf%lf%lf%lf",&x1,&y1,&x2,&y2,&x0,&y0,&len) != EOF){
		double len1,len2,len3; // 两个端点距离 + 一个点到直线的距离
		double vec1,vec2; // 存在两个角，所以两个向量积
		double minlen,maxlen;
		len1 = lengthline(x1,y1,x0,y0);
		len2 = lengthline(x2,y2,x0,y0); 
		vec1 = vector(x1-x2,y1-y2,x0-x2,y0-y2);
		vec2 = vector(x2-x1,y2-y1,x0-x1,y0-y1);
		if (vec1 > MIN && vec2 > MIN){ // 存在钝角 或直角
			len3 = length(x1,y1,x2,y2,x0,y0);
			minlen = min(len1,len2,len3);
			maxlen = max(len1,len2,len3);
		}else{
			if (len1 - len2 > MIN){
				minlen = len2;
				maxlen = len1;
			}else{
				minlen = len1;
				maxlen = len2;
			}
		}
		double resa,resb;
		if ( minlen - len < MINF )	
			resa = 0;
		else 	resa = minlen - len;
		if ( maxlen - len < MINF )	
			resb = 0;
		else	resb = maxlen - len;
		printf("%.2lf\n%.2lf\n",resa,resb);
	}
	return 0;
}