login
Pentanacci pi function: a(1)=a(2)=a(3)=a(4)=a(5)=1; for n>5, a(n) = pi(Sum_{j=1..5} a(n-j)) where pi = A000720.
2

%I #17 Apr 06 2023 02:21:11

%S 1,1,1,1,1,3,4,4,6,7,9,10,11,14,15,17,19,21,23,24,27,30,30,32,34,36,

%T 37,39,40,42,44,46,47,47,48,50,51,53,53,54,55,56,58,58,60,61,62,62,62,

%U 63,63,64,65,66,66,66,66,66,66,66,66,66,66,66,66,66,66,66,66,66,66,66,66

%N Pentanacci pi function: a(1)=a(2)=a(3)=a(4)=a(5)=1; for n>5, a(n) = pi(Sum_{j=1..5} a(n-j)) where pi = A000720.

%C Starting with other values of a(1), a(2), a(3), a(4), a(5) what behaviors are possible? Does the sequence always stick at a single integer after some point, or can it go into a loop, or is there a third pattern?

%C a(n) is equal to 66 for 54 <= n <= 10^7. - _G. C. Greubel_, Apr 06 2023

%H G. C. Greubel, <a href="/A100478/b100478.txt">Table of n, a(n) for n = 1..10000</a>

%H Andrew Booker, <a href="https://t5k.org/nthprime/index.php#piofx">The Nth Prime Page</a>.

%H I. Flores, <a href="http://www.fq.math.ca/Scanned/5-3/flores.pdf">k-Generalized Fibonacci numbers</a>, Fib. Quart., 5 (1967), 258-266.

%H V. E. Hoggatt, Jr. and M. Bicknell, <a href="http://www.fq.math.ca/Scanned/7-4/hoggatt-a.pdf">Diagonal sums of generalized Pascal triangles</a>, Fib. Quart., 7 (1969), 341-358, 393.

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PrimeCountingFunction.html">Prime Counting Function</a>

%F a(n) = pi(a(n-1) + a(n-2) + a(n-3) + a(n-4) + a(n-5)) with a(1) = a(2) = a(3) = a(4) = a(5) = 1.

%e a(6) = pi(a(1)+a(2)+a(3)+a(4)+a(5)) = pi(1+1+1+1+1) = pi(5) = 3.

%e a(7) = pi(a(2)+a(3)+a(4)+a(5)+a(6)) = pi(1+1+1+1+3) = pi(7) = 4.

%e a(8) = pi(a(3)+a(4)+a(5)+a(6)+a(7)) = pi(1+1+1+3+4) = pi(10) = 4.

%e a(9) = pi(a(4)+a(5)+a(6)+a(7)+a(8)) = pi(1+1+3+4+4) = pi(13) = 6.

%e a(10) = pi(a(5)+a(6)+a(7)+a(8)+a(9)) = pi(1+3+4+4+6) = pi(18) = 7.

%t a[n_]:= a[n]= If[n<6,1,PrimePi[Sum[a[n-j], {j,5}]]];

%t Table[a[n], {n,80}] (* _Robert G. Wilson v_, Dec 03 2004 *)

%o (SageMath)

%o @CachedFunction

%o def a(n): # a = A100478

%o if (n<6): return 1

%o else: return prime_pi(sum(a(n-j) for j in range(1,6)))

%o [a(n) for n in range(1, 81)] # _G. C. Greubel_, Apr 06 2023

%Y Cf. A001591, A038607.

%K easy,nonn

%O 1,6

%A _Jonathan Vos Post_, Nov 22 2004

%E Edited and extended by _Robert G. Wilson v_, Dec 03 2004