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.

A188496 a(n) = A188492(n+1) - A188495(n) - A002527(n).

Original entry on oeis.org

0, 0, 2, 4, 10, 28, 96, 304, 928, 2784, 8504, 26124, 80228, 245544, 751168, 2299184, 7040986, 21561028, 66015398, 202114264, 618817376, 1894692160, 5801169248, 17761879056, 54382725520, 166507388264, 509808051944, 1560917463152, 4779176035680
Offset: 0

Views

Author

N. J. A. Sloane, Apr 01 2011

Keywords

Comments

For n >= 2, a(n) is the number of permutations p on the set [n] with the properties that abs(p(i)-i) <= 3 for all i, p(1) <= 2, p(2) <= 4, and p(4) >= 2.
For n >= 2, a(n) is also the permanent of the n X n matrix that has ones on its diagonal, ones on its three superdiagonals (with the exception of a zero in the (1,4)-entry), ones on its three subdiagonals (with the exception of zeros in the (3,1), (4,1), and (5,2)-entries), and is zero elsewhere.
This is row 11 of Kløve's Table 3.

Programs

  • Maple
    with(LinearAlgebra):
    A188496:= n-> `if`(n<=1, 0, Permanent(Matrix(n, (i, j)->
                  `if`(abs(j-i)<4 and [i, j]<>[1, 4] and [i, j]<>[3, 1] and [i, j]<>[4, 1] and [i, j]<>[5, 2], 1, 0)))):
    seq(A188496(n), n=0..20);
  • Mathematica
    a[n_] := Permanent[Table[If[Abs[j - i] < 4 && {i, j} != {1, 4} && {i, j} != {3, 1} && {i, j} != {4, 1} && {i, j} != {5, 2}, 1, 0], {i, 1, n}, {j, 1, n}] ]; a[1] = 0; Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 0, 20}] (* Jean-François Alcover, Jan 06 2016, adapted from Maple *)
    LinearRecurrence[{1,3,3,13,21,19,3,-7,-9,-5,-3,-3,-1},{0,0,2,4,10,28,96,304,928,2784,8504,26124,80228},30] (* Harvey P. Dale, Aug 31 2016 *)
  • PARI
    concat([0,0], Vec(x^2*(2*x +2)/(x^13 +3*x^12 +3*x^11 +5*x^10 +9*x^9 +7*x^8 -3*x^7 -19*x^6 -21*x^5 -13*x^4 -3*x^3 -3*x^2 -x +1)  + O(x^100))) \\ Colin Barker, Dec 13 2014

Formula

a(n) = A002527(n-1) + A188495(n-1). - Nathaniel Johnston, Apr 11 2011
G.f.: x^2*(2*x +2)/(x^13 +3*x^12 +3*x^11 +5*x^10 +9*x^9 +7*x^8 -3*x^7 -19*x^6 -21*x^5 -13*x^4 -3*x^3 -3*x^2 -x +1). - Colin Barker, Dec 13 2014

Extensions

Name and comments edited, and a(12)-a(28) from Nathaniel Johnston, Apr 11 2011

A002529 a(n) = A002527(n+1) - A002527(n) - A002526(n).

Original entry on oeis.org

0, 0, 2, 6, 18, 46, 146, 460, 1436, 4352, 13252, 40532, 124396, 381140, 1166708, 3570684, 10932274, 33475170, 102499334, 313825690, 960844358, 2941873064, 9007393480, 27578681888, 84439657768, 258534813320, 791574775192, 2423623112104, 7420586212184, 22720153701768, 69563959091138
Offset: 0

Views

Author

Keywords

Comments

