OFFSET
1,2
COMMENTS
Primes are counted with multiplicity.
Closed under multiplication.
MAPLE
isA202237 := proc(n)
if type(n, 'odd') then
else
false;
end if;
end proc:
for n from 1 to 200 do
if isA202237(n) then
printf("%d, ", n);
end if;
end do: # R. J. Mathar, Dec 16 2011
MATHEMATICA
fQ[n_]:=Plus@@((Mod[#[[1]], 4]-2)*#[[2]]&/@If[==1, {}, FactorInteger[n]]==0 && OddQ[n]; Select[Range[600], fQ] (* Ray Chandler, Dec 20 2011 *)
PROG
(PARI) netprime(n)=local(fm=factor(n)); sum(k=1, matsize(fm)[1], if(fm[k, 1]==2, 0, if(fm[k, 1]%4==1, fm[k, 2], -fm[k, 2])))
ap(n)=forstep(k=1, n, 2, if(netprime(k)==0, print1(k", ")))
CROSSREFS
KEYWORD
nonn
AUTHOR
Franklin T. Adams-Watters, Dec 16 2011
STATUS
approved