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-6 of 6 results.

A078757 Values of A028470(n)/A000045(n+1).

Original entry on oeis.org

1, 17, 51, 449, 1853, 12853, 61557, 382024, 1971559, 11585969, 62088471, 355111613, 1939427729, 10943439733, 60338602299, 338172377293, 1873494595072, 10464657396101, 58113694771149, 324052035315389, 1801727076022631, 10038214290617749, 55845947547948897
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Mathematica
    CoefficientList[Series[-(x^7 + x^6 - 24 x^5 + 11 x^4 + 38 x^3 - 9 x^2 - 16 x - 1) / (x^8 + x^7 - 25 x^6 + 11 x^5 + 47 x^4 - 11 x^3 - 25 x^2 - x + 1), {x, 0, 30}], x] (* Vincenzo Librandi, Nov 13 2014 *)
    LinearRecurrence[{1,25,11,-47,-11,25,-1,-1},{1,17,51,449,1853,12853,61557,382024},30] (* Harvey P. Dale, Jul 15 2017 *)

Formula

a(n) = a(n-1)+25a(n-2)+11a(n-3)-47a(n-4)-11a(n-5)+25a(n-6)-a(n-7)-a(n-8).
G.f.: -x*(x^7+x^6-24*x^5+11*x^4+38*x^3-9*x^2-16*x-1)/(x^8+x^7-25*x^6+11*x^5+47*x^4-11*x^3-25*x^2-x+1). [Colin Barker, Jun 22 2012]

Extensions

More terms from Vincenzo Librandi, Nov 13 2014

A099390 Array T(m,n) read by antidiagonals: number of domino tilings (or dimer tilings) of the m X n grid (or m X n rectangle), for m>=1, n>=1.

Original entry on oeis.org

0, 1, 1, 0, 2, 0, 1, 3, 3, 1, 0, 5, 0, 5, 0, 1, 8, 11, 11, 8, 1, 0, 13, 0, 36, 0, 13, 0, 1, 21, 41, 95, 95, 41, 21, 1, 0, 34, 0, 281, 0, 281, 0, 34, 0, 1, 55, 153, 781, 1183, 1183, 781, 153, 55, 1, 0, 89, 0, 2245, 0, 6728, 0, 2245, 0, 89, 0, 1, 144, 571, 6336, 14824, 31529, 31529, 14824, 6336, 571, 144, 1
Offset: 1

Views

Author

Ralf Stephan, Oct 16 2004

Keywords

Comments

There are many versions of this array (or triangle) in the OEIS. This is the main entry, which ideally collects together all the references to the literature and to other versions in the OEIS. But see A004003 for further information. - N. J. A. Sloane, Mar 14 2015

Examples

			0,  1,  0,   1,    0,    1, ...
1,  2,  3,   5,    8,   13, ...
0,  3,  0,  11,    0,   41, ...
1,  5, 11,  36,   95,  281, ...
0,  8,  0,  95,    0, 1183, ...
1, 13, 41, 281, 1183, 6728, ...
		

References

  • S. R. Finch, Mathematical Constants, Cambridge, 2003, pp. 406-412.
  • P. E. John, H. Sachs, and H. Zernitz, Problem 5. Domino covers in square chessboards, Zastosowania Matematyki (Applicationes Mathematicae) XIX 3-4 (1987), 635-641.
  • R. P. Stanley, Enumerative Combinatorics, Vol. 1, Cambridge University Press, 2nd ed., pp. 547 and 570.
  • Darko Veljan, Kombinatorika: s teorijom grafova (Croatian) (Combinatorics with Graph Theory) mentions the value 12988816 = 2^4*901^2 for the 8 X 8 case on page 4.

Crossrefs

See A187596 for another version (with m >= 0, n >= 0). See A187616 for a triangular version. See also A187617, A187618.
See also A004003 for more literature on the dimer problem.
Main diagonal is A004003.

