OFFSET
1,1
COMMENTS
PROG
(PARI) a039654(n) = if(n<2, 0, my(m); while((m=sigma(n)-1)!=n, n=m); n);
a(n) = my(k=0, p=prime(n)); while(a039654(k) != p, k++); k; \\ Michel Marcus, Sep 26 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Sep 25 2017
STATUS
approved