pub fn hellinger_distance(p: &[f64], q: &[f64]) -> f64
Hellinger distance: H(p, q) = √(1 - BC(p, q)) / √2.
Closely related to Fisher-Rao but maps to [0, 1]. Often used when a [0, 1] range is preferred over [0, π].