UVA - 11796 - Dog Distance (计算几何~)

不得不感叹,计算几何真是太美丽了!!

AC代码:

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

struct Point {
	double x, y;
	Point(double x = 0, double y = 0) : x(x) , y(y) { }
};

typedef Point Vector;  

Vector operator + (Vector A, Vector B) { return Vector(A.x+B.x, A.y+B.y); }
Vector operator - (Vector A, Vector B) { return Vector(A.x-B.x, A.y-B.y); }
Vector operator * (Vector A, double p) { return Vector(A.x*p, A.y*p); }
Vector operator / (Vector A, double p) { return Vector(A.x/p, A.y/p); } 

bool operator < (const Point& a, const Point& b) {
	return a.x < b.x || (a.x == b.x && a.y < b.y);
} 

const double eps = 1e-10;
int dcmp(double x) {
	if(fabs(x) < eps) return 0; else return x < 0 ? -1 : 1;
}

bool operator == (const Point& a, const Point& b) {
	return dcmp(a.x - b.x) == 0 && dcmp(a.y - b.y) == 0;
}

double Dot(Vector A, Vector B) { return A.x*B.x + A.y*B.y; }
double Length(Vector A) { return sqrt(Dot(A, A)); }
double Angle(Vector A, Vector B) { return acos(Dot(A, B) / Length(A) / Length(B)); } 

double Cross(Vector A, Vector B) { return A.x*B.y - A.y*B.x; }
double Area2(Point A, Point B, Point C) { return Cross(B-A, C-A); }

Vector Rotate(Vector A, double rad) {
	return Vector(A.x*cos(rad) - A.y*sin(rad), A.x*sin(rad)+A.y*cos(rad) );
} 

Vector Normal(Vector A) {
    double L = Length(A);
    return Vector(-A.y/L, A.x/L);
}

Point GetLineIntersection(Point P, Vector v, Point Q, Vector w) {
	Vector u = P - Q;
	double t = Cross(w, u) / Cross(v, w);
	return P + v * t;
} 

double DistanceToLine(Point P, Point A, Point B) {
    Vector v1 = B-A, v2 = P - A;
    return fabs(Cross(v1,v2) / Length(v1));
}  

double DistanceToSegment(Point P, Point A, Point B) {
    if(A==B) return Length(P-A);
    Vector v1 = B - A, v2 = P - A, v3 = P - B;
    if(dcmp(Dot(v1, v2)) < 0) return Length(v2);
    else if(dcmp(Dot(v1, v3)) > 0) return Length(v3);
    else return fabs(Cross(v1, v2)) / Length(v1);
}  

Point GetLineProjection(Point P, Point A, Point B) {
	Vector v = B - A;
	return A + v * ( Dot(v, P-A) / Dot(v, v) );
}  

bool SegmentProperIntersection(Point a1, Point a2, Point b1, Point b2) {
	double c1 = Cross(a2 - a1, b1 - a1), c2 = Cross(a2 - a1, b2 - a1),
			c3 = Cross(b2 - b1, a1 - b1), c4 = Cross(b2 - b1, a2 - b1);
	return dcmp(c1) * dcmp(c2) < 0 && dcmp(c3) * dcmp(c4) < 0;
} 

bool OnSegment(Point p, Point a1, Point a2) {
	return dcmp(Cross(a1 - p, a2 - p)) == 0 && dcmp(Dot(a1 - p, a2 - p)) < 0;
} 

double ConvexPolygonArea(Point* p, int n) {
    double area = 0;
    for(int i = 1; i < n-1; i++)
        area += Cross(p[i] - p[0], p[i + 1] - p[0]);
    return area / 2;
} 

const int maxn = 60;
int T, A, B;
Point P[maxn], Q[maxn];
double Min, Max;

void update(Point P, Point A, Point B) {
	Min = min(Min, DistanceToSegment(P, A, B));
	Max = max(Max, Length(P-A));
	Max = max(Max, Length(P-B));
}

int main() {
	scanf("%d", &T);
	for(int cas = 1; cas <= T; cas++) {
		scanf("%d %d", &A, &B);
		for(int i = 0; i < A; i++) scanf("%lf %lf", &P[i].x, &P[i].y);
		for(int i = 0; i < B; i++) scanf("%lf %lf", &Q[i].x, &Q[i].y);

		double lenA = 0, lenB = 0;
		for(int i = 0; i < A-1; i++) lenA += Length(P[i+1]-P[i]);
		for(int i = 0; i < B-1; i++) lenB += Length(Q[i+1]-Q[i]);

		int Sa = 0, Sb = 0;//Sa,Sb为刚经过的拐点的编号
		Point Pa = P[0], Pb = Q[0];//Pa,Pb为甲乙现在的位置
		Min = 1e9, Max = -1e9;

		while(Sa < A-1 && Sb < B-1) {
			double La = Length(P[Sa+1] - Pa);//甲到下一个拐点的距离
			double Lb = Length(Q[Sb+1] - Pb);//乙到下一个拐点的距离
			double T = min(La/lenA, Lb/lenB);//可以先让甲和乙的速度分别是lenA和lenB

			Vector Va = (P[Sa + 1] - Pa) / La * T * lenA;//甲的位移向量,(P[Sa+1]-Pa)/La为单位向量
			Vector Vb = (Q[Sb + 1] - Pb) / Lb * T * lenB;//乙的位移向量,(Q[Sb+1]-Pb)/Lb为单位向量

			update(Pa, Pb, Pb + Vb - Va);				 //更新最大最小距离 

			Pa = Pa + Va;
			Pb = Pb + Vb;
			if(Pa == P[Sa+1]) Sa++;
			if(Pb == Q[Sb+1]) Sb++;
		}
		printf("Case %d: %.0lf\n", cas, Max - Min);//%.0lf可以表达四舍五入
	}
	return 0;
} 
时间: 2024-10-14 02:09:37

