OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..5000
Index entries for linear recurrences with constant coefficients, signature (3, -3, 1).
FORMULA
From G. C. Greubel, Aug 28 2017: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), n >= 3.
a(n) = n*(n+4), n >= 3.
G.f.: x*(7 - 7*x + 4*x^3 - 2*x^4)/(1 - x)^3.
E.g.f.: x*(x+5)*exp(x) + 2*x + x^2. (End)
MAPLE
A033676 := proc(n) local a, d; a := 0 ; for d in numtheory[divisors](n) do if d^2 <= n then a := max(a, d) ; fi; od: a; end: A163280 := proc(n, k) local r, T ; r := 0 ; for T from k^2 by k do if A033676(T) = k then r := r+1 ; if r = n then RETURN(T) ; fi; fi; od: end: A164005 := proc(n) if n = 0 then 0; else A163280(5, n) ; fi; end: seq(A164005(n), n=0..80) ; # R. J. Mathar, Aug 09 2009
MATHEMATICA
Join[{0, 7, 14}, Table[n*(n + 4), {n, 3, 50}]] (* G. C. Greubel, Aug 28 2017 *)
PROG
(PARI) x='x+O('x^50); concat([0], Vec(x*(7 - 7*x + 4*x^3 - 2*x^4)/(1 - x)^3)) \\ G. C. Greubel, Aug 28 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Omar E. Pol, Aug 08 2009
EXTENSIONS
Extended by R. J. Mathar, Aug 09 2009
STATUS
approved