OFFSET
1,1
COMMENTS
EXAMPLE
751 is a term because psi(x) = phi(751) = 750 has a solution that is x = 625 while there is no solution for sigma(y) = phi(751) = 750.
PROG
(PARI) is1(n) = my(N=eulerphi(n)); for(k=1, N, if(sigma(k)==N, return(1))); 0;
a001615(n) = my(f=factor(n)); prod(i=1, #f~, f[i, 1]^f[i, 2] + f[i, 1]^(f[i, 2]-1));
is2(n) = my(N=eulerphi(n)); for(k=1, N, if(a001615(k)==N, return(1))); 0;
isok(n) = !is1(n) && is2(n); \\ after Charles R Greathouse IV at A001615
CROSSREFS
KEYWORD
nonn
AUTHOR
Altug Alkan, Aug 25 2017
STATUS
approved