OFFSET
1,2
COMMENTS
Indices of records of A385416.
The corresponding record values are 1, 2, 4, 5, 6, 7, 11, 15, 16, 21, 22, 25, 37, 59, 62, 83, ... (see the link for more values).
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..86
Amiram Eldar, Table of n, a(n), A385416(a(n)) for n = 1..86.
FORMULA
The factorizations of the first 5 terms:
n | a(n) | A385416(a(n)) | factorizations
--+-------+---------------|------------------------------------------------------------------
1 | 1 | 1 | 1
2 | 6 | 2 | 2*3, 6
3 | 24 | 4 | 2 * 2 * 2 * 3, 2 * 2 * 6, 3 * 8, 24
4 | 30 | 5 | 2 * 3 * 5, 2 * 15, 3 * 10, 5 * 6, 30
5 | 60 | 6 | 2 * 2 * 3 * 5, 2 * 2 * 15, 2 * 3 * 10, 2 * 5 * 6, 2 * 30, 6 * 10
MATHEMATICA
T[n_, k_] := T[n, k] = If[n == 1 || k == 1, {{}}, Flatten[Table[Map[Prepend[#, d] &, T[d, k/d]], {d, Select[Divisors[k], 1 < # <= n &]}], 1]]; f[n_] := T[n, n];
expOdd[n_] := AllTrue[FactorInteger[n][[;; , 2]], OddQ]; s[n_] := Count[f[n], _?(And @@ (expOdd /@ #) &)];
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, 300}]; seq]
CROSSREFS
KEYWORD
nonn
AUTHOR
Amiram Eldar, Jun 28 2025
STATUS
approved
