automaton

automaton.py

from functools import reduce
from operator import getitem

def bst(n, i):
    if n % 2 == 0:
        h = 2 ** (n / 2)
        remainder, quotient = i % h, i // h
        return bst(n / 2, remainder), bst(n / 2, quotient)
    else:
        return i

class Rule:
    def __init__(self, rule):
        self.bst = bst(8, rule)
        
    def check(self, state):
        return reduce(getitem, state, self.bst)

class Layer(list):
    def __getitem__(self, key):
        if isinstance(key, slice):
            return [self.__getitem__(i) for i in range(key.start, key.stop)]
        elif key >= 0 and key < len(self):
            return list.__getitem__(self, key)
        else:
            return 0
    
first_layer = Layer([1])