Programs

  • Maple
    (Maple code for the even-numbered rows from N. J. A. Sloane, Mar 15 2015. This is not totally satisfactory since it uses floating point. However, it is useful for getting the initial values quickly.)
    Digits:=100;
    p:=evalf(Pi);
    z:=proc(h,d) global p; evalf(cos( h*p/(2*d+1) )); end;
    T:=proc(m,n) global z; round(mul( mul( 4*z(h,m)^2+4*z(k,n)^2, k=1..n), h=1..m)); end;
    [seq(T(1,n),n=0..10)]; # A001519
    [seq(T(2,n),n=0..10)]; # A188899
    [seq(T(3,n),n=0..10)]; # A256044
    [seq(T(n,n),n=0..10)]; # A004003
  • Mathematica
    T[?OddQ, ?OddQ] = 0;
    T[m_, n_] := Product[2*(2+Cos[2j*Pi/(m+1)]+Cos[2k*Pi/(n+1)]), {k, 1, n/2}, {j, 1, m/2}];
    Flatten[Table[Round[T[m-n+1, n]], {m, 1, 12}, {n, 1, m}]] (* Jean-François Alcover, Nov 25 2011, updated May 28 2022 *)
  • PARI
    {T(n, k) = sqrtint(abs(polresultant(polchebyshev(n, 2, x/2), polchebyshev(k, 2, I*x/2))))} \\ Seiichi Manyama, Apr 13 2020

Formula

T(m, n) = Product_{j=1..ceiling(m/2)} Product_{k=1..ceiling(n/2)} (4*cos(j*Pi/(m+1))^2 + 4*cos(k*Pi/(n+1))^2).

Extensions

Old link fixed and new link added by Frans J. Faase, Feb 04 2009
Entry edited by N. J. A. Sloane, Mar 15 2015

A005178 Number of domino tilings of 4 X (n-1) board.

Original entry on oeis.org

0, 1, 1, 5, 11, 36, 95, 281, 781, 2245, 6336, 18061, 51205, 145601, 413351, 1174500, 3335651, 9475901, 26915305, 76455961, 217172736, 616891945, 1752296281, 4977472781, 14138673395, 40161441636, 114079985111, 324048393905
Offset: 0

Views

Author

N. J. A. Sloane, David Singmaster, Frans J. Faase

Keywords

Comments

Or, number of perfect matchings in graph P_4 X P_{n-1}.
a(0) = 0, a(1) = 1 by convention.
It is easy to see that the g.f. for indecomposable tilings, i.e., those that cannot be split vertically into smaller tilings, is g = x + 4x^2 + 2x^3 + 3x^4 + 2x^5 + 3x^6 + 2x^7 + 3x^8 + ... = x + 4x^2 + x^3*(2+3x)/(1-x^2); then g.f. = 1/(1-g) = (1-x^2)/(1-x-5x^2-x^3+x^4). - Emeric Deutsch, Oct 16 2006
This is a divisibility sequence; that is, if n divides m, then a(n) divides a(m). - T. D. Noe, Dec 22 2008
From Artur Jasinski, Dec 20 2008: (Start)
All numbers in this sequence are:
congruent to 0 mod 100 if n is congruent to 14 or 29 mod 30
congruent to 1 mod 100 if n is congruent to 0 or 1 or 12 or 16 or 27 or 28 mod 30
congruent to 5 mod 100 if n is congruent to 2 or 11 or 17 or 26 mod 30
congruent to 11 mod 100 if n is congruent to 3 or 25 mod 30
congruent to 36 mod 100 if n is congruent to 4 or 9 or 19 or 24 mod 30
congruent to 45 mod 100 if n is congruent to 8 or 20 mod 30
congruent to 51 mod 100 if n is congruent to 13 or 15 mod 30
congruent to 61 mod 100 if n is congruent to 10 or 18 mod 30
congruent to 81 mod 100 if n is congruent to 6 or 7 or 21 or 22 mod 30
congruent to 95 mod 100 if n is congruent to 5 or 23 mod 30
(End)
This is the case P1 = 1, P2 = -7, Q = 1 of the 3 parameter family of 4th-order linear divisibility sequences found by Williams and Guy. - Peter Bala, Mar 31 2014

