OFFSET
1,2
COMMENTS
FORMULA
a(n) = phi(A103521(n)-1).
PROG
(PARI) b(n) = if(n==1, 2, floor(exp(Euler)*n*log(log(n^2))+2.5*n/log(log(n^2))));
f(p) = my(i=0); forprime(q=p+1, b(eulerphi(p-1))+1, i+=(eulerphi(q-1)<=eulerphi(p-1))); i;
forprime(p=2, 2e4, if(f(p)==0, print1(eulerphi(p-1), ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Jianing Song, Feb 11 2019
STATUS
approved