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.

A085810 Number of three-choice paths along a corridor of height 5, starting from the lower side.

Original entry on oeis.org

1, 2, 5, 13, 35, 96, 266, 741, 2070, 5791, 16213, 45409, 127206, 356384, 998509, 2797678, 7838801, 21963661, 61540563, 172432468, 483144522, 1353740121, 3793094450, 10628012915, 29779028189, 83438979561, 233790820762, 655067316176, 1835457822857, 5142838522138, 14409913303805
Offset: 1

Views

Author

Philippe Deléham, Jul 25 2003

Keywords

Comments

From Svjetlan Feretic, Jun 01 2013: (Start)
A three-choice path is a path whose steps lie in the set {(1,1), (1,0), (1,-1)}.
The paths under consideration "live" in a corridor like 0<=y<=5. Thus, the ordinate of a vertex of a path can take six values (0,1,2,3,4,5), but the height of the corridor is five.
a(1)=1 is the number of paths with zero steps, a(2)=2 is the number of paths with one step, a(3)=5 is the number of paths with two steps, ...
Narrower corridors produce A000012, A000079, A000129, A001519, A057960. An infinitely wide corridor would produce A005773.
(End)
Diagonal sums of A114164. - Paul Barry, Nov 15 2005
C(n):= a(n)*(-1)^n appears in the following formula for the nonpositive powers of rho*sigma, where rho:=2*cos(Pi/7) and sigma:=sin(3*Pi/7)/sin(Pi/7) = rho^2-1 are the ratios of the smaller and larger diagonal length to the side length in a regular 7-gon (heptagon). See the Steinbach reference where the basis <1,rho,sigma> is used in an extension of the rational field. (rho*sigma)^(-n) = C(n) + B(n)*rho + A(n)*sigma,n>=0, with B(n)= A181880(n-2)*(-1)^n, and A(n)= A116423(n+1)*(-1)^(n+1). For the nonnegative powers see A120757(n), |A122600(n-1)| and A181879(n), respectively. See also a comment under A052547.
a(n) is also the number of bi-wall directed polygons with n cells. (The definition of bi-wall directed polygons is given in the article on A122737.)

Crossrefs

Programs

  • Magma
    I:=[1,2,5]; [n le 3 select I[n] else 4*Self(n-1)-3*Self(n-2)-Self(n-3): n in [1..35]]; // Vincenzo Librandi, Sep 18 2015
    
  • Mathematica
    LinearRecurrence[{4,-3,-1}, {1,2,5}, 50] (* Roman Witula, Aug 09 2012 *)
    CoefficientList[Series[(1 - 2 x)/(1 - 4 x + 3 x^2 + x^3), {x, 0, 40}], x] (* Vincenzo Librandi, Sep 18 2015 *)
  • PARI
    x='x+O('x^30); Vec((1-2*x)/(1-4*x+3*x^2+x^3)) \\ G. C. Greubel, Apr 19 2018

Formula

a(n) = 4*a(n-1) - 3*a(n-2) - a(n-3).
From Paul Barry, Nov 15 2005: (Start)
G.f.: (1-2*x)/(1-4*x+3*x^2+x^3).
a(n) = Sum_{k=0..floor(n/2)} (Sum_{j=0..n-k} C(n-k, j)*C(j+k, 2k));
a(n) = Sum_{k=0..floor(n/2)} (Sum_{j=0..n-k} C(n-k, k+j)*C(k, k-j)*2^(n-2k-j));
a(n) = Sum_{k=0..floor(n/2)} (Sum_{j=0..n-2*k} C(n-j, n-2*k-j)*C(k, j)(-1)^j*2^(n-2*k-j)). (End)
a(n-1) = -B(n;-1) = (1/7)*((c(4)-c(1))*(1-c(1))^n + (c(1)-c(2))*(1-c(2))^n + (c(2)-c(4))*(1-c(4))^n), where a(-1):=0, c(j):=2*cos(2*Pi*j/7). Moreover, B(n;d), n in N, d in C, denotes the respective quasi-Fibonacci number defined in comments to A121449 or in Witula-Slota-Warzynski's paper (see also A077998, A006054, A052975, A094789, A121442). - Roman Witula, Aug 09 2012

