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

A126473 Number of strings over a 5 symbol alphabet with adjacent symbols differing by three or less.

Original entry on oeis.org

1, 5, 23, 107, 497, 2309, 10727, 49835, 231521, 1075589, 4996919, 23214443, 107848529, 501037445, 2327695367, 10813893803, 50238661313, 233396326661, 1084301290583, 5037394142315, 23402480441009, 108722104190981, 505095858086951, 2346549744920747
Offset: 0

Views

Author

R. H. Hardin, Dec 27 2006

Keywords

Comments

[Empirical] a(base,n) = a(base-1,n) + 7^(n-1) for base >= 3n-2; a(base,n) = a(base-1,n) + 7^(n-1)-2 when base = 3n-3.
From Johannes W. Meijer, Aug 01 2010: (Start)
The a(n) represent the number of n-move routes of a fairy chess piece starting in a given side square (m = 2, 4, 6 or 8) on a 3 X 3 chessboard. This fairy chess piece behaves like a king on the eight side and corner squares but on the central square the king goes crazy and turns into a red king, see A179596.
For the side squares the 512 red kings lead to 47 different red king sequences, see the cross-references for some examples.
The sequence above corresponds to four A[5] vectors with the decimal [binary] values 367 [1,0,1,1,0,1,1,1,1], 463 [1,1,1,0,0,1,1,1,1], 487 [1,1,1,1,0,0,1,1,1] and 493 [1,1,1,1,0,1,1,0,1]. These vectors lead for the corner squares to A179596 and for the central square to A179597.
This sequence belongs to a family of sequences with g.f. (1+x)/(1-4*x-k*x^2). Red king sequences that are members of this family are A003947 (k=0), A015448 (k=1), A123347 (k=2), A126473 (k=3; this sequence) and A086347 (k=4). Other members of this family are A000351 (k=5), A001834 (k=-1), A111567 (k=-2), A048473 (k=-3) and A053220 (k=-4)
Inverse binomial transform of A154244. (End)
Equals the INVERT transform of A055099: (1, 4, 14, 50, 178, ...). - Gary W. Adamson, Aug 14 2010
Number of one-sided n-step walks taking steps from {E, W, N, NE, NW}. - Shanzhen Gao, May 10 2011
For n>=1, a(n) equals the numbers of words of length n-1 on alphabet {0,1,2,3,4} containing no subwords 00 and 11. - Milan Janjic, Jan 31 2015

Crossrefs

Cf. 5 symbol differing by two or less A126392, one or less A057960.
Cf. Red king sequences side squares [numerical value A[5]]: A086347 [495], A179598 [239], A126473 [367], A123347 [335], A179602 [95], A154964 [31], A015448 [327], A152187 [27], A003947 [325], A108981 [11], A007483 [2]. - Johannes W. Meijer, Aug 01 2010
Cf. A055099.

Programs

  • Maple
    with(LinearAlgebra): nmax:=19; m:=2; A[5]:= [1,0,1,1,0,1,1,1,1]: A:=Matrix([[0,1,0,1,1,0,0,0,0],[1,0,1,1,1,1,0,0,0],[0,1,0,0,1,1,0,0,0],[1,1,0,0,1,0,1,1,0],A[5],[0,1,1,0,1,0,0,1,1],[0,0,0,1,1,0,0,1,0],[0,0,0,1,1,1,1,0,1],[0,0,0,0,1,1,0,1,0]]): for n from 0 to nmax do B(n):=A^n: a(n):= add(B(n)[m,k],k=1..9): od: seq(a(n), n=0..nmax); # Johannes W. Meijer, Aug 01 2010
    # second Maple program:
    a:= n-> (M-> M[1,2]+M[2,2])(<<0|1>, <3|4>>^n):
    seq(a(n), n=0..24);  # Alois P. Heinz, Jun 28 2021
  • Mathematica
    LinearRecurrence[{4, 3}, {1, 5}, 24] (* Jean-François Alcover, Dec 10 2024 *)
  • PARI
    a(n)=([0,1; 3,4]^n*[1;5])[1,1] \\ Charles R Greathouse IV, May 10 2016

