A371702 a(n) is the least positive k such that Product_{i=1..k} (1 + 1/(2*i-1)) >= n.

1, 1, 1, 3, 5, 8, 12, 16, 21, 26, 32, 39, 46, 54, 63, 72, 82, 92, 103, 115, 128, 141, 154, 169, 184, 199, 215, 232, 250, 268, 287, 306, 326, 347, 368, 390, 413, 436, 460, 484, 510, 535, 562, 589, 616, 645, 674, 703, 734, 765, 796, 828, 861, 894, 928, 963

Offset

0,4

Links

Table of n, a(n) for n=0..55.

Formula

a(n) ~ n^2 / Pi. - Vaclav Kotesovec, Apr 03 2024

Example

a(4) = 5: (2/1) * (4/3) * (6/5) * (8/7) * (10/9) = 256 / 63 = 4.063492... > 4.

Mathematica

a[n_] := For[k = 1, True, k++, If[(2 k)!!/(2 k - 1)!! >= n, Return[k]]]; Table[a[n], {n, 0, 55}]

Prog

(PARI) a(n) = my(k=1); while (prod(i=1, k, 1 + 1/(2*i-1)) < n, k++); k; \\ Michel Marcus, Apr 03 2024

Crossrefs

Cf. A000165, A001147, A001803, A046161, A072986, A371701.

Sequence in context: A236538 A194176 A186494 * A194180 A000212 A183139

Adjacent sequences: A371699 A371700 A371701 * A371703 A371704 A371705

Keyword

nonn

Author

Ilya Gutkovskiy, Apr 03 2024

Status

approved