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

A188493 a(n) = A188491(n-1) + A188495(n-1) + A188497(n-1).

Original entry on oeis.org

0, 0, 2, 6, 14, 31, 104, 344, 1084, 3236, 9784, 29964, 92241, 282780, 865064, 2646292, 8102454, 24813838, 75982346, 232630527, 712230076, 2180675264, 6676819512, 20443032008, 62591840320, 191641545768, 586762729889, 1796535598952, 5500587026592
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(j) <= 2+j for j = 1,2, and p(j) >= j-2 for j = 4,5.
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 zeros in the (1,4) and (2,5)-entries), 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 7 of Kløve's Table 3.

Programs

  • Maple
    with (LinearAlgebra):
    A188493:= 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] and [i, j]<>[1, 4] and [i, j]<>[2, 5], 1, 0)))):
    seq (A188493(n), n=0..20);
  • Mathematica
    a[n_] := Permanent[Table[If[Abs[j - i] < 4 && {i, j} != {4, 1} && {i, j} != {5, 2} && {i, j} != {1, 4} && {i, j} != {2, 5}, 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 *)

Formula

G.f.: -(x^10+2*x^9+2*x^7 +4*x^6-2*x^5-6*x^4 -9*x^3-2*x^2+2*x+2) *x^2 / (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 a(12)-a(28) from Nathaniel Johnston, Apr 08 2011

A188497 a(n) = A188493(n+1) - A188491(n) - A188495(n).

Original entry on oeis.org

0, 0, 2, 4, 7, 20, 72, 240, 722, 2140, 6508, 20077, 61776, 189056, 577856, 1768380, 5416230, 16587984, 50788707, 155489884, 476058864, 1457605616, 4462928950, 13664497400, 41837412392, 128096408137, 392202398144, 1200835918016, 3676688064688
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(j) <= 1+j for j=1,2, 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), (4,2), and (5,2)-entries), and is zero elsewhere.
This is row 12 of Kløve's Table 3.

Programs

  • Maple
    with(LinearAlgebra):
    A188497:= 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]<>[4, 2] and [i, j]<>[5, 2], 1, 0)))):
    seq (A188497(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} != {4, 2} && {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 07 2016, adapted from Maple *)
    CoefficientList[Series[x^2 (x^9 + 2 x^8 - 2 x^4 - 2 x^3 - 5 x^2 + 2) / ((1 - x) (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)), {x, 0, 33}], x] (* Vincenzo Librandi, Jan 07 2016 *)
  • PARI
    concat([0,0], Vec(-x^2*(x^9 +2*x^8 -2*x^4 -2*x^3 -5*x^2 +2) / ((x -1)*(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) = A188494(n-1) + A188498(n-1). - Nathaniel Johnston, Apr 11 2011
G.f.: -x^2*(x^9 +2*x^8 -2*x^4 -2*x^3 -5*x^2 +2) / ((x -1)*(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

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

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

A188492 a(n) = A002526(n+2) + A002526(n) - A002527(n+2) - A002527(n+1) + A002527(n) - A188493(n).

Original entry on oeis.org

0, 0, 2, 6, 14, 38, 124, 400, 1232, 3712, 11288, 34628, 106352, 325772, 996712, 3050352, 9340170, 28602014, 87576426, 268129662, 820931640, 2513509536, 7695861408, 23563048304, 72144604576, 220890113784, 676315440208, 2070725515096
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(j) <= 2+j for j = 1,2, 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 (4,1) and (5,2)-entries), and is zero elsewhere.
This is row 6 of Kløve's Table 3.

Programs

  • Maple
    with (LinearAlgebra):
    A188492:= 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] and [i, j]<>[1, 4], 1, 0)))):
    seq (A188492(n), n=0..20);
  • Mathematica
    a[n_] := Permanent[Table[If[Abs[j-i] < 4 && {i, j} != {4, 1} && {i, j} != {5, 2} && {i, j} != {1, 4}, 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 *)

Formula

a(n) = A002527(n-1) + A188495(n-1) + A188496(n-1). - Nathaniel Johnston, Apr 08 2011
G.f.: 2*x^2 * (x^2+2*x+1) / (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). - Alois P. Heinz, Apr 09 2011

Extensions

Name and comments edited, and a(12)-a(27) from Nathaniel Johnston, Apr 08 2011
Showing 1-5 of 5 results.