OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
Prime(4) = 7, prime(5) = 11. Sum of prime factors of 8 is 2+2+2 = 6, sum of prime factors of 9 is 3+3 = 6, sum of prime factors of 10 is 2+5= 7; 6+6+7 = 19 is prime, hence prime(4) = 7 is a term.
Prime(19) = 67, prime(20) = 71. Sum of prime factors of 68, 69, 70 is resp. 2+2+17 = 21, 3+23 = 26, 2+5+7 = 14; 21+26+14 = 61 is prime, hence prime(19) = 67 is a term.
Prime(26) = 101, prime(27) = 103. Sum of prime factors of 102 = 2*3*17 is 22, which is composite. Hence prime(26) = 101 is not in the sequence.
MAPLE
p:= 2: count:= 0: R:= NULL:
while count < 100 do
q:= nextprime(p);
v:= add(add(t[1]*t[2], t=ifactors(m)[2]), m=p+1..q-1);
if isprime(v) then count:= count+1; R:= R, p; fi;
p:= q;
od:
R; # Robert Israel, Oct 08 2020
MATHEMATICA
spfQ[n_]:=Module[{strt=First[n]+1, end=Last[n]-1}, PrimeQ[Total[Times@@@ Flatten[ FactorInteger[ Range[ strt, end]], 1]]]]; Transpose[Select[ Partition[ Prime[ Range[300]], 2, 1], spfQ]][[1]] (* Harvey P. Dale, May 06 2011 *)
PROG
(Magma) [ p: p in [ NthPrime(k): k in [2..258] ] | IsPrime(&+[ &+[ k[1]*k[2]: k in Factorization(c) ]: c in [p+1..NextPrime(p)-1] ] ) ]; /* Klaus Brockhaus, Mar 29 2007 */
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot, Mar 27 2007
EXTENSIONS
Edited, corrected and extended by Klaus Brockhaus, Mar 29 2007
STATUS
approved