Extensions

Name corrected and clarified, and offset 1 from Svjetlan Feretic, Jun 01 2013

A094789 Number of (s(0), s(1), ..., s(2n+1)) such that 0 < s(i) < 7 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n+1, s(0) = 1, s(2n+1) = 4.

Original entry on oeis.org

1, 4, 14, 47, 155, 507, 1652, 5373, 17460, 56714, 184183, 598091, 1942071, 6305992, 20475625, 66484244, 215873462, 700937471, 2275930827, 7389902771, 23994866364, 77910846021, 252974934692, 821404463698, 2667083556359
Offset: 1

Views

Author

Herbert Kociemba, Jun 11 2004

Keywords

Comments

In general, a(n) = (2/m)*Sum_{r=1..m-1} sin(r*j*Pi/m)*sin(r*k*Pi/m)*(2*cos(r*Pi/m))^(2n+1) counts (s(0), s(1), ..., s(2n+1)) such that 0 < s(i) < m and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n+1, s(0) = j, s(2n+1) = k.
With interpolated zeros (0,0,0,1,0,4,0,14,...) counts walks of length n between the start and fourth nodes on P_6. - Paul Barry, Jan 26 2005
The Hankel transforms of this sequence or of this sequence with the first term omitted give 1, -2, 1, 1, -2, 1, ... . - Wathek Chammam, Nov 16 2011
Diagonal of the square array A216201. - Philippe Deléham, Mar 28 2013

Crossrefs

Programs

  • Magma
    I:=[1,4,14]; [n le 3 select I[n] else 5*Self(n-1)-6*Self(n-2)+Self(n-3): n in [1..45]]; // Vincenzo Librandi, Nov 10 2014
    
  • Mathematica
    f[n_] := FullSimplify[ TrigToExp[(2/7)Sum[ Sin[Pi*k/7]Sin[4Pi*k/7](2Cos[Pi*k/7])^(2n + 1), {k, 1, 6}]]]; Table[ f[n], {n, 25}] (* Robert G. Wilson v, Jun 18 2004 *)
    LinearRecurrence[{5,-6,1}, {1,4,14}, 50] (* Roman Witula, Aug 09 2012 *)
    CoefficientList[Series[(x - 1) / (- 1 + 5 x - 6 x^2 + x^3), {x, 0, 40}], x] (* Vincenzo Librandi, Nov 10 2014 *)
  • PARI
    Vec(x*(x-1)/(-1 + 5*x - 6*x^2 + x^3) + O(x^40)) \\ Michel Marcus, Nov 10 2014

Formula

a(n) = (2/7)*Sum_{k = 1..6} sin(Pi*k/7)*sin(4*Pi*k/7)*(2*cos(Pi*k/7))^(2n + 1).
a(n) = 5*a(n-1) - 6*a(n-2) + a(n-3).
G.f.: x*(x-1)/(-1 + 5*x - 6*x^2 + x^3). - Corrected by Vincenzo Librandi, Nov 10 2014
a(n) = 2^n*B(n; 1/2) = (1/7)*((c(1) - c(4))*(c(4))^(2n) + (c(2) - c(1))*(c(1))^(2n) + (c(4) - c(2))*(c(2))^(2n)), where c(j) := 2*cos(2*Pi*j/7). Here B(n; d), n in N, d in C denotes the respective quasi-Fibonacci number - see A121449 and Witula-Slota-Warzynski paper for details (see also A052975, A085810, A077998, A006054, A121442). - Roman Witula, Aug 09 2012
a(n+1) = A216201(n,n+2) = A216201(n,n+3). - Philippe Deléham, Mar 28 2013

A052975 Expansion of (1-2*x)*(1-x)/(1-5*x+6*x^2-x^3).

Original entry on oeis.org

