OFFSET
1,1
COMMENTS
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000 (terms 1..560 from R. J. Mathar)
FORMULA
A138853={ p(i)^3+p(j)^3+p(k)^3 ; i>j>k>1 }
PROG
(PARI) isA138853(n)= local( c, d); n>494 && forprime( p=floor( sqrtn( n\3+1, 3))+1, floor( sqrtn( n-151, 3)), d=n-p^3; forprime( q=floor( sqrtn( d\2+1, 3))+1, min( p-1, floor( sqrtn( d-26, 3))), round( sqrtn( c=d-q^3, 3 ))^3==c || next; isprime( round( sqrtn( c, 3 ))) && return(1)))
forstep(n=3^3+5^3+7^3, 10^5, 2, isA138853(n)&print1(n", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
M. F. Hasler, Apr 13 2008
STATUS
approved