Say $X_1,X_2,\ldots$ are independent identically distributed Bernoulli random variables with mean $p$. This paper builds a new estimate $\hat p$ of $p$ that has the property that the relative error, $\hat p /p - 1$, of the estimate does not depend in any way on the value of $p$. This allows the construction of exact confidence intervals for $p$ of any desired level without needing any sort of limit or approximation. In addition, $\hat p$ is unbiased. For $ε$ and $δ$ in $(0,1)$, to obtain an estimate where $\mathbb{P}(|\hat p/p - 1| > ε) \leq δ$, the new algorithm takes on average at most $2ε^{-2} p^{-1}\ln(2δ^{-1})(1 - (14/3) ε)^{-1}$ samples. It is also shown that any such algorithm that applies whenever $p \leq 1/2$ requires at least $0.2ε^{-2} p^{-1}\ln((2-δ)δ^{-1})(1 + 2 ε)$ samples. The same algorithm can also be applied to estimate the mean of any random variable that falls in $[0,1]$.