OFFSET
1,3
REFERENCES
R. Schumacher, Explicit formulas for sums involving the squares of the first n Tribonacci numbers, Fib. Q., 58:3 (2020), 194-202. (Note that this paper uses an offset for the tribonacci numbers that is different from that used in A000073.)
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..1000
Index entries for linear recurrences with constant coefficients, signature (7,-9,-7,-56,96,108,252,-162,-114,-318,126,-16,136,-36,12,-21,3,-1,1).
FORMULA
From G. C. Greubel, Nov 22 2021: (Start)
a(n) = Sum_{j=0..n-1} j^2*A000073(j+1)^2.
G.f.: x^2*(1 -2*x +15*x^2 +62*x^3 -97*x^4 +96*x^5 +73*x^6 -64*x^7 -57*x^8 -194*x^9 -127*x^10 -138*x^11 -55*x^12 -12*x^13 -9*x^14 -4*x^15)/((1-x)*(1 +x +x^2 -x^3)^3*(1 -3*x -x^2 -x^3)^3). (End)
MATHEMATICA
T[n_]:= T[n]= If[n<2, 0, If[n==2, 1, T[n-1] +T[n-2] +T[n-3]]]; (* A000073 *)
A337285[n_]:= Sum[j^2*T[j+1]^2, {j, 0, n-1}];
Table[A337285[n], {n, 40}] (* G. C. Greubel, Nov 22 2021 *)
PROG
(Magma) R<x>:=PowerSeriesRing(Integers(), 40); [0] cat Coefficients(R!( x^2*(1 -2*x+15*x^2+62*x^3-97*x^4+96*x^5+73*x^6-64*x^7-57*x^8-194*x^9-127*x^10-138*x^11 -55*x^12-12*x^13-9*x^14-4*x^15)/((1-x)*(1+x+x^2-x^3)^3*(1-3*x-x^2 -x^3)^3) )); // G. C. Greubel, Nov 22 2021
(Sage)
@CachedFunction
def T(n): # A000073
if (n<2): return 0
elif (n==2): return 1
else: return T(n-1) +T(n-2) +T(n-3)
def A337285(n): return sum( j^2*T(j+1)^2 for j in (0..n-1) )
[A337285(n) for n in (1..40)] # G. C. Greubel, Nov 22 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, Sep 12 2020
STATUS
approved