OFFSET
1,1
COMMENTS
There are 88 such numbers up to 2*10^9, the last one in this range being 1945428.
LINKS
David W. Wilson, Table of n, a(n) for n = 1..88
PROG
(PARI) sumset_lim(a, b, lim)=my(v=[], u, t); if(a==b, for(i=1, #a, u=List(); for(j=i, #b, t=a[i]+b[j]; if(t>lim, break); listput(u, t)); v=vecsort(concat(v, Vec(u)), , 8)), for(i=1, #a, u=List(); for(j=1, #b, t=a[i]+b[j]; if(t>lim, break); listput(u, t)); v=vecsort(concat(v, Vec(u)), , 8))); v
makev(lim)=my(n=floor(sqrt(2*lim)-1/2), v); sumset_lim(v=vector(n, k, k*(k-1)/2), v, lim)
is(n)=for(i=1, #v, if(ispower(n-v[i], 4), return(0)); if(v[i]>n, return(1)))
v=makev(1e5);
for(n=1, 1e5, if(is(n), print1(n", "))) \\ Charles R Greathouse IV, Aug 17 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Giovanni Resta, Jan 15 2006
STATUS
approved