OFFSET
1,3
LINKS
Paul D. Hanna, Table of n, a(n) for n = 1..1000
EXAMPLE
G.f.: A(x) = x + x^2 - 3*x^3 + 13*x^4 - 65*x^5 + 351*x^6 - 1989*x^7 + 11650*x^8 - 69900*x^9 + 427168*x^10 - 2648444*x^11 + 16612986*x^12 + ...
A(x)^7 = x^7 + 7*x^8 + x^14 + 7*x^16 - 3*x^21 - 21*x^24 + 13*x^28 + 91*x^32 - 65*x^35 - 455*x^40 + 351*x^42 + 2457*x^48 - 1989*x^49 - 2273*x^56 - 69900*x^63 + ...
where A(x)^7 = A(x^7) + 7*A(x^8).
SPECIFIC VALUES.
A(1/7) = 0.15772708901641335150527781783735289338509817552641...
where A(1/7)^7 = A(1/7^7) + 7*A(1/7^8).
A(1/8) = 0.13674460851389170274849954118414318282336460429579...
where A(1/8)^7 = A(1/8^7) + 7*A(1/8^8).
A(1/9) = 0.12062970906064325012650139044513082950987815459592...
where A(1/9)^7 = A(1/9^7) + 7*A(1/9^8).
PROG
(PARI) {a(n) = my(A, V=[0, 1]); for(i=1, n+1, V=concat(V, 0); A=Ser(V);
V[#V] = polcoef( (subst(A, x, x^7) + 7*subst(A, x, x^8) - A^7 )/7, #V+5) ); polcoef(A, n)}
for(n=1, 50, print1(a(n), ", "))
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, Dec 01 2024
STATUS
approved