OFFSET
1,6
LINKS
Iain Fox, Table of n, a(n) for n = 1..400 (first 200 terms from T. D. Noe)
MATHEMATICA
ad[1] = 0; ad[n_] := Module[{f = FactorInteger[n]}, Total[n*f[[All, 2]]/ f[[All, 1]]]]; a[n_] := ad[Fibonacci[n]]; Array[a, 40] (* Jean-François Alcover, Feb 22 2018 *)
PROG
(Magma) Ad:=func<h | h*(&+[Factorisation(h)[i][2]/Factorisation(h)[i][1]: i in [1..#Factorisation(h)]])>; [n le 2 select 0 else Ad(Fibonacci(n)): n in [1..40]]; // Bruno Berselli, Oct 22 2013
(PARI) a(n) = my(f = factor(n=fibonacci(n))~); sum(i=1, #f, n/f[1, i]*f[2, i]) \\ Iain Fox, Oct 29 2018
(GAP) a:=Concatenation([0, 0], List(List([3..40], n->Factors(Fibonacci(n))), i->Product(i)*Sum(i, j->1/j))); # Muniru A Asiru, Oct 31 2018
(Python)
from sympy import fibonacci, factorint
def A068329(n):
f = fibonacci(n)
return sum((f*e//p for p, e in factorint(f).items())) if n > 2 else 0 # Chai Wah Wu, Jun 12 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Feb 27 2002
STATUS
approved