%I #18 May 15 2019 04:49:42
%S 1,1,0,1,0,1,1,0,0,0,0,0,0,1,0,1,1,0,1,0,0,0,0,0,0,1,1,0,1,0,1,1,0,0,
%T 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,1,0,
%U 1,0,0,0,0,0,0,1,1,0,0,0,1,0,1,1,0,1,1,1,1,0,1,0,1,0,0,0,1,1,0
%N Unique sequence of 1's and 0's such that (Sum_{n >= 0} a(n)*x^n)^2 mod 4 has coefficients which are all 1's and 2's (A083952).
%C Equals A084202 read mod 2.
%H Robert Israel, <a href="/A108336/b108336.txt">Table of n, a(n) for n = 0..10000</a>
%H N. Heninger, E. M. Rains and N. J. A. Sloane, <a href="https://arxiv.org/abs/math/0509316">On the Integrality of n-th Roots of Generating Functions</a>, arXiv:math/0509316 [math.NT], 2005-2006.
%H N. Heninger, E. M. Rains and N. J. A. Sloane, <a href="https://doi.org/10.1016/j.jcta.2006.03.018">On the Integrality of n-th Roots of Generating Functions</a>, J. Combinatorial Theory, Series A, 113 (2006), 1732-1745.
%p S:= 0: SS:= 0:
%p for i from 0 to 100 do
%p s:= coeff(SS,x,i);
%p if s = 0 or s = 3 then
%p SS:= SS + 2*expand(S*x^i)+x^(2*i) mod 4; S:= S + x^i;
%p fi
%p od:
%p seq(coeff(S,x,i),i=0..100); # _Robert Israel_, May 14 2019
%t max = 98; (* a = A084202 *) 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, max}]; A108336 = CoefficientList[ Series[Sqrt[Sum[a[i]*x^i, {i, 0, max}]], {x, 0, max}], x] // Mod[#, 2]& (* _Jean-François Alcover_, Apr 01 2016, after _Robert G. Wilson v_ *)
%Y Cf. A083952, A084202, A108335, A108337, A108340.
%K nonn,easy,nice
%O 0,1
%A _N. J. A. Sloane_ and _Nadia Heninger_, Jul 02 2005
|