|
%I #19 Aug 28 2023 11:31:35
%S 1,10,10,1100,10110100,10101000,1010011101011000,
%T 11001010101011110011010101010000,10111010011010101101010101010000,
%U 1010110010110101010110011001111101010011010010101010011001100000,1010101110101100101010010110101011010010101101010110010110100000
%N a(n) is the periodic part on the n-th diagonal from the right of rule-30 1-D cellular automaton, when started from a single ON cell.
%C As found by Brunnbauer (2019), if a period doubling occurs at n, then a(n) is of the form AB, where B is the inverse of A. Additionally, the number of trailing zeros of a(n) increases by one when n is even.
%H Paolo Xausa, <a href="/A364773/b364773.txt">Table of n, a(n) for n = 1..25</a>
%H Michael Brunnbauer, <a href="https://brunni.de/findings30/">Diagonals in elementary cellular automaton 30</a>, 2019 (<a href="/A364773/a364773.pdf">local PDF copy</a>, with author's permission).
%H Eric S. Rowland, <a href="https://wpmedia.wolfram.com/uploads/sites/13/2018/02/16-3-4.pdf">Local Nested Structure in Rule 30</a>, Complex Systems 16 (2006), pp. 239-258.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/Rule30.html">Rule 30</a>.
%H Stephen Wolfram, <a href="https://www.wolframscience.com/nks/notes-2-1--rule-30/">Notes on chapter 2, Rule 30</a>, from A new kind of science online, Wolfram Media, 2002.
%H <a href="/index/Ce#cell">Index entries for sequences related to cellular automata</a>
%e In the following diagram, showing the first 20 evolution steps of the CA, two diagonals are highlighted (the rest of the CA is represented by hyphens, for better visualization).
%e .
%e 2nd diagonal
%e a(2) = 10 __
%e \-
%e 7th diagonal __ -1-
%e a(7) = 1010011101011000 \ ---0-
%e 1----1-
%e --0----0-
%e ----1----1-
%e ------0----0-
%e --------0----1-
%e ----------1----0-
%e ------------1----1-
%e --------------1----0-
%e ----------------0----1-
%e ------------------1----0-
%e --------------------0----1-
%e ----------------------1----0-
%e ------------------------1----1-
%e --------------------------0----0-
%e ----------------------------0----1-
%e ------------------------------0----0-
%e --------------------------------1----1-
%e ----------------------------------0----0-
%e .
%t A364773list[steps_]:=Module[{d=2Ceiling[Log2[steps]],ca,n=1,p,plen,a={1}},ca=CellularAutomaton[30,{{1},0},{steps,{1-d,steps}}];While[++n<=2(d-1)&&(plen=Length[p=FindRepeat[Flatten[Rest[Split[Diagonal[ca,d-n]]]]]])>=IntegerLength[Last[a]]&&IntegerQ[Log2[plen]],AppendTo[a,FromDigits[p]]];a];
%t A364773list[80] (* Analyzes 80 evolution steps *)
%Y Cf. A070950, A094605 (periods), A363343 (diagonals), A364774 (base 10).
%K nonn
%O 1,2
%A _Paolo Xausa_, Aug 06 2023
| 