OFFSET
1,2
COMMENTS
Most of the elements of this sequence are prime. The "pseudoprimes" of these sequence are part of A244626.
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Jonas Kaiser, On the relationship between the Collatz conjecture and Mersenne prime numbers, arXiv:1608.00862 [math.GM], 2016.
MATHEMATICA
Select[Range[1, 6001, 6], # == 1 || PowerMod[2, (#-1)/3, #] == 1 && Mod[-PowerMod[2, (#-1)/6, #], #] == 1&] (* Jean-François Alcover, Nov 18 2018 *)
PROG
(PARI) is(n)=n%6==1 && Mod(2, n)^(n\3)==1 && (2*n-1)*Mod(2, n)^(n\6)==1 \\ Charles R Greathouse IV, Nov 08 2017
CROSSREFS
KEYWORD
nonn
AUTHOR
Jonas Kaiser, Nov 07 2017
STATUS
approved