1, 2, 6, 19, 61, 197, 638, 2069, 6714, 21794, 70755, 229725, 745889, 2421850, 7863641, 25532994, 82904974, 269190547, 874055885, 2838041117, 9215060822, 29921113293, 97153242650, 315454594314, 1024274628963, 3325798821581, 10798800928441, 35063486341682
Offset: 0

Views

Author

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

Keywords

Comments

Number of (s(0), s(1), ..., s(2n)) such that 0 < s(i) < 7 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n, s(0) = 3, s(2n) = 3. - Herbert Kociemba, Jun 11 2004
Counts all paths of length (2*n), n>=0, starting at the initial node and ending on the nodes 1, 2, 3, 4 and 5 on the path graph P_6, see the second Maple program. - Johannes W. Meijer, May 29 2010

Crossrefs

Programs

  • Magma
    I:=[1,2,6]; [n le 3 select I[n] else 5*Self(n-1)-6*Self(n-2)+Self(n-3): n in [1..30]]; // Vincenzo Librandi, Sep 18 2015
    
  • Maple
    spec := [S,{S=Sequence(Prod(Union(Sequence(Prod(Sequence(Z),Z)),Sequence(Z)),Z))},unlabeled ]: seq(combstruct[count ](spec,size=n), n=0..20);
    with(GraphTheory):G:=PathGraph(6): A:= AdjacencyMatrix(G): nmax:=25; n2:=2*nmax+1: for n from 0 to n2 do B(n):=A^n; a(n):=add(B(n)[k,1],k=1..5); od: seq(a(2*n),n=0..nmax); # Johannes W. Meijer, May 29 2010
  • Mathematica
    LinearRecurrence[{5,-6,1}, {1,2,6}, 50] (* Roman Witula, Aug 09 2012 *)
    CoefficientList[Series[(1 - 2 x) (1 - x)/(1 - 5 x + 6 x^2 - x^3), {x, 0, 40}], x] (* Vincenzo Librandi, Sep 18 2015 *)
  • PARI
    x='x+O('x^30); Vec((1-2*x)*(1-x)/(1-5*x+6*x^2-x^3)) \\ G. C. Greubel, Apr 19 2018

Formula

G.f.: (1-2*x)*(1-x)/(1-5*x+6*x^2-x^3).
a(n) = A028495(2*n). - Floor van Lamoen, Nov 02 2005
a(n) = Sum (1/7*(2-3*_alpha+_alpha^2)*_alpha^(-1-n), _alpha=RootOf(-1+5*_Z-6*_Z^2+_Z^3))
From Herbert Kociemba, Jun 11 2004: (Start)
a(n) = (2/7)*Sum_{r=1..6} sin(r*3*Pi/7)^2*(2*cos(r*Pi/7))^(2*n).
a(n) = 5*a(n-1) - 6*a(n-2) + a(n-3). (End)
a(n) = 2^n*A(n;1/2) = (1/7)*(s(2)^2*c(4)^(2n) + s(4)^2*c(1)^(2n) + s(1)^2*c(2)^(2n)), where c(j):=2*cos(2Pi*j/7) and s(j):=2*sin(2*Pi*j/7). Here A(n;d), n in N, d in C denotes the respective quasi-Fibonacci number - see A121449 and Witula-Slota-Warzynski paper for details (see also A094789, A085810, A077998, A006054, A121442). I note that my and the respective Herbert Kociemba's formulas are "compatible". - Roman Witula, Aug 09 2012
a(n) = A005021(n)-3*A005021(n-1)+2*A005021(n-2). - R. J. Mathar, Feb 27 2019

A121449 Expansion of (1 - 3*x + 2*x^2)/(1 - 4*x + 3*x^2 + x^3).

Original entry on oeis.org

1, 1, 3, 8, 22, 61, 170, 475, 1329, 3721, 10422, 29196, 81797, 229178, 642125, 1799169, 5041123, 14124860, 39576902, 110891905, 310712054, 870595599, 2439354329, 6834918465, 19151015274, 53659951372, 150351841201, 421276495414, 1180390506681, 3307380699281
Offset: 0

Views

Author

Philippe Deléham, Sep 06 2006

