|
%I #10 Mar 05 2018 04:56:57
%S 1,1,1,1,1,2,1,2,2,1,2,2,1,2,2,1,2,3,2,2,1,2,3,2,2,2,1,2,3,2,3,2,3,2,
%T 1,2,3,2,3,2,3,2,1,2,2,3,3,2,3,3,2,3,2,1,2,2,3,4,3,2,3,3,2,3,2,2,1,2,
%U 3,2,3,4,3,2,3,3,3,2,3,3,2,2,1,2,3,2
%N Number of factors L(i) > 1 of A274280(n), where L = A000032 (Lucas numbers, 1,3,4,...)
%H G. C. Greubel, <a href="/A274353/b274353.txt">Table of n, a(n) for n = 1..1000</a>
%e The products of distinct Lucas numbers, arranged in increasing order, comprise A274280. The list begins with 1, 3, 4, 7, 11, 12 = 3*4, so that a(6) = 2.
%t r[1] := 1; r[2] := 3; r[n_] := r[n] = r[n - 1] + r[n - 2];
%t s = {1}; z = 40; f = Map[r, Range[z]]; Take[f, 10]
%t Do[s = Union[s, Select[s*f[[i]], # <= f[[z]] &]], {i, z}];
%t infQ[n_] := MemberQ[f, n];
%t a = Table[#[[Flatten[Position[Map[Apply[Times, #] &, #], s[[n]]]][[1]]]] &[
%t Rest[Subsets[Map[#[[1]] &, Select[Map[{#, infQ[#]} &,
%t Divisors[s[[n]]]], #[[2]] && #[[1]] > 1 &]]]]], {n, 2, 200}];
%t ans = Join[{{1}}, a]; Take[ans, 8]
%t w = Map[Length, ans] (* A274353 *)
%t Map[Apply[Times, #] &, Select[ans, Length[#] == 2 &]] (* A274347 *)
%t Map[Apply[Times, #] &, Select[ans, Length[#] == 3 &]] (* A274348 *)
%t (* _Peter J. C. Moses_, Jun 17 2016 *)
%Y Cf. A000032, A274280, A274347, A274348, A160009, A272947, A274354.
%K nonn,easy
%O 1,6
%A _Clark Kimberling_, Jun 18 2016
|
