pub fn bhattacharyya_coefficient(p: &[f64], q: &[f64]) -> f64
Bhattacharyya coefficient: BC(p, q) = Σ √(p_i × q_i).
BC = 1 for identical distributions, BC = 0 for disjoint support. Related to Fisher-Rao: d_FR = 2 × arccos(BC).