Keywords

Comments

From Roman Witula, Aug 07 2012: (Start)
In the cited Witula-Slota-Warzynski paper three so-called quasi-Fibonacci numbers A(n;d), B(n;d) and C(n;d), where n = 0,1,...,d \in C are discussed. These numbers are created by each of the following relations:
(1+d*c(j))^n = A(n;d) + B(n;d)*c(j) + C(n;d)*c(2*j), for every j=1,2,4, where c(j):=2*cos(2*Pi*j/7).
In fact all these "numbers" are integer polynomials of the argument d.
In the sequel for d=-1 we obtain A(n;-1)=a(n), B(n+1;-1)=-A085810(n).
Moreover, we have A(n;1)=A077998(n), B(n;1)=A006054(n+1), C(n;1)=A006054(n), and A(n;2)=A121442(n).
We note that the elements of the sequences A(n;-1), B(n;-1), and C(n;-1) satisfy the following system of recurrence equations:
A(0;-1)=1, B(0;-1)=C(0;-1)=0,
A(n+1;-1)=A(n;-1)-2*B(n;-1)+C(n;-1),
B(n+1;-1)=-A(n;-1)+B(n;-1), C(n+1;-1)=-B(n;-1)+2*C(n;-1).
It is proved that binomial transforms of the sequences: A(n;1), B(n;1), and C(n;1) are equal to the following sequences:
A(n;1)*(A(n;-1)-C(n;-1))-B(n;1)*(B(n;-1)+C(n;-1))+C(n;1)*B(n;-1), -A(n;1)*C(n;-1)+B(n;1)*(A(n;-1)-C(n;-1))-C(n;1)*(B(n;-1)-C(n;-1)), and
A(n;1)*(B(n;-1)-C(n;-1))-B(n;1)*B(n;-1)+C(n;1)*(A(n;-1)-B(n;-1)+C(n;-1)), respectively, whereas we have
A(n;-1) = Sum_{k=0..n} binomial(n,k)*(A(k;1)*A(n-k;1)-A(k;1)*B(n-k;1)-B(k;1)*C(n-k;1)-A(n-k;1)*C(k;1)+2*B(n-k;1)*C(k;1)-C(k;1)*C(n-k;1)),
B(n;-1) = Sum_{k=0..n} binomial(n,k)*(-A(k;1)*B(n-k;1)+A(k;1)*C(n-k;1)+B(k;1)*B(n-k;1)-A(n-k;1)*C(k;1)+B(n-k;1)*C(k;1)-C(k;1)*C(n-k;1)), and
C(n;-1) = Sum_{k=0..n} binomial(n,k)*(-A(k;1)*B(n-k;1)+A(n-k;1)*B(k;1)+B(k;1)*B(n-k;1)-B(k;1)*C(n-k;1)-A(n-k;1)*C(k;1)) (see identities (3.50-52) and (3.61-63) in the Witula-Slota-Warzynski paper).
(End)

Crossrefs

Programs

  • Magma
    I:=[1,1,3]; [n le 3 select I[n] else 4*Self(n-1)-3*Self(n-2)-Self(n-3): n in [1..30]]; // Vincenzo Librandi, Sep 18 2015
    
  • Mathematica
    CoefficientList[Series[(1 - 3*x + 2*x^2)/(1-4*x + 3*x^2 + x^3), {x, 0, 200}], x] (* Stefan Steinerberger, Sep 11 2006 *)
    LinearRecurrence[{4,-3,-1},{1,1,3},50] (* Roman Witula, Aug 07 2012 *)
  • PARI
    x='x+O('x^30); Vec((1-3*x+2*x^2)/(1-4*x+3*x^2+x^3)) \\ G. C. Greubel, Apr 19 2018

Formula

