OFFSET
1,1
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000 (n = 1..70 from Colin Barker)
EXAMPLE
12160 is in the sequence because 12160+12545+12936 (consecutive octagonal numbers) = 37641 = 12541+12547+12553 (consecutive primes).
PROG
(PARI) L=List(); forprime(p=2, 20000000, q=nextprime(p+1); r=nextprime(q+1); t=p+q+r; if(issquare(36*t-180, &sq) && (sq-12)%18==0, u=(sq-12)\18; listput(L, 3*u^2-2*u))); Vec(L)
(Python)
from __future__ import division
from sympy import prevprime, nextprime
A298312_list, n, m = [], 1, 30
while len(A298312_list) < 10000:
k = prevprime(m//3)
k2 = nextprime(k)
if prevprime(k) + k + k2 == m or k + k2 + nextprime(k2) == m:
A298312_list.append(n*(3*n-2))
n += 1
m += 18*n + 3 # Chai Wah Wu, Jan 22 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Colin Barker, Jan 17 2018
STATUS
approved