OFFSET
1,1
COMMENTS
Numbers k such that 2k - 1 and 4k^2 + 2k + 1 are both prime.
LINKS
K. D. Bajpai, Table of n, a(n) for n = 1..9443
EXAMPLE
From K. D. Bajpai, Nov 16 2019: (Start)
a(3) = 6 is a term because (2*6)^3 - 1 = 1727 = 11*157, which is a semiprime.
a(4) = 7 is a term because (2*7)^3 - 1 = 2743 = 13*211, which is a semiprime.
9 is not in the sequence because (2*9)^3 - 1 = 5831 = 7*7*7*17, which is not semiprime.
(End)
MAPLE
issp:= n-> not isprime(n) and numtheory[bigomega](n)=2:
select( n-> issp((2*n)^3-1), [seq(n, n=1..200)]); # K. D. Bajpai, Nov 16 2019
MATHEMATICA
Select[Range@ 750, PrimeOmega[(2 #)^3 - 1] == 2 &] (* Michael De Vlieger, Aug 02 2018 *)
PROG
(PARI) for(k=1, 500, if(bigomega((2*k)^3-1)==2, print1(k, ", ")))
(Magma) IsSemiprime:=func<i | &+[d[2]: d in Factorization(i)] eq 2>; [n: n in [2..800] | IsSemiprime(s) where s is (2*n)^3-1]; // Vincenzo Librandi, Aug 04 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Jianing Song, Aug 01 2018
STATUS
approved