Conjugate Bayes for probit regression via unified skew-normal distributions