A376956 a(n) = least k such that n^(2k)/(2 k)! < 1.

1, 1, 2, 4, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 21, 22, 24, 25, 26, 28, 29, 30, 32, 33, 34, 36, 37, 38, 40, 41, 42, 44, 45, 46, 48, 49, 51, 52, 53, 55, 56, 57, 59, 60, 61, 63, 64, 65, 67, 68, 69, 71, 72, 74, 75, 76, 78, 79, 80, 82, 83, 84, 86, 87, 88

Offset

0,3

Comments

The numbers n^(2k)/(2 k)! are the coefficients in the Maclaurin series for cos x when x = 1. If m>a(n), then n^(2k)/(2 k)! < 1.

Links

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

Formula

a(n) ~ exp(1)*n/2 - log(n)/4. - Vaclav Kotesovec, Oct 13 2024

Mathematica

a[n_] := Select[Range[200], n^(2 #)/(2 #)! < 1 &, 1];
Flatten[Table[a[n], {n, 0, 200}]]

Prog

(Python)
from itertools import count
from math import gcd
def A376956(n):
    a, b = 1, 1
    for k in count(1):
        a *= n**2
        b *= (m:=k<<1)*(m-1)
        if a < b: return k
        c = gcd(a, b)
        a, b = a//c, b//c # Chai Wah Wu, Oct 16 2024

Crossrefs

Cf. A370507, A376284, A376952, A376953, A376954, A376955, A376957, A376958, A376959, A376960.

Sequence in context: A154112 A039101 A140102 * A288866 A084437 A037083

Adjacent sequences: A376953 A376954 A376955 * A376957 A376958 A376959

Keyword

nonn

Author

Clark Kimberling, Oct 12 2024

Status

approved