A318464 Additive with a(p^e) = A007895(e), where A007895(n) gives the number of terms in Zeckendorf representation of n.

0, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 3, 1, 1, 2, 2, 2, 2, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 3, 1, 2, 2, 2, 2, 3, 1, 2, 2, 3, 1, 2, 1, 2, 2, 2, 2, 3, 1, 3, 2, 2, 1, 3, 2, 2, 2, 2, 1, 3, 2, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 3, 1, 2, 3

Offset

1,6

Links

Antti Karttunen, Table of n, a(n) for n = 1..65537

Index entries for sequences computed from exponents in factorization of n.

Formula

a(n) = A007814(A318465(n)).

a(n) = A001222(A318469(n)).

Sum_{k=1..n} a(k) ~ n * (log(log(n)) + B + C), where B is Mertens's constant (A077761) and C = Sum_{k>=2} (A007895(k)-A007895(k-1)) * P(k) = 0.05631817952062180045..., where P(s) is the prime zeta function. - Amiram Eldar, Oct 09 2023

Mathematica

z[n_] := Length[DeleteCases[NestWhileList[# - Fibonacci[Floor[Log[Sqrt[5]*# + 3/2]/Log[GoldenRatio]]] &, n, # > 1 &], 0]]; a[n_] := Total[z /@ FactorInteger[n][[;; , 2]]]; a[1] = 0; Array[a, 100] (* Amiram Eldar, May 15 2023 *)

Prog

(PARI)
A072649(n) = { my(m); if(n<1, 0, m=0; until(fibonacci(m)>n, m++); m-2); }; \\ From A072649
A007895(n) = { my(s=0); while(n>0, s++; n -= fibonacci(1+A072649(n))); (s); }
A318464(n) = vecsum(apply(e -> A007895(e), factor(n)[, 2]));

Crossrefs

Cf. A007895, A077761, A318465, A318469.

Sequence in context: A351202 A333145 A335449 * A051265 A008647 A036475

Adjacent sequences: A318461 A318462 A318463 * A318465 A318466 A318467

Keyword

nonn

Author

Antti Karttunen, Aug 30 2018

Status

approved