|
%I #71 Aug 06 2024 22:02:11
%S 1,2,3,4,5,4,6,8,9,4,6,8,10,8,12,16,17,4,6,8,10,8,12,16,18,8,12,16,20,
%T 16,24,32,33,4,6,8,10,8,12,16,18,8,12,16,20,16,24,32,34,8,12,16,20,16,
%U 24,32,36,16,24,32,40,32,48,64,65,4,6,8,10,8,12
%N Number of vertices in the n-th row of the Sierpinski gasket (cf. A047999).
%C Partial sums give A233774.
%C The subsequence of odd terms is A083318. - _Gary W. Adamson_, Jan 13 2014
%C Equivalently, this is the coordination sequence for the Sierpinski gasket with respect to the apex. - _N. J. A. Sloane_, Sep 19 2020
%H N. J. A. Sloane, <a href="/A233775/b233775.txt">Table of n, a(n) for n = 0..10000</a>
%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>
%H Eric Weisstein's World of Mathematics, <a href="https://mathworld.wolfram.com/SierpinskiGasketGraph.html">Sierpiński Gasket Graph</a>.
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Sierpinski_triangle">Sierpinski triangle</a>.
%H <a href="/index/Con#coordination_sequences">Index entries for coordination sequences</a>
%F a(0)=1, a(n) = (2^t(n) + 1) * 2^(c(n) - 1) where t(n) = A007814(n) is the number of trailing zeros in the binary representation of n and c(n) = A000120(n) is the total number of ones in the binary representation of n. - _Johan Falk_, Jun 24 2020
%e Illustration of initial terms:
%e --------------------------------------------------------
%e Diagram n a(n) A233774(n)
%e --------------------------------------------------------
%e * 0 1 1
%e /T\
%e *---* 1 2 3
%e /T\ /T\
%e *---*---* 2 3 6
%e /T\ /T\
%e *---* *---* 3 4 10
%e /T\ /T\ /T\ /T\
%e *---*---*---*---* 4 5 15
%e /T\ /T\
%e *---* *---* 5 4 19
%e --------------------------------------------------------
%e After five stages the number of "black" triangles in the structure is A006046(5) = 11 and the number of "black" triangles in row 5 is A001316(5-1) = 2. The number of vertices in row 5 is equal to 4, so a(5) = 4.
%e Written as an irregular triangle the sequence begins:
%e 1;
%e 2;
%e 3;
%e 4,5;
%e 4,6,8,9;
%e 4,6,8,10,8,12,16,17;
%e 4,6,8,10,8,12,16,18,8,12,16,20,16,24,32,33;
%e ...
%p A000120 := n -> add(i, i=convert(n, base, 2)):
%p A007814 := n -> padic[ordp](n, 2):
%p A233775 := n->(2^A007814(n)+1)*(2^(A000120(n)-1); # _N. J. A. Sloane_, Sep 19 2020
%t A233775[n_] := If[n == 0, 1, (2^IntegerExponent[n, 2]+1)*2^(DigitSum[n, 2]-1)];
%t Array[A233775, 100, 0] (* _Paolo Xausa_, Aug 05 2024 *)
%o (PARI) print1("1, "); for(k=1,70, print1((2^valuation(k,2)+1) *2^(hammingweight(k)-1),", ")) \\ _Hugo Pfoertner_, Jul 27 2020
%Y Right border gives A094373.
%Y Cf. A001316, A006046, A047999, A233774, A083318.
%K nonn,tabf
%O 0,2
%A _Omar E. Pol_, Dec 16 2013
|
