AK F.*ing leetcode 流浪计划之半平面求交

已于 2023-11-14 19:04:23 修改
阅读量214 收藏
点赞数
分类专栏： 高阶算法 数学 图形学 文章标签： 半平面求交 凸包 极角排序
于 2023-11-14 17:21:38 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/chenbb1989/article/details/134274442
版权
图形学 同时被 3 个专栏收录
61 篇文章 39 订阅
订阅专栏
高阶算法
54 篇文章 7 订阅
订阅专栏
数学
38 篇文章 0 订阅
订阅专栏

欢迎关注更多精彩
关注我，学习常用算法与数据结构，一题多解，降维打击。

本期话题：半平面求交

背景知识

学习资料

视频讲解
https://www.bilibili.com/video/BV1jL411C7Ct/?spm_id_from=333.1007.top_right_bar_window_history.content.click&vd_source=fb27f95f25902a2cc94d4d8e49f5f777

文本资料
https://oi-wiki.org//geometry/half-plane/

基本问题转化

在很多题目中，给定的线段是没有方向的。此时，我们需要先把所有的线段都转化成点加向量的方式。使得向量的左边为有效区域。这样就可以使用模板求解了。

要注意的问题

  1. 主要的问题是浮点型的判断大小问题。在排序和判断点与线的关系时都用到浮点型判断。有些题型会卡精度，能用整数判断尽量不要使用浮点判断。
  2. atan2计算比较耗时，可以事先保存。

代码模板

求多边形的核

题目链接：https://vjudge.net/problem/UVA-1571

多边形的核就是取核区域内任意一点，站在该可以观察到多边形内任意一点。

利用半平面求交可以得到多边形的核

#include<stdio.h>
#include<cmath>
#include <algorithm>
#include <vector>
#include <list>
#include <cstring>
#include <set>


using namespace std;
const double EPS = 1e-14;

const int N = 2e6 + 10;

int cmp(double d) {
	if (abs(d) < EPS)return 0;
	if (d > 0)return 1;
	return -1;
}

class Point {
public:
	double x, y;
	int id;

	Point() {}
	Point(double a, double b) :x(a), y(b) {}
	Point(const Point& p) :x(p.x), y(p.y), id(p.id) {}

	void in() {
		scanf("%lf %lf", &x, &y);
	}
	void out() {
		printf("%.16f %.16f\n", x, y);
	}

	double dis() {
		return sqrt(x * x + y * y);
	}

	double dis2() {
		return x * x + y * y;
	}

	Point operator -() const {
		return Point(-x, -y);
	}

	Point operator -(const Point& p) const {
		return Point(x - p.x, y - p.y);
	}

	Point operator +(const Point& p) const {
		return Point(x + p.x, y + p.y);
	}
	Point operator *(double d)const {
		return Point(x * d, y * d);
	}

	Point operator /(double d)const {
		return Point(x / d, y / d);
	}


	void operator -=(Point& p) {
		x -= p.x;
		y -= p.y;
	}

	void operator +=(Point& p) {
		x += p.x;
		y += p.y;
	}
	void operator *=(double d) {
		x *= d;
		y *= d;
	}

	void operator /=(double d) {
		this ->operator*= (1 / d);
	}

	bool operator<(const Point& a) const {
		return x < a.x || (abs(x - a.x) < EPS && y < a.y);
	}

	bool operator==(const Point& a) const {
		return abs(x - a.x) < EPS && abs(y - a.y) < EPS;
	}
};

// 向量操作

double cross(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;
}

class Line {
public:
	Point front, tail;
	double ang;
	int u, v;
	Line() {}
	Line(const Point& a, const Point& b) :front(a), tail(b) {
		ang = atan2(front.y - tail.y, front.x - tail.x);
	}
};

int cmp(const Line& a, const Line& b) {
	//if (a.u == b.u && a.v == b.v)return 0;
	return cmp(a.ang - b.ang);

}


