cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A107069 Number of self-avoiding walks of length n on an infinite triangular prism starting at the origin.

Original entry on oeis.org

1, 4, 12, 34, 90, 222, 542, 1302, 3058, 7186, 16714, 38670, 89358, 205710, 472906, 1086138, 2491666, 5713318, 13094950, 30003190, 68731010, 157423986, 360530346, 825626942, 1890615518, 4329196974, 9912914314, 22698017834, 51972012258, 119000208806
Offset: 0

Views

Author

Jonathan Vos Post, May 10 2005

Keywords

Comments

The discrete space in which the walking happens is a triangular prism infinite in both directions along the x-axis. One vertex is the root, the origin. The basis is the set of single-step vectors, which we abbreviate as l (left), r (right), c (one step "clockwise" around the triangle) and c- (one step counterclockwise, more properly denoted c^-1).

Examples

			a(0) = 1, as there is one self-avoiding walk of length 0, namely the null-walk (the walk whose steps are the null set).
a(1) = 4 because (using the terminology in the Comment), the 4 possible 1-step walks are W_1 = {l,r,c,c-}.
a(2) = 12 because the set of legal 2-step walks are {l^2, lc, lc-, r^2, rc, rc-, c^2, cl, cr, c^-2, c-l, c-r}.
a(3) = 34 because we have every W_2 concatenated with {l,r,c,c-} except for those with immediate violations (lr etc.) and those two which go in a triangle {c^3, c^-3}; hence a(3) = 3*a(2) - 2 = 3*12 - 2 = 36 - 2 = 34.
		

Crossrefs

Programs

  • Python
    w = [[[(0, 0)]]]
    for n in range(1, 15):
        nw = []
        for walk in w[-1]:
            (x, t) = walk[-1]
            nss = [(x-1, t), (x+1, t), (x, (t+1)%3), (x, (t-1)%3)]
            for ns in nss:
                if ns not in walk:
                    nw.append(walk[:] + [ns])
        w.append(nw)
    print([len(x) for x in w])
    # Andrey Zabolotskiy, Sep 19 2019

Extensions

a(4) and a(5) corrected, a(6)-a(14) added by Andrey Zabolotskiy, Sep 19 2019
More terms from Andrey Zabolotskiy, Dec 04 2023