cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A191583 Sum of the distinct prime divisors of prime(n) + prime(n+1).

This page as a plain text file.
%I A191583 #15 Nov 21 2013 12:50:18
%S A191583 5,2,5,5,5,10,5,12,15,10,19,18,12,10,7,9,10,2,28,5,21,5,45,36,16,22,
%T A191583 17,5,42,10,48,69,28,5,10,20,7,21,24,13,10,36,5,23,16,48,40,10,24,23,
%U A191583 61,10,46,129,20,28,10,139,36,52,5,10,108,18,17,5,169,24
%N A191583 Sum of the distinct prime divisors of prime(n) + prime(n+1).
%C A191583 a(n) = A008472(A001043(n)). [_Reinhard Zumkeller_, Jun 28 2011]
%H A191583 Reinhard Zumkeller, <a href="/A191583/b191583.txt">Table of n, a(n) for n = 1..10000</a>
%e A191583 a(6) = 10 because prime(6) + prime(7) = 13+17 = 30 = 2*3*5 and 2+3+5 = 10.
%p A191583 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:
%t A191583 sdpd[n_]:=Total[Transpose[FactorInteger[n]][[1]]]; sdpd/@(Total/@ Partition[ Prime[Range[70]],2,1]) (* _Harvey P. Dale_, Mar 18 2012 *)
%o A191583 (PARI) vecsum(v)=sum(i=1,#v,v[i])
%o A191583 p=2;forprime(q=3,1e3,print1(vecsum(factor(p+q)[,1])", ");p=q)
%o A191583 \\ _Charles R Greathouse IV_, Jun 12 2011
%K A191583 nonn
%O A191583 1,1
%A A191583 _Michel Lagneau_, Jun 07 2011