login
Position of 1 + n^3 in A024670 (distinct sums of cubes of distinct positive integers).
1

%I #11 Jan 22 2024 00:03:37

%S 1,2,4,7,11,16,21,28,35,43,52,61,73,84,97,109,124,140,155,173,190,210,

%T 228,250,271,292,316,338,365,391,417,444,472,501,531,560,593,626,659,

%U 694,728,765,802,842,880,919,960,1001,1045,1086,1131,1176,1222,1268,1313

%N Position of 1 + n^3 in A024670 (distinct sums of cubes of distinct positive integers).

%H Robert Israel, <a href="/A024673/b024673.txt">Table of n, a(n) for n = 2..10000</a>

%p with(priqueue):

%p initialize(PQ):

%p insert([-9,1,2],PQ): v:= 0:

%p R:= NULL: count:= 0: i:= 0:

%p while count < 100 do

%p t:= extract(PQ);

%p w:= -t[1];

%p if floor(w) > v then i:= i+1; v:= floor(w) fi;

%p if t[2] = 1 then R:= R, i; count:= count+1 fi;

%p wp:= w + (t[3]+1)^3 - t[3]^3;

%p insert([-wp,t[2],t[3]+1], PQ);

%p if t[3] = t[2]+1 then

%p wp:= t[3]^3 + (t[3]+1)^3 + 1/10;

%p insert([-wp, t[3], t[3]+1],PQ);

%p fi

%p od:

%p R; # _Robert Israel_, Jan 21 2024

%K nonn

%O 2,2

%A _Clark Kimberling_

%E Offset corrected by _Sean A. Irvine_, Jul 19 2019