OFFSET
1,1
EXAMPLE
For n = 8, a(8) is 13, the sum of 2 and 11 which are the distinct prime factors of 44, given by a(8) in A393096.
MATHEMATICA
sopf[k_]:=Total[First/@FactorInteger[k]]; q[k_]:=PrimeQ[sopf[k]]&&PrimeQ[k-sopf[k]]; sopf/@Select[Range[2850], q] (* James C. McMahon, Mar 16 2026 *)
PROG
(Python)
from sympy import isprime, primefactors
def ok(n):return isprime(sum(primefactors(n))) and isprime(n-sum(primefactors(n)))
print(list(map(lambda n: sum(primefactors(n)), filter(ok, range(1, 2000)))))
(PARI) sopf(k) = vecsum(factor(k)[, 1]);
isok(k) = my(s=sopf(k)); isprime(s) && isprime(k-s);
apply(sopf, select(isok, [1..4000])) \\ Michel Marcus, Mar 09 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Guy Siviour, Mar 09 2026
STATUS
approved
