OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
4889 = 2^6 + 3^6 + 4^6 = 64 + 729 + 4096.
66377 = 4^6 + 5^6 + 6^6 = 4096 + 15625 + 46656.
MAPLE
N:= 10^10; # to find all terms <= N
A := {}:
for a from 1 to iroot(N, 6) do
for b from 1 to a-1 while a^6 + b^6 < N do
for c from (a+b) mod 2 + 1 to b-1 by 2 do
r:= a^6 + b^6 + c^6;
if r > N then break fi;
if isprime(r) then A:= A union {r} fi;
od od od:
sort(convert(A, list)); # Robert Israel, Dec 15 2015
MATHEMATICA
Union[Select[Total/@Subsets[Range[20]^6, {3}], PrimeQ]] (* Harvey P. Dale, Apr 20 2013 *)
PROG
(PARI) {m=16; p=m^6; w=[]; for(i=1, m-2, for(j=i+1, m-1, for(k=j+1, m, if((n=i^6+j^6+k^6)<p&&isprime(n), w=concat(w, n))))); w=vecsort(w); for(h=1, #w, print1(w[h], ", "))} /* Klaus Brockhaus, Feb 16 2007 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Tomas Xordan, Feb 11 2007
EXTENSIONS
Edited, corrected and extended by Klaus Brockhaus, Feb 16 2007
STATUS
approved