Formula

From Johannes W. Meijer, Aug 01 2010: (Start)
G.f.: (1+x)/(1-4*x-3*x^2).
a(n) = 4*a(n-1) + 3*a(n-2) with a(0) = 1 and a(1) = 5.
a(n) = ((1+3/sqrt(7))/2)*(A)^(-n) + ((1-3/sqrt(7))/2)*(B)^(-n) with A = (-2 + sqrt(7))/3 and B = (-2-sqrt(7))/3.
Lim_{k->oo} a(n+k)/a(k) = (-1)^(n+1)*A000244(n)/(A015530(n)*sqrt(7)-A108851(n))
(End)
a(n) = A015330(n)+A015330(n+1). - R. J. Mathar, May 09 2023

Extensions

Edited by Johannes W. Meijer, Aug 10 2010

A214992 Power ceiling-floor sequence of (golden ratio)^4.

Original entry on oeis.org

7, 47, 323, 2213, 15169, 103969, 712615, 4884335, 33477731, 229459781, 1572740737, 10779725377, 73885336903, 506417632943, 3471038093699, 23790849022949, 163064905066945, 1117663486445665, 7660579500052711
Offset: 0

Views

Author

Clark Kimberling, Nov 08 2012, Jan 24 2013

Keywords

Comments

Let f = floor and c = ceiling. For x > 1, define four sequences as functions of x, as follows:
p1(0) = f(x), p1(n) = f(x*p1(n-1));
p2(0) = f(x), p2(n) = c(x*p2(n-1)) if n is odd and p2(n) = f(x*p1(n-1)) if n is even;
p3(0) = c(x), p3(n) = f(x*p3(n-1)) if n is odd and p3(n) = c(x*p3(n-1)) if n is even;
p4(0) = c(x), p4(n) = c(x*p4(n-1)).
The present sequence is given by a(n) = p3(n).
Following the terminology at A214986, call the four sequences power floor, power floor-ceiling, power ceiling-floor, and power ceiling sequences. In the table below, a sequence is identified with an A-numbered sequence if they appear to agree except possibly for initial terms. Notation: S(t)=sqrt(t), r = (1+S(5))/2 = golden ratio, and Limit = limit of p3(n)/p2(n).
x ......p1..... p2..... p3..... p4.......Limit
r^2.....A001519 A001654 A061646 A001906..-1+S(5)
r^3.....A024551 A001076 A015448 A049652..-1+S(5)
r^4.....A049685 A157335 A214992 A004187..-19+9*S(5)
r^5.....A214993 A049666 A015457 A214994...(-9+5*S(5))/2
r^6.....A007805 A156085 A214995 A049660..-151+68*S(5)
2+S(2)..A007052 A214996 A214997 A007070..(1+S(2))/2
1+S(3)..A057960 A002605 A028859 A077846..(1+S(3))/2
2+S(3)..A001835 A109437 A214998 A001353..-4+3*S(3)
S(5)....A214999 A215091 A218982 A218983..1.26879683...
2+S(5)..A024551 A001076 A015448 A049652..-1+S(5)
2+S(6)..A218984 A090017 A123347 A218985..S(3/2)
2+S(7)..A218986 A015530 A126473 A218987..(1+S(7))/3
2+S(8)..A218988 A057087 A086347 A218989..(1+S(2))/2
3+S(8)..A001653 A084158 A218990 A001109..-13+10*S(2)
3+S(10).A218991 A005668 A015451 A218992..-2+S(10)
...
Properties of p1, p2, p3, p4:
(1) If x > 2, the terms of p2 and p3 interlace: p2(0) < p3(0) < p2(1) < p3(1) < p2(2) < p3(2)... Also, p1(n) <= p2(n) <= p3(n) <= p4(n) <= p1(n+1) for all x>0 and n>=0.
(2) If x > 2, the limits L(x) = limit(p/x^n) exist for the four functions p(x), and L1(x) <= L2(x) <= L3(x) <= L4 (x). See the Mathematica programs for plots of the four functions; one of them also occurs in the Odlyzko and Wilf article, along with a discussion of the special case x = 3/2.
(3) Suppose that x = u + sqrt(v) where v is a nonsquare positive integer. If u = f(x) or u = c(x), then p1, p2, p3, p4 are linear recurrence sequences. Is this true for sequences p1, p2, p3, p4 obtained from x = (u + sqrt(v))^q for every positive integer q?
(4) Suppose that x is a Pisot-Vijayaraghavan number. Must p1, p2, p3, p4 then be linearly recurrent? If x is also a quadratic irrational b + c*sqrt(d), must the four limits L(x) be in the field Q(sqrt(d))?
(5) The Odlyzko and Wilf article (page 239) raises three interesting questions about the power ceiling function; it appears that they remain open.

Examples

			a(0) = ceiling(r) = 7, where r = ((1+sqrt(5))/2)^4 = 6.8...; a(1) = floor(7*r) = 47; a(2) = ceiling(47) = 323.
		

Crossrefs

Programs

  • Mathematica
    (* Program 1.  A214992 and related sequences *)
    x = GoldenRatio^4; z = 30; (* z = # terms in sequences *)
    z1 = 100; (* z1 = # digits in approximations *)
    f[x_] := Floor[x]; c[x_] := Ceiling[x];
    p1[0] = f[x]; p2[0] = f[x]; p3[0] = c[x]; p4[0] = c[x];
    p1[n_] := f[x*p1[n - 1]]
    p2[n_] := If[Mod[n, 2] == 1, c[x*p2[n - 1]], f[x*p2[n - 1]]]
    p3[n_] := If[Mod[n, 2] == 1, f[x*p3[n - 1]], c[x*p3[n - 1]]]
    p4[n_] := c[x*p4[n - 1]]
    Table[p1[n], {n, 0, z}]  (* A049685 *)
    Table[p2[n], {n, 0, z}]  (* A157335 *)
    Table[p3[n], {n, 0, z}]  (* A214992 *)
    Table[p4[n], {n, 0, z}]  (* A004187 *)
    Table[p4[n] - p1[n], {n, 0, z}]  (* A004187 *)
    Table[p3[n] - p2[n], {n, 0, z}]  (* A098305 *)
    (* Program 2.  Plot of power floor and power ceiling functions, p1(x) and p4(x) *)
    f[x_] := f[x] = Floor[x]; c[x_] := c[x] = Ceiling[x];
    p1[x_, 0] := f[x]; p1[x_, n_] := f[x*p1[x, n - 1]];
    p4[x_, 0] := c[x]; p4[x_, n_] := c[x*p4[x, n - 1]];
    Plot[Evaluate[{p1[x, 10]/x^10, p4[x, 10]/x^10}], {x, 2, 3}, PlotRange -> {0, 4}]
    (* Program 3. Plot of power floor-ceiling and power ceiling-floor functions, p2(x) and p3(x) *)
    f[x_] := f[x] = Floor[x]; c[x_] := c[x] = Ceiling[x];
    p2[x_, 0] := f[x]; p3[x_, 0] := c[x];
    p2[x_, n_] := If[Mod[n, 2] == 1, c[x*p2[x, n - 1]], f[x*p2[x, n - 1]]]
    p3[x_, n_] := If[Mod[n, 2] == 1, f[x*p3[x, n - 1]], c[x*p3[x, n - 1]]]
    Plot[Evaluate[{p2[x, 10]/x^10, p3[x, 10]/x^10}], {x, 2, 3}, PlotRange -> {0, 4}]

Formula

a(n) = floor(r*a(n-1)) if n is odd and a(n) = ceiling(r*a(n-1)) if n is even, where a(0) = ceiling(r), r = (golden ratio)^4 = (7 + sqrt(5))/2.
a(n) = 6*a(n-1) + 6*a(n-2) - a(n-3).
G.f.: (7 + 5*x - x^2)/((1 + x)*(1 - 7*x + x^2)).
a(n) = (10*(-2)^n+(10+3*sqrt(5))*(7-3*sqrt(5))^(n+2)+(10-3*sqrt(5))*(7+3*sqrt(5))^(n+2))/(90*2^n). - Bruno Berselli, Nov 14 2012
a(n) = 7*A157335(n) + 5*A157335(n-1) - A157335(n-2). - R. J. Mathar, Feb 05 2020
E.g.f.: exp(-x)*(5 + 2*exp(9*x/2)*(155*cosh(3*sqrt(5)*x/2) + 69*sqrt(5)*sinh(3*sqrt(5)*x/2)))/45. - Stefano Spezia, Oct 28 2024

A180140 Eight rooks and one berserker on a 3 X 3 chessboard. G.f.: (1+x+x^2)/(1-3*x-5*x^2).

Original entry on oeis.org

1, 4, 18, 74, 312, 1306, 5478, 22964, 96282, 403666, 1692408, 7095554, 29748702, 124723876, 522915138, 2192364794, 9191670072, 38536834186, 161568852918, 677390729684, 2840016453642, 11907003009346, 49921091296248
Offset: 0

Views

Author

Johannes W. Meijer, Aug 13 2010, Jun 15 2013

Keywords

Comments

a(n) gives the number of n-move routes of a fairy chess piece starting in a given side square (m = 2, 4, 6, 8) on a 3 X 3 chessboard. This fairy chess piece behaves like a rook on the four side and four corner (m = 1, 3, 7, 9) squares but on the center square (m = 5) it goes berserk and turns into a berserker. For this sequence, the berserker can move to three of the side squares and three of the corners from the center.
The berserker is one of the Lewis chessmen which were discovered in 1831 on the Isle of Lewis. They are carved from walrus ivory in Scandinavian style of the 12th century. The pawns look like decorated tombstones. The pieces have all human representations with facial expressions varying from gloom to anger. Some of the rooks show men biting their shield in the manner of berserkers. According to Hooper and Whyld none looks happy.
Let A be the adjacency matrix of the graph G, where V(G) = {v1, v2, v3, v4, v5, v6, v7, v8, v9}. Then the (m, k) entry of A^n is the number of different vm-vk walks of length n in G, see the Chartrand reference. In the adjacency matrix A, see the Maple program, the A[1], A[3], A[7] and A[9] vectors represent the rook moves on the corner squares, the A[2], A[4], A[6] and A[8] vectors represent the rook moves on the side squares and the A[5] vector represents the moves of the berserker. On a 3 X 3 chessboard there are 2^9 = 512 ways a berserker could move from the center square (off the center the berserker behaves like a rook) so there are 512 different berserkers.
For the side squares the 512 berserker vectors lead to 42 different sequences, see the overview of berserker sequences. There are 16 berserker vectors that lead to the sequence given above. Their decimal [binary] values are: 111 [001 101 111] , 207 [011 001 111], 231 [011 100 111], 237 [011 101 101], 303 [100 101 111], 363 [101 101 011], 366 [101 101 110], 399 [110 001 111], 423 [110 100 111], 429 [110 101 101], 459 [111 001 011], 462 [111 001 110], 483 [111 100 011], 486 [111 100 110], 489 [111 101 001] and 492 [111 101 100]. These berserker vectors lead for the corner squares to sequence 4*A179606 (with leading term 1 added) and for the central square to sequence 6*A179606 (with leading term 1 added).
This sequence belongs to a family of sequences with GF(x)=(1+x-k*x^2)/(1-3*x+(k-4)*x^2), see A180142.

References

  • Gary Chartrand, Introductory Graph Theory, pp. 217-221, 1984.
  • David Hooper and Kenneth Whyld, The Oxford Companion to Chess, pp. 131, 225, 1992.

Crossrefs

Cf. A180141 (corner squares) and A180147 (central square).
Cf. Berserker sequences side squares: 4*A007482 (with leading 1 added), A180144, A003500 (n>=1 and a(0)=1), A180142, A000302, A180140 (this sequence), 2*A001077 (n>=1 and a(0)=1), A180146, 4*A154964 (n>=1 and a(0)=1), 4*A123347 (with leading 1 added).

Programs

  • Maple
    nmax:=22; m:=2; A[1]:=[0, 1, 1, 1, 0, 0, 1, 0, 0]: A[2]:=[1, 0, 1, 0, 1, 0, 0, 1, 0]: A[3]:= [1, 1, 0, 0, 0, 1, 0, 0, 1]: A[4]:= [1, 0, 0, 0, 1, 1, 1, 0, 0]: A[5]:=[0, 0, 1, 1, 0, 1, 1, 1, 1]: A[6]:=[0, 0, 1, 1, 1, 0, 0, 0, 1]: A[7]:=[1, 0, 0, 1, 0, 0, 0, 1, 1]: A[8]:=[0, 1, 0, 0, 1, 0, 1, 0, 1]: A[9]:=[0, 0, 1, 0, 0, 1, 1, 1, 0]: A:= Matrix([A[1], A[2], A[3], A[4], A[5], A[6], A[7], A[8], A[9]]): for n from 0 to nmax do B(n):=A^n: a(n):= add(B(n)[m, k], k=1..9): od: seq(a(n), n=0..nmax);
  • Mathematica
    CoefficientList[Series[(1+x+x^2)/(1-3*x-5*x^2), {x, 0, 22}],x] (* or *) LinearRecurrence[{3,5,0},{1,4,18},23] (* Indranil Ghosh, Mar 05 2017 *)
  • PARI
    print(Vec((1 + x + x^2)/(1- 3*x - 5*x^2) + O(x^23))); \\ Indranil Ghosh, Mar 05 2017

Formula

G.f.: (1+x+x^2)/(1-3*x-5*x^2).
a(n) = 3*a(n-1) + 5*a(n-2) for n>=3 with a(0)=1, a(1)=4 and a(2)=18.
a(n) = ((22+54*A)*A^(-n-1) + (22+54*B)*B^(-n-1))/145 with A=(-3+sqrt(29))/10 and B=(-3-sqrt(29))/10 for n>=1 with a(0)=1.
5*a(n) = 2*( A015523(n) + 3*A015523(n+1)), n>0 - R. J. Mathar, May 11 2013

A116412 Riordan array ((1+x)/(1-2x),x(1+x)/(1-2x)).

Original entry on oeis.org

1, 3, 1, 6, 6, 1, 12, 21, 9, 1, 24, 60, 45, 12, 1, 48, 156, 171, 78, 15, 1, 96, 384, 558, 372, 120, 18, 1, 192, 912, 1656, 1473, 690, 171, 21, 1, 384, 2112, 4608, 5160, 3225, 1152, 231, 24, 1, 768, 4800, 12240, 16584, 13083, 6219, 1785, 300, 27, 1, 1536, 10752
Offset: 0

Views

Author

Paul Barry, Feb 13 2006

Keywords

Comments

Row sums are A003688. Diagonal sums are A116413. Product of A007318 and A116413 is A116414. Product of A007318 and A105475.
Subtriangle of triangle given by (0, 3, -1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Jan 18 2012

Examples

			Triangle begins
1,
3, 1,
6, 6, 1,
12, 21, 9, 1,
24, 60, 45, 12, 1,
48, 156, 171, 78, 15, 1
Triangle T(n,k), 0<=k<=n, given by (0, 3, -1, 0, 0, 0, ...) DELTA (1, 0, 0, 0, 0, ...) begins :
1
0, 1
0, 3, 1
0, 6, 6, 1
0, 12, 21, 9, 1
0, 24, 60, 45, 12, 1
0, 48, 156, 171, 78, 15, 1
... - _Philippe Deléham_, Jan 18 2012
		

Crossrefs

Programs

  • Mathematica
    With[{n = 10}, DeleteCases[#, 0] & /@ CoefficientList[Series[(1 + x)/(1 - (y + 2) x - y x^2), {x, 0, n}, {y, 0, n}], {x, y}]] // Flatten (* Michael De Vlieger, Apr 25 2018 *)

Formula

Number triangle T(n,k)=sum{j=0..n, C(k+1,j)*C(n-j,k)2^(n-k-j)}
From Vladimir Kruchinin, Mar 17 2011: (Start)
T((m+1)*n+r-1, m*n+r-1) * r/(m*n+r) = sum(k=1..n, k/n * T((m+1)*n-k-1, m*n-1) * T(r+k-1,r-1)), n>=m>1.
T(n-1,m-1) = m/n * sum(k=1..n-m+1, k*A003945(k-1)*T(n-k-1,m-2)), n>=m>1. (End)
G.f.: (1+x)/(1-(y+2)*x -y*x^2). - Philippe Deléham, Jan 18 2012
Sum_{k, 0<=k<=n} T(n,k)*x^k = A104537(n), A110523(n), (-2)^floor(n/2), A057079(n), A003945(n), A003688(n+1), A123347(n), A180035(n) for x = -4, -3, -2, -1, 0, 1, 2, 3 respectively. - Philippe Deléham, Jan 18 2012
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) + T(n-2,k-1), T(0,0) = 1, T(1,0) = 3, T(1,1) = 1, T(2,0) = T(2,1) = 6, T(2,2) = 1, T(n,k) = 0 if k>n or if k<0. - Philippe Deléham, Oct 31 2013

A218984 Power floor sequence of 2+sqrt(6).

Original entry on oeis.org

4, 17, 75, 333, 1481, 6589, 29317, 130445, 580413, 2582541, 11490989, 51129037, 227498125, 1012250573, 4503998541, 20040495309, 89169978317, 396760903885, 1765383572173, 7855056096461, 34950991530189, 155514078313677, 691958296315085, 3078861341887693
Offset: 0

Views

Author

Clark Kimberling, Nov 11 2012

Keywords

Comments

See A214992 for a discussion of power floor sequence and the power floor function, p1(x) = limit of a(n,x)/x^n. The present sequence is a(n,r), where r = 2+sqrt(6), and the limit p1(r) = 3.77794213613376987528458445727451673384055973517...

Examples

			a(0) = [r] = 4, where r = 2+sqrt(6); a(1) = [4*r] = 17; a(2) = [17*r] = 75.
		

Crossrefs

Programs

  • Mathematica
    x = 2 + Sqrt[6]; z = 30; (* z = # terms in sequences *)
    f[x_] := Floor[x]; c[x_] := Ceiling[x];
    p1[0] = f[x]; p2[0] = f[x]; p3[0] = c[x]; p4[0] = c[x];
    p1[n_] := f[x*p1[n - 1]]
    p2[n_] := If[Mod[n, 2] == 1, c[x*p2[n - 1]], f[x*p2[n - 1]]]
    p3[n_] := If[Mod[n, 2] == 1, f[x*p3[n - 1]], c[x*p3[n - 1]]]
    p4[n_] := c[x*p4[n - 1]]
    t1 = Table[p1[n], {n, 0, z}]  (* A218984 *)
    t2 = Table[p2[n], {n, 0, z}]  (* A090017 *)
    t3 = Table[p3[n], {n, 0, z}]  (* A123347 *)
    t4 = Table[p4[n], {n, 0, z}]  (* A218985 *)
  • PARI
    Vec((4 - 3*x - 2*x^2) / ((1 - x)*(1 - 4*x - 2*x^2)) + O(x^40)) \\ Colin Barker, Nov 13 2017

Formula

a(n) = [x*a(n-1)], where x=2+sqrt(6), a(0) = [x].
a(n) = 5*a(n-1) - 2*a(n-2) - 2*a(n-3).
G.f.: (4 - 3*x - 2*x^2)/(1 - 5*x + 2*x^2 + 2*x^3).
a(n) = (1/30)*(6 + (57-23*sqrt(6))*(2-sqrt(6))^n + (2+sqrt(6))^n*(57+23*sqrt(6))). - Colin Barker, Nov 13 2017

A218985 Power ceiling sequence of 2+sqrt(6).

Original entry on oeis.org

5, 23, 103, 459, 2043, 9091, 40451, 179987, 800851, 3563379, 15855219, 70547635, 313900979, 1396699187, 6214598707, 27651793203, 123036370227, 547449067315, 2435869009715, 10838374173491, 48225234713395, 214577687200563, 954761218229043, 4248200247317299
Offset: 0

Views

Author

Clark Kimberling, Nov 11 2012

Keywords

Comments

See A214992 for a discussion of power ceiling sequence and the power ceiling function, p4(x) = limit of a(n,x)/x^n. The present sequence is a(n,r), where r = 2+sqrt(6), and the limit p4(r) = 5.2127890589687233047437696796862841514303439...
See A218984 for the power floor function, p1(x). For comparison of p4 and p1, limit(p4(r)/p1(r)) = 2*(1+sqrt(6))/5 = 1.379795897113271239278913629882356556786378...

Examples

			a(0) = ceiling(r) = 5, where r = 2+sqrt(6).
a(1) = ceiling(5*r) = 23; a(2) = ceiling(23*r) = 103.
		

Crossrefs

Programs

  • Mathematica
    (See A218984.)
  • PARI
    Vec((5 - 2*x - 2*x^2) / ((1 - x)*(1 - 4*x - 2*x^2)) + O(x^40)) \\ Colin Barker, Nov 13 2017

Formula

a(n) = [x*a(n-1)], where x=2+sqrt(6), a(0) = [x].
a(n) = 5*a(n-1) - 2*a(n-2) - 2*a(n-3).
G.f.: (5 - 2*x - 2*x^2)/(1 - 5*x + 2*x^2 + 2*x^3).
a(n) = (1/15)*(-3 + (39-16*sqrt(6))*(2-sqrt(6))^n + (2+sqrt(6))^n*(39+16*sqrt(6))). - Colin Barker, Nov 13 2017

A179600 Eight white kings and one red king on a 3 X 3 chessboard. G.f.: (1 + x)/(1 - 2*x - 10*x^2 - 4*x^3).

Original entry on oeis.org

1, 3, 16, 66, 304, 1332, 5968, 26472, 117952, 524496, 2334400, 10385568, 46213120, 205619520, 914912512, 4070872704, 18113348608, 80595074304, 358607125504, 1595618388480, 7099688329216, 31589989045248, 140559334936576
Offset: 0

Views

Author

Johannes W. Meijer, Jul 28 2010

Keywords

Comments

The a(n) represent the number of n-move routes of a fairy chess piece starting in a given corner square (m = 1, 3, 7 or 9) on a 3 X 3 chessboard. This fairy chess piece behaves like a king on the eight side and corner squares but on the central square the king goes crazy and turns into a red king, see A179596.
The sequence above corresponds to 6 red king vectors, i.e., A[5] vectors, with decimal values 335, 359, 365, 455, 461 and 485. These vectors lead for the side squares to A123347 and for the central square to A179601.

Programs

  • Maple
    with(LinearAlgebra): nmax:=24; m:=1; A[1]:= [0,1,0,1,1,0,0,0,0]: A[2]:= [1,0,1,1,1,1,0,0,0]: A[3]:= [0,1,0,0,1,1,0,0,0]: A[4]:=[1,1,0,0,1,0,1,1,0]: A[5]:= [1,1,1,0,0,0,1,1,1]: A[6]:= [0,1,1,0,1,0,0,1,1]: A[7]:= [0,0,0,1,1,0,0,1,0]: A[8]:= [0,0,0,1,1,1,1,0,1]: A[9]:= [0,0,0,0,1,1,0,1,0]: A:=Matrix([A[1],A[2],A[3],A[4],A[5],A[6],A[7],A[8],A[9]]): for n from 0 to nmax do B(n):=A^n: a(n):= add(B(n)[m,k],k=1..9): od: seq(a(n), n=0..nmax);
  • PARI
    Vec((1+x)/(1 - 2*x - 10*x^2 - 4*x^3) + O(x^40)) \\ Jinyuan Wang, Mar 10 2020

Formula

G.f.: (1+x)/(1 - 2*x - 10*x^2 - 4*x^3).
a(n) = 2*a(n-1) + 10*a(n-2) + 4*a(n-3) with a(0)=1, a(1)=3 and a(2)=16.
a(n) = (4*(-1/2)^(-n) + (1+sqrt(6))*A^(-n-1) + (1-sqrt(6))*B^(-n-1))/20 with A = (-1+sqrt(6)/2) and B = (-1-sqrt(6)/2).
Lim_{k->infinity} a(n+k)/a(k) = (-1)^(n+1)*(A016116(n+1)/(A041007(n-1)*sqrt(6) - A041006(n-1))) for n => 1.

A179601 Eight white kings and one red king on a 3 X 3 chessboard. G.f.: (1+4*x)/(1 - 2*x - 10*x^2 - 4*x^3).

Original entry on oeis.org

1, 6, 22, 108, 460, 2088, 9208, 41136, 182704, 813600, 3618784, 16104384, 71651008, 318820992, 1418569600, 6311953152, 28084886272, 124963582464, 556023840256, 2474023050240, 11008138832896, 48980603529216, 217938687588352
Offset: 0

Views

Author

Johannes W. Meijer, Jul 28 2010

Keywords

Comments

The a(n) represent the number of n-move routes of a fairy chess piece starting in the central square (m = 5) on a 3 X 3 chessboard. This fairy chess piece behaves like a king on the eight side and corner squares but on the central square the king goes crazy and turns into a red king, see A179596.
The sequence above corresponds to 6 red king vectors, i.e., A[5] vectors, with decimal values 335, 359, 365, 455, 461 and 485. These vectors lead for the corner squares to A179600 and for the side squares to A123347.

Crossrefs

Cf. A041006, A041007, A123347, A179596, A179597 (central square), A179600.

Programs

  • Maple
    with(LinearAlgebra): nmax:=22; m:=5; A[1]:= [0,1,0,1,1,0,0,0,0]: A[2]:= [1,0,1,1,1,1,0,0,0]: A[3]:= [0,1,0,0,1,1,0,0,0]: A[4]:= [1,1,0,0,1,0,1,1,0]: A[5]:= [1,1,1,0,0,0,1,1,1]: A[6]:= [0,1,1,0,1,0,0,1,1]: A[7]:= [0,0,0,1,1,0,0,1,0]: A[8]:= [0,0,0,1,1,1,1,0,1]: A[9]:= [0,0,0,0,1,1,0,1,0]: A:=Matrix([A[1],A[2],A[3],A[4],A[5],A[6],A[7],A[8],A[9]]): for n from 0 to nmax do B(n):=A^n: a(n):= add(B(n)[m,k],k=1..9): od: seq(a(n), n=0..nmax);

Formula

G.f.: ( -1-4*x ) / ( (2*x+1)*(2*x^2 + 4*x - 1) ).
a(n) = 2*a(n-1) + 10*a(n-2) + 4*a(n-3) with a(0)=1, a(1)=6 and a(2)=22.
a(n) = (-2/5)*(-1/2)^(-n) + ((2+3*A)*A^(-n-1) + (2+3*B)*B^(-n-1))/10 with A = (-1+sqrt(6)/2) and B = (-1-sqrt(6)/2).
Limit_{k->oo} a(n+k)/a(k) = (-1)^(n+1)*A016116(n+1)/(A041007(n-1)*sqrt(6) - A041006(n-1)) for n => 1.
Showing 1-8 of 8 results.