|
MATHEMATICA
|
f[p_, e_] := (-1)^e; q[1] = True; q[n_] := Plus @@ f @@@ FactorInteger[n] == 0; seq[len_, nmax_] := Module[{s = Table[0, {len}], v = {1}, n = 2, c = 0, m}, While[c <= len && n <= nmax, If[q[n], v = Join[v, {n}], m = Length[v]; v = {}; If[0 <= m <= len && s[[m]] == 0, c++; s[[m]] = n - m]]; n++]; s]; seq[4, 10^6]
|