OFFSET
1,1
COMMENTS
34 is the smallest term in common with A108605.
LINKS
J.W.L. (Jan) Eerland, Table of n, a(n) for n = 1..10000
EXAMPLE
34 = 2*17 (semiprime) and 2 + 17 = 19 is prime.
MATHEMATICA
A108606=Select[Range[1000], Plus@@(Transpose[FactorInteger[ # ]])[[2]]==2&& PrimeQ[Plus@@IntegerDigits[ # ]]&]
DeleteCases[ParallelTable[If[PrimeOmega[n]==2&&PrimeQ[Total[IntegerDigits[n]]], n, a], {n, 0, 126181}], a] (* J.W.L. (Jan) Eerland, Dec 21 2021 *)
PROG
(PARI) select(isA108606(n)={bigomega(n)==2&&isprime(sumdigits(n))}, [1..1000]) \\ J.W.L. (Jan) Eerland, Dec 23 2021
(Python)
from sympy import isprime, factorint
def ok(n): return isprime(sum(map(int, str(n)))) and sum(factorint(n).values()) == 2
print([k for k in range(455) if ok(k)]) # Michael S. Branicky, Aug 22 2022
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Zak Seidov, Jun 12 2005
STATUS
approved