Me gustaría generar números aleatorios cuya distribución se imita a la de los números primos. Por lo que el número de genera números aleatorios menos de $n$ debe crecer como $n / \log n$, la mayoría de los intervalos de $[n,n+n^\epsilon]$ debe contener aproximadamente $n^\epsilon / \log n$ los números que se generan (Selberg de intervalos cortos de tiempo), etc. ¿Alguien sabe de un computacionalmente factible método para la generación de este tipo de "mirar-uno-como el de los números primos"?
Addendum. He implementado Henry y Xoff sugerencias.
Aquí hay varias instancias de los diez primeros "pseudoprimes":
$$ 4, 5, 9, 10, 17, 23, 27, 28, 31, 44 $$
$$ 7, 8, 9, 10, 12, 15, 18, 19, 27, 34 $$
$$ 6, 11, 15, 16, 23, 26, 27, 29, 45, 49 $$
Y aquí está la distribución acumulativa pseudoprimes a a $10^6$ (rojo), junto con
una parcela de $n / \log n$ (púrpura), para una ejecución aleatorio: