nonzeroElements p = [1..p]
inverse a p = head [ b | b <- nonzeroElements p, b*a `mod` p == 1]
inverses p = [ (a, inverse a p) | a <- nonzeroElements p ]