// 点在直线哪一边>0 左边，<0边
double SideJudge(const Line& a, const Point& b) {
	//return cmp(cross(a.front - a.tail, b - a.tail));
	return cross(a.front - a.tail, b - a.tail);
}


int LineSort(const Line& a, const Line& b) {
	int c = cmp(a, b);
	if (c)return c < 0;
	return	cross(b.front - b.tail, a.front - b.tail) > 0;
}

/*
点p 到 p+r 表示线段1
点q 到 q+s 表示线段2
线段1 上1点用 p' = p+t*r (0<=t<=1)
线段2 上1点用 q' = q+u*s (0<=u<=1)
让两式相等求交点 p+t*r = q+u*s
两边都叉乘s
(p+t*r)Xs = (q+u*s)Xs
pXs + t*rXs = qXs
t = (q-p)Xs/(rXs)
同理，
u = (p-q)Xr/(sXr) -> u = (q-p)Xr/(rXs)

以下分4种情况：
1. 共线，sXr==0 && (q-p)Xr==0, 计算 (q-p)在r上的投影在r长度上的占比t0，
计算(q+s-p)在r上的投影在r长度上的占比t1，查看[t0, t1]是否与范围[0，1]有交集。
如果t0>t1, 则比较[t1, t0]是否与范围[0，1]有交集。
t0 = (q-p)*r/(r*r)
t1 = (q+s-p)*r/(r*r) = t0 + s · r / (r · r)
2. 平行sXr==0 && (q-p)Xr!=0
3. 0<=u<=1 && 0<=t<=1 有交点
4. 其他u, t不在0到范围内，没有交点。
*/
pair<double, double> intersection(const Point& q, const Point& s, const Point& p, const Point& r, bool &oneline) {
	// 计算 (q-p)Xr
	auto qpr = cross(q - p, r);
	auto qps = cross(q - p, s);

	auto rXs = cross(r, s);
	if (cmp(rXs) == 0) {
		oneline = true;
		return { -1, -1 }; // 平行或共线
	}
	// 求解t, u
	// t = (q-p)Xs/(rXs)
	auto t = qps / rXs;

	// u = (q-p)Xr/(rXs)
	auto u = qpr / rXs;

	return { u, t };
}

Point LineCross(const Line& a, const Line& b, bool &f) {
	Point dira = a.front - a.tail;
	Point dirb = b.front - b.tail;
	bool oneline=false;
	auto p = intersection(a.tail, dira, b.tail, dirb, oneline);
	if (oneline)f = false;
	return a.tail + dira * p.first;
}

class HalfPlane {
public:
	vector<Line> lines;
	vector<int> q;
	vector<Point> t;
	int len;

	HalfPlane() {
		lines.resize(N);
		q.resize(N);
		t.resize(N);
	}

	void reset() {
		len = 0;
	}

	void addLine(const Line& a) {
		lines[len++] = a;
	}

	bool run() {
		sort(lines.begin(), lines.begin() + len, LineSort);

		int l = -1, r = 0;
		q[0] = 0;
		for (int i = 1; i < len; ++i) {
			if (cmp(lines[i], lines[i - 1]) == 0)continue;
			while (r - l > 1 && SideJudge(lines[i], t[r]) < 0)r--;
			while (r - l > 1 && SideJudge(lines[i], t[l + 2]) < 0)l++;
			q[++r] = i;
			bool f=true;
			t[r] = LineCross(lines[q[r]], lines[q[r - 1]], f);
		}
		while (r - l > 1 && SideJudge(lines[q[l + 1]], t[r]) < 0)r--;

		//if (r - l > 1) {
		//	bool f = true;
		//	t[r + 1] = LineCross(lines[q[l + 1]], lines[q[r]], f);
		//	r++;
		//	if (!f)r -= 2;
		//}

		 统计交点
		//l++;
		//vector<Point> ans(r - l);
		//for (int i = 0; i < ans.size(); ++i) {
		//	ans[i] = t[i + l + 1];
		//}


		return r-l>2;
	}
};

Point oiPs[N * 2];
pair<int, int> ori[N * 2];
HalfPlane hp;

