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.

Showing 1-2 of 2 results.

A116994 Prime partial sums of triangular numbers with prime indices.

Original entry on oeis.org

3, 1759, 3323, 469303, 605113, 641969, 1110587, 1426669, 11148289, 18352349, 20473721, 21820391, 24710753, 30048589, 36690923, 40785301, 97060681, 155135369, 160593239, 168132247, 361391623, 377965069, 416572171, 645803201
Offset: 1

Views

Author

Jonathan Vos Post, Apr 02 2006

Keywords

Examples

			a(1) = Sum_{i=1..1} prime(i)*(prime(i)+1)/2 = T(2) = 3.
a(2) = Sum_{i=1..11} prime(i)*(prime(i)+1)/2 = T(2)+T(3)+T(5)+T(7)+T(11)+T(13)+T(17)+T(19)+T(23)+T(29)+T(31) = 1759.
a(3) = Sum_{i=1..13} prime(i)*(prime(i)+1)/2 = 3323.
a(4) = Sum_{i=1..53} prime(i)*(prime(i)+1)/2 = T(2) + ... + T(241) = 469303.
a(5) = Sum_{i=1..57} prime(i)*(prime(i)+1)/2 = T(2) + ... + T(269) = 605113.
a(6) = Sum_{i=1..58} prime(i)*(prime(i)+1)/2 = T(2) + ... + T(271) = 641969.
a(7) = Sum_{i=1..68} prime(i)*(prime(i)+1)/2 = T(2) + ... + T(337) = 1110587.
		

Crossrefs

Programs

  • Maple
    T:=n->n*(n+1)/2: a:=proc(n): if isprime(sum(T(ithprime(j)),j=1..n))=true then sum(T(ithprime(j)),j=1..n) else fi end: seq(a(n),n=1..500); # Emeric Deutsch, Apr 06 2006
  • Mathematica
    Select[Accumulate[Table[(n(n+1))/2,{n,Prime[Range[500]]}]],PrimeQ] (* Harvey P. Dale, Jan 25 2015 *)

Formula

A000040 INTERSECTION A085739. Primes in A085739.

Extensions

More terms from Emeric Deutsch, Apr 06 2006

A116911 Prime partial sums of pentagonal numbers with prime indices.

Original entry on oeis.org

5, 17, 4957, 129277, 2826443, 3861083, 5126483, 9451573, 19811083, 53751743, 68136617, 98729003, 264616831, 388771421, 498157871, 608312141, 682548511, 779346653, 918754301, 1174179079, 1700023891, 2056298683, 2149703411
Offset: 1

Views

Author

Jonathan Vos Post, Apr 03 2006

Keywords

Comments

See also: A116994 Prime partial sums of triangular numbers with prime indices. A116995 Pentagonal numbers with prime indices.

Examples

			a(1) = Sum_{i=1..1} prime(i)*(3*prime(i)-1)/2 = P(2) = 5.
a(2) = Sum_{i=1..2} prime(i)*(3*prime(i)-1)/2 = P(2) + P(3) = 17.
a(3) = Sum_{i=1..11} prime(i)*(3*prime(i)-1)/2 = P(2) + P(3) + P(5) + P(7) + P(11) + P(13) + P(17) + P(19) + P(23) + P(29) + P(31) = 4957.
a(4) = P(2) + ... + P(103) = 129277.
		

Crossrefs

Programs

  • Maple
    P:=n->n*(3*n-1)/2: seq(P(n),n=0..10): a:=proc(n) if isprime(sum(P(ithprime(j)),j=1..n))=true then sum(P(ithprime(j)),j=1..n) else fi end: seq(a(n),n=1..600); # Emeric Deutsch, Apr 15 2006
  • Mathematica
    Module[{nn=4000,pn,pr},pn=PolygonalNumber[5,Range[nn]];pr=Table[If[ PrimeQ[ n],1,0],{n,nn}];Select[Accumulate[Pick[pn,pr,1]],PrimeQ]] (* Harvey P. Dale, Jan 27 2020 *)

Formula

A000040 INTERSECTION {Partial sums of A116995(n)}. (Sum_{i=1..k} A000326(A000040(i))) iff in A000040. (Sum_{i=1..k} prime(i)*(3*prime(i)-1)/2) iff in A000040.

Extensions

More terms from Emeric Deutsch, Apr 15 2006
Showing 1-2 of 2 results.