OFFSET
1,1
COMMENTS
a(1)=491525857=A052239(12).
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..100
PROG
(Python)
from gmpy2 import is_prime, next_prime
A287550_list, p = [], 2
q, r, s = p+72, p+144, p+216
while s <= 10**10:
np = next_prime(p)
if np == q and is_prime(r) and is_prime(s) and next_prime(q) == r and next_prime(r) == s:
A287550_list.append(p)
p, q, r, s = np, np+72, np+144, np+216 # Chai Wah Wu, Jun 03 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, May 26 2017
EXTENSIONS
a(8)-a(21) from Chai Wah Wu, Jun 03 2017
STATUS
approved