OFFSET
1,2
COMMENTS
The least term that is not a power of 2 is a(16) = 2^24 * 3^8.
Indices of records of the multiplicative function f(n) with f(p^e) = A008680(e).
All the terms are cubefull numbers since f(1) = 1 and f(n) = 0 if n is a noncubefull number.
The corresponding record values are 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 20, ... (see the link for more values).
Every exponent must be the index of the first occurrence of A008680(e) in A008680. So possible exponents of prime factors of terms are 0, 8, 12, 15, 18, 20, ... - David A. Corneth, Jun 30 2025
LINKS
David A. Corneth, Table of n, a(n) for n = 1..810 (first 216 terms from Amiram Eldar, terms <= 10^120)
EXAMPLE
256 in the sequence as 256 = 1^3 * 4^4 * 1^5 = 2^3 * 1^4 * 2^5 so there are two ways to write 256 as b^3 * c^4 * d^5, with b, c and d >= 1 and no smaller positive integer can be written in at least two ways like that. - David A. Corneth, Jun 30 2025
MATHEMATICA
f[p_, e_] := Floor[(1+(-1)^e)*(-1)^Floor[e/2]/8 + (e^2 + 12*e + 90)/120]; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; With[{lps = Cases[Import["https://oeis.org/A025487/b025487.txt", "Table"], {_, _}][[;; , 2]]}, sm = -1; seq = {}; Do[s1 = s[lps[[i]]]; If[s1 > sm, sm = s1; AppendTo[seq, lps[[i]]]], {i, 1, Length[lps]}]; seq]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jun 10 2025
STATUS
approved
