[go: up one dir, main page]

login
A269175
a(n) = number of distinct k for which A269174(k) = n; number of finite predecessors for pattern encoded in the binary expansion of n in Wolfram's Rule 124 cellular automaton.
6
1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 1, 2, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 2, 1, 0, 2, 2, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 0, 2, 0, 0, 2, 1, 1
OFFSET
0,32
COMMENTS
At positions A000225 seems to occur the record values of this sequence: 1, 0, 1, 1, 1, 2, 2, 3, 4, 5, 7, 9, 12, 16, 21, 28, 37, 49, 65, 86, 114, 151, 200, 265, 351, 465, ... which seem to match with A000931 (Padovan sequence), or more exactly, with A182097 (Number of compositions (ordered partitions) into parts 2 and 3). Note that these values give also the number of predecessors for each "repunit-pattern" (2^n)-1 in Rule 110 cellular automaton, as rules 110 and 124 are mirror images of each other.
PROG
(Scheme)
(definec (A269175 n) (let loop ((p 0) (s 0)) (cond ((> p n) s) (else (loop (+ 1 p) (+ s (if (= n (A269174 p)) 1 0))))))) ;; Very straightforward and very slow.
;; Somewhat optimized version:
(definec (A269175 n) (if (zero? n) 1 (let ((nwid-1 (- (A000523 n) 1))) (let loop ((p (if (< n 2) 0 (A000079 nwid-1))) (s 0)) (cond ((> (A000523 p) nwid-1) s) (else (loop (+ 1 p) (+ s (if (= n (A269174 p)) 1 0)))))))))
CROSSREFS
Cf. A269176 (indices of zeros), A269177 (of nonzeros), A269178 (of ones).
Sequence in context: A299173 A351726 A285720 * A086076 A334348 A316717
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 22 2016
STATUS
approved