login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A038575 Number of prime factors of n-th Fibonacci number, counted with multiplicity. 18

%I #35 Feb 02 2022 15:11:29

%S 0,0,0,1,1,1,3,1,2,2,2,1,6,1,2,3,3,1,5,2,4,3,2,1,9,3,2,4,4,1,7,2,4,3,

%T 2,3,10,3,3,3,6,2,7,1,5,5,3,1,12,3,6,3,4,2,8,4,7,5,3,2,12,2,3,5,6,3,7,

%U 3,5,5,7,2,14,2,4,6,5,4,8,2,9,7,3,1,13,4,3,4,9,2,12,5,6,4,2,6,16,4,5,6,10,2,8

%N Number of prime factors of n-th Fibonacci number, counted with multiplicity.

%C Row lengths of table A060441. - _Reinhard Zumkeller_, Aug 30 2014

%H Amiram Eldar, <a href="/A038575/b038575.txt">Table of n, a(n) for n = 0..1408</a> (terms 0..1000 from T. D. Noe derived from Kelly's data)

%H Blair Kelly, <a href="http://mersennus.net/fibonacci/">Fibonacci and Lucas Factorizations</a>.

%H Douglas Lind, <a href="https://fq.math.ca/Scanned/6-6/advanced6-6.pdf">Problem H-145</a>, Advanced Problems and Solutions, The Fibonacci Quarterly, Vol. 6, No. 6 (1968), p. 351; <a href="https://www.fq.math.ca/Scanned/8-4/advanced8-4.pdf">Factor Analysis</a>, Solution to Problem H-145 by the proposer, ibid., Vol. 8, No. 4 (1970), pp. 386-387.

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/FibonacciNumber.html">Fibonacci Number</a>.

%F For n > 0: a(n) = A001222(A000045(n)). - _Reinhard Zumkeller_, Aug 30 2014

%F a(n) >= A001222(n) - 1 (Lind, 1968). - _Amiram Eldar_, Feb 02 2022

%e a(12) = 6 because Fibonacci(12) = 144 = 2^4 * 3^2 has 6 prime factors.

%p with(numtheory):with(combinat):a:=proc(n) if n=0 then 0 else bigomega(fibonacci(n)) fi end: seq(a(n), n=0..102); # _Zerinvary Lajos_, Apr 11 2008

%t Join[{0, 0}, Table[Plus@@(Transpose[FactorInteger[Fibonacci[n]]][[2]]), {n, 3, 102}]]

%t Join[{0},PrimeOmega[Fibonacci[Range[110]]]] (* _Harvey P. Dale_, Apr 14 2018 *)

%o (Haskell)

%o a038575 n = if n == 0 then 0 else a001222 $ a000045 n

%o -- _Reinhard Zumkeller_, Aug 30 2014

%o (PARI) a(n)=bigomega(fibonacci(n)) \\ _Charles R Greathouse IV_, Sep 14 2015

%o (Python)

%o from sympy import primeomega, fibonacci

%o def a(n): return 0 if n == 0 else primeomega(fibonacci(n))

%o print([a(n) for n in range(103)]) # _Michael S. Branicky_, Feb 02 2022

%Y Cf. A022307 (number of distinct prime factors), A086597 (number of primitive prime factors).

%Y Cf. also A001222, A000045, A060441.

%K nonn

%O 0,7

%A _Jeff Burch_

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 25 04:42 EDT 2024. Contains 371964 sequences. (Running on oeis4.)