OFFSET
1,1
LINKS
Eric Weisstein's World of Mathematics, Pratt Certificate.
MATHEMATICA
a[1] = 1; a[n_] := a[n] = 1 + Plus @@ (a@ PrimePi@ # & /@ First /@ FactorInteger[ Prime@ n - 1]); t = Table[ 0, {25}]; k = 2; While[k < 23420001, b = a[k]/2; If[b < 1001 && t[[b]] == 0, t[[b]] = Prime@ k; Print[{b, Prime@ k }]]; k++]; t
CROSSREFS
KEYWORD
nonn,more
AUTHOR
STATUS
approved