OFFSET
1,4
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1000
FORMULA
a(n) = Sum_{j=1..floor((n+1)/2)} A023531(j)*Fibonacci(n-j+1). - G. C. Greubel, Jan 19 2022
MATHEMATICA
Table[t=0; m=3; p=BitShiftRight[n]; n--; While[n>p, t += Fibonacci[n+1]; n -= m++]; t, {n, 120}] (* G. C. Greubel, Jan 19 2022 *)
PROG
(Magma)
b:= func< n, j | IsIntegral((Sqrt(8*j+9) -3)/2) select Fibonacci(n-j+1) else 0 >;
A024318:= func< n | (&+[b(n, j): j in [1..Floor((n+1)/2)]]) >;
[A024318(n) : n in [1..80]]; // G. C. Greubel, Jan 19 2022
(Sage)
def b(n, j): return fibonacci(n-j+1) if ((sqrt(8*j+9) -3)/2).is_integer() else 0
def A024318(n): return sum( b(n, j) for j in (1..floor((n+1)/2)) )
[A024318(n) for n in (1..120)] # G. C. Greubel, Jan 19 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved