editing
approved
editing
approved
(PARI) is(n)=if(!isprime((n+3)\4), return(0)); for(k=n, n+3, if(bigomega(k)!=3, return(0))); 1 \\ Charles R Greathouse IV, Feb 05 2017
(PARI) list(lim)=my(v=List(), u=v, t); forprime(p=2, lim\4, forprime(q=2, min(lim\(2*p), p), t=p*q; forprime(r=2, min(lim\t, q), listput(u, t*r)))); u=Set(u); for(i=4, #u, if(u[i]-u[i-3]==3, listput(v, u[i-3]))); Vec(v) \\ Charles R Greathouse IV, Feb 05 2017
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed
okQ[{a_, b_, c_, d_}]:=Union[{a, b, c, d}]=={3}; Flatten[Position[Partition[ PrimeOmega[ Range[11000]], 4, 1], _?(okQ)]] (* Harvey P. Dale, Sep 23 2012 *)
A subset of A045940. - _Zak Seidov, _, Nov 05 2006
with(numtheory): a:=proc(n) if bigomega(n)=3 and bigomega(n+1)=3 and bigomega(n+2)=3 and bigomega(n+3)=3 then n else fi end: seq(a(n), n=1..15000); - _# _Emeric Deutsch_, Nov 07 2006
More terms from _Zak Seidov, _, Nov 05 2006
proposed
editing
editing
proposed
A subset of A045940. - Zak Seidov, Nov 05 2006 - Zak Seidov, Nov 05 2006
approved
editing
editing
approved