OFFSET
0,8
COMMENTS
Limit a(n)/a(n+1) -> r = -0.530852489019085 where A(r)=0.
Let A_n(x) be the power series formed from the first n terms of this sequence. Then a(0) = 1, a(n) = floor(1 - [x^n] (A_(n-1)(x))^2/2). Replacing 2 with a larger integer k generates the related sequences A084203-A084212. - Charlie Neder, Jan 16 2019
LINKS
Paul D. Hanna, Table of n, a(n) for n = 0..1024
N. Heninger, E. M. Rains and N. J. A. Sloane, On the Integrality of n-th Roots of Generating Functions, arXiv:math/0509316 [math.NT], 2005-2006; J. Combinatorial Theory, Series A, 113 (2006), 1732-1745.
MATHEMATICA
a[n_] := a[n] = Block[{s = Sum[a[i]*x^i, {i, 0, n - 1}]}, If[ IntegerQ@ Last@ CoefficientList[ Series[ Sqrt[s + x^n], {x, 0, n}], x], 1, 2]]; Table[a[n], {n, 0, 42}]; CoefficientList[ Series[ Sqrt[ Sum[ a[i]*x^i, {i, 0, 42}]], {x, 0, 42}], x] (* Robert G. Wilson v, Nov 11 2007 *)
PROG
(PARI) /* Using Charlie Neder's formula */
{a(n) = my(A=[1]); for(i=0, n, A=concat(A, 0); A[#A] = floor(1 - polcoeff( Ser(A)^2, #A-1)/2) ); A[n+1]}
for(n=0, 50, print1(a(n), ", ")) \\ Paul D. Hanna, Jan 17 2019
CROSSREFS
KEYWORD
sign
AUTHOR
Paul D. Hanna, May 19 2003
STATUS
approved