tratando de relacionar $P(x) = \sum_{p \in \mathcal{P}} \frac{x^p}{p!}$ al logaritmo de la de Riemann zeta función, me sale :
$$P(x) = \sum_{p \in \mathcal{P}} \frac{x^p}{p!} = \frac{x}{(2 i \pi)^2}\int_{|u|=r} \frac{1}{(u-2x)(u-x)}\left( \int_{c-i\infty}^{c+i\infty} (-\ln u)^{-s} \Gamma(s) \sum_{m=1}^\infty \frac{\mu(m)}{m} \ln \zeta(ms) ds\right) du$$
que espero que deben ser simplifiable.
vamos a :
$$f(z) = \sum_{n=1}^\infty a_n z^n, \qquad\qquad g(z) = \sum_{n=1}^\infty \frac{a_n}{n!} z^n, \qquad \qquad F(s) = \sum_{n=1}^\infty a_n n^{-s}$$
(aquí se $a_n = 1$ fib $n$ es primo, por lo tanto $F(s) = \sum_{m=1}^\infty \frac{\mu(m)}{m} \ln \zeta(ms)$)
por la integral de Cauchy fórmula, para cualquier $r < R$ el radio de convergencia de $f$ :
$$a_n = \frac{n!}{2 i \pi}\int_{|u|=r} \frac{f(u)}{(u-z)^{n+1}} du$$
por tanto, para cualquier $|z| < r/2$ : $$g(z) = \frac{1}{2 i \pi}\sum_{n=1}^\infty z^n \int_{|u|=r} \frac{f(u)}{(u-z)^{n+1}} du = \frac{1}{2 i \pi}\int_{|u|=r} f(u) \sum_{n=1}^\infty \frac{z^n}{(u-z)^{n+1}} du$$ $$ =\frac{z}{2 i \pi}\int_{|u|=r} \frac{f(u)}{(u-2z)(u-z)} du$$
mientras que a partir de $n^{-s} \Gamma(s) = \int_0^\infty x^{s-1} e^{-nx} dx$ (que se extiende a $f(z)$ $[0,1[$ por la continuación analítica si es necesario) : $$\Gamma(s) F(s) = \int_0^\infty x^{s-1} f(e^{-x})dx$$
por lo tanto, por la inversa Mellin transformar : $$f(e^{-x}) = \frac{1}{2 i \pi}\int_{c-i\infty}^{c+i\infty} x^{-s} F(s) \Gamma(s) ds$$
y $$g(z) = \frac{z}{2 i \pi}\int_{|u|=r} \frac{f(u)}{(u-2z)(u-z)} du = \frac{z}{(2 i \pi)^2}\int_{|u|=r} \frac{1}{(u-2z)(u-z)}\left( \int_{c-i\infty}^{c+i\infty} (-\ln u)^{-s} F(s) \Gamma(s) ds\right) du = \ldots$$