OFFSET
1,1
COMMENTS
There are 670 semiprimes of form prime-1 below 10^5.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
FORMULA
MAPLE
q:= n-> (n::even) and andmap(isprime, [n+1, n/2]):
select(q, [$1..5000])[]; # Alois P. Heinz, Jul 19 2023
MATHEMATICA
Select[Range[6000], Plus@@Last/@FactorInteger[#]==2&&PrimeQ[#+1]&] (* Vladimir Joseph Stephan Orlovsky, May 08 2011 *)
Select[Range[3000], PrimeOmega[#]==2&&PrimeQ[#+1]&] (* Harvey P. Dale, Oct 16 2012 *)
Select[ Prime@ Range@ 430 - 1, PrimeOmega@# == 2 &] (* Robert G. Wilson v, Feb 18 2014 *)
PROG
(Magma) IsSemiprime:=func<n | &+[d[2]: d in Factorization(n)] eq 2>; [s: n in [2..500] | IsSemiprime(s) where s is NthPrime(n)-1]; // Vincenzo Librandi, Oct 17 2012
(Haskell)
a077065 n = a077065_list !! (n-1)
a077065_list = filter ((== 1) . a010051' . (`div` 2)) a006093_list
-- Reinhard Zumkeller, Nov 22 2013, Oct 27 2012
(PARI) [x-1|x<-primes(10^4), bigomega(x-1)==2] \\ Charles R Greathouse IV, Nov 22 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Oct 23 2002
STATUS
approved