OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..3614
MAPLE
N:= 1000: # to get all terms <= N
kmax:= floor(LambertW(2*ln(2)*N)/ln(2)):
S:= {};
for k from 2 to kmax do
bmax:= floor((N/k)^(1/(k-1)));
S:= S union select(`<=`, {seq(seq(b^k-a^k, a=0..b-1), b=1..bmax), seq(seq(b^k+a^k, a=0..b), b=1..floor(N^(1/k)))}, N);
od:
sort(convert({$1..N} minus S, list)); # Robert Israel, Jan 03 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
David W. Wilson, Feb 09 2003
STATUS
approved