DH: x <- [2..p-2], not [0..p-1]
DH: x <- [2..p-2], not [0..p-1]

g^0 = g^(p-1) = 1 mod p
g^1 = g^p     = g mod p

nasrallahmounir/ocaml-nocrypto,mirleft/ocaml-nocrypto,mirleft/ocaml-nocrypto,nasrallahmounir/ocaml-nocrypto