OFFSET
1,1
COMMENTS
EXAMPLE
(See the Example section of A207669.)
MATHEMATICA
t = Table[IntegerDigits[n, 3], {n, 1, 1000}];
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], Modulus -> 3],
AppendTo[u, n]], {n, 1, 400}]
u (* A207669 *)
Complement[Range[200], %] (* A207670 *)
b[n_] := FromDigits[IntegerDigits[u, 3][[n]]]
Table[b[n], {n, 1, 50}] (* A207671 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Clark Kimberling, Feb 26 2012
STATUS
approved