|
|
A096054
|
|
a(n) = (36^n/6)*B(2n,1/6)/B(2n) where B(n,x) is the n-th Bernoulli polynomial and B(k)=B(k,0) is the k-th Bernoulli number.
|
|
3
|
|
|
1, 91, 3751, 138811, 5028751, 181308931, 6529545751, 235085301451, 8463265086751, 304679288612371, 10968470088963751, 394865064451017691, 14215143591303768751, 511745180725868773411, 18422826609078989373751
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
FORMULA
|
a(n) = (1/12)*(36^n - 2*9^n - 3*4^n+6).
G.f.: x*(1 - 6*x)*(1 + 47*x + 36*x^2) / ((1 - x)*(1 - 4*x)*(1 - 9*x)*(1 - 36*x)).
a(n) = 50*a(n-1) - 553*a(n-2) + 1800*a(n-3) - 1296*a(n-4) for n>4.
(End)
|
|
PROG
|
(PARI) a(n)=(1/12)*36^n-(1/6)*9^n-(1/4)*4^n+1/2
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|