Es que ya se ha señalado que la convergencia de la siguiente manera fácilmente a partir de la
El primer Número Teorema (Adam Hughes), o incluso de menos precisa
estimaciones tales como Mertens (Daniel Fischer; el de Chebyshev obligado
$p_k \gg k \log k$ sería suficiente) $-$, pero también que la convergencia
es frustrantemente lento, con la suma de más de $p > x$ descomposición sólo como $1/\log x$.
He aquí otro enfoque, a través de Euler producto
$$
\sum_{n=1}^\infty \frac1{n^s} =: \zeta(s) = \prod_p \frac1{1-p^{s}}
$$
(el producto se extiende sobre todos los números primos $p$),
lo que hace que sea práctico para estimar el $\sum_p 1 / p \log p$ a de alta precisión.
El valor numérico que resulta ser 1.636616323351260868569658...;
en estos días, una vez que uno tiene una expansión decimal, Google encontrará a menudo
una referencia, y aquí el cálculo se indica en el arXiv preprint
Richard J. Mathar:
Veinte Dígitos de Algunas de las Integrales de la Primer Función Zeta
(preprint, 2008), arXiv: 0811.4738
(ver Tabla 2.4 en la parte inferior de la página 4).
Tomando logaritmos del producto de Euler nos encontramos
$$
\log \zeta(s) = \sum_p -\log (1-p^{s})
= \sum_p \frac1{p^s}
+ \sum_p \frac1{2p^{2}}
+ \sum_p \frac1{3p^{3}}
+ \cdots.
$$
Podemos aislar la primera contribución $\sum_p 1/p^s$, teniendo un adecuado
combinación lineal de $\log \zeta(s)$, $\log \zeta(2s)$, $\log \zeta(3s)$, etc.,
encontrar
$$
\sum_p \frac1{p^s} = \sum_{m=1}^\infty \frac{\mu(m)}{m} \log \zeta(ms),
$$
donde $\mu$ es la función de Möbius.
Ahora desde $1 / \log p = \int_1^\infty p^{-s} ds$, podemos integrar
la fórmula para $\sum_p \frac1{p^s}$ encontrar
$$
\sum_p \frac1{p \log p}
= \sum_{m=1}^\infty \frac{\mu(m)}{m} \int_1^\infty \log \zeta(ms) \, ds
= \sum_{m=1}^\infty \frac{\mu(m)}{m^2} \int_m^\infty \log \zeta(s) \, ds.
$$
Esto claramente converge, porque
$\int_m^\infty \log \zeta(ms)$ decae como $2^{-m}$ grandes $m$, mientras que
para $s=1+\epsilon$ sabemos que $\zeta(s)$ crece como $1/\epsilon$, por lo que
$\log\zeta(s)$ crece sólo como $\log(1/\epsilon)$ que es integrable.
Por otra parte, podemos evaluar la $\zeta(s)$ (e $(s-1)\zeta(s)$ cerca de $s=1$)
alta precisión en el uso de Euler-Maclaurin suma de $\sum_{n=1}^\infty 1/n^s$.
Esto hace que la fórmula susceptible de técnicas conocidas para la integración numérica.
Una de estas técnicas es implementado en gp, y el comando
sum(m=1,199,moebius(m)*intnum(s=m,200,log(zeta(s)))/m^2)
sólo toma unos pocos segundos para devolver el valor numérico que se informó anteriormente.