OFFSET
0,2
COMMENTS
B(n,p)=sum(i=0,n,p^i*sum(j=0,i,binomial(n,j)*B(j))) where B(k)=k-th Bernoulli number
FORMULA
a(n)=(1/10)*(11*36^n+2*9^n+3*4^n-6); a(0)=1, a(1)=42, a(2)=1446, a(3)=51486 and a(n)=50*a(n-1)-553*a(n-2)+1800*a(n-3)-1296*a(n-4)
PROG
(PARI) a(n)=sum(i=0, 2*n, 6^i*sum(j=0, i, binomial(2*n, j)*bernfrac(j)))/bernfrac(2*n)
(Maxima) a[0]:1$ a[1]:42$ a[2]:1446$ a[3]:51486$ a[n]:=(1/10)*(11*36^n+2*9^n+3*4^n-6)$ A096048(n):=a[n]$ makelist(A096048(n), n, 0, 30); /* Martin Ettl, Nov 13 2012 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Jun 17 2004
STATUS
approved