int bigDevid(int a, int b) {
	for (int i = max(abs(a), abs(b)); i >= 1; i--) {
		if (a % i == 0 && b % i == 0)return i;
	}
	return 1;
}


void  solve() {
	int n, m = 1;
	//FILE* fp = fopen("ans.txt", "w");
	while (scanf("%d", &n) != EOF && n) {
		int a, b;
		for (int i = 0; i < n; ++i) {
			scanf("%d %d", &a, &b);
			oiPs[i] = Point(a, b);
			ori[i] = { a,b };
		}
		oiPs[n] = oiPs[0];
		ori[n] = ori[0];

		hp.reset();
		for (int i = 0; i < n; ++i) {
			hp.addLine(Line(oiPs[i+1], oiPs[i]));
			hp.lines[i].u = ori[i+1].first - ori[i].first;
			hp.lines[i].v = ori[i+1].second - ori[i].second;

			int bd = bigDevid(hp.lines[i].u, hp.lines[i].v);
			hp.lines[i].u /= bd;
			hp.lines[i].v /= bd;
		}

		auto ps = hp.run();
		if (ps)puts("1");
		else puts("0");
		m++;
	}
}


int main() {
	solve();
	return 0;

}

/*
4
0 0
0 1
1 1
1 0
8
0 0
3 0
4 3
2 2
3 4
4 4
4 5
0 5
0


8
0 0
0 1
1 1
1 2
0 2
0 3
3 3
3 0
*/

练习一

链接：https://www.luogu.com.cn/problem/P4196

求多个凸多边形的交面积。

对每条边进行半平面求交，再利用三角形求多边形面积。


#include<stdio.h>
#include<cmath>
#include <algorithm>
#include <vector>
#include <list>
#include <cstring>
#include <set>


using namespace std;
const double EPS = 1e-14;

const int N = 2e6 + 10;

int cmp(double d) {
	if (abs(d) < EPS)return 0;
	if (d > 0)return 1;
	return -1;
}

class Point {
public:
	double x, y;
	int id;

	Point() {}
	Point(double a, double b) :x(a), y(b) {}
	Point(const Point& p) :x(p.x), y(p.y), id(p.id) {}

	void in() {
		scanf("%lf %lf", &x, &y);
	}
	void out() {
		printf("%.16f %.16f\n", x, y);
	}

	double dis() {
		return sqrt(x * x + y * y);
	}

	double dis2() {
		return x * x + y * y;
	}

	Point operator -() const {
		return Point(-x, -y);
	}

	Point operator -(const Point& p) const {
		return Point(x - p.x, y - p.y);
	}

	Point operator +(const Point& p) const {
		return Point(x + p.x, y + p.y);
	}
	Point operator *(double d)const {
		return Point(x * d, y * d);
	}

	Point operator /(double d)const {
		return Point(x / d, y / d);
	}


	void operator -=(Point& p) {
		x -= p.x;
		y -= p.y;
	}

	void operator +=(Point& p) {
		x += p.x;
		y += p.y;
	}
	void operator *=(double d) {
		x *= d;
		y *= d;
	}

	void operator /=(double d) {
		this ->operator*= (1 / d);
	}

	bool operator<(const Point& a) const {
		return x < a.x || (abs(x - a.x) < EPS && y < a.y);
	}

	bool operator==(const Point& a) const {
		return abs(x - a.x) < EPS && abs(y - a.y) < EPS;
	}
};

// 向量操作

double cross(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;
}

class Line {
public:
	Point front, tail;
	double ang;
	int u, v;
	Line() {}
	Line(const Point& a, const Point& b) :front(a), tail(b) {
		ang = atan2(front.y - tail.y, front.x - tail.x);
	}
};

int cmp(const Line& a, const Line& b) {
	//if (a.u == b.u && a.v == b.v)return 0;
	return cmp(a.ang - b.ang);

}


// 点在直线哪一边>0 左边，<0边
double SideJudge(const Line& a, const Point& b) {
	//return cmp(cross(a.front - a.tail, b - a.tail));
	return cross(a.front - a.tail, b - a.tail);
}


