login
a(n) = Product_{i in row n of A245562} Fibonacci(i+1).
3

%I #32 Sep 28 2021 01:21:44

%S 1,1,1,2,1,1,2,3,1,1,1,2,2,2,3,5,1,1,1,2,1,1,2,3,2,2,2,4,3,3,5,8,1,1,

%T 1,2,1,1,2,3,1,1,1,2,2,2,3,5,2,2,2,4,2,2,4,6,3,3,3,6,5,5,8,13,1,1,1,2,

%U 1,1,2,3,1,1,1,2,2,2,3,5,1,1,1,2,1,1,2,3,2,2,2,4,3,3,5,8,2,2,2,4,2

%N a(n) = Product_{i in row n of A245562} Fibonacci(i+1).

%C This is the Run Length Transform of S(n) = Fibonacci(n+1).

%C The Run Length Transform of a sequence {S(n), n>=0} is defined to be the sequence {T(n), n>=0} given by T(n) = Product_i S(i), where i runs through the lengths of runs of 1's in the binary expansion of n. E.g., 19 is 10011 in binary, which has two runs of 1's, of lengths 1 and 2. So T(19) = S(1)*S(2). T(0)=1 (the empty product).

%H Alois P. Heinz, <a href="/A246028/b246028.txt">Table of n, a(n) for n = 0..8191</a>

%H N. J. A. Sloane, On the No. of ON Cells in Cellular Automata, Video of talk in Doron Zeilberger's Experimental Math Seminar at Rutgers University, Feb 05 2015: <a href="https://vimeo.com/119073818">Part 1</a>, <a href="https://vimeo.com/119073819">Part 2</a>

%H N. J. A. Sloane, <a href="http://arxiv.org/abs/1503.01168">On the Number of ON Cells in Cellular Automata</a>, arXiv:1503.01168, 2015

%H Chai Wah Wu, <a href="https://arxiv.org/abs/1610.06166">Sums of products of binomial coefficients mod 2 and run length transforms of sequences</a>, arXiv:1610.06166 [math.CO], 2016.

%H <a href="/index/Ce#cell">Index entries for sequences related to cellular automata</a>

%F a(n) = Sum_{k=0..n} ((binomial(n-k,2k)*binomial(n,k)) mod 2). - _Chai Wah Wu_, Oct 19 2016

%p with(combinat); ans:=[];

%p for n from 0 to 100 do lis:=[]; t1:=convert(n,base,2); L1:=nops(t1); out1:=1; c:=0;

%p for i from 1 to L1 do

%p if out1 = 1 and t1[i] = 1 then out1:=0; c:=c+1;

%p elif out1 = 0 and t1[i] = 1 then c:=c+1;

%p elif out1 = 1 and t1[i] = 0 then c:=c;

%p elif out1 = 0 and t1[i] = 0 then lis:=[c,op(lis)]; out1:=1; c:=0;

%p fi;

%p if i = L1 and c>0 then lis:=[c,op(lis)]; fi;

%p od:

%p a:=mul(fibonacci(i+1), i in lis);

%p ans:=[op(ans),a];

%p od:

%p ans;

%t a[n_] := Sum[Mod[Binomial[n-k, 2k] Binomial[n, k], 2], {k, 0, n}];

%t a /@ Range[0, 100] (* _Jean-François Alcover_, Feb 28 2020, after _Chai Wah Wu_ *)

%o (PARI) a(n)=my(s=1,k); while(n, n>>=valuation(n,2); k=valuation(n+1,2); if(k>1, s*=fibonacci(k+1)); n>>=k); s \\ _Charles R Greathouse IV_, Oct 21 2016

%o (PARI) a(n)=sum(k=0,n, !bitand(n-3*k,2*k) && !bitand(n-k,k)) \\ _Charles R Greathouse IV_, Oct 21 2016

%o (Python)

%o def A246028(n): return sum(int(not (~(n-k) & 2*k) | (~n & k)) for k in range(n+1)) # _Chai Wah Wu_, Sep 27 2021

%Y Cf. A245562, A000045, A001045, A071053, A245565, A245564.

%K nonn,easy

%O 0,4

%A _N. J. A. Sloane_, Aug 15 2014; revised Sep 05 2014