a(0)=a(1)=1, a(2)=3, a(n+1) = 4*a(n) - 3*a(n-1) - a(n-2) for n>=2.
7*a(n) = (2-c(4))*(1-c(1))^n + (2-c(1))*(1-c(2))^n + (2-c(2))*(1-c(4))^n = (s(2))^2*(1-c(1))^n + (s(4))^2*(1-c(2))^n + (s(1))^2*(1-c(4))^n, where c(j):=2*cos(2*Pi*j/7) and s(j):=2*sin(2*Pi*j/7) -- it is the special case, for d=-1, of the Binet's formula for the respective quasi-Fibonacci number A(n;d) discussed in the Witula-Slota-Warzynski paper. - Roman Witula, Aug 07 2012

Extensions

More terms from Stefan Steinerberger, Sep 11 2006
a(27)-a(29) from Vincenzo Librandi, Sep 18 2015

A215404 a(n) = 4*a(n-1) - 3*a(n-2) - a(n-3), with a(0)=0, a(1)=0 and a(2)=1.

Original entry on oeis.org

0, 0, 1, 4, 13, 39, 113, 322, 910, 2561, 7192, 20175, 56563, 158535, 444276, 1244936, 3488381, 9774440, 27387681, 76739023, 215018609, 602469686, 1688083894, 4729907909, 13252910268, 37133833451, 104046695091, 291532369743, 816855560248, 2288778436672, 6413014696201
Offset: 0

Views

Author

Roman Witula, Aug 09 2012

Keywords

Comments

We have a(n)=C(n;-1), A121449(n)=A(n;-1), A085810(n+1)=-B(n+1;-1), where A(n;d), B(n;d), and C(n;d), n in N, d in C, are so-called quasi-Fibonacci numbers defined and discussed in the comments to A121449 and in Witula-Slota-Warzynski's paper. It follows from formulas (3.47-49) in this paper that the value of A(n;1/3), B(n;1/3) and C(n;1/3) could be obtained from special convolution type identities for sequences a(n), A121449, and A085810.

Crossrefs

Programs

  • Magma
    I:=[0,0,1]; [n le 3 select I[n] else 4*Self(n-1)-3*Self(n-2)-Self(n-3): n in [1..35]]; // Vincenzo Librandi, Sep 18 2015
  • Mathematica
    LinearRecurrence[{4,-3,-1}, {0,0,1}, 50]
    CoefficientList[Series[x^2/(1 - 4 x + 3 x^2 + x^3), {x, 0, 40}], x] (* Vincenzo Librandi, Sep 18 2015 *)
  • PARI
    Vec(x^2/(1-4*x+3*x^2+x^3)+O(x^99)) \\ Charles R Greathouse IV, Oct 01 2012
    

Formula

G.f.: x^2/(1-4*x+3*x^2+x^3).
a(n) = (1/7)*((c(2)-c(4))*(1-c(1))^n + (c(4)-c(1))*(1-c(2))^n + (c(1)-c(2))*(1-c(4))^n), where c(j):=2*cos(2*Pi*j/7) - this formula is the Binet formula for a(n) (see the Binet formula (3.17) for the respective quasi-Fibonacci number C(n;d) for value d=-1 in the Witula-Slota-Warzynski paper).

A215492 a(n) = 21*a(n-2) + 7*a(n-3), with a(0)=0, a(1)=3, and a(2)=6.

Original entry on oeis.org

0, 3, 6, 63, 147, 1365, 3528, 29694, 83643, 648270, 1964361, 14199171, 45789471, 311933118, 1060973088, 6871121775, 24463966674, 151720368891, 561841152579, 3357375513429, 12860706786396, 74437773850062, 293576471108319, 1653218198356074, 6686170310225133
Offset: 0

Views

Author

Roman Witula, Aug 13 2012

Keywords

Comments

We have a(n)=B(n;3), where B(n;d), n=1,2,..., d \in C, denote one of the quasi-Fibonacci numbers defined in the comments to A121449 and in the Witula-Slota-Warzynski paper. Its conjugate sequences A(n;3) and C(n;3) are discussed in A121458 and A215484 respectively. Similarly as in A121458 we deduce that each of the following elements a(3*n), a(3*n+1), a(3*n+2) is divided by 3*7^n for every n=0,1,... .

Crossrefs

