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.

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

A164073 a(n) = 2*a(n-2) for n > 2; a(1) = 1, a(2) = 3.

Original entry on oeis.org

1, 3, 2, 6, 4, 12, 8, 24, 16, 48, 32, 96, 64, 192, 128, 384, 256, 768, 512, 1536, 1024, 3072, 2048, 6144, 4096, 12288, 8192, 24576, 16384, 49152, 32768, 98304, 65536, 196608, 131072, 393216, 262144, 786432, 524288, 1572864, 1048576, 3145728, 2097152, 6291456
Offset: 1

Views

Author

Klaus Brockhaus, Aug 09 2009

Keywords

Comments

Interleaving of A000079 and A007283.
Binomial transform is A048654. Second binomial transform is A111567. Third binomial transform is A081179 without initial 0. Fourth binomial transform is A164072. Fifth binomial transform is A164031.
Absolute second differences are the sequence itself. - Eric Angelini, Jul 30 2013
Least number having n - 1 Gaussian prime factors, counted with multiplicity, excluding units. See A239628 for a similar sequence. - T. D. Noe, Mar 31 2014
Writing the prime factorizations of the terms of this sequence, the exponents of prime factor 2 give the integers repeated (A004526), while the exponents of prime factor 3 give the sequence of alternating 0's and 1's (A000035). - Alonso del Arte, Nov 30 2016

Crossrefs

Programs

  • Magma
    [ n le 2 select 2*n-1 else 2*Self(n-2): n in [1..42] ];
    
  • Mathematica
    terms = 50; CoefficientList[Series[x * (1 + 3 * x)/(1 - 2 * x^2), {x, 0, terms}], x] (* T. D. Noe, Mar 31 2014 *)
    Flatten[Table[{2^n, 3 * 2^n}, {n, 0, 31}]] (* Alonso del Arte, Nov 30 2016 *)
    CoefficientList[Series[x (1 + 3 x)/(1 - 2 x^2), {x, 0, 44}], x] (* Michael De Vlieger, Dec 13 2016 *)
  • PARI
    a(n) = (5 + (-1)^n) * 2^((2*n-9)\/4)
    
  • PARI
    Vec(x*(1+3*x)/(1-2*x^2)+O(x^99)) \\ Charles R Greathouse IV, Dec 13 2016

Formula

a(n) = (5 + (-1)^n) * 2^(1/4 * (2*n - 1 + (-1)^n))/4.
G.f.: x*(1 + 3 * x)/(1 - 2 * x^2).
a(n) = A074323(n), n>=1.
a(n) = A162255(n-1), n>=2.
a(n) = A072946(n-2), n > 2. - R. J. Mathar, Aug 17 2009
a(n+3) = a(n + 2) * a(n + 1)/a(n). - Reinhard Zumkeller, Mar 04 2011
a(n) = (2/3)a(n - 1) for odd n > 1, a(n) = 3a(n - 1) for even n. - Alonso del Arte, Nov 30 2016

A111566 a(n) = ((1+sqrt(8))*(2+sqrt(2))^n + (1-sqrt(8))*(2-sqrt(2))^n)/2.

Original entry on oeis.org

1, 6, 22, 76, 260, 888, 3032, 10352, 35344, 120672, 412000, 1406656, 4802624, 16397184, 55983488, 191139584, 652591360, 2228086272, 7607162368, 25972476928, 88675582976, 302757378048, 1033678346240, 3529198628864, 12049437822976, 41139354034176, 140458540490752
Offset: 0

Views

Author

Creighton Dement, Aug 06 2005

Keywords

Comments

Binomial transform of A048655: generalized Pellian with second term equal to 5.
Floretion Algebra Multiplication Program, FAMP Code: 1vesseq[K*J] with K = + .5'i + .5'j + .5k' + .5'kk' and J = + .5i' + .5j' + 2'kk' + .5'ki' + .5'kj'.

Crossrefs

