OFFSET
1,1
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
4934 = 2*2467 is a term because psi(4934)*phi(4934) = (2^2 - 1)*(2467^2 - 1) = 2^3*3^3*137*617 is not a totient number (A002202).
PROG
(PARI) is(k) = if(isprime(k), 0, my(f = factor(k)); !istotient(prod(i = 1, #f~, (f[i, 1]^2 - 1) * f[i, 1]^(2*f[i, 2] - 2)))); \\ Amiram Eldar, Dec 05 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Altug Alkan, Sep 21 2017
STATUS
approved