|
%I #34 Feb 24 2021 02:48:18
%S 1,3,4,6,6,6,8,12,10,6,8,12,12,12,16,24,18,6,8,12,12,12,16,24,20,12,
%T 16,24,24,24,32,48,34,6,8,12,12,12,16,24,20,12,16,24,24,24,32,48,36,
%U 12,16,24,24,24,32,48,40,24,32,48,48,48,64,96,66,6,8,12,12,12,16,24,20,12,16
%N Convolution of binomial(1,n) and Gould's sequence A001316.
%C A universal function related to the spherical growth of repeated truncations of maps.
%C a(n) = (number of ones in row n of triangle A249133) = (number of odd terms in row n of triangle A249095) = A000120(A249184(n)). - _Reinhard Zumkeller_, Nov 14 2014
%H Reinhard Zumkeller, <a href="/A105321/b105321.txt">Table of n, a(n) for n = 0..10000</a>
%H David Applegate, Omar E. Pol and N. J. A. Sloane, <a href="/A000695/a000695_1.pdf">The Toothpick Sequence and Other Sequences from Cellular Automata</a>, Congressus Numerantium, Vol. 206 (2010), 157-191. [There is a typo in Theorem 6: (13) should read u(n) = 4.3^(wt(n-1)-1) for n >= 2.]
%H T. Pisanski and T. W. Tucker, <a href="http://www.imfm.si/preprinti/PDF/00696.pdf">Growth in Repeated Truncations of Maps</a>, Preprint series, Univ. of Ljubljana, Vol. 38 (2000), 696.
%H T. Pisanski and T. W. Tucker, <a href="http://www.imfm.si/preprinti/PDF/00696.pdf">Growth in Repeated Truncations of Maps</a>, Atti. Sem. Mat. Fis. Univ. Modena, Vol. 49 (2001), 167-176.
%H N. J. A. Sloane, <a href="/wiki/Catalog_of_Toothpick_and_CA_Sequences_in_OEIS">Catalog of Toothpick and Cellular Automata Sequences in the OEIS</a> - _Omar E. Pol_, May 29 2010
%F G.f. (1+x)*Product{k>=0, 1+2x^(2^k)};
%F a(n) = Sum_{k=0..n, binomial(1, n-k)*Sum_{j=0..k, binomial(k, j) mod 2}}.
%F a(n) = 2*A048460(n) for n>=2. - _Omar E. Pol_, Jan 02 2011
%F a((2*n-1)*2^p) = (2^p+2)*A001316(n-1), p >= 0 and n >= 1, with a(0) = 1. - _Johannes W. Meijer_, Jan 28 2013
%F a(n) = A001316(n) + A001316(n-1) for n > 0. - _Reinhard Zumkeller_, Nov 14 2014
%e Contribution from _Omar E. Pol_, May 29 2010: (Start)
%e If written as a triangle:
%e 1;
%e 3;
%e 4;
%e 6,6;
%e 6,8,12,10;
%e 6,8,12,12,12,16,24,18;
%e 6,8,12,12,12,16,24,20,12,16,24,24,24,32,48,34;
%e 6,8,12,12,12,16,24,20,12,16,24,24,24,32,48,36,12,16,24,24,24,32,48,40,24,32,48,48,48,64,96,66;
%e (End)
%p nmax := 74: A001316 := n -> if n <= -1 then 0 else 2^add(i, i=convert(n, base, 2)) fi: for p from 0 to ceil(log[2](nmax)) do for n from 1 to nmax/(p+2)+1 do a((2*n-1)*2^p) := (2^p+2) * A001316(n-1) od: od: a(0) :=1: seq(a(n), n=0..nmax); # _Johannes W. Meijer_, Jan 28 2013
%t f[n_] := Sum[Binomial[1, n - k]Mod[Binomial[k, j], 2], {k, 0, n}, {j, 0, k}]; Array[f, 75, 0] (* _Robert G. Wilson v_, Jun 28 2010 *)
%o (Haskell)
%o a105321 n = if n == 0 then 1 else a001316 n + a001316 (n - 1)
%o -- _Reinhard Zumkeller_, Nov 14 2014
%o (PARI) a(n) = sum(k=0, n, binomial(1, n-k)*sum(j=0, k, binomial(k, j) % 2)); \\ _Michel Marcus_, Apr 29 2018
%Y Cf. A139250, A139251, A173522, A048460.
%Y Cf. A001316, A000120, A249184, A249095, A249133.
%K easy,nonn
%O 0,2
%A _Paul Barry_, Apr 01 2005
|
