OFFSET
1,1
COMMENTS
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
PROG
(Python)
from __future__ import division
from sympy import nextprime, prevprime
A167807_list = []
for i in range(3, 10**6):
n = i*(i+1)*(2*i+1)//6
p2 = prevprime(n//3)
p1, p3 = prevprime(p2), nextprime(p2)
q = p1+p2+p3
while q <= n:
if q == n:
A167807_list.append(n)
p1, p2, p3 = p2, p3, nextprime(p3)
q = p1+p2+p3 # Chai Wah Wu, Dec 31 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Nov 12 2009
EXTENSIONS
a(6)-a(22) from Donovan Johnson, Nov 15 2009
STATUS
approved