|
|
A005086
|
|
Number of Fibonacci numbers 1,2,3,5,... dividing n.
|
|
23
|
|
|
1, 2, 2, 2, 2, 3, 1, 3, 2, 3, 1, 3, 2, 2, 3, 3, 1, 3, 1, 3, 3, 2, 1, 4, 2, 3, 2, 2, 1, 4, 1, 3, 2, 3, 2, 3, 1, 2, 3, 4, 1, 4, 1, 2, 3, 2, 1, 4, 1, 3, 2, 3, 1, 3, 3, 3, 2, 2, 1, 4, 1, 2, 3, 3, 3, 3, 1, 3, 2, 3, 1, 4, 1, 2, 3, 2, 1, 4, 1, 4, 2, 2, 1, 4, 2, 2, 2, 3, 2, 4, 2, 2, 2, 2, 2, 4, 1, 2, 2, 3, 1, 4, 1, 4, 4
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
|
|
LINKS
|
|
|
FORMULA
|
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A079586 - 1 = 2.359885... . - Amiram Eldar, Dec 31 2023
|
|
MAPLE
|
with(combinat): for n from 1 to 200 do printf(`%d, `, sum(floor(n/fibonacci(k))-floor((n-1)/fibonacci(k)), k=2..15)) od:
|
|
MATHEMATICA
|
f[n_] := Block[{k = 1}, While[Fibonacci[k] <= n, k++ ]; Count[ Mod[n, Array[ Fibonacci, k - 1]], 0] - 1]; Array[f, 105] (* Robert G. Wilson v, Dec 10 2006 *)
|
|
PROG
|
(PARI) isfib(n)=my(k=n^2); k+=(k+1)<<2; issquare(k) || (n>0 && issquare(k-8))
(Python)
from sympy import divisors
from sympy.ntheory.primetest import is_square
def A005086(n): return sum(1 for d in divisors(n, generator=True) if is_square(m:=5*d**2-4) or is_square(m+8)) # Chai Wah Wu, Mar 30 2023
(Python)
from itertools import count, takewhile
def F(f=1, g=1):
while True:
f, g = g, f+g;
yield f
def a(n):
return sum(1 for f in takewhile(lambda x: x<=n, F()) if n%f == 0)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|