login
A359738
a(n) = [x^n] (2*x^4 + 2*x^3 + 2*x^2 + x + 1)/(x^2 + 1).
1
1, 1, 1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1
OFFSET
0
LINKS
Madeline Beals-Reid, A Quadratic Relation in the Bernoulli Numbers, The Pump Journal of Undergraduate Research, 6 (2023), 29-39.
FORMULA
Let B(x) = x/(1 - exp(-x)), the e.g.f. of the Bernoulli numbers with B(1) = 1/2.
a(n) = signum([x^n] B(x)^2) = signum([x^n] z^2 / (exp(-z) - 1)^2).
a(n) = signum([x^n] (x + 1)*B(x) - x*B'(x)).
a(n) = A057077(n-3), n>2. - R. J. Mathar, Jan 27 2025
E.g.f.: x*(2 + x) + cos(x) - sin(x). - Stefano Spezia, Jan 27 2025
MAPLE
ogf := (2*z^4 + 2*z^3 + 2*z^2 + z + 1)/(z^2 + 1):
ser := series(ogf, z, 100): seq(coeff(ser, z, n), n = 0..74);
CROSSREFS
KEYWORD
sign,easy
AUTHOR
Peter Luschny, Jan 23 2023
STATUS
approved