cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A061446 Primitive part of Fibonacci(n).

Original entry on oeis.org

1, 1, 2, 3, 5, 4, 13, 7, 17, 11, 89, 6, 233, 29, 61, 47, 1597, 19, 4181, 41, 421, 199, 28657, 46, 15005, 521, 5777, 281, 514229, 31, 1346269, 2207, 19801, 3571, 141961, 321, 24157817, 9349, 135721, 2161, 165580141, 211, 433494437, 13201, 109441
Offset: 1

Views

Author

David Broadhurst, Jun 10 2001

Keywords

Comments

Fib(n) = A000045(n) = Product_{d|n} a(d), Lucas(n) = A000204(n) = Product_{d|2n and 2^m|d iff 2^m|2n} a(d) (e.g., Lucas(4) = 7 = a(8), Lucas(6) = 18 = a(12)*a(4)). - Len Smiley, Nov 11 2001
A 2001 Iranian Mathematical Olympiad question shows such a sequence exists whenever gcd(a(m),a(n)) = a(gcd(m,n)).
The problem of the characterization of the family of all GCD-morphic sequences F, i.e., F such that GCD(F(m),F(n)) = F(GCD(m,n)), was posed by A. K. Kwasniewski (GCD-morphic Problem). Dziemianczuk and Bajguz (2008) showed that any GCD-morphic sequence is coded by a certain natural number-valued sequence. - Maciej Dziemianczuk, Jan 15 2009
This is the LCM-transform of the Fibonacci numbers (cf. Nowicki). - N. J. A. Sloane, Jan 02 2016

Crossrefs

Cf. A000010 (comments on product formulas).

Programs

  • Maple
    N:= 200; # to get a(1) to a(N)
    L[0]:= 1:
    for i from 1 to N do L[i]:=ilcm(L[i-1],combinat:-fibonacci(i)) od:
    seq(L[i]/L[i-1],i=1..N); # Robert Israel, Aug 03 2015
  • Mathematica
    t={1}; Do[f=Fibonacci[n]; Do[f=f/GCD[f,t[[d]]], {d,Most[Divisors[n]]}]; AppendTo[t,f], {n,2,100}]; t
    (* Second program: *)
    a[n_] := Product[Fibonacci[d]^MoebiusMu[n/d], {d, Divisors[n]}];
    Array[a, 45] (* Jean-François Alcover, Jul 04 2019 *)
  • PARI
    a(n)=my(d=divisors(n)); fibonacci(n)/lcm(apply(fibonacci,d[1..#d-1])) \\ Charles R Greathouse IV, Oct 06 2016

Formula

Let r=(1+sqrt(5))/2. For n>2, the primitive part of F(n)=(r^n-(-1/r)^n)/sqrt(5) is Phi_n(-r^2)/r^phi(n) where Phi_n is n-th cyclotomic polynomial and phi is Euler's totient function A000010.
a(n) = Product_{d|n} Fib(d)^mu(n/d), where mu = A008683 (Bliss, Fulan, Lovett, Sommars, eq. (7)). - Jonathan Sondow, Jun 11 2013
a(n) = lcm(Fib(1),Fib(2),...,Fib(n))/lcm(Fib(1),Fib(2),...,Fib(n-1)). - Thomas Ordowski, Aug 03 2015
a(n) = Product_{k=1..n} Fib(gcd(n,k))^(mu(n/gcd(n,k))/phi(n/gcd(n,k))) = Product_{k=1..n} Fib(n/gcd(n,k))^(mu(gcd(n,k))/phi(n/gcd(n,k))) where mu = A008683, phi = A000010. - Richard L. Ollerton, Nov 08 2021

Extensions

More terms from Vladeta Jovovic, Nov 09 2001
Edited by N. J. A. Sloane at the suggestion of Andrew S. Plewe, May 29 2007
Edited by Charles R Greathouse IV, Oct 28 2009

A126015 Lucas numbers divided by their primitive part.

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 1, 1, 4, 3, 1, 7, 1, 3, 44, 1, 1, 18, 1, 7, 116, 3, 1, 47, 11, 3, 76, 7, 1, 738, 1, 1, 796, 3, 319, 322, 1, 3, 2084, 47, 1, 5058, 1, 7, 25916, 3, 1, 2207, 29, 123, 14284, 7, 1, 5778, 2189, 47, 37396, 3, 1, 695842, 1, 3, 465044, 1, 5731, 237618, 1, 7, 256316
Offset: 1

Views

Author

T. D. Noe, Dec 16 2006

Keywords

Comments

Note that a(n)=1 for prime n and for n=2^k for k=0,1,2,...

Crossrefs

Cf. A000032 (Lucas numbers), A105602 (Fibonacci numbers divided by their primitive part).

Formula

a(n)=L(n)/A061447(n)

A274584 Composite numbers n such that primitive part of Fibonacci(n) is prime.

Original entry on oeis.org

4, 8, 9, 10, 14, 15, 16, 18, 20, 21, 22, 26, 28, 30, 32, 33, 34, 35, 38, 39, 40, 42, 45, 51, 52, 54, 60, 62, 63, 65, 66, 72, 74, 75, 76, 82, 93, 94, 98, 105, 106, 111, 112, 119, 121, 122, 123, 124, 132, 135, 136, 140, 142, 144, 145, 152, 156, 158, 160, 172, 180
Offset: 1

Views

Author

Arkadiusz Wesolowski, Jun 29 2016

Keywords

Comments

For every n > 1, A105602(a(n)) > 1.

Crossrefs

Cf. A061446. Supersequence of A061254.

Programs

  • Magma
    lst:=[]; for n in [4..180] do if not IsPrime(n) then d:=Divisors(n); p:=Truncate(&*[Fibonacci(d[i])^MoebiusMu(Truncate(n/d[i])): i in [1..#d]]); if IsPrime(p) then Append(~lst, n); end if; end if; end for; lst;
Showing 1-3 of 3 results.