OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
{n: A025457(n) >= 2}. - R. J. Mathar, Jun 15 2018
MAPLE
N:= 2000: # for terms <= N
S2:= {}: S1:= {}:
for x from 1 while x^3 < N do
for y from 1 to x while x^3 + y^3 < N do
for z from 1 to y while x^3 + y^3 + z^3 < N do
for w from 1 to z do
v:= x^3 + y^3 + z^3 + w^3;
if v > N then break fi;
if member(v, S1) then S2:= S2 union {v}
else S1:= S1 union {v}
fi
od od od od:
sort(convert(S2, list)); # Robert Israel, Feb 24 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved