| DATA
|
1, 2, 3, 4, 5, 6, 8, 9, 10, 11, 14, 15, 16, 18, 20, 21, 22, 23, 24, 26, 27, 29, 30, 32, 33, 35, 36, 39, 40, 41, 42, 44, 45, 46, 48, 50, 51, 53, 54, 55, 56, 58, 60, 63, 64, 65, 66, 68, 69, 70, 72, 74, 75, 78, 81, 82, 83, 86, 87, 88, 89, 90, 95, 96, 98, 99, 102, 105, 106, 110, 111
|
| 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, if(eulerphi(k)==r && lcm(znstar(k)[2])==r/2, return(1)); if(k==N, return(0)))
for(n=1, 100, if(isA328413(n), print1(n, ", ")))
|