OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
FORMULA
Sum_{n>=1} 1/a(n) = ((Sum_{p prime} 1/(p*(p-1)))^2 - Sum_{p prime} 1/(p*(p-1))^2)/2 - ((Sum_{p prime} 1/p^2)^2 - Sum_{p prime} 1/p^4)/2 = 0.094618784518255407065... . - Amiram Eldar, Jan 06 2026
EXAMPLE
Table of n, a(n) for select n:
n a(n)
--------------------
1 72 = 2^3 * 3^2
2 108 = 2^2 * 3^3
3 144 = 2^4 * 3^2
4 200 = 2^3 * 5^2
5 216 = 2^3 * 3^3
6 288 = 2^5 * 3^2
7 324 = 2^2 * 3^4
8 392 = 2^3 * 7^2
9 400 = 2^4 * 5^2
10 432 = 2^4 * 3^3
11 500 = 2^2 * 5^3
14 675 = 3^3 * 5^2
MATHEMATICA
nn = 3600; s = Union@ Flatten@ Table[a^2*b^3, {b, Surd[nn, 3]}, {a, Sqrt[nn/b^3]}]; Select[s, And[Length[#] == 2, AnyTrue[#, # > 2 &]] &@ FactorInteger[#][[;; , -1]] &]
PROG
(PARI) isok(k) = {my(e = factor(k)[, 2]); #e == 2 && vecmin(e) > 1 && vecmax(e) > 2; } \\ Amiram Eldar, Jan 06 2026
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael De Vlieger, Dec 23 2025
STATUS
approved
