OFFSET
1,1
MATHEMATICA
(* first do *) Needs["Combinatorica`"] (* then *) f[n_] := Block[{d = Divisors@n, k, mx}, k = 1 + Length@d; mx = 2^Length[d]; While[k < mx && !PrimeQ[Plus @@ NthSubset[k, d]], k++ ]; If[k == mx, Length@d, 0]];
a[n_] := a[n] = Module[{p = If[n == 1, 2, a[n-1]]}, While[f[p^n] == 0, p = NextPrime[p]]; p]; Array[a, 13] (* second part of the program added by Amiram Eldar, Jul 30 2024 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Robert G. Wilson v, Jun 01 2009
EXTENSIONS
Offset changed to 1 and name corrected by Amiram Eldar, Jul 30 2024
STATUS
approved