|
|
A135072
|
|
Minimal values of m associated with A135061.
|
|
2
|
|
|
1, 3, 4, 6, 9, 13, 15, 19, 23, 28, 36, 37, 44, 50, 52, 57, 63, 73, 78, 87, 90, 96, 104, 109, 115, 123, 133, 139, 147, 157, 162, 169, 178, 189, 195, 202, 212, 224, 230, 251, 248, 260, 278, 284, 294, 310, 309, 316, 325, 337, 351, 371, 376, 385, 399, 401, 419, 427, 437, 451, 469, 472, 480, 490, 503, 519
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
If t = floor(2*n^(3/2))+1, then a(n) = 1 + floor((t-sqrt(t^2-4*n^3))/2). - Robert Israel, Mar 06 2017
|
|
MAPLE
|
f:= proc(n) local t; t:= floor(2*n^(3/2))+1; 1 + floor((t-sqrt(t^2-4*n^3))/2) end proc:
|
|
MATHEMATICA
|
a[n_] := With[{t = Floor[2n^(3/2)]+1}, 1 + Floor[(t-Sqrt[t^2-4n^3])/2]];
|
|
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); ); minm; }
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|