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.

A049996 a(n) is the index k such that F(k)=d(n), where d=A049874 (difference sequence of ordered products of Fibonacci numbers).

This page as a plain text file.
%I A049996 #13 May 28 2019 19:31:41
%S A049996 1,1,1,3,1,3,3,4,3,1,5,4,3,6,5,1,3,7,6,3,4,8,7,3,1,5,9,8,4,3,6,10,9,5,
%T A049996 1,3,7,11,10,6,3,4,8,12,11,7,3,1,5,9,13,12,8,4,3,6,10,14,13,9,5,1,3,7,
%U A049996 11,15,14,10,6,3,4,8,12,16,15,11,7,3,1,5,9,13,17
%N A049996 a(n) is the index k such that F(k)=d(n), where d=A049874 (difference sequence of ordered products of Fibonacci numbers).
%H A049996 Michel Marcus, <a href="/A049996/b049996.txt">Table of n, a(n) for n = 1..1010</a>
%t A049996 Block[{nn = 123, s, t}, s = Differences@ Take[#, nn] &@ Union@ Flatten[Table[Fibonacci[i]*Fibonacci[j], {i, 0, nn}, {j, i + 1, nn}]]; t = Fibonacci@ Range@ nn; Array[First@ FirstPosition[t, s[[#]] ] &, Length@ s]] (* _Michael De Vlieger_, May 27 2019 *)
%o A049996 (PARI) ifib(n) = if (n==1, 1, log(n*sqrt(5) + 1/2)\log((1+sqrt(5))/2));
%o A049996 lista(nn) = {my(out = List([0])); for (i=0, nn, for (j=i+1, nn, listput(out, fibonacci(i)*fibonacci(j)););); my(v = Vec(vecsort(select(x->(x < fibonacci(nn+1)), out), , 8))); my(w = vector(#v-1, k, v[k+1] - v[k])); vector(#w, k, ifib(w[k]));} \\ _Michel Marcus_, May 27 2019
%Y A049996 Cf. A049862, A049874.
%K A049996 nonn
%O A049996 1,4
%A A049996 _Clark Kimberling_
%E A049996 More terms from _Michel Marcus_, May 27 2019