cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A063640 Primes of form p*q*r + 1, where p, q and r are primes.

This page as a plain text file.
%I A063640 #24 Mar 08 2023 14:59:29
%S A063640 13,19,29,31,43,53,67,71,79,103,131,139,149,173,191,223,239,269,283,
%T A063640 293,311,317,367,389,419,431,439,443,499,509,557,599,607,619,643,647,
%U A063640 653,659,683,743,773,787,797,823,827,907,947,971,1031,1039,1087,1091
%N A063640 Primes of form p*q*r + 1, where p, q and r are primes.
%C A063640 Odd primes p such that (p-1)/2 is a semiprime. - _Robert G. Wilson v_, Sep 01 2007
%H A063640 Harry J. Smith, <a href="/A063640/b063640.txt">Table of n, a(n) for n = 1..1000</a>
%p A063640 q:= n-> isprime(n) and numtheory[bigomega](n-1)=3:
%p A063640 select(q, [$2..1100])[];  # _Alois P. Heinz_, Mar 08 2023
%t A063640 Take[ Select[ Union@ Flatten@ Table[ Prime@p*Prime@q*Prime@r + 1, {p, 48}, {q, p}, {r, q}], PrimeQ@ # &], 53] (* Or *)
%t A063640 semiPrimeQ[x_] := Plus @@ Last /@ FactorInteger[x] == 2; Select[Prime@ Range@ 182, semiPrimeQ[(# - 1)/2] &] (* _Robert G. Wilson v_, Sep 01 2007 *)
%t A063640 2#+1&/@Select[Table[(n-1)/2,{n,Prime[Range[200]]}],PrimeOmega[#]==2&] (* _Harvey P. Dale_, Oct 11 2018 *)
%o A063640 (PARI)  n=0; for (m=1, 10^9, p=prime(m); if (bigomega(p - 1) == 3, write("b063640.txt", n++, " ", p); if (n==1000, break)) )  \\ _Harry J. Smith_, Aug 26 2009
%Y A063640 A090866 is a subsequence.
%Y A063640 Cf. A014612, A063639.
%K A063640 nonn
%O A063640 1,1
%A A063640 _Reinhard Zumkeller_, Jul 21 2001