OFFSET
1,1
LINKS
Harry J. Smith, Table of n, a(n) for n = 1..500
EXAMPLE
If k = p(p+2) is a product of twin primes then phi(k) + 1 = p^2, sigma(k) + 1 = (p+2)^2, so k is in the sequence, A037074 a proper subset. There are many solutions not of this form, such as 72, 123, and 165.
PROG
(PARI) { n=0; for (a=1, 10^9, if (issquare(eulerphi(a) + 1) && issquare(sigma(a) + 1), write("b063532.txt", n++, " ", a); if (n==500, break)) ) } \\ Harry J. Smith, Aug 25 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Aug 02 2001
STATUS
approved