A053872 Primes arising in A053782.

37, 199, 277, 367, 997, 1049, 1103, 1451, 1709, 1777, 2137, 2213, 2953, 3041, 3407, 3889, 3989, 4091, 5309, 7193, 7883, 11113, 13757, 15083, 16073, 17093, 28643, 28909, 30259, 32203, 35089, 35977, 40879, 42157, 43451, 43777, 48119, 51949

Offset: 1

Author

G. L. Honaker, Jr., Mar 30 2000

Keywords

• nonn

Links

• Chai Wah Wu, Table of n, a(n) for n = 1..10000

Crossrefs

Cf. A000040, A002808, A053782.

Programs

Mathematica

s=0;lst={};Do[If[PrimeQ[n],NULL,s+=n;If[PrimeQ[s],AppendTo[lst,s]]],{n,2,7!}];lst (* Vladimir Joseph Stephan Orlovsky, Feb 01 2009 *)

PARI

lista(nn) = {my(s = 0); forcomposite(c=1, nn, s += c; if (isprime(s), print1(s, ", ")););} \\ Michel Marcus, May 13 2018

Python

from sympy import isprime
A053872_list, n, m, s = [], 1, 4, 4
while len(A053872_list) < 10000:
    if isprime(s):
        A053872_list.append(s)
    m += 1
    if isprime(m):
        m += 1
    n += 1
    s += m # Chai Wah Wu, May 13 2018

Extensions

More terms from Reiner Martin, Jul 17 2001