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.

A096163 Primes p of the form qrs + 1 where q, r and s are distinct primes.

This page as a plain text file.
%I A096163 #4 Jun 06 2021 11:38:27
%S A096163 31,43,67,71,79,103,131,139,191,223,239,283,311,367,419,431,439,443,
%T A096163 499,599,607,619,643,647,659,683,743,787,823,827,907,947,971,1031,
%U A096163 1039,1087,1091,1103,1163,1223,1259,1399,1427,1447,1499,1511,1543,1559,1571,1579
%N A096163 Primes p of the form qrs + 1 where q, r and s are distinct primes.
%C A096163 Each composite number qrs = a(n)-1 is a squarefree 3-almost prime. This sequence is a subsequence of A078330 which, besides having 3 as its first term, first differs by including 2311 = 2*3*5*7*11 + 1 (a squarefree 5-almost prime plus 1).
%t A096163 With[{nn=50},Take[Union[Select[Times@@@Subsets[Prime[Range[2nn]],{3}]+1,PrimeQ]],nn]] (* _Harvey P. Dale_, Jun 06 2021 *)
%o A096163 (PARI) /* Here are five equivalent PARI programs */ forprime(p=2,2400, if(moebius(p-1)==-1 && omega(p-1)==3, print1(p,","))) forprime(p=2,2400, if(moebius(p-1)==-1 && bigomega(p-1)==3, print1(p,","))) forprime(p=2,2400, if(bigomega(p-1)==3 && omega(p-1)==3, print1(p,","))) forprime(p=2,2400, if(omega(p-1)==3 && issquarefree(p-1), print1(p,","))) forprime(p=2,2400, if(bigomega(p-1)==3 && issquarefree(p-1), print1(p,",")))
%Y A096163 Cf. A078330 (primes p with mu(p-1) = -1).
%K A096163 nonn
%O A096163 1,1
%A A096163 _Rick L. Shepherd_, Jun 18 2004