login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

a(n) = floor( (4/5)*( (9/4)^(n+1)-1 ) ).
4

%I #20 Dec 03 2023 18:39:47

%S 1,3,8,19,45,102,232,524,1181,2659,5984,13466,30300,68177,153400,

%T 345151,776590,1747330,3931495,8845865,19903197,44782195,100759939,

%U 226709865,510097199,1147718699,2582367075,5810325919,13073233320,29414774972,66183243689,148912298302

%N a(n) = floor( (4/5)*( (9/4)^(n+1)-1 ) ).

%D N. Tokuda, An Improved Shellsort, IFIP Transactions, A-12 (1992) 449-457.

%H Marcin Ciura, <a href="https://doi.org/10.1007/3-540-44669-9_12">Best Increments for the Average Case of Shellsort</a>, in R. Freivalds, (ed.), Fundamentals of Computation Theory: 13th International Symposium, FCT 2001, Riga, Latvia, August 2001, Lecture Notes in Computer Science, vol. 2138, Springer, pp. 106-117.

%t A361506[n_]:=Floor[4/5((9/4)^(n+1)-1)];Array[A361506,50,0] (* _Paolo Xausa_, Dec 02 2023 *)

%Y Other sequences used for Shell sort: A003462, A033622, A036562, A036564, A036569, A055875, A055876, A108870, A361507.

%K nonn

%O 0,2

%A _N. J. A. Sloane_, Mar 20 2023, at the suggestion of _Don Knuth_.

%E Some terms corrected by _Paolo Xausa_, Dec 02 2023