OFFSET
1,1
COMMENTS
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
EXAMPLE
a(6) = 10 because prime(6) + prime(7) = 13+17 = 30 = 2*3*5 and 2+3+5 = 10.
MAPLE
with(numtheory):for n from 1 to 100 do:x:=ithprime(n)+ithprime(n+1):y:=factorset(x):n1:=nops(y):s:=0:for k from 1 to n1 do:s:=s+y[k]:od:printf(`%d, `, s):od:
MATHEMATICA
sdpd[n_]:=Total[Transpose[FactorInteger[n]][[1]]]; sdpd/@(Total/@ Partition[ Prime[Range[70]], 2, 1]) (* Harvey P. Dale, Mar 18 2012 *)
PROG
(PARI) vecsum(v)=sum(i=1, #v, v[i])
p=2; forprime(q=3, 1e3, print1(vecsum(factor(p+q)[, 1])", "); p=q)
\\ Charles R Greathouse IV, Jun 12 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Jun 07 2011
STATUS
approved