OFFSET
1,2
COMMENTS
This is sequence A090482(n) for prime n.
MATHEMATICA
nn = 25; t = Table[-1, {nn}]; t[[1]] = 0; t[[8]] = 0; cnt = 2; p = 1; While[cnt < nn, p = NextPrime[p]; s = DivisorSigma[0, p - 1] + DivisorSigma[0, p + 1]; If[PrimeQ[s], i = PrimePi[s]; If[i <= nn && t[[i]] == -1, t[[i]] = p; cnt++]]]; t (* T. D. Noe, Apr 28 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Apr 23 2011
STATUS
approved