|
|
A135061
|
|
a(n) = minimum (floor(n^3/m) + m) for any integer m >= 1.
|
|
2
|
|
|
2, 5, 10, 16, 22, 29, 37, 45, 54, 63, 72, 83, 93, 104, 116, 128, 140, 152, 165, 178, 192, 206, 220, 235, 250, 265, 280, 296, 312, 328, 345, 362, 379, 396, 414, 432, 450, 468, 487, 505, 525, 544, 563, 583, 603, 623, 644, 665, 686, 707, 728, 749, 771, 793, 815, 838, 860, 883, 906, 929, 952, 976, 1000
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
FORMULA
|
a(n) = floor(2*n^(3/2)), with the minimum occurring at m = ceiling(n^(3/2)). See link for proof. - Robert Israel, Mar 03 2017
|
|
MAPLE
|
f:= proc(n) local m, U, L;
U:= floor(2*n^(3/2));
m:=floor(n^(3/2));
L:= floor(min(n^3/m+m, n^3/(m+1)+m+1));
if L <> U then error("Conjecture fails") fi;
L
end proc:
|
|
PROG
|
(PARI) a(n)=local( minsum=0, cursum =0, minm=0, lastminsum=0); minsum = n^3 + 1; lastminsum= n^3 + 1; minm =1; for(m=1, n^3, cursum = floor(n^3/m + m); lastminsum = minsum; if(cursum < minsum, minsum = cursum); if(cursum < lastminsum, minm=m); ); print1("minm="minm" "); print1("minsum="minsum" ");
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
a(2) corrected and more terms added by Robert Israel, Mar 01 2017
|
|
STATUS
|
approved
|
|
|
|