int LineSort(const Line& a, const Line& b) {
	int c = cmp(a, b);
	if (c)return c < 0;
	return	cross(b.front - b.tail, a.front - b.tail) > 0;
}

/*
点p 到 p+r 表示线段1
点q 到 q+s 表示线段2
线段1 上1点用 p' = p+t*r (0<=t<=1)
线段2 上1点用 q' = q+u*s (0<=u<=1)
让两式相等求交点 p+t*r = q+u*s
两边都叉乘s
(p+t*r)Xs = (q+u*s)Xs
pXs + t*rXs = qXs
t = (q-p)Xs/(rXs)
同理，
u = (p-q)Xr/(sXr) -> u = (q-p)Xr/(rXs)

以下分4种情况：
1. 共线，sXr==0 && (q-p)Xr==0, 计算 (q-p)在r上的投影在r长度上的占比t0，
计算(q+s-p)在r上的投影在r长度上的占比t1，查看[t0, t1]是否与范围[0，1]有交集。
如果t0>t1, 则比较[t1, t0]是否与范围[0，1]有交集。
t0 = (q-p)*r/(r*r)
t1 = (q+s-p)*r/(r*r) = t0 + s · r / (r · r)
2. 平行sXr==0 && (q-p)Xr!=0
3. 0<=u<=1 && 0<=t<=1 有交点
4. 其他u, t不在0到范围内，没有交点。
*/
pair<double, double> intersection(const Point& q, const Point& s, const Point& p, const Point& r, bool &oneline) {
	// 计算 (q-p)Xr
	auto qpr = cross(q - p, r);
	auto qps = cross(q - p, s);

	auto rXs = cross(r, s);
	if (cmp(rXs) == 0) {
		oneline = true;
		return { -1, -1 }; // 平行或共线
	}
	// 求解t, u
	// t = (q-p)Xs/(rXs)
	auto t = qps / rXs;

	// u = (q-p)Xr/(rXs)
	auto u = qpr / rXs;

	return { u, t };
}

Point LineCross(const Line& a, const Line& b, bool &f) {
	Point dira = a.front - a.tail;
	Point dirb = b.front - b.tail;
	bool oneline=false;
	auto p = intersection(a.tail, dira, b.tail, dirb, oneline);
	if (oneline)f = false;
	return a.tail + dira * p.first;
}


class HalfPlane {
public:
	vector<Line> lines;

	void addLine(const Line& a) {
		lines.push_back(a);
	}

	vector<Point> run() {
		sort(lines.begin(), lines.end(), LineSort);
		vector<int> q(lines.size() + 10);
		vector<Point> t(lines.size() + 10);

		int l = -1, r = 0;
		q[0] = 0;
		for (int i = 1; i < lines.size(); ++i) {
			if (cmp(lines[i], lines[i - 1]) == 0)continue;
			while (r - l > 1 && SideJudge(lines[i], t[r]) < 0)r--;
			while (r - l > 1 && SideJudge(lines[i], t[l + 2]) < 0)l++;
			q[++r] = i;
			bool f = true;
			t[r] = LineCross(lines[q[r]], lines[q[r - 1]], f);
		}
		while (r - l > 1 && SideJudge(lines[q[l + 1]], t[r]) < 0)r--;
		if (r - l > 1) {
			bool f = true;
			t[r + 1] = LineCross(lines[q[l + 1]], lines[q[r]], f);
			r++;
		}

		// 统计交点
		l++;
		vector<Point> ans(r - l);
		for (int i = 0; i < ans.size(); ++i) {
			ans[i] = t[i + l + 1];
		}

		return ans;
	}
};

Point oiPs[N];

