OFFSET
1,3
COMMENTS
This sequence can be used as a filter for certain sequences involving Fibonacci numbers as it matches to any sequence that is obtained as f(A000045(n)), where f(n) is any function that depends only on the prime signature of n (see the index entry for "sequences computed from exponents in ...").
Matching in this context means that the sequence a matches with the sequence b iff for all i, j: a(i) = a(j) => b(i) = b(j). In other words, iff the sequence b partitions the natural numbers to the same or coarser equivalence classes (as/than the sequence a) by the distinct values it obtains.
LINKS
Antti Karttunen (terms 1..374) & Hans Havermann, Table of n, a(n) for n = 1..1300
EXAMPLE
From Michael De Vlieger, May 18 2017: (Start)
a(6) = 8 because Fibonacci(6) = 8, the multiplicity of the prime factor of 8 is 3; the smallest p^3 = 2^3 = 8.
a(7) = 2 because Fibonacci(7) = 13, the multiplicity of the prime factor of 13 is 1; the smallest p^1 = 2^1 = 2.
a(15) = 30 because Fibonacci(15) = 610. The multiplicities of the prime factors of 610, in order from greatest to least, are {1, 1, 1}, the smallest prime power product p^1 * q^1 * r^1 = 2 * 3 * 5 = 30.
a(18) = 120 because Fibonacci(18) = 2584 = 2^3 * 17 * 19 -> 2^3 * 3 * 5 = 120. (End)
MATHEMATICA
Table[If[# == 1, 1, Times @@ MapIndexed[Prime[First[#2]]^#1 &,
Sort[FactorInteger[#][[All, -1]], Greater]]] &@ Fibonacci@ n, {n, 79}] (* Michael De Vlieger, May 18 2017 *)
PROG
(PARI)
A046523(n) = my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]) \\ From Charles R Greathouse IV, Aug 17 2011
f0 = 0; f1 = 1; for(n=1, 10000, write("b278245.txt", n, " ", A046523(f1)); old_f0 = f0; f0 = f1; f1 = f1 + old_f0; );
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Nov 16 2016
STATUS
approved