login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A104589 a(1)=1. a(n) = a(n-1) + (sum of terms, from among terms a(1) through a(n-1), which are prime or 1). 4
1, 2, 5, 13, 34, 55, 76, 97, 215, 333, 451, 569, 1256, 1943, 2630, 3317, 4004, 4691, 10069, 25516, 40963, 56410, 71857, 87304, 102751, 118198, 133645, 149092, 164539, 179986, 195433, 210880, 226327, 241774, 257221, 529889, 802557, 1075225 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
By Dirichlet's Theorem there are an infinite number of primes in this sequence.
LINKS
FORMULA
a(n) = 3*a(n-1) - a(n-2) if a(n-1) is prime, else a(n) = 2*a(n-1) - a(n-2) for n>3. - John Tyler Rascoe, Jul 20 2022
EXAMPLE
The noncomposites among the first 8 terms of the sequence are 1, 2, 5, 13 and 97. The sum of these is 1+2+5+13+97 = 118. So a(9) = a(8) + 118 = 215.
MATHEMATICA
f[lst_] := Append[lst, Last@ lst + Plus @@ Select[lst, (PrimeQ@ # || # == 1) &]]; Nest[f, {1}, 38] (* Robert G. Wilson v, Jul 02 2007 *)
PROG
(PARI) lista(nn) = my(va = vector(nn), s = 1); va[1] = 1; for (n=2, nn, va[n] = va[n-1] + s; if (isprime(va[n]), s += va[n]); ); va; \\ Michel Marcus, Jul 21 2022
CROSSREFS
Cf. A008578.
Sequence in context: A271940 A273721 A112841 * A154101 A122024 A252932
KEYWORD
nonn
AUTHOR
Leroy Quet, Jun 12 2007
EXTENSIONS
More terms from Robert G. Wilson v, Jul 02 2007
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 18 22:56 EDT 2024. Contains 370952 sequences. (Running on oeis4.)