OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
38 is in the sequence because 36 = 2^2 * 3^2 and 40 = 2^3 * 5.
MAPLE
with(numtheory): a:=proc(n) if nops(factorset(n-2))=nops(factorset(n+2)) then n else fi end: seq(a(2*n), n=2..133); # Emeric Deutsch, Mar 12 2006
MATHEMATICA
Select[Range[2, 262, 2], PrimeNu[# - 2] == PrimeNu[# + 2] &] (* Amiram Eldar, Feb 18 2020 *)
Select[Mean/@SequencePosition[PrimeNu[Range[300]], {x_, _, _, _, x_}], EvenQ] (* Harvey P. Dale, Oct 11 2023 *)
PROG
(PARI) g(n) = forstep(x=4, n, 2, p1=omega(x-2); p2=omega(x+2); if(p1==p2, print(x", ")))
(Magma) [k:k in [4.. 270 by 2]| #PrimeDivisors(k-2) eq #PrimeDivisors(k+2)]; // Marius A. Burtea, Feb 18 2020
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Cino Hilliard, Mar 03 2006
STATUS
approved