Programs

  • Magma
    I:=[0,3,6]; [n le 3 select I[n] else 21*Self(n-2)+7*Self(n-3): n in [1..30]]; // Vincenzo Librandi, Sep 18 2015
  • Mathematica
    LinearRecurrence[{0,21,7}, {0,3,6}, 50]
    CoefficientList[Series[(3 x + 6 x^2)/(1 - 21 x^2 - 7 x^3), {x, 0, 33}], x] (* Vincenzo Librandi, Sep 18 2015 *)
  • PARI
    concat(0,Vec((3+6*x)/(1-21*x^2-7*x^3)+O(x^99))) \\ Charles R Greathouse IV, Oct 01 2012
    

Formula

a(n) = (1/7)*((c(1)-c(4))*(1+3*c(1))^n + (c(2)-c(1))*(1+3*c(2))^n + (c(4)-c(2))*(1+3*c(4))^n), where c(j):=2*cos(2*Pi*j/7) (for the proof see Witula-Slota-Warzynski paper).
G.f.: (3*x+6*x^2)/(1-21*x^2-7*x^3).

A271944 Expansion of 2*x*(1 + x)/(1 - x - 9*x^2 + x^3).

Original entry on oeis.org

0, 2, 4, 22, 56, 250, 732, 2926, 9264, 34866, 115316, 419846, 1422824, 5086122, 17471692, 61823966, 213983072, 752927074, 2616950756, 9179311350, 31978941080, 111975792474, 390606950844, 1366410142030, 4769896907152, 16676981234578, 58239643256916
Offset: 0

Views

Author

Vincenzo Librandi, Apr 18 2016

Keywords

Crossrefs

Cf. A121442.

Programs

  • Magma
    [n le 2 select 2*n else Self(n)+9*Self(n-1)-Self(n-2): n in [0..30]];
    
  • Mathematica
    CoefficientList[Series[2 x (1 + x)/(1 - x - 9 x^2 + x^3), {x, 0, 33}], x]
  • PARI
    x='x+O('x^99); concat(0, Vec(2*x*(1+x)/(1-x-9*x^2+x^3))) \\ Altug Alkan, Apr 18 2016
    
  • Sage
    gf = 2*x*(1+x)/(1-x-9*x^2+x^3); taylor(gf, x, 0, 40).list() # Bruno Berselli, Apr 18 2016

Formula

G.f.: 2*x*(1 + x)/(1 - x - 9*x^2 + x^3).
a(n) = a(n-1) + 9*a(n-2) - a(n-3).

A271945 Expansion of 4*x^2/(1 - x - 9*x^2 + x^3).

Original entry on oeis.org

0, 0, 4, 4, 40, 72, 428, 1036, 4816, 13712, 56020, 174612, 665080, 2180568, 7991676, 26951708, 96696224, 331269920, 1174584228, 4059317284, 14299305416, 49658576744, 174293008204, 606920893484, 2125899390576, 7413894423728, 25940068045428, 90539218468404
Offset: 0

Views

Author

Vincenzo Librandi, Apr 18 2016

Keywords

Crossrefs

Programs

  • Magma
    [n le 2 select 2*n*(n-1) else Self(n)+9*Self(n-1)-Self(n-2): n in [0..30]];
    
  • Mathematica
    CoefficientList[Series[4 x^2 /(1 - x - 9 x^2 + x^3), {x, 0, 30}], x]
    LinearRecurrence[{1,9,-1},{0,0,4},30] (* Harvey P. Dale, Jul 18 2021 *)
  • PARI
    x='x+O('x^99); concat([0, 0], Vec(4*x^2/(1-x-9*x^2+x^3))) \\ Altug Alkan, Apr 18 2016
    
  • Sage
    gf = 4*x^2/(1 - x - 9*x^2 + x^3); taylor(gf, x, 0, 30).list() # Bruno Berselli, Apr 18 2016

Formula

G.f.: 4*x^2/(1 - x - 9*x^2 + x^3).
a(n) = a(n-1) + 9*a(n-2) - a(n-3).
Showing 1-8 of 8 results.