pub fn add(a: u64, b: u64) -> u64
Add a + b mod p. Result may be non-canonical (in [0, 2^64)).
[0, 2^64)