Gracias a todos. Resulta que $p=e$ es exactamente correcto, que todo depende de usar el logaritmo, y luego transformar una suma infinita en una integral. Mis pasos fueron más o menos los mismos:
Sea $\tilde{N}$ denotan el producto después de $j$ iteraciones. Entonces, tenemos $$\frac{\tilde{N}}{N_0} = \left(p \cdot r_1\right) \cdot \left(p \cdot r_2\right) \cdot \left(p \cdot r_3\right) \cdot \cdots \:,$$ donde $p$ es constante, y cada $r_j$ es un número aleatorio comprendido entre $0$ y $1$ . A continuación, toma el logaritmo natural de ambos lados, convirtiendo el producto de la derecha en suma: $$\ln\left(\frac{\tilde{N}}{N_0}\right) = \ln\left(p \cdot r_1\right) + \ln\left(p \cdot r_2\right) + \ln\left(p \cdot r_3\right) + \cdots$$
Para satisfacer la pregunta tal como está planteada, tenemos que el lado izquierdo (en promedio) se resuelve en cero, y el número de términos de la derecha se aproxima al infinito. Esto hace que $r$ una variable continua: después de muchas iteraciones, hemos muestreado todos los números reales entre cero y uno. Es decir, podemos sustituir la suma por una integral: $$0 \approx \ln\left(\frac{\tilde{N}}{N_0}\right) = \int_0^1 \left(\ln p + \ln r\right) dr \:,$$ y la respuesta cae directamente. Sólo $p=e$ cumple lo anterior.