OFFSET
0,5
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..449
FORMULA
a(n) = Sum_{k=0..n} Stirling1(n,k)*A000045(k).
From Vaclav Kotesovec, Jan 21 2019: (Start)
a(n) = -(-1)^n * cos(sqrt(5)*Pi/2) * (Gamma((3 + sqrt(5))/2) * Gamma(n - (1 + sqrt(5))/2) - Gamma((3 - sqrt(5))/2) * Gamma(n + (sqrt(5) - 1)/2)) / (Pi*sqrt(5)).
a(n) ~ -(-1)^n * n! / (sqrt(5) * Gamma((sqrt(5)-1)/2) * n^((3 - sqrt(5))/2)).
a(n) = -2*(n-2)*a(n-1) - (n^2 - 5*n + 5)*a(n-2). (End)
MATHEMATICA
FullSimplify[nmax = 23; CoefficientList[Series[2 Sqrt[1 + x] Sinh[Sqrt[5] Log[1 + x]/2]/Sqrt[5], {x, 0, nmax}], x] Range[0, nmax]!]
Table[Sum[StirlingS1[n, k] Fibonacci[k], {k, 0, n}], {n, 0, 23}]
PROG
(PARI) {a(n) = sum(k=0, n, stirling(n, k, 1)*fibonacci(k))};
vector(25, n, n--; a(n)) \\ G. C. Greubel, Feb 07 2019
(Magma) [(&+[StirlingFirst(n, k)*Fibonacci(k): k in [0..n]]): n in [0..25]]; // G. C. Greubel, Feb 07 2019
(Sage) [sum((-1)^(k+n)*stirling_number1(n, k)*fibonacci(k) for k in (0..n)) for n in (0..25)] # G. C. Greubel, Feb 07 2019
CROSSREFS
KEYWORD
sign
AUTHOR
Ilya Gutkovskiy, Jan 20 2019
STATUS
approved