login
Complex cross sequence (see Comments lines for definition).
11

%I #91 Nov 11 2019 00:47:12

%S 0,1,3,11,19,27,35,67,83,99,115,163,179,195,211,275,323,355,387,467,

%T 483,499,515,579,627,675,707,787,803,819,835,899,947,995,1027,1107,

%U 1123,1139,1155,1219,1267,1315,1347,1427,1443,1459,1475,1539,1587,1635,1667,1747,1763,1779,1795,1859,1907,1955,1987,2067

%N Complex cross sequence (see Comments lines for definition).

%C The sequence arises from a "hybrid" cellular automaton on the infinite square grid, which consist of two successive generations using toothpicks of length 2 (cf. A139250) followed by two successive generations using the rules of the D-toothpick sequence A220500.

%C In other words (and more precisely) we have that:

%C 1) If n is congruent to 1 or 2 mod 4 (cf. A042963), for example: 1, 2, 5, 6, 9, 10, ..., the elements added to the structure at n-th stage must be toothpicks of length 2. These toothpicks are connected to the structure by their midpoints.

%C 2) If n is a positive integer of the form 4*k-1 (cf. A004767), for example: 3, 7, 11, 15, ..., the elements added to the structure at n-th stage must be D-toothpicks of length sqrt(2) and eventually D-toothpicks of length sqrt(2)/2, in both cases the D-toothpicks are connected to the structure by their endpoints, in the same way as in the even-indexed stages of A220500.

%C 3) If n is a positive multiple of 4 (cf. A008586) the elements added to the structure at n-th stage must be toothpicks of length 1 connected by their endpoints, in the same way as in the odd-indexed stages of A220500.

%C a(n) is the total number of elements in the structure after n generations.

%C A289841 (the first differences) gives the number of elements added at n-th stage.

%C Note that after 19 generations the structure is a 72-gon which essentially looks like a diamond (as a square that has been rotated 45 degrees).

%C The surprising fact is that from n = 20 up to 27 the structure is gradually transformed into a square cross.

%C The diamond mentioned above can be interpreted as the center of the cross. The diamond has an area equal to 384 and it contains 222 polygonal regions (or enclosures) of 11 distinct shapes. Missing two heptagonal shapes which are in the arms of the square cross only.

%C In total the complex cross contains 13 distinct shapes of polygonal regions. There are ten polygonal shapes that have an infinite number of copies. On the other hand, three of these polygonal shapes have a finite number of copies because they are in the center of the cross only. For example: there are only four copies of the concave 14-gon, which is also the largest polygon in the structure.

%C For n => 27 the shape of the square cross remains forever because its four arms grow indefinitely.

%C Every arm has a minimum width equal to 8, and a maximum width equal to 12.

%C Every arm also has a periodic structure which can be dissected in infinitely many clusters of area equal to 64. Every cluster is a 30-gon that contains 40 polygonal regions of nine distinct shapes.

%C If n is a number of the form 8*k-3 (cf. A017101) and greater than 19, for example: 27, 35, 43, 51, ..., then at n-th stage a new cluster is finished in every arm of the cross.

%C The behavior is similar to A290220 and A294020 in the sense that these three cellular automata have the property of self-limiting their growth only in some directions of the square grid. - _Omar E. Pol_, Oct 29 2017

%H Colin Barker, <a href="/A289840/b289840.txt">Table of n, a(n) for n = 0..1000</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 <a href="/index/Ce#cell">Index entries for sequences related to cellular automata</a>

%H <a href="/index/Rec#order_09">Index entries for linear recurrences with constant coefficients</a>, signature (1,0,0,0,0,0,0,1,-1).

%F From _Colin Barker_, Nov 11 2017: (Start)

%F G.f.: x*(1 + 2*x + 8*x^2 + 8*x^3 + 8*x^4 + 8*x^5 + 32*x^6 + 16*x^7 + 15*x^8 + 14*x^9 + 40*x^10 + 8*x^11 + 8*x^12 + 8*x^13 + 32*x^14 + 32*x^15 + 16*x^16 + 16*x^17 + 32*x^18 + 16*x^24) / ((1 - x)^2*(1 + x)*(1 + x^2)*(1 + x^4)).

%F a(n) = a(n-1) + a(n-8) - a(n-9) for n>19.

%F (End)

%o (PARI) concat(0, Vec(x*(1 + 2*x + 8*x^2 + 8*x^3 + 8*x^4 + 8*x^5 + 32*x^6 + 16*x^7 + 15*x^8 + 14*x^9 + 40*x^10 + 8*x^11 + 8*x^12 + 8*x^13 + 32*x^14 + 32*x^15 + 16*x^16 + 16*x^17 + 32*x^18 + 16*x^24) / ((1 - x)^2*(1 + x)*(1 + x^2)*(1 + x^4))+ O(x^50))) \\ _Colin Barker_, Nov 12 2017

%Y Cf. A004767, A008586, A017101, A042963, A139250, A220500, A289841, A290220 (a simpler cross), A294020.

%K nonn,easy

%O 0,3

%A _Omar E. Pol_, Jul 14 2017