OFFSET
1,2
PROG
(PARI) isA328413(n) = my(r=4*n, N=floor(exp(Euler)*r*log(log(r^2))+2.5*r/log(log(r^2)))); for(k=r+1, N+1, if(eulerphi(k)==r && lcm(znstar(k)[2])==r/2, return(1)); if(k==N+1, return(0)))
for(n=1, 100, if(isA328413(n), print1(n, ", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Jianing Song, Oct 14 2019
STATUS
editing