|
MATHEMATICA
|
abQ[n_] := DivisorSigma[1, n] > 2*n; nobAbQ[n_] := And @@ abQ /@ DivisorSigma[{0, 1}, n]; seq[len_, nmax_] := Module[{s = Table[0, {len}], count = 0, n = 1, n1, d}, While[n < nmax && count < len, n1 = n; If[nobAbQ[n], While[nobAbQ[++n1]]; d = n1 - n; If[d <= len && s[[d]] == 0, count++; s[[d]] = n]]; n = n1 + 1]; TakeWhile[s, # > 0 &]]; seq[6, 10^6]
|