目录

【几何】个人练习-Leetcode-1453. Maximum Number of Darts Inside of a Circular Dartboard

题目链接:https://leetcode.cn/problems/maximum-number-of-darts-inside-of-a-circular-dartboard/description/

题目大意:给出一系列点和一个圆的半径,(寻找一个圆心)求这个半径的圆最多能覆盖多少个点。

思路:几何上,如果一个圆能够覆盖N个点,那么在这N个点中,一定存在两个点,使得这个圆移动一下使得这两个点在圆上后,依然能够覆盖这原来N个点(详细的证明看网站上的题解,感觉还是比较intuitive的)。因此只需要遍历点对,寻找过这两个点,半径为r的圆的圆心,再计算这个圆覆盖的点数,求最大即可。注意两个点一个半径并无法确定圆心,因为这个圆心可能有两个,对称的,在纸上画画就能看出来。

然而代码写起来是有点繁杂,好多地方忘了用浮点数,debug了挺久。并且在判点是否在圆内圆外的函数中,我本地IDE上只需要>=0就行了,但这样子在leetcode网站上总有case过不了,跑出来答案不一样。于是修改了一下boundary,才通过。

完整代码

cpp 复制代码
class Solution {
public:
    inline int calD(vector<vector<int>>& darts, int r2, double cx, double cy) {
        int num = 0;
        for (auto d : darts) {
            double dis2 = (d[0] - cx) * (d[0] - cx) + (d[1] - cy) * (d[1] - cy);
            if (r2 - dis2 >= -1e-5)
                num++; 
        }
        return num;
    }


    int numPoints(vector<vector<int>>& darts, int r) {
        int n = darts.size();
        int ans = 1;
        int r2 = r*r;
        for (int i = 0; i < n; i++) {
            for (int j = i+1; j < n; j++) {
                double midx = 1.0*(darts[i][0] + darts[j][0]) / 2;
                double midy = 1.0*(darts[i][1] + darts[j][1]) / 2;
                double half = sqrt((darts[i][0] - darts[j][0]) * (darts[i][0] - darts[j][0]) + (darts[i][1] - darts[j][1]) * (darts[i][1] - darts[j][1]))/2;
                double p = sqrt(r*r - half*half);

                if (darts[i][0] == darts[j][0]) {
                    ans = max(ans, calD(darts, r2, darts[i][0] + p, midy));
                    ans = max(ans, calD(darts, r2, darts[i][0] - p, midy));
                }
                else if (darts[i][1] == darts[j][1]) {
                    ans = max(ans, calD(darts, r2, midx, darts[i][1] + p));
                    ans = max(ans, calD(darts, r2, midx, darts[i][1] - p));
                }
                else {
                    double k = 1.0*(darts[i][1] - darts[j][1]) / (darts[i][0] - darts[j][0]);
                    k = -1.0 / k;
                    ans = max(ans, calD(darts, r2, midx + p * 1 / sqrt(1 + k*k), midy + p * k / sqrt(1 + k*k)));
                    ans = max(ans, calD(darts, r2, midx - p * 1 / sqrt(1 + k*k), midy - p * k / sqrt(1 + k*k)));
                }
            }
        }
        return ans;
    }
};
本文是转载文章,点击查看原文
如有侵权,请联系 xyy@jishuzhan.net 删除
相关推荐
地平线开发者20 分钟前
【征程 6】工具链 VP 示例中 Cmakelists 解读
算法·自动驾驶
邪神与厨二病22 分钟前
2025蓝桥杯python A组题解
数据结构·c++·python·算法·蓝桥杯·单调栈·反悔贪心
明月看潮生1 小时前
青少年编程与数学 02-016 Python数据结构与算法 13课题、回溯
数据结构·python·算法·青少年编程
炫友呀1 小时前
python求π近似值
python·学习·算法
拓端研究室TRL1 小时前
Python与R语言用XGBOOST、NLTK、LASSO、决策树、聚类分析电商平台评论信息数据集
开发语言·python·算法·决策树·r语言
小伯宝宝么么2 小时前
关于编译原理——递归下降分析器的设计
java·开发语言·算法
梭七y3 小时前
【力扣hot100题】(089)最长有效括号
算法·leetcode·职场和发展
珂朵莉MM4 小时前
2025第四届大学生算法挑战赛 赛前测试赛 题解
算法
LabVIEW开发4 小时前
LabVIEW 图像处理中常见的边缘检测算法
图像处理·算法·labview
nuo5342025 小时前
第十六届蓝桥杯 省赛C/C++ 大学B组
c语言·c++·算法·蓝桥杯·图论