OFFSET
1,1
COMMENTS
a(9) > 10^11, if it exists.
MATHEMATICA
d[n_] := DivisorSigma[0, n] - 2^PrimeNu[n]; seq[len_, nmax_] := Module[{s = Table[0, {len}], dprev = -1, n = 1, c = 0, k = 0}, 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[6, 10^6]
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Nov 14 2021
STATUS
approved