cp's OEIS Frontend

A207187 Numbers matching polynomials y(k,x) that have x+1 as a factor; see Comments.

%I A207187 #6 Mar 30 2012 18:58:12
%S A207187 3,5,9,13,19,22,25,27,30,33,39,43,49,52,55,59,65,68,71,83,89,92,95,
%T A207187 101,104,107,110,113,116,119,121,124,127,133,136,139,142,145,148,151,
%U A207187 157,169,172,175,181,185,191,194,197,209,215,218,221,224,227,230
%N A207187 Numbers matching polynomials y(k,x) that have x+1 as a factor; see Comments.
%C A207187 The polynomials y(k,x) range through all monic polynomials with coefficients in {-1,0,1}, ordered as at A206821.
%e A207187 The first 13 polynomials:
%e A207187 1 .... 1
%e A207187 2 .... x
%e A207187 3 .... x + 1
%e A207187 4 .... x^2
%e A207187 5 .... x^2 - 1
%e A207187 6 .... x^2 - x
%e A207187 7 .... x^2 - x - 1
%e A207187 8 .... x^2 + 1
%e A207187 9 .... x^2 + x
%e A207187 10 ... x^2 + x + 1
%e A207187 11 ... x^3
%e A207187 12 ... x^3 - 1
%e A207187 13 ... x^3 - x
%e A207187 The list exemplifies these sequences:
%e A207187 A207187 (multiples of x + 1): 3,5,9,13,...
%e A207187 A207188 (multiples of x): 2,4,6,9,11,13,...
%e A207187 A207189 (multiples of x - 1): 5,6,12,13,...
%e A207187 A207190 (multiples of x^2 + 1): 8,20,25,27,...
%t A207187 t = Table[IntegerDigits[n, 2], {n, 1, 2000}];
%t A207187 b[n_] := Reverse[Table[x^k, {k, 0, n}]]
%t A207187 p[n_] := p[n] = t[[n]].b[-1 + Length[t[[n]]]]
%t A207187 TableForm[Table[{n, p[n], Factor[p[n]]}, {n, 1, 6}]]
%t A207187 f[k_] := 2^k - k; g[k_] := 2^k - 2 + f[k - 1];
%t A207187 q1[n_] := p[2^(k - 1)] - p[n + 1 - f[k]]
%t A207187 q2[n_] := p[n - f[k] + 2]
%t A207187 y1 = Table[p[n], {n, 1, 4}];
%t A207187 Do[AppendTo[y1,
%t A207187   Join[Table[q1[n], {n, f[k], g[k] - 1}],
%t A207187    Table[q2[n], {n, g[k], f[k + 1] - 1}]]], {k, 3, 10}]
%t A207187 y = Flatten[y1]; (* monic polynomials over {-1,0,1} *)
%t A207187 TableForm[Table[{n, y[[n]], Factor[y[[n]]]}, {n, 1, 10}]]
%t A207187 Table[y[[n]] /. x -> -1, {n, 1, 300}];
%t A207187 Flatten[Position[%, 0]]  (* A207187 *)
%t A207187 Table[y[[n]] /. x -> 0, {n, 1, 300}] ;
%t A207187 Flatten[Position[%, 0]]  (* A207188 *)
%t A207187 Table[y[[n]] /. x -> 1, {n, 1, 1200}] ;
%t A207187 Flatten[Position[%, 0]]  (* A207189 *)
%t A207187 Table[y[[n]] /. x -> I, {n, 1, 600}] ;
%t A207187 Flatten[Position[%, 0]]  (* A207190 *)
%Y A207187 Cf. A206821.
%K A207187 nonn
%O A207187 1,1
%A A207187 _Clark Kimberling_, Feb 16 2012