%I #9 Mar 17 2019 12:21:46

%S 1,2,3,4,5,7,9,12,16,21,27,35,44,56,70,87,108,133,163,199,242,292,352,

%T 421,504,599,712,841,994,1167,1371,1602,1873,2179,2535,2936,3401,3924,

%U 4528,5206,5985,6858,7857,8976,10252,11679,13299,15109,17159,19446,22028

%N Positions of Fibonacci numbers in ordered sequence A160009 of all products of Fibonacci numbers.

%H Rémy Sigrist, <a href="/A272948/b272948.txt">Table of n, a(n) for n = 1..121</a>

%H <a href="http://mathoverflow.net/questions/238505/distinctness-of-products-of-fibonacci-numbers">Distinctness of Fibonacci products</a>

%H Rémy Sigrist, <a href="/A272948/a272948.gp.txt">PARI program for A272948</a>

%e A160009 = (0,1,2,3,5,6,8,10,13,15,16,21,...), so that a = (1,2,3,4,5,7,9,12,...).

%t s = {1}; nn = 60; f = Fibonacci[2 + Range[nn]]; Do[s = Union[s, Select[s*f[[i]], # <= f[[nn]] &]], {i, nn}]; s = Prepend[s, 0]; Take[s, 100] (* A160009 *)

%t isFibonacciQ[n_] := Apply[Or, Map[IntegerQ, Sqrt[{# + 4, # - 4} &[5 n^2]]]];

%t ans = Join[{{0}}, {{1}}, Table[#[[Flatten[Position[Map[Apply[Times, #] &, #], s[[n]]]][[1]]]] &[Rest[Subsets[Rest[Map[#[[1]] &, Select[Map[{#, isFibonacciQ[#]} &, Divisors[s[[n]]]], #[[2]] &]]]]]], {n, 3, 500}]]

%t Map[Length, ans] (* A272947 *)

%t Flatten[Position[Map[Length, ans], 1]] (* A272948 *)

%t Map[Apply[Times, #] &, Select[ans, Length[#] == 1 &]] (* A000045 *)

%t Map[Apply[Times, #] &, Select[ans, Length[#] == 2 &]] (* A271354 *)

%t Map[Apply[Times, #] &, Select[ans, Length[#] == 3 &]] (* A272949 *)

%t Map[Apply[Times, #] &, Select[ans, Length[#] == 4 &]] (* A272950 *)

%t (* _Peter J. C. Moses_, May 11 2016 *)

%Y Cf. A000045, A160009, A272947, A271354, A272949, A273950.

%K nonn

%O 1,2

%A _Clark Kimberling_, May 13 2016

%E More terms from _Rémy Sigrist_, Mar 17 2019