OFFSET
1,1
EXAMPLE
30 = 2*3*5 is a term because psi(phi(30)) = phi(psi(30)-30).
60 = 30*2 is not a term because it is not a squarefree number.
PROG
(PARI) a001615(n) = my(f=factor(n)); prod(i=1, #f~, f[i, 1]^f[i, 2] + f[i, 1]^(f[i, 2]-1));
isok(n) = issquarefree(n) && a001615(eulerphi(n))==eulerphi(a001615(n)-n); \\ after Charles R Greathouse IV at A001615
CROSSREFS
KEYWORD
nonn
AUTHOR
Altug Alkan, Sep 08 2017
STATUS
approved