login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A078465 Primonacci numbers: a(n)=a(n-2)+a(n-3)+a(n-5)+a(n-7)+a(n-11)+...+a(n-p(k))+... until n > p(k), where p(k) is the k-th prime. a(1)=a(2)=1. 2
1, 1, 1, 2, 2, 4, 5, 8, 12, 16, 26, 36, 55, 81, 118, 177, 257, 384, 564, 833, 1233, 1813, 2685, 3956, 5845, 8629, 12731, 18807, 27746, 40976, 60481, 89282, 131816, 194562, 287253, 424018, 625968, 924077 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,4

COMMENTS

a(n)/a(n-1) -> 1.476229...=1/x, where x satisfies the Sum x^p(n)=1 equation, i.e. x^2+x^3+x^5+x^7+x^11+... =1. (What constant is it?)

LINKS

T. D. Noe, Table of n, a(n) for n=1..500

EXAMPLE

a(12) = 36 = a(12-2)+a(12-3)+a(12-5)+a(12-7)+a(12-11) = a(10)+a(9)+a(7)+a(5)+a(1) = 16+12+5+2+1 = 36.

MATHEMATICA

a[1] = a[2] = 1; a[n_] := a[n] = Sum[a[n - Prime[k]], {k, 1, PrimePi[n]}]; Table[a[n], {n, 1, 38}] (* Jean-Fran├žois Alcover, Mar 22 2011 *)

PROG

(Haskell)

import Data.List (genericIndex)

a078465 n = a078465_list `genericIndex` (n-1)

a078465_list = 1 : 1 : f 3 where

   f x = (sum $ map (a078465 . (x -)) $

         takeWhile (< x) a000040_list) : f (x + 1)

-- Reinhard Zumkeller, Jul 20 2012

CROSSREFS

Cf. A078974 (the constant 1.47622...), A084256 (the constant 1/1.47622...)

Sequence in context: A095719 A153952 A050364 * A094992 A172128 A274154

Adjacent sequences:  A078462 A078463 A078464 * A078466 A078467 A078468

KEYWORD

easy,nice,nonn

AUTHOR

Miklos Kristof, Jan 02 2003

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 15 22:25 EDT 2019. Contains 328038 sequences. (Running on oeis4.)