login

Reminder: The OEIS is hiring a new managing editor, and the application deadline is January 26.

A091207
Indices of primes that are also irreducible GF(2)[X]-polynomials.
3
1, 2, 4, 5, 6, 8, 11, 12, 13, 15, 17, 18, 19, 21, 25, 27, 29, 32, 33, 37, 39, 43, 44, 47, 50, 52, 53, 61, 65, 75, 78, 81, 84, 90, 93, 95, 102, 103, 107, 110, 111, 112, 113, 115, 118, 121, 123, 126, 128, 134, 135, 136, 138, 144, 149, 151, 153, 156, 158, 162, 163, 164
OFFSET
1,2
FORMULA
a(n) = A049084(A091206(n)).
MAPLE
filter:= proc(p) local L, i, x;
L:= convert(p, base, 2);
Irreduc(add(L[i]*x^(i-1), i=1..nops(L))) mod 2
end proc:
select(t -> filter(ithprime(t)), [$1..300]); # Robert Israel, Jun 12 2018
MATHEMATICA
okQ[p_] := PrimeQ[p] && Module[{id, pol}, id = IntegerDigits[p, 2] // Reverse; pol = id.x^Range[0, Length[id] - 1]; IrreduciblePolynomialQ[pol, Modulus -> 2]];
Select[Range[300], okQ[Prime[#]]&] (* Jean-François Alcover, Feb 07 2023 *)
PROG
(PARI) isok(n) = polisirreducible(Mod(1, 2)*Pol(binary(prime(n)))); \\ Michel Marcus, Jun 13 2018
CROSSREFS
Complement of A091210.
Sequence in context: A099747 A249053 A082851 * A284525 A353187 A099247
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 03 2004
STATUS
approved