OFFSET
0,2
COMMENTS
a(23) = 80200, a(26) = 112015, and there are no other terms below 158000.
MATHEMATICA
f[n_] := Log2[NestWhile[DivisorSigma[0, #] &, n!, # != 2^IntegerExponent[#, 2] &]]; (* A380802 *) seq[len_] := Module[{v = Table[0, {len}], c = 0, k = 1, i}, While[c < len, i = f[k] + 1; If[i <= len && v[[i]] == 0, c++; v[[i]] = k]; k++]; v]; seq[15]
PROG
(PARI) f(n) = {my(m = n!); while(1 << valuation(m, 2) != m, m = numdiv(m)); valuation(m, 2); } \\ A380802
list(len) = {my(v = vector(len), c = 0, k = 1, i); while(c < len, i = f(k) + 1; if(i <= len && v[i] == 0, c++; v[i] = k); k++); v; }
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Feb 04 2025
STATUS
approved
