A108335 -id:A108335 - cp's OEIS frontend

A083952 Integer coefficients a(n) of A(x), where a(n) = 1 or 2 for all n, such that A(x)^(1/2) has only integer coefficients.

1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2

Offset: 0

Paul D. Hanna, May 09 2003

More generally, the sequence "integer coefficients of A(x), where 1<=a(n)<=m, such that A(x)^(1/m) consists entirely of integer coefficients", appears to have a unique solution for all m. [That is true - see Theorem 17 of Heninger-Rains-Sloane (2006). - N. J. A. Sloane, Aug 27 2015]

Is this sequence periodic? [It is not periodic for m = 2 or 3. Larger cases remain open. - N. J. A. Sloane, Aug 27 2015]

Robert G. Wilson v, Table of n, a(n) for n = 0..10000
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.

Cf. A084202 (A(x)^(1/2)), A108335 (A084202 mod 4), A108336 (A084202 mod 2), A108340 (a(n) mod 2). Positions of 1's: A108783.

Cf. A083953, A083954, A083945, A083946.

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, 104}] (* Robert G. Wilson v, Nov 25 2006 *)
s = 0; a[n_] := a[n] = Block[{}, If[IntegerQ@ Last@ CoefficientList[ Series[ Sqrt[s + x^n], {x, 0, n}], x], s = s + x^n; 1, s = s + 2 x^n; 2]]; Table[ a@n, {n, 0, 104}] (* Robert G. Wilson v, Sep 08 2007 *)

A083952_upto(N=99)=vector(N+1, n, if(n>1, (denominator(polcoeff(sqrt(O(x^n)+N+=x^(n-1)),n-1))>1 && N+=x^(n-1))+1, N=1)) \\ M. F. Hasler, Jan 27 2025

More terms from N. J. A. Sloane, Jul 02 2005

A108336 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).

Original entry on oeis.org

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, 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, 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

Offset: 0

N. J. A. Sloane and Nadia Heninger, Jul 02 2005

Equals A084202 read mod 2.

Robert Israel, Table of n, a(n) for n = 0..10000
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.
N. Heninger, E. M. Rains and N. J. A. Sloane, On the Integrality of n-th Roots of Generating Functions, J. Combinatorial Theory, Series A, 113 (2006), 1732-1745.

Cf. A083952, A084202, A108335, A108337, A108340.

S:= 0: SS:= 0:
for i from 0 to 100 do
  s:= coeff(SS,x,i);
  if s = 0 or s = 3 then
     SS:= SS + 2*expand(S*x^i)+x^(2*i) mod 4; S:= S + x^i;
  fi
od:
seq(coeff(S,x,i),i=0..100); # Robert Israel, May 14 2019

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 *)

cp's OEIS Frontend

A083952 Integer coefficients a(n) of A(x), where a(n) = 1 or 2 for all n, such that A(x)^(1/2) has only integer coefficients.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Extensions