UVA - 11796 - Dog Distance (计算几何~)的相关文章

uva 11796 Dog Distance (计算几何-点和直线)

C Dog Distance Input Standard Input Output Standard Output Two dogs, Ranga and Banga, are running randomly following two different paths. They both run for T seconds with different speeds. Ranga runs with a constant speed of R m/s, whereas Banga runs

相对运动模拟:UVa 11796 Dog Distance

11796 - Dog Distance 题目大意 甲和乙两条狗分别沿着给定的折线匀速奔跑,同时出发.同时到达,试计算甲和乙在奔跑过程中的最远距离和最近距离的差. 解题思路 模拟整个过程,把总过程分割成一个个子过程:甲和乙的路线都是一条线段,因为运动是相对的,因此也可以认为甲静止不动,乙自己沿着某个相对运动的直线(用相对位移向量表示)走,因此问题转化为点到线段的最小最大距离,注意子过程的最小最大值不代表整个过程的最小最大值,所以每处理完一个子过程就要不断更新. 假设现在甲的位置是Pa,刚经过编号

●UVA 11796 Dog Distance

题链: https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2896 题解: 计算几何,骚操作 对于简单情况,即两只狗的路线均为一条线段, 可以从相对运动的角度出发,考虑一直狗不动,另一只狗在运动. 而由于两只狗此时都在做匀速直线运动,所以在那只不懂的狗看来,另一直狗也在匀速直线前行.(物理老师:速度的矢量和!) 所以这个简单情况下,问题就变为

UVa 11796 Dog Distance

题意: 有甲乙两条狗分别沿着一条折线奔跑,已知它们同时从各自的起点出发,同时到达各自的终点.求整个过程中两条狗的最大距离Max与最小距离Min的差值. 分析: 假设甲乙的路线都是一条线段的简单情况.运动是相对的,我们假定甲不动,乙相对甲的运动也是匀速直线运动.所以就将问题转化成了点到直线的最小和最大距离. 甲或乙每到达一个拐点所对应的时刻称作“关键点”,那么没两个关键点之间的运动都可看做上面分析的简单情况.我们只要及时更新甲乙的位置即可. LenA和LenB分别是两条路线的总长,因为运动时间相同

简单几何(相对运动距离最值) UVA 11796 Dog Distance

题目传送门 题意:两只狗在折线上跑,速度未知,同时出发,同时达到.问跑的过程中,两狗的最大距离和最小距离的差 分析:训练指南P261,考虑相对运动,设A静止不动,B相对A运动,相对的运动向量:Vb - Va(可以理解为速度矢量),那么就是pa到线段pb-pb+Vb-Va的距离最值 /************************************************ * Author :Running_Time * Created Time :2015/10/22 星期四 10:21

【UVA】11796 - Dog Distance(相对运动)

因为a,b是同时出发,同时达到,但是他们的速度不一定一样,所以我们可以设他们的速度为La,Lb(La,Lb为a,b狗的总路程) 那么,如何a,b都是沿着直线运动的时候如何求他们之间的最短最长距离呢? 因为运动都是相对的,所以我们可以把a看成不懂的,而b相对于a的移动方向就是Vb - Va,因此就可以看成a和线段 b + (Vb -Va)之间的关系了 至于方向Va,Vb向量怎么求,我们可以利用单位向量 X 移动时间 X 移动速度 得到. 经典的一道题,顺便还修复了一个模板的BUG: 1412818

UVA11796 Dog Distance 计算几何

计算几何: 题解转自:点击打开链接 首先我们想一下如果甲乙都仅沿着两条线段向前跑,那么他们之间的最短和最长距离怎么算? 假设甲的速度向量为v1(速度向量指甲单位时间所走的位移向量),乙的速度向量为v2. 那么可以把甲看成是静止不动的,而让乙一个人以v2-v1的速度向量去运动(画图验证下). 且假设甲和乙都运动了t秒时间,那么我们可以知道上面的过程类似于甲不动,乙从自己的起点运动了t*(v2-v1)的位移向量. 而乙已知起点和位移向量,那么它就是在一条线段上. 最终我们要求的就是甲这个不动的点到乙

UVA 10140 - Prime Distance(数论)

10140 - Prime Distance 题目链接 题意:求[l,r]区间内最近和最远的素数对. 思路:素数打表,打到sqrt(Max)即可,然后利用大的表去筛素数,由于[l, r]最多100W,所以可以去遍历一遍,找出答案.注意1的情况,一开始没判断1,结果WA了 代码: #include <stdio.h> #include <string.h> #include <algorithm> using namespace std; #define INF 0x3f

UVA 11178-Morley&#39;s Theorem(计算几何_莫雷定理)

Problem D Morley's Theorem Input: Standard Input Output: Standard Output Morley's theorem states that that the lines trisecting the angles of an arbitrary plane triangle meet at the vertices of an equilateral triangle. For example in the figure below