fn table_index(p: i32, m: i32) -> usize
Position p * p + p + m of (p, m) in a harmonic table.
p * p + p + m
(p, m)