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.

Showing 1-4 of 4 results.

A210662 Triangle read by rows: T(n,k) (1 <= k <= n) = number of monomer-dimer tilings of an n X k board.

Original entry on oeis.org

1, 2, 7, 3, 22, 131, 5, 71, 823, 10012, 8, 228, 5096, 120465, 2810694, 13, 733, 31687, 1453535, 65805403, 2989126727, 21, 2356, 196785, 17525619, 1539222016, 135658637925, 11945257052321, 34, 7573, 1222550, 211351945, 36012826776, 6158217253688, 1052091957273408, 179788343101980135
Offset: 1

Views

Author

N. J. A. Sloane, Mar 28 2012

Keywords

Comments

The triangle is half of a symmetric square array, since T(n,k) = T(k,n).
Equivalently, ways of paving n X k grid cells using only singletons and dominoes. Also, the number of tilings of an n X k chessboard with the two polyominoes (0,0) and (0,0)+(0,1).
Also, matchings of the n X k grid graph. The n X k grid graph is also denoted P_m X P_n. For k=2, this is called the ladder graph L_n.
In statistical mechanics, this is a special case of the Monomer-Dimer Problem, which deals with monomer-dimer coverings of a finite patch of a lattice.
Right hand diagonal is A028420. Left hand diagonal is A000045.
Taken as a full square array, columns (and rows) 1-7 respectively match A000045(n+1), A030186, A033506(n-1), A033507(n-1), A033508(n-1), A033509(n-1), A033510(n-1), and have recurrences of order 2 3 6 9 20 36 72. - R. H. Hardin, Dec 11 2012

Examples

			Triangle begins:
1
2 7
3 22 131
5 71 823 10012
8 228 5096 120465 2810694
13 733 31687 1453535 65805403 2989126727
21 2356 196785 17525619 1539222016 135658637925 11945257052321
34 7573 1222550 211351945 36012826776 6158217253688 1052091957273408 179788343101980135...
The 7 matchings of the P(2) X P(2)-graph are:
  . .   .-.   . .   . .   . .   . .   .-.
              |       |         | |
  . .   . .   . .   . .   .-.   . .   .-.
		

References

  • Steven R. Finch, Mathematical Constants, Cambridge, 2003, pp. 406-412.
  • Per Hakan Lundow, "Computation of matching polynomials and the number of 1-factors in polygraphs", Research reports, No 12, 1996, Department of Mathematics, Umea University.

Programs

  • Sage
    from sage.combinat.tiling import TilingSolver, Polyomino
    def T(n, k):
        p = Polyomino([(0, 0)])
        q = Polyomino([(0, 0), (0, 1)])
        T = TilingSolver([p, q], box=[n, k], reusable=True)
        return T.number_of_solutions()
    # Ralf Stephan, May 22 2014

Formula

T(1,n) = A000045(n+1), T(2,n) = A030186(n), T(3,n) = A033506(n), T(4,n) = A033507(n), T(5,n) = A033508(n), T(6,n) = A033509(n), T(7,n) = A033510(n), T(8,n) = A033511(n), T(9,n) = A033512(n), T(10,n) = A033513(n), T(11,n) = A033514(n), T(n,n) = A028420(n).

Extensions

Typo in term 27 corrected by Alois P. Heinz, Dec 03 2012
Reviewed by Ralf Stephan, May 22 2014

A033507 Number of matchings in graph P_{4} X P_{n}.

Original entry on oeis.org

1, 5, 71, 823, 10012, 120465, 1453535, 17525619, 211351945, 2548684656, 30734932553, 370635224561, 4469527322891, 53898461609719, 649966808093412, 7838012982224913, 94519361817920403, 1139818186429110279, 13745178487929574337, 165754445655292452448
Offset: 0

Views

Author

Keywords

Examples

			a(1) = 5: the graph is
. o-o-o-o
and the five matchings are
. o o o o
. o-o o o
. o o-o o
. o o o-o
. o-o o-o
		

References

  • H. Hosoya and A. Motoyama, An effective algorithm for obtaining polynomials for dimer statistics. Application of operator technique on the topological index to two- and three-dimensional rectangular and torus lattices, J. Math. Phys., 26(1985), 157-167.

Crossrefs

Column 4 of triangle A210662. Row sums of A100265.
For perfect matchings see A005178.
Bisection (even part) gives A260034.

Programs

  • GAP
    a:=[1,5,71,823,10012,120465, 1453535,17525619,211351945];; for n in [10..30] do a[n]:=9*a[n-1]+41*a[n-2]-41*a[n-3]-111*a[n-4]+91*a[n-5] +29*a[n-6]-23*a[n-7]-a[n-8]+a[n-9]; od; a; # G. C. Greubel, Oct 26 2019
  • Magma
    R:=PowerSeriesRing(Integers(), 30); Coefficients(R!( (1-x)*(1 -3*x-18*x^2+2*x^3+12*x^4+x^5-x^6)/(1-9*x-41*x^2+41*x^3+111*x^4-91*x^5 -29*x^6+23*x^7+x^8-x^9) )); // G. C. Greubel, Oct 26 2019
    
  • Maple
    a:=array(0..20,[1, 5, 71, 823, 10012, 120465, 1453535, 17525619, 211351945]):
    for j from 9 to 20 do
      a[j]:=9*a[j-1]+41*a[j-2]-41*a[j-3]-111*a[j-4]+91*a[j-5]+
            29*a[j-6]-23*a[j-7]-a[j-8]+a[j-9]
    od:
    convert(a,list);
    # Sergey Perepechko, Apr 24 2013
  • Mathematica
    LinearRecurrence[{9,41,-41,-111,91,29,-23,-1,1},{1,5,71,823,10012,120465, 1453535,17525619,211351945},30] (* Harvey P. Dale, Mar 27 2015 *)
  • PARI
    my(x='x+O('x^30)); Vec((1-x)*(1 -3*x-18*x^2+2*x^3+12*x^4+x^5-x^6)/(1-9*x-41*x^2+41*x^3+111*x^4-91*x^5 -29*x^6+23*x^7+x^8-x^9)) \\ G. C. Greubel, Oct 26 2019
    
  • Sage
    def A033507_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P( (1-x)*(1 -3*x-18*x^2+2*x^3+12*x^4+x^5-x^6)/(1-9*x-41*x^2+41*x^3+111*x^4-91*x^5 -29*x^6+23*x^7+x^8-x^9) ).list()
    A033507_list(30) # G. C. Greubel, Oct 26 2019
    

