OFFSET
1,1
MATHEMATICA
s = {4, 6, 9}; Do[a = s[[-1]] + s[[-2]] + s[[-3]]; n = 10; While[MemberQ[s, n] || 2 != PrimeOmega[n] || 2 != PrimeOmega[a + n], n++]; AppendTo[s, n], {120}]; s
PROG
(PARI) issp(k) = bigomega(k)==2; \\ A001358
lista(nn) = my(va = vector(nn)); va[1]=4; va[2]=6; va[3]=9; my(vs = vecsort(va)); my(s=sum(k=1, 3, va[k])); for (n=4, nn, my(k=1); while (!(issp(k) && issp(k+s) && !vecsearch(vs, k)), k++); va[n]=k; vs = vecsort(va); s += k - va[n-3]; ); va; \\ Michel Marcus, Aug 04 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Jun 13 2022
STATUS
approved