OFFSET
1,2
COMMENTS
There are 57 terms less than 200000. - David Wasserman, Dec 20 2002
LINKS
Robert G. Wilson v, Table of n, a(n) for n = 1..1276 (first 161 terms from Robert Israel)
MAPLE
N:= 80: # to get a(1)..a(N)
A067513:= n -> nops(select(t -> isprime(t+1), numtheory:-divisors(n))):
A:= Vector(N):
count:= 1: A[1]:= 1:
for n from 2 by 2 while count < N do
v:= A067513(n);
if v <= N and A[v]=0 then A[v]:= n; count:= count+1; fi
od:
convert(A, list); # Robert Israel, Jun 08 2018
PROG
(PARI) a067513(n) = sumdiv(n, d, isprime(d+1));
a(n) = {my(k=1); while(a067513(k) != n, k++); k; } \\ Michel Marcus, Jun 09 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladeta Jovovic, Feb 14 2002
EXTENSIONS
More terms from David Wasserman, Dec 20 2002
STATUS
approved