For n >= 2, a(n) is the number of permutations p on the set [n] with the properties that abs(p(i)-i) <= 3 for all i and p(j) <= 2+j for j = 1,2.
For n >= 2, a(n) is also the permanent of the n X n matrix that has ones on its diagonal, ones on its three superdiagonals, ones on its three subdiagonals (with the exception of zeros in the (4,1) and (5,2)-entries), and is zero elsewhere.
This is row 3 of Kløve's Table 3.

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Programs

  • Maple
    with (LinearAlgebra):
    A002529:= n-> `if` (n<=1, 0, Permanent (Matrix (n, (i, j)->
                  `if` (abs(j-i)<4 and [i, j]<>[4, 1] and [i, j]<>[5, 2], 1, 0)))):
    seq (A002529(n), n=0..20);
  • Mathematica
    a[n_] := If [n <= 1, 0, Permanent[Table[If[Abs[j-i]<4 && {i, j} != {4, 1} && {i, j} != {5, 2}, 1, 0], {i, 1, n}, {j, 1, n}]]]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Mar 12 2014, after Maple *)

Formula

a(n) = A188379(n+1) - A188492(n) - A188493(n). - Nathaniel Johnston, Apr 08 2011
G.f.: 2*x^2 * (x^4+x^3-x^2-x-1) / (x^14+2*x^13+2*x^11 +4*x^10 -2*x^9 -10*x^8 -16*x^7-2*x^6 +8*x^5+10*x^4 +2*x^2+2*x-1). - Alois P. Heinz, Apr 08 2011

Extensions

Name and comments edited and terms a(12)-a(30) from Nathaniel Johnston, Apr 08 2011

A002526 Number of permutations of length n within distance 3 of a fixed permutation.

Original entry on oeis.org

1, 1, 2, 6, 24, 78, 230, 675, 2069, 6404, 19708, 60216, 183988, 563172, 1725349, 5284109, 16177694, 49526506, 151635752, 464286962, 1421566698, 4352505527, 13326304313, 40802053896, 124926806216, 382497958000, 1171122069784, 3585709284968, 10978628154457
Offset: 0

Views

Author

Keywords

Comments

For positive n, a(n) equals the permanent of the n X n matrix with 1's along the seven central diagonals, and 0's everywhere else. - John M. Campbell, Jul 09 2011

References

  • D. H. Lehmer, Permutations with strongly restricted displacements. Combinatorial theory and its applications, II (Proc. Colloq., Balatonfured, 1969), pp. 755-770. North-Holland, Amsterdam, 1970.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

The 14 sequences in Kløve's Table 3 are A002526, A002527, A002529, A188379, A188491, A188492, A188493, A188494, A002528, A188495, A188496, A188497, A188498, A002526.
Cf. A002524.
Column k=3 of A306209.

Programs

  • Magma
    R:=PowerSeriesRing(Integers(), 50); Coefficients(R!( (1-x-2*x^2-2*x^4+x^7+x^8)/(1-2*x-2*x^2-10*x^4-8*x^5+2*x^6+16*x^7+10*x^8 +2*x^9-4*x^10-2*x^11-2*x^13-x^14) )); // G. C. Greubel, Jan 22 2022
    
  • Mathematica
    CoefficientList[Series[(1-x-2x^2-2x^4+x^7+x^8)/(1-2x-2x^2-10x^4-8x^5+ 2x^6+ 16x^7+10x^8+2x^9-4x^10-2x^11-2x^13-x^14),{x,0,50}],x] (* or *) LinearRecurrence[{2,2,0,10,8,-2,-16,-10,-2,4,2,0,2,1},{1,1,2,6,24,78, 230, 675,2069,6404,19708,60216,183988,563172},51] (* Harvey P. Dale, Jun 22 2011 *)
  • PARI
    Vec((1-x-2*x^2-2*x^4+x^7+x^8)/(1-2*x-2*x^2-10*x^4-8*x^5+2*x^6+16*x^7+10*x^8+2*x^9-4*x^10-2*x^11-2*x^13-x^14)+O(x^99)) \\ Charles R Greathouse IV, Jul 16 2011
    
  • Sage
    [( (1-x-2*x^2-2*x^4+x^7+x^8)/(1-2*x-2*x^2-10*x^4-8*x^5+2*x^6+16*x^7+10*x^8 +2*x^9-4*x^10-2*x^11-2*x^13-x^14) ).series(x,n+1).list()[n] for n in (0..40)] # G. C. Greubel, Jan 22 2022

