OFFSET
1,1
COMMENTS
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..600
MAPLE
isA025475not1 := proc(n)
if n <= 1 then
false;
elif isprime(n) then
false;
elif nops(numtheory[factorset](n)) = 1 then
true;
else
false;
end if;
end proc:
A025475not1 := proc(n)
option remember;
local a;
if n = 1 then
4;
else
for a from procname(n-1)+1 do
if isA025475not1(a) then
return a;
end if;
end do:
end if;
end proc:
A225104w := proc(n)
local a, i, ppi, ppj ;
a := 0 ;
for i from 1 do
ppi := A025475not1(i) ;
if ppi >= n/2 then
break;
end if;
ppj := n-ppi ;
if isA025475not1(ppj) then
a := a+1 ;
end if;
end do:
a ;
end proc:
for n from 1 do
if A225104w(n) >= 3 then
print(n) ;
end if;
end do: # R. J. Mathar, Jun 13 2013
MATHEMATICA
nn = 36490; p = Sort[Flatten[Table[Prime[n]^i, {n, PrimePi[Sqrt[nn]]}, {i, 2, Log[Prime[n], nn]}]]]; Transpose[Sort[Select[Tally[Flatten[Table[p[[i]] + p[[j]], {i, Length[p] - 1}, {j, i + 1, Length[p]}]]], #[[1]] <= nn && #[[2]] > 2 &]]][[1]] (* T. D. Noe, Apr 29 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alex Ratushnyak, Apr 28 2013
STATUS
approved
