OFFSET
1,1
LINKS
Index entries for linear recurrences with constant coefficients, signature (30, -273, 820, -576).
FORMULA
For n > 0, a(n) = 5^(2*n+1)/(2*n+1)*sum(k = 0, 2*n + 1, (1/5)^k*C(2*n + 1, k)*B(k)) where B(k) is the k-th Bernoulli number.
G.f.: x*(16/(1 - 16*x) + 9/(1 - 9*x) + 4/(1 - 4*x) + 1/(1 - x)). - Harvey P. Dale, May 04 2011
a(1) = 30, a(2) = 354, a(3) = 4890, a(4) = 72354, a(n) = 30*a(n-1) - 273*a(n-2) + 820*a(n-3) - 576*a(n-4). - Harvey P. Dale, May 04 2011
MAPLE
MATHEMATICA
Table[1 + 4^n + 9^n + 16^n, {n, 20}] (* or *) LinearRecurrence[ {30, -273, 820, -576}, {30, 354, 4890, 72354}, 20] (* Harvey P. Dale, May 04 2011 *)
PROG
(Magma) [1+4^n+9^n+16^n : n in [1..20]]; // Wesley Ivan Hurt, Nov 26 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Benoit Cloitre, Mar 06 2004
EXTENSIONS
Corrected and extended by Harvey P. Dale, May 04 2011
STATUS
approved