Examples

			For n=2 the graph is
  o-o-o-o
and there is one perfect tiling:
  o-o o-o
For n=3 the graph is
  o-o-o-o
  | | | |
  o-o-o-o
and there are five perfect tilings:
  o o o o
  | | | |
  o o o o
two like:
  o o o-o
  | | ...
  o o o-o
and this
  o-o o-o
  .......
  o-o o-o
and this
  o o-o o
  | ... |
  o o-o o
a(n+1)=r(n)-r(n-2), r(n)=if n=0 then 1 else sum(sum(binomial(k,j)*sum(binomial(j,i-j)*5^(i-j)*binomial(k-j,n-i-3*(k-j))*(-1)^(n-i-3*(k-j)),i,j,n-k+j),j,0,k),k,1,n), n>1. - _Vladimir Kruchinin_, Sep 08 2010
		

References

  • F. Faase, On the number of specific spanning subgraphs of the graphs G X P_n, Ars Combin. 49 (1998), 129-154.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • R. P. Stanley, Enumerative Combinatorics I, p. 292.

Crossrefs

Row 4 of array A099390.
For all matchings see A033507.
Cf. A003757. - T. D. Noe, Dec 22 2008
Bisection (odd part) gives A188899. - Alois P. Heinz, Oct 28 2012
Column k=2 of A250662.

Programs

  • Maple
    a[0]:=1: a[1]:=1: a[2]:=5: a[3]:=11: for n from 4 to 26 do a[n]:=a[n-1]+5*a[n-2]+a[n-3]-a[n-4] od: seq(a[n],n=0..26); # Emeric Deutsch, Oct 16 2006
    A005178:=-(-1-4*z-z**2+z**3)/(1-z-5*z**2-z**3+z**4) # conjectured (correctly) by Simon Plouffe in his 1992 dissertation; gives sequence apart from an initial 1
  • Mathematica
    CoefficientList[Series[x(1-x^2)/(1-x-5x^2-x^3+x^4), {x,0,30}], x] (* T. D. Noe, Dec 22 2008 *)
    LinearRecurrence[{1, 5, 1, -1}, {0, 1, 1, 5}, 28] (* Robert G. Wilson v, Aug 08 2011 *)
    a[0] = 0; a[n_] := Product[2(2+Cos[2j Pi/5]+Cos[2k Pi/n]), {k, 1, (n-1)/2}, {j, 1, 2}] // Round;
    Table[a[n], {n, 0, 27}] (* Jean-François Alcover, Aug 20 2018 *)
  • Maxima
    r(n):=if n=0 then 1 else sum(sum(binomial(k,j)*sum(binomial(j,i-j)*5^(i-j)*binomial(k-j,n-i-3*(k-j))*(-1)^(n-i-3*(k-j)),i,j,n-k+j),j,0,k),k,1,n); a(n):=r(n)-r(n-2); /* Vladimir Kruchinin, Sep 08 2010 */

Formula

