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}]
Complement[Range[200], %] (*
A207670 *)
b[n_] := FromDigits[IntegerDigits[u, 3][[n]]]
Table[b[n], {n, 1, 50}] (*
A207671 *)