fn full_shr<const N: usize>(a: &ArrayPlusOne<u64, N>, shift: u32) -> [u64; N]
Narrowing right shift of an (N+1)-digit array by at most 63 bits