OFFSET
0,3
COMMENTS
Limit a(n)/a(n+1) -> r = -0.166670835025545 where A(r)=0.
MATHEMATICA
kmax = 25;
A[x_] = Sum[a[k] x^k, {k, 0, kmax}];
coes = CoefficientList[A[x]^7 + O[x]^(kmax + 1), x];
r = {a[0] -> 1, a[1] -> 1};
coes = coes /. r;
Do[r = Flatten @ Append[r, Reduce[1 <= coes[[k]] <= 7, a[k-1], Integers] // ToRules];
coes = coes /. r, {k, 3, kmax + 1}];
Table[a[k], {k, 0, kmax}] /. r (* Jean-François Alcover, Jul 26 2018 *)
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, May 20 2003
STATUS
approved