如何用海伦公式快速判断点在直线的哪一侧

一、海伦公式的定义与推导

1. 海伦公式的定义

海伦公式（Heron’s Formula）是用于计算三角形面积的一种方法，适用于已知三角形三边长度的情况。公式如下：

$$S=\sqrt{s(s-a)(s-b)(s-c)}$$

其中：

• $ a $、$ b $、$ c $是三角形的三边长度。
• $ s $ 是三角形的半周长，计算公式为：
  $$s=\frac{a+b+c}{2}$$

2. 海伦公式的推导

海伦公式的推导可以通过三角形的面积公式和余弦定理来完成。以下是简要的推导过程：

1. 三角形面积公式：
   $$S=\frac{1}{2}ab\sin C$$
   其中 $ C$ 是边 $ a $ 和边 $ b $ 之间的夹角。

2. 余弦定理：
   $$c^2=a^2+b^2-2ab\cos C$$

3. 结合上述公式：
   通过代数运算和三角恒等式，可以将面积公式转换为仅用三边 $ a $、$ b$、$ c $ 表示的形式，最终得到海伦公式。

二、用海伦公式快速判断点在直线的哪一侧的原理

1. 行列式公式

行列式公式可以快速判断点$ p $ 是否在由两点 $ a $ 和 $ b $构成的直线的某一侧。行列式公式如下：

$$D=(x_b-x_a)(y_p-y_a)-(y_b-y_a)(x_p-x_a)$$

• 如果 $ D > 0 $，点 $p $ 在直线 $ ab $ 的左侧。
• 如果 $ D < 0 $，点$ p $ 在直线$ ab $ 的右侧。
• 如果 $ D \approx 0 $，点$ p $ 在直线 $ ab $ 上。

2. 海伦公式的应用

如果行列式公式的结果接近零（即点 $ p $ 可能在直线上），可以用海伦公式进一步验证。通过计算三角形 ( abp ) 的面积：

• 如果面积为零，则点 $p $ 在直线$ ab $ 上。
• 如果面积不为零，则点 $ p $ 不在直线上。
  

三、代码实现

以下是用 C++ 实现的代码，结合行列式公式和海伦公式来判断点在直线的哪一侧：

#include <iostream>
#include <cmath>
#include <vector>

using namespace std;

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

// 行列式公式判断点 p 在直线 ab 的哪一侧
int determinantSign(const Point& a, const Point& b, const Point& p) {
    double d = (b.x - a.x) * (p.y - a.y) - (b.y - a.y) * (p.x - a.x);
    if (d > 1e-9) return 1;  // 左侧
    if (d < -1e-9) return -1; // 右侧
    return 0;                 // 在直线上
}

// 海伦公式计算三角形面积
double heronArea(const Point& a, const Point& b, const Point& p) {
    double ab = sqrt(pow(b.x - a.x, 2) + pow(b.y - a.y, 2));
    double ap = sqrt(pow(p.x - a.x, 2) + pow(p.y - a.y, 2));
    double bp = sqrt(pow(p.x - b.x, 2) + pow(p.y - b.y, 2));
    double s = (ab + ap + bp) / 2.0;
    return sqrt(s * (s - ab) * (s - ap) * (s - bp));
}

int main() {
    Point a(0, 0), b(4, 0);
    vector<Point> points = {
        {2, 3},  // 在左侧
        {2, -3}, // 在右侧
        {2, 0},  // 在直线上
        {0, 0},  // 在直线上
        {4, 0}   // 在直线上
    };

    for (const auto& p : points) {
        int side = determinantSign(a, b, p);
        if (side == 0) {
            // 如果行列式结果接近零，进一步用海伦公式验证
            double area = heronArea(a, b, p);
            if (area < 1e-9) {
                cout << "Point (" << p.x << ", " << p.y << ") is on the line AB." << endl;
            }
        } else if (side == 1) {
            cout << "Point (" << p.x << ", " << p.y << ") is on the left side of line AB." << endl;
        } else {
            cout << "Point (" << p.x << ", " << p.y << ") is on the right side of line AB." << endl;
        }
    }

    return 0;
}

代码说明

1. 行列式公式：

   • 通过行列式公式快速判断点 ( p ) 的位置关系。
   • 如果行列式结果大于 ( 1e-9 )，点在左侧；小于 ( -1e-9 )，点在右侧；否则认为点在直线上。

2. 海伦公式：

   • 如果行列式结果接近零（即点可能在直线上），进一步用海伦公式计算三角形面积。
   • 如果面积小于 ( 1e-9 )，则确认点在直线上。

3. 输出结果：

   • 根据行列式公式和海伦公式的结合结果，输出点的位置关系。

四、总结

通过行列式公式和海伦公式的结合，可以快速且准确地判断点在直线的哪一侧。行列式公式提供了快速判断，而海伦公式用于进一步验证点是否在直线上。这种方法既高效又可靠，适用于多种场景下的几何计算。希望本文能帮助你更好地理解和应用海伦公式。