2021 ICPC 昆明 I Mr Main and Windmills(直线与线段的交点)

2021 ICPC 昆明 I Mr. Main and Windmills(直线与线段的交点)

I Mr. Main and Windmills

大意:给出一条线段 , 一个人从线段的起点走到线段的终点 , 线段的一侧有若干风车 , 当前的人在线段上的每一个位置观察风车都会得到一个顺序 。多次询问第 i 号风车被观察的位置第k次改变时人在线段上的位置。

思路:不难发现 , 两个风车交换位置当且仅当人走过 两风车所在直线与线段交点的时候 , 两两枚举风车求直线与线段交点 , 然后根据和起始点的距离排序后根据要求输出即可。

易错点:这里线段与直线求交会有一个易错点。

如果先求 线段所在直线与风车直线的交点(line_make_point) , 然后再判断交点是否在线段上(point_on_segment) , 这样误差会巨大。因为直线求交会有除法 , 求出的交点存在误差 ,然后判断点在线段上时会用到叉积 , 叉积的几何意义就是形成三角形的面积 , 如果线段特别特别长 , 叉积就会很大 , 从而在这里产生错误。

cpp 复制代码
if(!line_make_point(l , r , now)) continue;
if(!point_on_segment(now , st , ed)) continue;

解决方法:

1. 对于求交问题 , 先判断在求交

对应在这里 , 就可以先判断线段和直线是否相交(toleft) , 相交求交点即可 , 这样是不会有判断误差的产生的。

cpp 复制代码
if(toleft(st , p[i] , p[j]) * toleft(ed , p[i] , p[j]) > 0) continue;
line_make_point(l , r , now);

2. double 换成 long double , 容限(eps) 调大即可

这里推荐第一种 , 第一种更规范

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define IOS std::ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
#define int long long
const int N = 1e3 + 10;
const int mod = 1e9 + 7;
typedef pair<int,int>PII;

//--------------------------------------------------------------
const double eps = 1e-9;
const double pi = acos(-1);
inline double sqr(double x) {return x * x;} //平方
int sign(double x){
	if(fabs(x) < eps) return 0;
	if(x > 0) return 1;
	return -1;
}//符号
struct point{
	double x , y;
	point(){}
	point(double a , double b) : x(a) , y(b){}
	friend point operator + (const point &a , const point &b){
		return point(a.x + b.x , a.y + b.y);
	}
	friend point operator - (const point &a , const point &b){
		return point(a.x - b.x , a.y - b.y);
	}
	friend bool operator == (const point &a , const point &b){
		return !sign(a.x - b.x) && !sign(a.y - b.y);
	}
	friend point operator * (const point &a , const double &b){
		return point(a.x * b , a.y * b);
	}
	friend point operator * (const double &a , const point &b){
		return point(a * b.x , a * b.y);
	}
	friend point operator / (const point &a , const double &b){
		return point(a.x / b , a.y / b);
	}
	//向量模长 
	double norm(){ 
		return sqrt(sqr(x) + sqr(y));
	}
}; 

struct line{
	point a , b;
	line(){}
	line(point x , point y) : a(x) , b(y) {}
};

double det(const point &a , const point &b){
	return a.x * b.y - a.y * b.x;
}//叉积 判断两点共线 

double dot(const point &a , const point &b){
	return a.x * b.x + a.y * b.y;
}//点积

double dist(const point &a , const point &b){
	return (a - b).norm();
}//两点距离
 
point rotate_point(const point &a , const point &p , double A){
	double tx = p.x - a.x , ty = p.y - a.y;
	return point(a.x + tx * cos(A) - ty * sin(A) , a.y + tx * sin(A) + ty * cos(A));
}// p 点 绕 a 点逆时针旋转 A 弧度

int toleft(const point &p , const point &a , const point &b) {
	return sign(det(b - a , p - a));
	// 1 左 0 上 -1 右
}//只适用凸多边形

//判断点 p 是否在线段 st 上(包括端点)
bool point_on_segment(point p , point s , point t){
	return sign(det(p - s , t - s)) == 0 && sign(dot(p - s , p - t)) <= 0;
}

bool parallel(line a , line b){
	return !sign(det(a.a - a.b , b.a - b.b));
}

bool line_make_point(line a , line b , point &res){
	if(parallel(a , b)) return 0;
	double s1 = det(a.a - b.a , b.b - b.a);
	double s2 = det(a.b - b.a , b.b - b.a);
	res = (s1 * a.b - s2 * a.a) / (s1 - s2);
	return 1;
}
//--------------------------------------------------------------
//--------------------------------------------------------------

int n , m;
point st , ed , p[N] , now;
double x , y;
int h , k;

vector<tuple<double , double , double>>ans[N];

signed main(){

	IOS
	cout << fixed << setprecision(10);
	cin >> n >> m;
	
	cin >> x >> y;
	st = point{x , y};
	cin >> x >> y;
	ed = point{x , y};
	
	for(int i = 1 ; i <= n ; i ++){
		cin >> x >> y;
		p[i] = point{x , y};
	}
	
	line l = line{st , ed};
	
	for(int i = 1 ; i <= n ; i ++){
		for(int j = i + 1 ; j <= n ; j ++){
			line r = line{p[i] , p[j]};
			if(toleft(st , p[i] , p[j]) * toleft(ed , p[i] , p[j]) > 0) continue;
			line_make_point(l , r , now);
			ans[i].emplace_back(now.x , now.y , dist(now , st));
			ans[j].emplace_back(now.x , now.y , dist(now , st));
		} 
	}
	
	for(int i = 1 ; i <= n ; i ++) sort(ans[i].begin() , ans[i].end() , 
		[&](tuple<double , double , double> a , tuple<double , double , double> b){
			return get<2>(a) < get<2>(b);
		}
	);
	
	for(int i = 1 ; i <= m ; i ++){
		cin >> h >> k;
		if(ans[h].size() < k){
			cout << "-1\n";
		}else{
			auto [x , y , z] = ans[h][k - 1];
			cout << x << " " << y << "\n";
		}
	}
	
	return 0;
}
//freopen("文件名.in","r",stdin);
//freopen("文件名.out","w",stdout);
相关推荐
shymoy12 分钟前
Radix Sorts
数据结构·算法·排序算法
风影小子21 分钟前
注册登录学生管理系统小项目
算法
黑龙江亿林等保23 分钟前
深入探索哈尔滨二级等保下的负载均衡SLB及其核心算法
运维·算法·负载均衡
lucy1530275107926 分钟前
【青牛科技】GC5931:工业风扇驱动芯片的卓越替代者
人工智能·科技·单片机·嵌入式硬件·算法·机器学习
杜杜的man42 分钟前
【go从零单排】迭代器(Iterators)
开发语言·算法·golang
小沈熬夜秃头中୧⍤⃝1 小时前
【贪心算法】No.1---贪心算法(1)
算法·贪心算法
木向2 小时前
leetcode92:反转链表||
数据结构·c++·算法·leetcode·链表
阿阿越2 小时前
算法每日练 -- 双指针篇(持续更新中)
数据结构·c++·算法
skaiuijing2 小时前
Sparrow系列拓展篇:对调度层进行抽象并引入IPC机制信号量
c语言·算法·操作系统·调度算法·操作系统内核
Star Patrick2 小时前
算法训练(leetcode)二刷第十九天 | *39. 组合总和、*40. 组合总和 II、*131. 分割回文串
python·算法·leetcode