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

Author

Paul D. Hanna, May 09 2003

Keywords

• nonn
• nice

Comments

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]

Links

• 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.

Crossrefs

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.

Programs

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

PARI

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

Extensions

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