OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
FORMULA
Sum_{n>=1} 1/a(n) = Product_{p prime}(1 + 1/(p^2*(p-1))) - Sum_{p prime} 1/(p^2*(p-1)) - 1 - ((Sum_{p prime} (1/(p^2*(p-1))))^2 - Sum_{p prime} (1/(p^4*(p-1)^2)))/2 = 0.00023303003383679282415... . - Amiram Eldar, Dec 23 2025
EXAMPLE
Table of n, a(n) for select n:
n a(n)
-------------------------------------
1 27000 = 2^3 * 3^3 * 5^3
2 54000 = 2^4 * 3^3 * 5^3
3 74088 = 2^3 * 3^3 * 7^3
4 81000 = 2^3 * 3^4 * 5^3
5 108000 = 2^5 * 3^3 * 5^3
6 135000 = 2^3 * 3^3 * 5^4
7 148176 = 2^4 * 3^3 * 7^3
8 162000 = 2^4 * 3^4 * 5^3
9 216000 = 2^6 * 3^3 * 5^3
10 222264 = 2^3 * 3^4 * 7^3
41 1157625 = 3^3 * 5^3 * 7^3
161 9261000 = 2^3 * 3^3 * 5^3 * 7^3
MATHEMATICA
nn = 2^20; s = Union@ Flatten@ Table[a^5*b^4*c^3, {c, Surd[nn, 3]}, {b, Surd[nn/(c^3), 4]}, {a, Surd[nn/(b^4*c^3), 5]}]; Select[s, PrimeNu[#] > 2 &]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael De Vlieger, Dec 22 2025
STATUS
approved
