|
%I #14 Mar 15 2014 06:50:06
%S 1,0,1,1,0,2,0,1,2,0,2,1,0,3,0,2,2,0,3,0,1,3,0,3,2,0,4,0,2,3,0,3,1,0,
%T 4,0,3,3,0,5,0,2,4,0,4,2,0,5,0,3,3,0,4,0,1,4,0,4,3,0,6,0,3,5,0,5,2,0,
%U 6,0,4,4,0,6,0,2,5,0,5,3,0,6,0,3,4,0,4
%N Number of partitions of n into distinct parts all of which are Fibonacci numbers greater than 1.
%C a(n) > 0 if n+1 is a term of the lower Wythoff sequence, A000201; a(n) = 0 if n+1 is a term of the upper Wythoff sequence, A001950.
%H Alois P. Heinz, <a href="/A239002/b239002.txt">Table of n, a(n) for n = 0..10946</a>
%F G.f.: Product_{i>=3} (1+x^Fibonacci(i)). - _Alois P. Heinz_, Mar 15 2014
%p F:= combinat[fibonacci]:
%p b:= proc(n, i) option remember; `if`(n=0, 1, `if`(i<3, 0,
%p b(n, i-1)+`if`(F(i)>n, 0, b(n-F(i), i-1))))
%p end:
%p a:= proc(n) local j; for j from ilog[(1+sqrt(5))/2](n+1)
%p while F(j+1)<=n do od; b(n, j)
%p end:
%p seq(a(n), n=0..100); # _Alois P. Heinz_, Mar 15 2014
%t f = Table[Fibonacci[n], {n, 3, 75}]; b[n_] := SeriesCoefficient[Product[1 + x^f[[k]], {k, n}], {x, 0, n}]; u = Table[b[n], {n, 0, 60}] (* A239002 *)
%t Flatten[Position[u, 0]] (* A001950 *)
%Y Cf. A000201, A001950, A000045, A000119, A239003, A000009.
%K nonn,easy,look
%O 0,6
%A _Clark Kimberling_, Mar 08 2014
|