a(n) = a(n-1) + 5*a(n-2) + a(n-3) - a(n-4).
G.f.: x*(1 - x^2)/(1 - x - 5*x^2 - x^3 + x^4).
Limit_{n->oo} a(n)/a(n-1) = (1 + sqrt(29) + sqrt(14 + 2*sqrt(29)))/4 = 2.84053619409... - Philippe Deléham, Jun 12 2005
a(n) = (5*sqrt(29)/145)*(((1+sqrt(29)+sqrt(14+2*sqrt(29)))/4)^n+((1+sqrt(29)-sqrt(14+2*sqrt(29)))/4)^n-((1-sqrt(29)+sqrt(14-2*sqrt(29)))/4)^n-((1-sqrt(29)-sqrt(14-2*sqrt(29)))/4)^n). - Tim Monahan, Jul 30 2011
From Peter Bala, Mar 31 2014: (Start)
a(n) = ( T(n,alpha) - T(n,beta) )/(alpha - beta), where alpha = (1 + sqrt(29))/4 and beta = (1 - sqrt(29))/4 and T(n,x) denotes the Chebyshev polynomial of the first kind.
a(n) = the bottom left entry of the 2 X 2 matrix T(n, M), where M is the 2 X 2 matrix [0, 7/4; 1, 1/2].
a(n) = U(n-1,i*(1 + sqrt(5))/4)*U(n-1,i*(1 - sqrt(5))/4), where U(n,x) denotes the Chebyshev polynomial of the second kind.
See the remarks in A100047 for the general connection between Chebyshev polynomials and 4th-order linear divisibility sequences. (End)
a(n) = A129113(n+2) - A129113(n). - R. J. Mathar, May 03 2021

Extensions

Amalgamated with (former) A003692, Dec 30 1995
Name changed and 0 prepended by T. D. Noe, Dec 22 2008
Edited by N. J. A. Sloane, Nov 15 2009

A189006 Array A(m,n) read by antidiagonals: number of domino tilings of the m X n grid with upper left corner removed iff m*n is odd, (m>=0, n>=0).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 3, 3, 1, 1, 1, 1, 5, 4, 5, 1, 1, 1, 1, 8, 11, 11, 8, 1, 1, 1, 1, 13, 15, 36, 15, 13, 1, 1, 1, 1, 21, 41, 95, 95, 41, 21, 1, 1, 1, 1, 34, 56, 281, 192, 281, 56, 34, 1, 1, 1, 1, 55, 153, 781, 1183, 1183, 781, 153, 55, 1, 1, 1, 1, 89, 209, 2245, 2415, 6728, 2415, 2245, 209, 89, 1, 1
Offset: 0

Views

Author

Alois P. Heinz, Apr 15 2011

Keywords

Examples

			A(3,3) = 4, because there are 4 domino tilings of the 3 X 3 grid with upper left corner removed:
  . .___. . .___. . .___. . .___.
  ._|___| ._|___| ._| | | ._|___|
  | |___| | | | | | |_|_| |___| |
  |_|___| |_|_|_| |_|___| |___|_|
Array begins:
  1, 1,  1,  1,   1,    1,    1, ...
  1, 1,  1,  1,   1,    1,    1, ...
  1, 1,  2,  3,   5,    8,   13, ...
  1, 1,  3,  4,  11,   15,   41, ...
  1, 1,  5, 11,  36,   95,  281, ...
  1, 1,  8, 15,  95,  192, 1183, ...
  1, 1, 13, 41, 281, 1183, 6728, ...
		

Crossrefs

Rows m=0+1, 2-12 give: A000012, A000045(n+1), A002530(n+1), A005178(n+1), A189003, A028468, A189004, A028470, A189005, A028472, A210724, A028474.
Main diagonal gives: A189002.

Programs

  • Maple
    with(LinearAlgebra):
    A:= proc(m, n) option remember; local i, j, s, t, M;
          if m=0 or n=0 then 1
        elif m1 or j>1 or s=0 then
                   if j
    				
  • Mathematica
    A[1, 1] = 1; A[m_, n_] := A[m, n] = Module[{i, j, s, t, M}, Which[m == 0 || n == 0, 1, m < n, A[n, m], True, s = Mod[n*m, 2];M[i_, j_] /; j < i := -M[j, i]; M[, ] = 0; For[i = 1, i <= n, i++, For[j = 1, j <= m, j++, t = (i-1)*m+j-s; If[i > 1 || j > 1 || s == 0, If[j < m, M[t, t+1] = 1]; If[i < n, M[t, t+m] = 1-2*Mod[j, 2]]]]]; Sqrt[Det[Array[M, {n*m-s, n*m-s}]]]]]; Table[Table[A[m, d-m], {m, 0, d}], {d, 0, 15}] // Flatten (* Jean-François Alcover, Dec 26 2013, translated from Maple *)

