OFFSET
1,2
COMMENTS
a(17) > 10^11, if it exists.
MATHEMATICA
f[p_, e_] := 2^DigitCount[e, 2, 1] - 1; d[1] = 0; d[n_] := Plus @@ f @@@ FactorInteger[n]; seq[len_, nmax_] := Module[{s = Table[0, {len}], dprev = 0, n = 2, c = 1, k = 1}, s[[1]] = 1; While[k < len && n < nmax, d1 = d[n]; If[d1 == dprev, c++, If[c > 0 && c <= len && s[[c]] == 0, k++; s[[c]] = n - c]; c = 1]; n++; dprev = d1]; TakeWhile[s, # > 0 &]]; seq[8, 10^4]
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Nov 12 2021
STATUS
approved