Su suma contiene la suma de todos los recíprocos de potencias de números naturales no unitarios de primos, es decir, cada número de la forma $\frac{1}{p^k}$ para el primo $p$, y $k\in \{2,3,...\}=\mathbb{N}\setminus\{1\}$. Dado que cada número de esta forma está determinado de forma única por la elección de $p$ y $k$, obtenemos
$\sum_{p}\sum_{i=2}^{\infty}\dfrac{1}{p^i}=\sum_p\dfrac{1}{p(p-1)}$ donde $p$ es cada primo.
Dado que $\sum_{i=1}^{\infty}\dfrac{1}{p^i}=\dfrac{1}{p-1}$ (esto se sigue de una teoría numérica básica; una forma fácil de verlo es considerar el decimal repetitivo, $0.111...$ en base $p$), hemos eliminado el término principal, $\dfrac{1}{p}$, por lo tanto
$\sum_{p}\sum_{i=2}^{\infty}\dfrac{1}{p^i}=\dfrac{1}{p-1}-\dfrac{1}{p}=\dfrac{p}{p(p-1)}-\dfrac{p-1}{p(p-1)}=\dfrac{p-(p-1)}{p(p-1)}=\dfrac{1}{p(p-1)}$
En el caso que nos ocupa, ya no nos estamos concentrando en los primos, pero el resultado de conteo anterior se mantiene para los números no primos, $n$, el problema ahora es: ¿cómo evitar contar dos veces números como 4^2=4^2=16? De nuevo, la solución es elegir primos, el primo aquí es el $2$, por lo tanto lo contaremos como el $2^4$ y esta táctica maneja efectivamente todas las potencias de 4. Pero eso no resuelve el problema planteado por números como $36=2^2\cdot 3^2$. Aquí tenemos primos distintos que son parte de la factorización, en este caso, podemos contar esto como $6^2$. También tenemos que lidiar con casos como
$4^2\cdot 3^2=2^4\cdot 3^2=12^2=144$
Para esto, lo consideramos como $12^2$. En este punto, tal vez esté surgiendo el patrón: Queremos tomar nuestra suma de recíprocos de potencias (no unitarias) sobre cada número cuyos exponentes de su factorización prima tienen mcd de 1. Es decir, estamos interesados en números:
$b=p_1^{k_1}\cdot p_2^{k_2}\cdot ...\cdot p_f^{k_f}$
para los cuales $\gcd(k_1,...,k_f)=1$. La razón por la que este es el criterio que requerimos es que buscamos eliminar recuentos de números que ya tenemos, y cualquier número que tenga el anterior $\gcd(k_1,...,k_f)=r\neq 1$, entonces podemos tomar la raíz $r$ de ese número, y obtener un número que ciertamente tiene $\gcd(k_1,...,k_f)=1, por lo tanto tenemos una forma única de representar cada número que nos interesa como una potencia de un número en la forma anterior con $\gcd(k_1,...,k_f)=1$.
Entonces, nuestra suma es, en términos de la notación anterior:
$\sum_b\sum_{i=2}^{\infty}\dfrac{1}{b^i}=\sum_b\dfrac{1}{b(b-1)}$.
No estoy seguro si hay una forma agradable de caracterizar los números relevantes $b$ o no, pero siempre que se pueda encontrar una caracterización lo suficientemente buena, esto podría resultar útil.
0 votos
La primera suma debería tener potencias solo una vez, mis disculpas si eso no está claro
0 votos
La suma parcial para el límite $n = 10^9$ es 0.8744322328535789186668554166 (Calculado con PARI/GP)
0 votos
Sí, lo he intentado usando el programa mencionado anteriormente. Aunque, el mío solo puede hacerlo realistamente hasta el orden de 10^5 términos...
0 votos
@gebra Toma un conjunto de bits para filtrar primero las potencias perfectas. Un array de un GB te lleva a $2^{33}$, el filtrado no lleva mucho tiempo. Suma "hacia abajo", comenzando desde los recíprocos de las potencias más grandes para reducir los efectos del redondeo de los puntos flotantes.
0 votos
El problema con el tamizado es que no se puede ejecutar una sola vez y producir valores que se acerquen al límite; primero tendrías que elegir el límite.
0 votos
Permitir $1$ como un número perfecto solo añade $1$ a la suma