Formula

G.f.: (1-x-2*x^2-2*x^4+x^7+x^8)/(1-2*x-2*x^2-10*x^4-8*x^5+2*x^6+16*x^7+10*x^8 +2*x^9-4*x^10-2*x^11-2*x^13-x^14).
a(0)=1, a(1)=1, a(2)=2, a(3)=6, a(4)=24, a(5)=78, a(6)=230, a(7)=675, a(8)=2069, a(9)=6404, a(10)=19708, a(11)=60216, a(12)=183988, a(13)=563172, a(n) = 2*a(n-1) +2*a(n-2) +10*a(n-4) +8*a(n-5) -2*a(n-6) -16*a(n-7) -10*a(n-8) -2*a(n-9) +4*a(n-10) +2*a(n-11) +2*a(n-13) +a(n-14). - Harvey P. Dale, Jun 22 2011

A188379 a(n) = A002526(n+1) - A002527(n+1).

Original entry on oeis.org

0, 0, 0, 6, 18, 46, 115, 374, 1204, 3752, 11300, 34324, 105124, 322989, 989692, 3028484, 9267328, 28374898, 86891022, 266058106, 814585879, 2494006074, 7636057864, 23380074400, 71584762200, 219176102664, 671066472872, 2054652945289
Offset: 0

Views

Author

N. J. A. Sloane, Apr 01 2011

Keywords

Comments

For n >= 3, a(n) is the number of permutations p on the set [n] with the properties that abs(p(i)-i) <= 3 for all i and p(j) <= 2+j for j = 1,2,3.
For n >= 3, a(n) is also the permanent of the n X n matrix that has ones on its diagonal, ones on its three superdiagonals, ones on its three subdiagonals (with the exception of zeros in the (4,1),(5,2), and (6,3)-entries), and is zero elsewhere.
This is row 4 of Kløve's Table 3.

Programs

  • Maple
    with (LinearAlgebra):
    A188379:= n-> `if` (n<=2, 0, Permanent (Matrix (n, (i, j)->
                  `if` (abs(j-i)<4 and [i, j]<>[4, 1] and [i, j]<>[5, 2] and [i, j]<>[6, 3], 1, 0)))):
    seq (A188379(n), n=0..20);
  • Mathematica
    a[n_] := Permanent[Table[If[Abs[j - i] < 4 && {i, j} != {4, 1} && {i, j} != {5, 2} && {i, j} != {6, 3}, 1, 0], {i, 1, n}, {j, 1, n}]]; a[1] = 0; Table[an = a[n]; Print["a(", n, ") = ", an]; an, {n, 0, 20}] (* Jean-François Alcover, Jan 07 2016, adapted from Maple *)
    CoefficientList[Series[-(x^10 + 2 x^9 + 2 x^7 + 4 x^6 - 2 x^5 - 8 x^4 - 13 x^3 - 2 x^2 + 6 x+6) x^3 / (x^14 + 2 x^13 + 2 x^11 + 4 x^10 - 2 x^9 - 10 x^8 - 16 x^7 - 2 x^6 + 8 x^5 + 10 x^4 + 2 x^2 + 2 x - 1), {x, 0, 33}], x] (* Vincenzo Librandi, Jan 07 2016 *)

Formula

a(n) = A002529(n-1) + A188492(n-1) + A188493(n-1). - Nathaniel Johnston, Apr 08 2011
G.f.: -(x^10 +2*x^9 +2*x^7 +4*x^6 -2*x^5 -8*x^4 -13*x^3 -2*x^2 +6*x+6) * x^3 / (x^14 +2*x^13 +2*x^11 +4*x^10 -2*x^9 -10*x^8 -16*x^7 -2*x^6 +8*x^5 +10*x^4 +2*x^2 +2*x-1). - Alois P. Heinz, Apr 07 2011

Extensions

Name and comments edited by Nathaniel Johnston, Apr 08 2011
Showing 1-4 of 4 results.