点到线段的最短矩离 及垂足的计算

过P做MN的垂线，垂足为Q，若Q在线段MN以内(包括与点M点N重合)，则最短距离为垂线段长度，若垂足在MN以外，则最短距离为PM，PN中的较小者。（若P与MN共线，垂线长度为零，同样适用）
所以可由 Q位置判断

一、代数法

若Q在线段MN上 则
$$|MQ|+|QN|=|MN|$$
或不开方的
$$|MQ{|}^2+|QN{|}^2\le |MN{|}^2$$

所以最小距离可表示为
$$\{\begin{array}{lcr}|PQ|& & \text{|}MQ|+|QN|=|MN|\\ & & \\ MIN(|PN|,|PM|)& & |MQ|+|QN|>|MN|\end{array}$$
或
$$\{\begin{array}{lcr}|PQ|& & \text{|}MQ{|}^2+|QN{|}^2\le |MN{|}^2\\ & & \\ MIN(|PN|,|PM|)& & |MQ{|}^2+|QN{|}^2>|MN{|}^2\end{array}$$

计算垂足
$$\angle b=\angle a+90°$$
$$\tan (b)=\tan (a+\frac{\pi }{2})=-\cot (a)$$
MN斜率为k PQ斜率为 -1/k

$$\{\begin{array}{lcr}\frac{y_q-y_p}{x_q-x_p}=-\frac{1}{k}& & \\ & & \\ y=kx+b& & \end{array}$$

点Q在MN上求解

$$\{\begin{array}{l}{x}_q=\frac{x_p+k{y}_p-kb}{1+k^2}\\ \\ y_q=\frac{k{x}_p+k^2{y}_p+b}{1+k^2}\end{array}$$

其中
$$\{\begin{array}{lcr}b=y_m-\frac{y_n-y_m}{x_n-x_m}{x}_m=\frac{x_n{y}_m-x_m{y}_n}{x_n-x_m}& & \\ & & \\ k=\frac{y_n-y_m}{x_n-x_m}& & \end{array}$$

代入 M N P 坐标 可算出 垂足坐标

水平与坚直另分一种情况直接求解

二、向量法

$$\overrightarrow{MN}\dot{}\overrightarrow{MP}=|MN||MP|cos\alpha$$
$\angle PMN>90°$ 点积结果为负
$\angle PMN=90°$ 点积结果为零
$\angle PMN<90°$ 点积结果为正

$$|MP|cos\alpha =\frac{\overrightarrow{MN}\dot{}\overrightarrow{MP}}{|MN|}$$
结果为MQ的长度 带符号 另比一次

$$\frac{|MP|cos\alpha }{|MN|}=\frac{\overrightarrow{MN}\dot{}\overrightarrow{MP}}{|MN{|}^2}$$

结果为MQ与MN长度比带符号 设为$\lambda$
若Q在线段MN上 则
$0\le \lambda \le 1$
其它情况Q在线段外
垂足坐标
$$\begin{gathered} x_q=\lambda (x_n-x_m)+x_m \\ {y}_q=\lambda (y_n-y_m)+y_m \end{gathered}$$

 typedef double real;
 typedef struct {
 	real x;
 	real y;
} PT;

 real dist(PT* a, PT* b, PT* c,PT* foot)
{
  real ab = length(a,b);
  if (ab < EPS){  
    return length(a,c);
  }   
  
  real lambda = dot(a,b,a,c)/(ab*ab);

  real dx = b->x - a->x;
  real dy = b->y - a->y;

  foot->x = lambda*dx + a->x;
  foot->y = lambda*dy + a->y;

  if( lambda >= 0 && lambda <= 1) {
    return  length(c,foot); // or cross(a,b,a,c)/ab;
  }else {
    real ac = length(a,c);
    real bc = length(b,c);
    return MIN(ac,bc);
  }
}