Programs

  • Magma
    Z:=PolynomialRing(Integers()); N:=NumberField(x^2-2); S:=[ ((1+2*r)*(2+r)^n+(1-2*r)*(2-r)^n)/2: n in [0..23] ]; [ Integers()!S[j]: j in [1..#S] ]; // Klaus Brockhaus, Jul 27 2009
    
  • Mathematica
    LinearRecurrence[{4,-2},{1,6},30] (* Harvey P. Dale, Jan 31 2015 *)
  • PARI
    x='x+O('x^30); Vec((1+2*x)/(1-4*x+2*x^2)) \\ G. C. Greubel, Jan 27 2018

Formula

a(n) = 4*a(n-1) - 2*a(n-2), a(0) = 1, a(1) = 6.
Program "FAMP" returns: a(n) = A007052(n) - A006012(n) + A111567(n).
From R. J. Mathar, Apr 02 2008: (Start)
O.g.f.: (1+2*x)/(1-4*x+2*x^2).
a(n) = A007070(n) + 2*A007070(n-1). (End)
a(n) = Sum_{k=0..n} A207543(n,k)*2^k. - Philippe Deléham, Feb 25 2012
a(n) = 4*A007070(n) - A007052(n+1). - Yuriy Sibirmovsky, Sep 13 2016
E.g.f.: exp(2*x)*(cosh(sqrt(2)*x) + 2*sqrt(2)*sinh(sqrt(2)*x)). - Stefano Spezia, May 26 2024

Extensions

Edited by N. J. A. Sloane, Jul 27 2009 using new definition from Al Hakanson (hawkuu(AT)gmail.com)

A108012 a(n)= 8*a(n-1) -16*a(n-2) +4*a(n-4).

Original entry on oeis.org

0, 4, 18, 76, 320, 1360, 5832, 25200, 109568, 478784, 2100512, 9244352, 40784896, 180284672, 798121088, 3537391360, 15692333056, 69661541376, 309407486464, 1374824795136, 6110847909888, 27168232722432, 120809925167104
Offset: 0

Views

Author

Roger L. Bagula, May 30 2005

Keywords

Programs

  • Mathematica
    M = {{0, 0, 0, 2}, {1, 4, 0, 0}, {0, 2, 0, 0}, {0, 0, 1, 4}} v[1] = {0, 1, 1, 2}; v[n_] := v[n] = M.v[n - 1]; digits = 50; a = Table[v[n][[1]], {n, 1, digits}]
    LinearRecurrence[{8,-16,0,4},{0,4,18,76},30] (* Harvey P. Dale, Aug 11 2017 *)

Formula

G.f.: 2*x*(-2+7*x+2*x^2)/( (2*x^2+4*x-1) * (2*x^2-4*x+1)). [Sep 28 2009]
a(n) = A111567(n-1) + A090017(n+1) - A090017(n) - A090017(n-1), n>0.

Extensions

Definition replaced by recurrence by the Associate Editors of the OEIS, Sep 28 2009

A113859 Expansion of (7-14*x+6*x^2)/((1-x)*(2*x^2-4*x+1)); related to the binomial transform of Pell numbers A000129 (see formula and comment for A007070).

Original entry on oeis.org

7, 21, 69, 233, 793, 2705, 9233, 31521, 107617, 367425, 1254465, 4283009, 14623105, 49926401, 170459393, 581984769, 1987020289, 6784111617, 23162405889, 79081400321, 270000789505, 921840357377, 3147359850497, 10745758687233
Offset: 0

Views

Author

Creighton Dement, Jan 25 2006

Keywords

Comments

If g.f. (x^6+5*x^4+6*x^2+1)/(x^7+6*x^5+10*x^3+4*x) is expanded, where (x^6+5*x^4+6*x^2+1) and (x^7+6*x^5+10*x^3+4*x) are the 7th and 8th Fibonacci polynomials, respectively, the sequence: [0, 7/8, 0, -21/16, 0, 69/32, 0, -233/64, 0, 793/128, 0, -2705/256, ] is returned. (a(n)) is seen to be the numerators of the bisection of this sequences, apart from signs.

Crossrefs

Programs

  • Maple
    with(combinat, fibonacci): seq(fibonacci(i, x), i=1..15); [[generates sequence of Fibonacci polynomials]]

Formula

a(n+1) - a(n) = A007070(n+2), a(n) - 2*a(n+1) + a(n+2) = A007052(n+3) (Number of order consecutive partitions of n), a(n+3) - 3*a(n+2) + 3*a(n+1) - a(n) = A003480(n+4), a(n+2) - a(n) = A111567(n+3)
Showing 1-5 of 5 results.