login
Positions of 1's in A083952, where A083952 gives the 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.
4

%I #15 Feb 10 2024 09:24:19

%S 0,2,6,10,12,26,30,32,36,50,52,56,60,62,126,130,132,136,150,152,160,

%T 164,166,170,172,174,176,180,184,192,194,198,200,202,214,216,220,226,

%U 228,230,234,236,240,242,244,260,262,264,272,274,278,282,286

%N Positions of 1's in A083952, where A083952 gives the 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.

%H Robert G. Wilson v, <a href="/A108783/b108783.txt">Table of n, a(n) for n = 1..1356</a>

%t 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]]; Union@ Table[ If[ a[n] == 1, n, 0], {n, 0, 300}] (* _Robert G. Wilson v_, Nov 25 2006 *)

%Y Cf. A083952, A108337, A108338. See A111363 for another version.

%K nonn

%O 1,2

%A _N. J. A. Sloane_, following a suggestion from _Paul D. Hanna_, Jun 30 2005