OFFSET
1,2
COMMENTS
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..21951
FORMULA
Other identities and observations. For all n >= 1:
A255573(a(n)) = n.
EXAMPLE
The reducible polynomials matching the first four terms:
1 = 1(base 2) matches 1
4 = 100(base 2) matches x^2
6 = 110(base 2) matches x^2 + x
8 = 1000(base 2) matches x^3
9 = 1001(base 2) matches x^3 + 1
MATHEMATICA
t = Table[IntegerDigits[n, 2], {n, 1, 850}];
b[n_] := Reverse[Table[x^k, {k, 0, n}]]
p[n_, x_] := t[[n]].b[-1 + Length[t[[n]]]]
Table[p[n, x], {n, 1, 15}]
u = {}; Do[n++; If[IrreduciblePolynomialQ[p[n, x]],
AppendTo[u, n]], {n, 300}];
u (* A206074 *)
Complement[Range[200], u] (* A205783 *)
b[n_] := FromDigits[IntegerDigits[u, 2][[n]]]
Table[b[n], {n, 1, 40}] (* A206073 *)
PROG
(PARI)
isA205783(n) = ((n > 0) && !polisirreducible(Pol(binary(n))));
n = 0; i = 0; while(n < 32768, n++; if(isA205783(n), i++; write("b205783.txt", i, " ", n)));
CROSSREFS
After 1 a subsequence of A091212 (69 is the first term missing from here).
KEYWORD
nonn
AUTHOR
Clark Kimberling, Feb 03 2012
STATUS
approved