OFFSET
1,3
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Binomial Number.
EXAMPLE
10 = 3^2+1^2; 13 = 7^2-6^2; 54 = 3^3+3^3.
MAPLE
N:= 200: # for all terms <= N
kmax:= floor(LambertW(2*ln(2)*N)/ln(2)):
S:= {0};
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(S, list)); # Robert Israel, Jan 03 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
David W. Wilson, Feb 09 2003
STATUS
approved