[CF772B Volatile Kite] 题解

最新推荐文章于 2024-09-29 16:01:26 发布

toothless03

最新推荐文章于 2024-09-29 16:01:26 发布

阅读量65

点赞数 1

分类专栏：题解文章标签：算法

本文链接：https://blog.csdn.net/2301_76305180/article/details/128749557

版权

题解专栏收录该内容

10 篇文章 0 订阅

订阅专栏

[CF772B Volatile Kite] 题解

先来双倍经验：P3744 李彬的几何，真的是一毛一样的

显然对于三个点 $A, B, C$ ，如果 $B$ 移动到 $A C$ 上时为最大值，如果 $A C$ 和 $B$ 相距 $d$ ，辣么 $B$ 可以向 $A C$ 移动 $\displaystyle\frac{d}{2}$ ， $A C$ 可以整体向 $B$ 移动 $\displaystyle\frac{d}{2}$ ，所以最终答案即为 $\displaystyle\frac{d}{2}$ 。我们显然可以看出来， $B$ 到其他对角线的距离小于 $B$ 到 $B$ 两侧点的连线。所以我们就可以暴力枚举连续的 $3$ 个点即可，复杂度 $O (n)$ 。

那如何计算点到直线的距离呢？

对于 $B (n, m)$ 到 $A C$ 的距离当然可以先将 $A C$ 的直线方程计算成形如 $a x + b y = c$ 的形式，距离就是 $\displaystyle\frac{|an+bm+c|}{\sqrt{a^2+b^2}}$ 但是实在是有亿点麻烦。。。。。。

所以我们可以使用向量叉积 $|\vec{BA}\times\vec{BC}|$ 除以 $∣ A C ∣$ 即可。

附上 AC 代码

#include <iostream>
#include <algorithm>
#include <cmath>
#include <queue>
#include <stack>
#include <climits>
#include <cstring>
#include <map>
#include <bitset>
using namespace std;
#define int long long
#define double long double
#define pr pair<int, int>
const int INF = LONG_LONG_MAX;
const int mod = 1000000007;
int gcd(int x, int y) {
	if (y == 0)
		return x;
	return gcd(y, x % y);
}
struct point {
	double x, y;
	friend void operator>>(istream& in, point& p) {
		in >> p.x >> p.y;
	}
	friend double operator*(point a, point b) {
		return abs(a.x * b.y - a.y * b.x);
	}
	friend double len(point p) {
		return sqrt(p.x * p.x + p.y * p.y);
	}
	friend point operator-(point a, point b) {
		return {a.x - b.x, a.y - b.y};
	}
};
signed main() {
	ios::sync_with_stdio(false);
	cin.tie(0);
	int n; cin >> n;
	vector<point> pts(n + 2);
	for (int i = 0; i < n; i++)
		cin >> pts[i];
	pts[n] = pts[0], pts[n + 1] = pts[1];
	double ans = INF;
	for (int i = 1; i <= n; i++)
		ans = min((pts[i - 1] - pts[i]) * (pts[i + 1] - pts[i]) / len(pts[i + 1] - pts[i - 1]), ans);
	ans /= 2.0;
	printf("%.8Lf\n", ans);
}