OFFSET
1,1
COMMENTS
Except term 7, same as A181602.
MAPLE
isA181669 := proc(n::integer)
if isprime(n) then
if modp(n, 6) = 5 then
if numtheory[bigomega](n-1) =2 then
if isprime(n+2) or issqr(n+2) and isprime(sqrt(n+2)) then
true ;
else
false;
end if;
else
false;
end if;
else
false;
end if;
else
false;
end if;
end proc:
for n from 2 to 100000 do
if isA181669(n) then
printf("%d, ", n) ;
end if;
end do: # R. J. Mathar, Oct 17 2025
MATHEMATICA
Select[6*Range[3000]-1, PrimeQ[#]&&PrimeOmega[#-1]==2&&AnyTrue[ {#+2, Sqrt[ #+2]}, PrimeQ]&] (* Harvey P. Dale, Jul 01 2022 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Giovanni Teofilatto, Nov 18 2010
STATUS
approved
