OFFSET
1,1
COMMENTS
a(15) > 10^12. - Giovanni Resta, Aug 21 2017
EXAMPLE
Prime number 2477 is a term since phi(2477*2478/2) = 1856*1857/2.
PROG
(PARI) isok(n) = isprime(n) && ispolygonal(eulerphi(n*(n+1)/2), 3);
(PARI) is(n) = ispolygonal(eulerphi(n\2+1)*(n-1), 3) && isprime(n) \\ Charles R Greathouse IV, Aug 22 2017
(Python)
from __future__ import division
from sympy.ntheory.primetest import is_square
from sympy import totient, nextprime
A291199_list, p = [], 3
while p < 10**8:
if is_square(8*(p-1)*totient((p+1)//2)+1):
A291199_list.append(p)
p = nextprime(p) # Chai Wah Wu, Aug 22 2017
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Altug Alkan, Aug 20 2017
EXTENSIONS
a(5)-a(14) from Giovanni Resta, Aug 21 2017
STATUS
approved