Formula

From Sergey Perepechko, Apr 24 2013: (Start)
a(n) = 9*a(n-1) +41*a(n-2) -41*a(n-3) -111*a(n-4) +91*a(n-5) +29*a(n-6) -23*a(n-7) -a(n-8) +a(n-9).
G.f.: (1-x) * (1 -3*x -18*x^2 +2*x^3 +12*x^4 +x^5 -x^6) / (1 -9*x -41*x^2 +41*x^3 +111*x^4 -91*x^5 -29*x^6 +23*x^7 +x^8 -x^9). (End)

Extensions

Edited by N. J. A. Sloane, Nov 15 2009

A033506 Number of matchings in graph P_{3} X P_{n}.

Original entry on oeis.org

1, 3, 22, 131, 823, 5096, 31687, 196785, 1222550, 7594361, 47177097, 293066688, 1820552297, 11309395995, 70254767718, 436427542283, 2711118571311, 16841658983944, 104621568809247, 649916534985369, 4037327172325542
Offset: 0

Views

Author

Keywords

References

  • Miklos Bona, editor, Handbook of Enumerative Combinatorics, CRC Press, 2015, pages 50, 999.
  • Per Hakan Lundow, "Computation of matching polynomials and the number of 1-factors in polygraphs", Research reports, No 12, 1996, Department of Mathematics, Umea University.

Crossrefs

Column 3 of triangle A210662. Row sums of A100245.

Programs

  • GAP
    a:=[1,3,22,131,823,5096];; for n in [7..30] do a[n]:=4*a[n-1] +14*a[n-2]-10*a[n-4]+a[n-6]; od; a; # G. C. Greubel, Oct 26 2019
  • Magma
    R:=PowerSeriesRing(Integers(), 30); Coefficients(R!( (1-2*x-x^2)*(1+x-x^2)/((1+x)*(1-5*x-9*x^2+9*x^3+x^4-x^5)) )); // G. C. Greubel, Oct 26 2019
    
  • Maple
    seq(coeff(series((1-2*x-x^2)*(1+x-x^2)/((1+x)*(1-5*x-9*x^2+9*x^3+x^4-x^5 )), x, n+1), x, n), n = 0..30); # G. C. Greubel, Oct 26 2019
  • Mathematica
    CoefficientList[Series[(1-2x-x^2)(1+x-x^2)/((1+x)(1-5x-9x^2+9x^3+x^4-x^5) ), {x, 0, 30}], x] (* Harvey P. Dale, Dec 05 2014 *)
    LinearRecurrence[{4, 14, 0, -10, 0, 1}, {1, 3, 22, 131, 823, 5096}, 30] (* Harvey P. Dale, Dec 05 2014 *)
    Table[RootSum[-1 +# +9#^2 -9#^3 -5#^4 +#^5 &, 1436541#^n + 3941068#^(n+1) -6086452#^(n+2) -2800519#^(n+3) +591744#^(n+4) &]/10204570 -(-1)^n/5, {n, 20}] (* Eric W. Weisstein, Oct 02 2017 *)
  • PARI
    my(x='x+O('x^30)); Vec((1-2*x-x^2)*(1+x-x^2)/((1+x)*(1-5*x-9*x^2 +9*x^3+x^4-x^5))) \\ G. C. Greubel, Oct 26 2019
    
  • Sage
    def A033506_list(prec):
        P. = PowerSeriesRing(ZZ, prec)
        return P((1-2*x-x^2)*(1+x-x^2)/((1+x)*(1-5*x-9*x^2+9*x^3+x^4-x^5)) ).list()
    A033506_list(30) # G. C. Greubel, Oct 26 2019
    

Formula

G.f.: (1-2*x-x^2)*(1+x-x^2)/((1+x)*(1-5*x-9*x^2+9*x^3+x^4-x^5)). - Sergey Perepechko, Apr 19 2013

A260036 Number of configurations of the general monomer-dimer model for an 8 X 2n square lattice.

Original entry on oeis.org

1, 7573, 211351945, 6158217253688, 179788343101980135, 5249581929453966097649, 153282416794739031814924079, 4475693398084103779028604877129, 130685784089846859975392644067590259, 3815894599555857747953551843550514992237
Offset: 0

Views

Author

N. J. A. Sloane, Jul 19 2015

Keywords

Crossrefs

Bisection (even part) of A033511.

Programs

Extensions

a(0), a(5)-a(9) from Alois P. Heinz, Mar 08 2016
Showing 1-4 of 4 results.