void  solve() {
	int n, m;
	scanf("%d", &n);
	HalfPlane hp;
	int a, b;
	while (n--) {
		scanf("%d", &m);
		for (int i = 0; i < m; ++i) {
			scanf("%d%d", &a, &b);
			oiPs[i].x = a;
			oiPs[i].y = b;
		}

		oiPs[m] = oiPs[0];
		for (int i = 0; i < m; ++i) {
			hp.addLine(Line(oiPs[i + 1], oiPs[i]));
		}
	}

	auto keyPoints = hp.run();

	double ans = 0;
	for (int i = 2; i < keyPoints.size(); ++i) {
		ans += cross(keyPoints[i - 1] - keyPoints[0], keyPoints[i] - keyPoints[0]);
	}

	printf("%.3f\n", ans / 2);
}


int main() {
	solve();
	return 0;

}

/*
3
3
-1 2
-2 1
-1 1

3
1 1
2 1
1 2

3
1 1
3 0
2 2

*/

本人码农，希望通过自己的分享，让大家更容易学懂计算机知识。创作不易，帮忙点击公众号的链接。

闪电彬彬
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
vscode安装leetcode-LeetCode.Practice:LeetCode联系题
07-07
【标题】：在VSCode中安装与使用LeetCode.Practice插件进行LeetCode练习 【描述】：Visual Studio Code（VSCode）是一款流行的、免费且开源的代码编辑器，广泛用于编写各种编程语言的代码。为了帮助开发者更方便地...
LeetCode 1. 两数之和 Two Sum
Umbrella Corporation
08-29 2万+
给定一个整数数组和一个目标值，找出数组中和为目标值的两个数。 你可以假设每个输入只对应一种答案，且同样的元素不能被重复利用。 示例: 给定 nums = [2, 7, 11, 15], target = 9 因为 nums[0] + nums[1] = 2 + 7 = 9 所以返回 [0, 1]   解决方案 方法一：暴力法 暴力法很简单。遍历每个元素 xxx，并查找是否...
【leetcode】001.两数之和（C语言/C++，超详细）
慕雪华年的博客
10-31 1万+
文章目录两数之和 leetcode-001题目要求实现步骤1.模板中四个值的意义2.在数组中找到两个连续的元素3.判断相加是否等于target4.返回元素下标代码示例1static的作用*returnSize是什么？代码示例2什么是malloc函数？malloc使用后free释放总结 两数之和 leetcode-001 题目来源leetcode 如下图所示 右侧给出了题目的基本模板 /** * Note: The returned array must be malloced, assume calle
Leetcode算法Java全解答--18. 四数之和
cmqwan的博客
11-20 1万+
Leetcode算法Java全解答–18. 四数之和 文章目录Leetcode算法Java全解答--18. 四数之和题目想法结果总结代码我的答案大佬们的答案测试用例其他 题目 给定一个包含 n 个整数的数组 nums 和一个目标值 target，判断 nums 中是否存在四个元素 a，b，c 和 d ，使得 a + b + c + d 的值与 target 相等？找出所有满足条件且不重复的四元组。...
IDEA上班摸鱼神器之LeetCode刷题插件
Randy_Wang_的博客
12-21 4902
上班无聊的时候想刷题，网页操作怕被老板发现？刷题的时候先在ide敲代码再复制到网页，切换太频繁？你需要一款神器插件——LeetCode Editor。 插件github地址： https://github.com/shuzijun/leetcode-editor 安装 与安装其他插件步骤一样，打开Plugins界面，搜索关键字，搜索结果第一个就是我们需要安装的插件，点击Install即可安装（这里我已经安装过了）。 ...
Error:no matching function for call to......&leetcode2
qq_40831429的博客
04-11 6181
刷leetcode，写成如下的时候这样报错 struct ListNode *result = new ListNode(); 一开始连结构体都没有初始化，然后就改成了这样。太久不用C++忘了很多。 结构体指针在使用之前必须通过初始化或赋值操作，把某个结构体变量的首地址赋给它，使它指向该结构体变量。 现在的报错内容为 error: no matching function for c...
LeetCode解题汇总目录
热门推荐
Michael是个半路程序员
09-06 7万+
此篇为学习完《数据结构与算法之美》后，在LeetCode刷题的汇总目录，方便大家查找（Ctrl+Find），一起刷题，一起PK交流！如果本文对你有帮助，可以给我点赞加油！ Updated on 2022.6.12刷题可以按照不同的专题进行，便于加强某个知识点的理解。另有《剑指Offer》面试题（更新完毕）、《程序员面试金典》面试题（更新完毕）、LintCode代码能力测试CAT。LeetCode MySQL 数据库解题目录点击访问 我的 LeetCode 主页我的 GitHub 主页注：本文超过字数限制，
leetcode 224. Basic Calculator | 224. 基本计算器（中缀表达式求值）
寒泉
09-11 4万+
题目 https://leetcode.com/problems/basic-calculator/ 题解 class Solution { public int calculate(String s) { if (s.trim().startsWith("-")) s = "0" + s; Stack<Integer> numStack = new Stack<>(); Stack<Character
VScode配置Leetcode环境
Sxueyuan的博客
08-19 4184
在本地VScode进行Leetcode刷题的配置方法
[环境配置] vscode 配置 leetcode 插件
Chris_zhangrx的博客
12-28 3246
文章目录配置步骤可能会遇到的问题 配置步骤 Extension 商店搜索： leetcode，选择第一个下载量最高的下载就可以了 安装完成后左侧就会多一个 leetcode 的图标栏 点击上面的图标，会有 Sign in to LeetCode 字样 上图中： 红框标记的是选择中文还是 us 版的 leetcode 网站。 黄框标记的是登录方式。 若依赖没有问题直接选择相应的网站，登录自己的帐号就可以了。 可能会遇到的问题 这里需要的依赖是： VS Code 1.30.1+ Node.
LeetCode 15. 三数之和
y4ung
08-17 5411
题目描述 15. 三数之和 思路 思路1 比较容易想到的就是，求三数之和等于0，可以等价于求两个数的和，然后看这个和的相反数是否在nums里面。 但是 T_T这样的话复杂度太高了,会超时，捂脸，最后三个case，怎么改都超时…… bool find(vector&lt;int&gt; nums, int nums_sz, int target, int m, i...
leetcode-java.zip_leetcode
09-24
《LeetCode Java版：算法刷题与学习指南》 LeetCode是一个广受欢迎的在线平台，致力于帮助程序员提升技能，特别是算法和数据结构方面的能力。它提供了丰富的编程挑战，支持多种编程语言，包括Java。本资源...
LeetCode.zip_leetcode
07-15
《LeetCode题解：全面提升编程能力的黄金宝典》 LeetCode，作为全球知名的在线编程挑战平台，汇聚了各类算法题目，旨在提升程序员的编程思维和解决实际问题的能力。这个名为"LeetCode.zip_leetcode"的压缩包，显然...
LeetCode_No80.rar_leetcode_leetcode c++代码
09-24
LeetCode 是一个非常知名的在线编程挑战平台，专为程序员提升技能和准备面试而设计。它提供了各种难度级别的问题，涵盖多种编程语言，包括C++。在这个名为 "LeetCode_No80.rar" 的压缩包文件中，我们可以推测它包含...
leetcode详解.zip_leetcode_leetcode详解
07-15
程序员招工组必备，leetcode刷题中文详解，迅速提升代码能力
udp分片报文pcap文件
07-18
udp分片报文pcap文件：64kudp大包分成每片658字节总计100片
计算机教学导案.doc
07-18
计算机
摩根：2024年大市前瞻.pdf
最新发布
07-18
摩根：2024年大市前瞻
计算机控制专业技术专题实习任务书温度控制.doc
07-18
计算机
regexp '^[a-zA-Z]+[a-zA-Z0-9_\\.\\/\\-]*@leetcode\\.com$'什么意思
07-31
这个正则表达式的含义是匹配以字母开头，后面可以是字母（大写或小写）、数字、下划线 '_'、句点 '.' 和/或横杠 '-' 的字符串，最后以'@leetcode.com'结尾。换句话说，这个正则表达式用于匹配符合特定格式的电子邮件...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值