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.

A034345 Number of binary [ n,4 ] codes without 0 columns.

Original entry on oeis.org

0, 0, 0, 1, 4, 11, 27, 63, 134, 276, 544, 1048, 1956, 3577, 6395, 11217, 19307, 32685, 54413, 89225, 144144, 229647, 360975, 560259, 858967, 1301757, 1950955, 2893102, 4246868, 6174084, 8892966, 12696295, 17973092, 25237467, 35163431, 48629902, 66774760, 91063984
Offset: 1

Views

Author

Keywords

Comments

"We say that the sequence (a_n) is quasi-polynomial in n if there exist polynomials P_0, ..., P_{s-1} and an integer n_0 such that, for all n >= n_0, a_n = P_i(n) where i == n (mod s)." [This is from the abstract of Lisonek (2007), and he states that the condition "n >= n_0" makes his definition broader than the one in Stanley's book. From Section 5 of his paper, we conclude that (a(n): n >= 1) is a quasi-polynomial in n.] - Petros Hadjicostas, Oct 02 2019

Crossrefs

Column k=4 of A034253 and first differences of A034358.

Programs

  • Sage
    # Fripertinger's method to find the g.f. of column k >= 2 of A034253 (for small k):
    def A034253col(k, length):
        G1 = PSL(k, GF(2))
        G2 = PSL(k-1, GF(2))
        D1 = G1.cycle_index()
        D2 = G2.cycle_index()
        f1 = sum(i[1]*prod(1/(1-x^j) for j in i[0]) for i in D1)
        f2 = sum(i[1]*prod(1/(1-x^j) for j in i[0]) for i in D2)
        f = f1 - f2
        return f.taylor(x, 0, length).list()
    # For instance the Taylor expansion for column k = 4 (this sequence) gives
    print(A034253col(4, 30)) #

Extensions

More terms by Petros Hadjicostas, Oct 02 2019