OFFSET
1,2
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
phi(57) = 36 = 8*9/2, a triangular number; so 57 is a term of the sequence.
MAPLE
filter := m -> issqr(1 + 8*numtheory:-phi(m)) : select(filter, [$(1 .. 700)]);
MATHEMATICA
Select[Range[700], IntegerQ[Sqrt[8 * EulerPhi[#] + 1]] &] (* Amiram Eldar, Feb 27 2023 *)
PROG
(PARI) isok(m) = ispolygonal(eulerphi(m), 3); \\ Michel Marcus, Feb 27 2023
(Python)
from itertools import islice, count
from sympy.ntheory.primetest import is_square
from sympy import totient
def A360944_gen(startvalue=1): # generator of terms >= startvalue
return filter(lambda n:is_square((totient(n)<<3)+1), count(max(1, startvalue)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Bernard Schott, Feb 26 2023
STATUS
approved