A028471 Number of perfect matchings (or domino tilings) in the graph P_9 X P_2n.

Original entry on oeis.org

1, 55, 6336, 817991, 108435745, 14479521761, 1937528668711, 259423766712000, 34741645659770711, 4652799879944138561, 623139489426439754945, 83456125990631342400791, 11177167872295392172767936, 1496943834332592837945956455, 200483802581126644843760725601
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

  • Mathematica
    T[?OddQ, ?OddQ] = 0;
    T[m_, n_] := Product[2(2+Cos[2 j Pi/(m+1)]+Cos[2 k Pi/(n+1)]), {k, 1, n/2}, {j, 1, m/2}];
    a[n_] := T[2n, 9] // Round;
    Table[a[n], {n, 0, 20}] (* Jean-François Alcover, May 28 2022 *)
  • PARI
    {a(n) = sqrtint(polresultant(polchebyshev(2*n, 2, x/2), polchebyshev(9, 2, I*x/2)))} \\ Seiichi Manyama, Apr 13 2020

Formula

a(n) = 209*a(n-1) - 11936*a(n-2) + 274208*a(n-3) - 3112032*a(n-4) + 19456019*a(n-5) - 70651107*a(n-6) + 152325888*a(n-7) - 196664896*a(n-8) + 152325888*a(n-9) - 70651107*a(n-10) + 19456019*a(n-11) - 3112032*a(n-12) + 274208*a(n-13) - 11936*a(n-14) + 209*a(n-15) - a(n-16). - Jay Anderson (horndude77(AT)gmail.com), Apr 07 2007
G.f.: (1 - 154x + 6777x^2 - 123961x^3 + 1132714x^4 - 5684515x^5 + 16401668x^6 - 27757938x^7 + 27757938*x^8 - 16401668x^9 + 5684515x^10 - 1132714x^11 + 123961x^12 -6777x^13 + 154x^14 - x^15)/(1 - 209x + 11936x^2 - 274208x^3 + 3112032x^4 - 19456019x^5 + 70651107x^6 - 152325888x^7 + 196664896x^8 - 152325888x^9 + 70651107x^10 -19456019x^11 + 3112032x^12 - 274208x^13 + 11936x^14 - 209x^15 + x^16). - Sergey Perepechko, Nov 23 2012

Extensions

Edited by N. J. A. Sloane, Jul 03 2008 at the suggestion of R. J. Mathar

A241908 Number of perfect matchings in graph P_{13} X P_{2n}.

Original entry on oeis.org

1, 377, 413351, 536948224, 731164253833, 1012747193318519, 1412218550274852671, 1974622635952709613247, 2764079753958605286860951, 3870940598132705729413670953, 5422065916132126528319352874496, 7595338059193606161156363370300487, 10640045682768766172108553992086690201
Offset: 0

Views

Author

Sergey Perepechko, May 01 2014

Keywords

Comments

In Karavaev and Perepechko generating functions G_m(x) for P_m X P_n graphs were found for all values of m up to 27.

References

  • A. M. Karavaev and S. N. Perepechko, Generating functions for dimer problem on rectangular lattices (in Russian), Information Processes, 13(2013), No4, 374-400.

Crossrefs

Row 13 of array A099390.

Programs

  • PARI
    {a(n) = sqrtint(polresultant(polchebyshev(2*n, 2, x/2), polchebyshev(13, 2, I*x/2)))} \\ Seiichi Manyama, Apr 13 2020
Showing 1-6 of 6 results.