login
A107048
Denominators of coefficients that satisfy: 2^n = Sum_{k=0..n} c(k)*x^k for n>=0, where c(k) = A107047(k)/a(k).
11
1, 1, 4, 108, 6912, 21600000, 2332800000, 1921161110400000, 31476303632793600000, 16727798278915463577600000, 209097478486443294720000000000
OFFSET
0,3
FORMULA
A107047(n)/a(n) = Sum_{k=0..n} T(n, k)*2^k where T(n, k) = A107045(n, k)/A107046(n, k) = [A079901^-1](n, k) (matrix inverse of A079901).
EXAMPLE
2^0 = 1;
2^1 = 1 + 1;
2^2 = 1 + 1*2 + (1/4)*2^2;
2^3 = 1 + 1*3 + (1/4)*3^2 + (7/108)*3^3;
2^4 = 1 + 1*4 + (1/4)*4^2 + (7/108)*4^3 + (77/6912)*4^4.
Initial fractional coefficients are:
A107047/A107048 = {1, 1, 1/4, 7/108, 77/6912, 32387/21600000,
395159/2332800000, 31824093937/1921161110400000,
44855117331581/31476303632793600000, ... }.
PROG
(PARI) {a(n)=denominator(sum(k=0, n, 2^k*(matrix(n+1, n+1, r, c, if(r>=c, (r-1)^(c-1)))^-1)[n+1, k+1]))}
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
Paul D. Hanna, May 10 2005
STATUS
approved