OFFSET
1,2
COMMENTS
Numbers k that divide A119564(k).
A term is prime if and only if it is in A023394.
If k is in A307843, then it is a term of this sequence.
The terms that are not divisors of Fermat numbers are 1605, 4369, 32113, 94945, ... (they are all composite). Are there infinitely many of them?
Note that 2^(2^k) - 2^k + 1 = (2^(2^k) - 1) - (2^k - 2).
MATHEMATICA
q[k_] := Mod[PowerMod[2, 2^k, k] - PowerMod[2, k, k] + 1, k] == 0; Select[Range[1, 10^5, 2], q] (* Amiram Eldar, Jun 10 2024 *)
PROG
(PARI) isok(k) = Mod(Mod(2, k)^(2^k) - Mod(2, k)^k + 1, k) == 0; \\ Michel Marcus, Jun 12 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Thomas Ordowski, Jun 10 2024
EXTENSIONS
More terms from Amiram Eldar, Jun 10 2024
STATUS
approved