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-10 of 126 results. Next

A080923 First differences of A003946.

Original entry on oeis.org

1, 3, 8, 24, 72, 216, 648, 1944, 5832, 17496, 52488, 157464, 472392, 1417176, 4251528, 12754584, 38263752, 114791256, 344373768, 1033121304, 3099363912, 9298091736, 27894275208, 83682825624, 251048476872, 753145430616
Offset: 0

Views

Author

Paul Barry, Feb 26 2003

Keywords

Comments

Sum of consecutive pairs of elements of A025192.
The alternating sign sequence with g.f. (1-x^2)/(1+3x) gives the diagonal sums of A110168. - Paul Barry, Jul 14 2005
Let M = an infinite lower triangular matrix with the odd integers (1,3,5,...) in every column, with the leftmost column shifted up one row. Then A080923 = lim_{n->inf} M^n. - Gary W. Adamson, Feb 18 2010
a(n+1), n >= 0, with o.g.f. ((1-x^2)/(1-3*x)-1)/x = (3-x)/(1-3*x) provides the coefficients in the formal power series for tan(3*x)/tan(x) = (3-z)/(1-3*z) = Sum_{n>=0} a(n+1)*z^n, with z = tan(x)^2. Convergence holds for 0 <= z < 1/3, i.e., |x| < Pi/6, approximately 0.5235987758. For the numerator and denominator of this o.g.f. see A034867 and A034839, respectively. - Wolfdieter Lang, Jan 18 2013

Crossrefs

Essentially the same as A005051, A026097 and A083583.

Programs

  • Mathematica
    CoefficientList[Series[(1 - x^2) / (1 - 3 x), {x, 0, 20}], x] (* Vincenzo Librandi, Aug 05 2013 *)

Formula

G.f.: (1-x^2)/(1-3*x).
G.f.: 1/(1 - 3*x + x^2 - 3*x^3 + x^4 - 3*x^5 + ...). - Gary W. Adamson, Jan 06 2011
a(n) = 2^3*3^(n-2), n >= 2, a(0) = 1, a(1) = 3. - Wolfdieter Lang, Jan 18 2013

A003945 Expansion of g.f. (1+x)/(1-2*x).

Original entry on oeis.org

1, 3, 6, 12, 24, 48, 96, 192, 384, 768, 1536, 3072, 6144, 12288, 24576, 49152, 98304, 196608, 393216, 786432, 1572864, 3145728, 6291456, 12582912, 25165824, 50331648, 100663296, 201326592, 402653184, 805306368, 1610612736, 3221225472, 6442450944, 12884901888
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 3.
Number of Hamiltonian cycles in K_3 X P_n.
Number of ternary words of length n avoiding aa, bb, cc.
For n > 0, row sums of A029635. - Paul Barry, Jan 30 2005
Binomial transform is {1, 4, 13, 40, 121, 364, ...}, see A003462. - Philippe Deléham, Jul 23 2005
Convolved with the Jacobsthal sequence A001045 = A001786: (1, 4, 12, 32, 80, ...). - Gary W. Adamson, May 23 2009
Equals (n+1)-th row sums of triangle A161175. - Gary W. Adamson, Jun 05 2009
a(n) written in base 2: a(0) = 1, a(n) for n >= 1: 11, 110, 11000, 110000, ..., i.e.: 2 times 1, (n-1) times 0 (see A003953(n)). - Jaroslav Krizek, Aug 17 2009
INVERTi transform of A003688. - Gary W. Adamson, Aug 05 2010
An elephant sequence, see A175655. For the central square four A[5] vectors, with decimal values 42, 138, 162 and 168, lead to this sequence. For the corner squares these vectors lead to the companion sequence A083329. - Johannes W. Meijer, Aug 15 2010
A216022(a(n)) != 2 and A216059(a(n)) != 3. - Reinhard Zumkeller, Sep 01 2012
Number of length-n strings of 3 letters with no two adjacent letters identical. The general case (strings of r letters) is the sequence with g.f. (1+x)/(1-(r-1)*x). - Joerg Arndt, Oct 11 2012
Sums of pairs of rows of Pascal's triangle A007318, T(2n,k)+T(2n+1,k); Sum_{n>=1} A000290(n)/a(n) = 4. - John Molokach, Sep 26 2013

Crossrefs

Essentially same as A007283 (3*2^n) and A042950.
Generating functions of the form (1+x)/(1-k*x) for k=1 to 12: A040000, A003945, A003946, A003947, A003948, A003949, A003950, A003951, A003952.
Generating functions of the form (1+x)/(1-k*x) for k=13 to 30: A170732, A170733, A170734, A170735, A170736, A170737, A170738, A170739, A170740, A170741, A170742, A170743, A170744, A170745, A170746, A170747, A170748.
Generating functions of the form (1+x)/(1-k*x) for k=31 to 50: A170749, A170750, A170751, A170752, A170753, A170754, A170755, A170756, A170757, A170758, A170759, A170760, A170761, A170762, A170763, A170764, A170765, A170766, A170767, A170768, A170769.
Cf. A003688.

Programs

  • Maple
    k := 3; if n = 0 then 1 else k*(k-1)^(n-1); fi;
  • Mathematica
    Join[{1}, 3*2^Range[0, 60]] (* Vladimir Joseph Stephan Orlovsky, Jun 09 2011 *)
    Table[2^n+Floor[2^(n-1)], {n,0,30}] (* Martin Grymel, Oct 17 2012 *)
    CoefficientList[Series[(1+x)/(1-2x),{x,0,40}],x] (* or *) LinearRecurrence[ {2},{1,3},40] (* Harvey P. Dale, May 04 2017 *)
  • PARI
    a(n)=if(n,3<Charles R Greathouse IV, Jan 12 2012

Formula

a(0) = 1; for n > 0, a(n) = 3*2^(n-1).
a(n) = 2*a(n-1), n > 1; a(0)=1, a(1)=3.
More generally, the g.f. (1+x)/(1-k*x) produces the sequence [1, 1 + k, (1 + k)*k, (1 + k)*k^2, ..., (1+k)*k^(n-1), ...], with a(0) = 1, a(n) = (1+k)*k^(n-1) for n >= 1. Also a(n+1) = k*a(n) for n >= 1. - Zak Seidov and N. J. A. Sloane, Dec 05 2009
The g.f. (1+x)/(1-k*x) produces the sequence with closed form (in PARI notation) a(n)=(n>=0)*k^n+(n>=1)*k^(n-1). - Jaume Oliver Lafont, Dec 05 2009
Binomial transform of A000034. a(n) = (3*2^n - 0^n)/2. - Paul Barry, Apr 29 2003
a(n) = Sum_{k=0..n} (n+k)*binomial(n, k)/n. - Paul Barry, Jan 30 2005
a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 1. - Philippe Deléham, Jul 10 2005
Binomial transform of A000034. Hankel transform is {1,-3,0,0,0,...}. - Paul Barry, Aug 29 2006
a(0) = 1, a(n) = 2 + Sum_{k=0..n-1} a(k) for n >= 1. - Joerg Arndt, Aug 15 2012
a(n) = 2^n + floor(2^(n-1)). - Martin Grymel, Oct 17 2012
E.g.f.: (3*exp(2*x) - 1)/2. - Stefano Spezia, Jan 31 2023

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A025192 a(0)=1; a(n) = 2*3^(n-1) for n >= 1.

Original entry on oeis.org

1, 2, 6, 18, 54, 162, 486, 1458, 4374, 13122, 39366, 118098, 354294, 1062882, 3188646, 9565938, 28697814, 86093442, 258280326, 774840978, 2324522934, 6973568802, 20920706406, 62762119218, 188286357654, 564859072962, 1694577218886, 5083731656658, 15251194969974
Offset: 0

Views

Author

Keywords

Comments

Warning: there is considerable overlap between this entry and the essentially identical A008776.
Shifts one place left when plus-convolved (PLUSCONV) with itself. a(n) = 2*Sum_{i=0..n-1} a(i). - Antti Karttunen, May 15 2001
Let M = { 0, 1, ..., 2^n-1 } be the set of all n-bit numbers. Consider two operations on this set: "sum modulo 2^n" (+) and "bitwise exclusive or" (XOR). The results of these operations are correlated.
To give a numerical measure, consider the equations over M: u = x + y, v = x XOR y and ask for how many pairs (u,v) is there a solution? The answer is exactly a(n) = 2*3^(n-1) for n >= 1. The fraction a(n)/4^n of such pairs vanishes as n goes to infinity. - Max Alekseyev, Feb 26 2003
Number of (s(0), s(1), ..., s(2n+2)) such that 0 < s(i) < 6 and |s(i) - s(i-1)| = 1 for i = 1,2,...,2n+2, s(0) = 3, s(2n+2) = 3. - Herbert Kociemba, Jun 10 2004
Number of compositions of n into parts of two kinds. For a string of n objects, before the first, choose first kind or second kind; before each subsequent object, choose continue, first kind, or second kind. For example, compositions of 3 are 3; 2,1; 1,2; and 1,1,1. Using parts of two kinds, these produce respectively 2, 4, 4 and 8 compositions, 2+4+4+8 = 18. - Franklin T. Adams-Watters, Aug 18 2006
In the compositions the kinds of parts are ordered inside a run of identical parts, see example. Replacing "ordered" by "unordered" gives A052945. - Joerg Arndt, Apr 28 2013
Number of permutations of {1, 2, ..., n+1} such that no term is more than 2 larger than its predecessor. For example, a(3) = 18 because all permutations of {1, 2, 3, 4} are valid except 1423, 1432, 2143, 3142, 2314, 3214, in which 1 is followed by 4. Proof: removing (n + 1) gives a still-valid sequence. For n >= 2, can insert (n + 1) either at the beginning or immediately following n or immediately following (n - 1), but nowhere else. Thus the number of such permutations triples when we increase the sequence length by 1. - Joel B. Lewis, Nov 14 2006
Antidiagonal sums of square array A081277. - Philippe Deléham, Dec 04 2006
Equals row sums of triangle A160760. - Gary W. Adamson, May 25 2009
Let M = a triangle with (1, 2, 4, 8, ...) as the left border and all other columns = (0, 1, 2, 4, 8, ...). A025192 = lim_{n->oo} M^n, the left-shifted vector considered as a sequence. - Gary W. Adamson, Jul 27 2010
Number of nonisomorphic graded posets with 0 and uniform hasse graph of rank n with no 3-element antichain. ("Uniform" used in the sense of Retakh, Serconek and Wilson. By "graded" we mean that all maximal chains have the same length n.) - David Nacin, Feb 13 2012
Equals partial sums of A003946 prefaced with a 1: (1, 1, 4, 12, 36, 108, ...). - Gary W. Adamson, Feb 15 2012
Number of vertices (or sides) of the (n-1)-th iteration of a Gosper island. - Arkadiusz Wesolowski, Feb 07 2013
Row sums of triangle in A035002. - Jon Perry, May 30 2013
a(n) counts walks (closed) on the graph G(1-vertex; 1-loop, 1-loop, 2-loop, 2-loop, 3-loop, 3-loop, ...). - David Neil McGrath, Jan 01 2015
From Tom Copeland, Dec 03 2015: (Start)
For n > 0, a(n) are the traces of the even powers of the adjacency matrix M of the simple Lie algebra B_3, tr(M^(2n)) where M = Matrix(row 1; row 2; row 3) = Matrix[0,1,0; 1,0,2; 0,1,0], same as the traces of Matrix[0,2,0; 1,0,1; 0,1,0] (cf. Damianou). The traces of the odd powers vanish.
The characteristic polynomial of M equals determinant(x*I - M) = x^3 - 3x = A127672(3,x), so 1 - 3*x^2 = det(I - x M) = exp(-Sum_{n>=1} tr(M^n) x^n / n), implying Sum_{n>=1} a(n+1) x^(2n) / (2n) = -log(1 - 3*x^2), giving a logarithmic generating function for the aerated sequence, excluding a(0) and a(1).
a(n+1) = tr(M^(2n)), where tr(M^n) = 3^(n/2) + (-1)^n * 3^(n/2) = 2^n*(cos(Pi/6)^n + cos(5*Pi/6)^n) = n-th power sum of the eigenvalues of M = n-th power sum of the zeros of the characteristic polynomial.
The relation det(I - x M) = exp(-Sum_{n>=1} tr(M^n) x^n / n) = Sum_{n>=0} P_n(-tr(M), -tr(M^2), ..., -tr(M^n)) x^n/n! = exp(P.(-tr(M), -tr(M^2), ...)x), where P_n(x(1), ..., x(n)) are the partition polynomials of A036039 implies that with x(2n) = -tr(M^(2n)) = -a(n+1) for n > 0 and x(n) = 0 otherwise, the partition polynomials evaluate to zero except for P_2(x(1), x(2)) = P_2(0,-6) = -6.
Because of the inverse relation between the partition polynomials of A036039 and the Faber polynomials F_k(b1,b2,...,bk) of A263916, F_k(0,-3,0,0,...) = tr(M^k) gives aerated a(n), excluding n=0,1. E.g., F_2(0,-3) = -2(-3) = 6, F_4(0,-3,0,0) = 2 (-3)^2 = 18, and F_6(0,-3,0,0,0,0) = -2(-3)^3 = 54. (Cf. A265185.)
(End)
Number of permutations of length n > 0 avoiding the partially ordered pattern (POP) {1>2, 1>3, 1>4} of length 4. That is, number of length n permutations having no subsequences of length 4 in which the first element is the largest. - Sergey Kitaev, Dec 08 2020
For n > 0, a(n) is the number of 3-colorings of the grid graph P_2 X P_(n-1). More generally, for q > 1, the number of q-colorings of the grid graph P_2 X P_n is given by q*(q - 1)*((q - 1)*(q - 2) + 1)^(n - 1). - Sela Fried, Sep 25 2023
For n > 1, a(n) is the largest solution to the equation phi(x) = a(n-1). - M. Farrokhi D. G., Oct 25 2023
Number of dotted compositions of degree n. - Diego Arcis, Feb 01 2024

Examples

			There are a(3)=18 compositions of 3 into 2 kinds of parts. Here p:s stands for "part p of sort s":
01:  [ 1:0  1:0  1:0  ]
02:  [ 1:0  1:0  1:1  ]
03:  [ 1:0  1:1  1:0  ]
04:  [ 1:0  1:1  1:1  ]
05:  [ 1:0  2:0  ]
06:  [ 1:0  2:1  ]
07:  [ 1:1  1:0  1:0  ]
08:  [ 1:1  1:0  1:1  ]
09:  [ 1:1  1:1  1:0  ]
10:  [ 1:1  1:1  1:1  ]
11:  [ 1:1  2:0  ]
12:  [ 1:1  2:1  ]
13:  [ 2:0  1:0  ]
14:  [ 2:0  1:1  ]
15:  [ 2:1  1:0  ]
16:  [ 2:1  1:1  ]
17:  [ 3:0  ]
18:  [ 3:1  ]
- _Joerg Arndt_, Apr 28 2013
G.f. = 1 + 2*x + 6*x^2 + 18*x^3 + 54*x^4 + 162*x^5 + 486*x^6 + 1458*x^7 + ...
		

References

  • Richard P. Stanley, Enumerative combinatorics, Vol. 1, Cambridge University Press, Cambridge, 1997, pp. 96-100.

Crossrefs

First differences of 3^n (A000244). Other self-convolved sequences: A000108, A007460, A007461, A007462, A007463, A007464, A061922.
Apart from initial term, same as A008776.

Programs

  • Haskell
    a025192 0 = 1
    a025192 n = 2 * 3 ^ (n -1)
    a025192_list = 1 : iterate (* 3) 2  -- Reinhard Zumkeller, Nov 27 2012
  • Maple
    A025192 := proc(n): if n=0 then 1 else 2*3^(n-1) fi: end: seq(A025192(n),n=0..26);
  • Mathematica
    Join[{1},2*3^(Range[30]-1)]  (* Harvey P. Dale, Mar 22 2011 *)
  • PARI
    a(n)=max(1,2*3^(n-1)) \\ Charles R Greathouse IV, Jul 25 2011
    
  • PARI
    Vec((1-x)/(1-3*x) + O(x^100)) \\ Altug Alkan, Dec 05 2015
    
  • Python
    [1]+[2*3**(n-1) for n in range(1,30)] # David Nacin, Mar 04 2012
    

Formula

G.f.: (1-x)/(1-3*x).
E.g.f.: (2*exp(3*x) + exp(0))/3. - Paul Barry, Apr 20 2003
a(n) = phi(3^n) = A000010(A000244(n)). - Labos Elemer, Apr 14 2003
a(0) = 1, a(n) = Sum_{k=0..n-1} (a(k) + a(n-k-1)). - Benoit Cloitre, Jun 24 2003
a(n) = A002326((3^n-1)/2). - Vladimir Shevelev, May 26 2008
a(1) = 2, a(n) = 3*a(n-1). - Vincenzo Librandi, Jan 01 2011
a(n) = lcm(a(n-1), Sum_{k=1..n-1} a(k)) for n >= 3. - David W. Wilson, Sep 27 2011
a(n) = ((2*n-1)*a(n-1) + (3*n-6)*a(n-2))/(n-1); a(0)=1, a(1)=2. - Sergei N. Gladkovskii, Jul 16 2012
From Sergei N. Gladkovskii, Jul 17 2012: (Start)
For the e.g.f. E(x) = (2/3)*exp(3*x) + exp(0)/3 we have
E(x) = 2*G(0)/3 where G(k) = 1 + k!/(3*(9*x)^k - 3*(9*x)^(2*k+1)/((9*x)^(k+1) + (k+1)!/G(k+1))); (continued fraction, 3rd kind, 3-step).
E(x) = 1+2*x/(G(0)-3*x) where G(k) = 3*x + 1 + k - 3*x*(k+1)/G(k+1); (continued fraction, Euler's 1st kind, 1-step). (End)
a(n) = A114283(0,0). - Reinhard Zumkeller, Nov 27 2012
G.f.: 1 + ((1/2)/G(0) - 1)/x where G(k) = 1 - 2^k/(2 - 4*x/(2*x - 2^k/G(k+1) )); (recursively defined continued fraction). - Sergei N. Gladkovskii, Dec 22 2012
G.f.: 1 + x*W(0), where W(k) = 1 + 1/(1 - x*(2*k+3)/(x*(2*k+4) + 1/W(k+1))); (continued fraction). - Sergei N. Gladkovskii, Aug 28 2013
G.f.: 1 / (1 - 2*x / (1 - x)). - Michael Somos, Apr 03 2014
Construct the power matrix T(n,j) = [A(n)^*j]*[S(n)^*(j-1)] where A(n)=(2,2,2,...) and S(n)=(0,1,0,0,...). (* is convolution operation.) Then a(n) = Sum_{j=1..n} T(n,j). - David Neil McGrath, Jan 01 2015
G.f.: 1 + 2*x/(1 + 2*x)*( 1 + 5*x/(1 + 5*x)*( 1 + 8*x/(1 + 8*x)*( 1 + 11*x/(1 + 11*x)*( 1 + .... - Peter Bala, May 27 2017
Sum_{n>=0} 1/a(n) = 7/4. - Bernard Schott, Oct 02 2021
From Amiram Eldar, May 08 2023: (Start)
Sum_{n>=0} (-1)^n/a(n) = 5/8.
Product_{n>=1} (1 - 1/a(n)) = A132019. (End)

Extensions

Additional comments from Barry E. Williams, May 27 2000
a(22) corrected by T. D. Noe, Feb 08 2008
Maple programs simplified by Johannes W. Meijer, Jun 02 2011

A015448 a(0) = 1, a(1) = 1, and a(n) = 4*a(n-1) + a(n-2) for n >= 2.

Original entry on oeis.org

1, 1, 5, 21, 89, 377, 1597, 6765, 28657, 121393, 514229, 2178309, 9227465, 39088169, 165580141, 701408733, 2971215073, 12586269025, 53316291173, 225851433717, 956722026041, 4052739537881, 17167680177565, 72723460248141, 308061521170129, 1304969544928657, 5527939700884757
Offset: 0

Views

Author

Keywords

Comments

If one deletes the leading 0 in A084326, takes the inverse binomial transform, and adds a(0)=1 in front, one obtains this sequence here. - Al Hakanson (hawkuu(AT)gmail.com), May 02 2009
For n >= 1, row sums of triangle
m |k=0 1 2 3 4 5 6 7
====+=============================================
0 | 1
1 | 1 4
2 | 1 4 16
3 | 1 8 16 64
4 | 1 8 48 64 256
5 | 1 12 48 256 256 1024
6 | 1 12 96 256 1280 1024 4096
7 | 1 16 96 640 1280 6144 4096 16384
which is triangle for numbers 4^k*C(m,k) with duplicated diagonals. - Vladimir Shevelev, Apr 12 2012
a(n) = a(n;-2) = 3^n*Sum_{k=0..n} binomial(n,k)*F(k+1)*(-2/3)^k, where a(n;d), n=0,1,...,d, denotes the delta-Fibonacci numbers defined in comments to A000045 (see also the papers of Witula et al.). We note that (see A033887) F(3n+1) = 3^n*a(n,2/3) = Sum_{k=0..n} binomial(n,k)*F(k-1)*(-2/3)^k, which implies F(3n+1) + 3^(-n)*a(n) = Sum_{k=0..n} binomial(n,k)*L(k)*(-2/3)^k, where L(k) denotes the k-th Lucas number. - Roman Witula, Jul 12 2012
a(n+1) is (for n >= 0) the number of length-n strings of 5 letters {0,1,2,3,4} with no two adjacent nonzero letters identical. The general case (strings of L letters) is the sequence with g.f. (1+x)/(1-(L-1)*x-x^2). - Joerg Arndt, Oct 11 2012
Starting with offset 1 the sequence is the INVERT transform of (1, 4, 4*3, 4*3^2, 4*3^3, ...); i.e., of A003946: (1, 4, 12, 36, 108, ...). - Gary W. Adamson, Aug 06 2016
a(n+1) equals the number of quinary sequences of length n such that no two consecutive terms differ by 3. - David Nacin, May 31 2017

Crossrefs

Cf. A001076, A147722 (INVERT transform), A109499 (INVERTi transform), A154626 (Binomial transform), A086344 (inverse binomial transform), A003946, A049310.

Programs

Formula

a(n) = Fibonacci(3n-1) = ( (1+sqrt(5))*(2-sqrt(5))^n - (1-sqrt(5))*(2+sqrt(5))^n )/ (2*sqrt(5)).
O.g.f.: (1-3*x)/(1-4*x-x^2). - Len Smiley, Dec 09 2001
a(n) = Sum_{k=0..n} 3^k*A055830(n,k). - Philippe Deléham, Oct 18 2006
a(n) = upper left term in the 2 X 2 matrix [1,2; 2,3]^n. - Gary W. Adamson, Mar 02 2008
[a(n), A001076(n)] = [1,4; 1,3]^n * [1,0]. - Gary W. Adamson, Mar 21 2008
a(n) = A167808(3*n-1) for n > 0. - Reinhard Zumkeller, Nov 12 2009
a(n) = Fibonacci(3n+1) mod Fibonacci(3n), n > 0.
a(n) = (A000032(3*n)-Fibonacci(3*n))/2 = (A014448(n)-A014445(n))/2.
For n >= 2, a(n) = F_n(4) + F_(n+1)(4), where F_n(x) is a Fibonacci polynomial (cf. A049310): F_n(x) = Sum_{i=0..floor((n-1)/2)} binomial(n-i-1,i)*x^(n-2*i-1). - Vladimir Shevelev, Apr 13 2012
a(n) = A001076(n+1) - 3*A001076(n). - R. J. Mathar, Jul 12 2012
From Gary Detlefs and Wolfdieter Lang, Aug 20 2012: (Start)
a(n) = (5*F(n)^3 + 5*F(n-1)^3 + 3*(-1)^n*F(n-2))/2,
a(n) = (F(n+1)^3 + 2*F(n)^3 - F(n-2)^3)/2, n >= 0, with F(-1) = 1 and F(-2) = -1. Second line from first one with 3*(-1)^n* F(n-2) = F(n-1)^3 - 4*F(n-2)^3 - F(n-3)^3 (in Koshy's book, p. 89, 32. (with a - sign) and 33. For the Koshy reference see A000045) and the F^3 recurrence (see row n=4 of A055870, or Koshy p. 87, 1.). First line from the preceding R. J. Mathar formula with F(3*n) = 5*F(n)^3 + 3*(-1)^n*F(n) (Koshy p. 89, 46.) and the above mentioned formula, Koshy's 32. and 33., with n -> n+2 in order to eliminate F(n+1)^3. (End)
For n > 0, a(n) = L(n-1)*L(n)*F(n) + F(n+1)*(-1)^n with L(n)=A000032(n). - J. M. Bergot, Dec 10 2015
For n > 1, a(n)^2 is the denominator of continued fraction [4,4,...,4, 6, 4,4,...4], which has n-1 4's before, and n-1 4's after, the middle 6. - Greg Dresden, Sep 18 2019
From Gary Detlefs and Wolfdieter Lang, Mar 06 2023: (Start)
a(n) = A001076(n) + A001076(n-1), with A001076(-1) = 1. See the R. J. Mathar formula above.
a(n+1) = i^n*(S(n-1,-4*i) - i*S(n-2,-4*i)), for n >= 0, with i = sqrt(-1), and the Chebyshev S-polynomials (see A049310) with S(n, -1) = 0. From the simplified Fibonacci trisection formula for {F(3*n+2)}_{n>=0}. (End)
a(n) = Sum_{k=0..n} A046854(n-1,k)*4^k. - R. J. Mathar, Feb 10 2024
E.g.f.: exp(2*x)*(5*cosh(sqrt(5)*x) - sqrt(5)*sinh(sqrt(5)*x))/5. - Stefano Spezia, Jun 03 2024

A003948 Expansion of (1+x)/(1-5*x).

Original entry on oeis.org

1, 6, 30, 150, 750, 3750, 18750, 93750, 468750, 2343750, 11718750, 58593750, 292968750, 1464843750, 7324218750, 36621093750, 183105468750, 915527343750, 4577636718750, 22888183593750, 114440917968750, 572204589843750, 2861022949218750, 14305114746093750
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 6.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954, m is 2, 3, 4, 5, 6. - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001
Hamiltonian path in S_4 X P_2n.
For n>=1, a(n+1) is equal to the number of functions f:{1,2,...,n+1}->{1,2,3,4,5,6} such that for fixed, different x_1, x_2,...,x_n in {1,2,...,n+1} and fixed y_1, y_2,...,y_n in {1,2,3,4,5,6} we have f(x_i)<>y_i, (i=1..n). - Milan Janjic, May 10 2007
For n>=1, a(n) equals the numbers of words of length n over the alphabet {0..5} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015 [Corrected by David Nacin, May 30 2017]
a(n) equals the numbers of sequences of length n on {0,...,5} where no two adjacent terms differ by three. - David Nacin, May 30 2017
It appears that these are the only n>1 for which alpha(n)=2n, where alpha(n) is the entry point of n in the Fibonacci sequence, see A001177. - Philippe Schnoebelen, Apr 11 2024

Crossrefs

Programs

  • GAP
    Concatenation([1], List([1..30], n-> 6*5^(n-1) )); # G. C. Greubel, Sep 24 2019
  • Magma
    [1] cat [6*5^(n-1): n in [1..30]]; // G. C. Greubel, Sep 24 2019
    
  • Maple
    k := 6; if n = 0 then 1 else k*(k-1)^(n-1); fi;
  • Mathematica
    q = 6; Join[{a = 1}, Table[If[n != 0, a = q*a - a, a = q*a], {n, 0, 25}]] (* and *) Join[{1}, 6*5^Range[0, 25]] (* Vladimir Joseph Stephan Orlovsky, Jul 11 2011 *)
    Join[{1},NestList[5#&,6,30]] (* Harvey P. Dale, Dec 31 2013 *)
    CoefficientList[Series[(1+x)/(1-5x), {x,0,30}], x] (* Michael De Vlieger, Dec 10 2016 *)
  • PARI
    Vec((1+x)/(1-5*x)+O(x^30)) \\ Charles R Greathouse IV, Nov 20 2012
    
  • Sage
    [1]+[6*5^(n-1) for n in (1..30)] # G. C. Greubel, Sep 24 2019
    

Formula

G.f.: (1+x)/(1-5*x).
a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 4. - Philippe Deléham, Jul 10 2005
The Hankel transform of this sequence is [1,-6,0,0,0,0,0,0,0,0,...]. - Philippe Deléham, Nov 21 2007
a(n) = 6*5^(n-1) for n>0, a(0)=1. - Vincenzo Librandi, Nov 18 2010
G.f.: 2/x - 5 - 8/(x*U(0)) where U(k)= 1 + 2/(3^k - 3^k/(2 + 1 - 12*x*3^k/(6*x*3^k + 1/U(k+1)))) ; (continued fraction, 4-step). - Sergei N. Gladkovskii, Oct 30 2012
E.g.f.: (6*exp(5*x) - 1)/5. - Ilya Gutkovskiy, Dec 10 2016
Sum_{n>=0} 1/a(n) = 29/24. - Bernard Schott, Oct 25 2021

Extensions

Definition corrected by Frans J. Faase, Feb 07 2009
Edited by N. J. A. Sloane, Dec 04 2009

A003950 Expansion of g.f.: (1+x)/(1-7*x).

Original entry on oeis.org

1, 8, 56, 392, 2744, 19208, 134456, 941192, 6588344, 46118408, 322828856, 2259801992, 15818613944, 110730297608, 775112083256, 5425784582792, 37980492079544, 265863444556808, 1861044111897656, 13027308783283592, 91191161482985144, 638338130380896008
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 8.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954 m is 2, 3, 4, 5, 6 . - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001.
For n>=1, a(n) equals the number of words of length n on the alphabet {0,1,...,7} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015 [Corrected by David Nacin, May 31 2017]
a(n) is the number of octonary sequences of length n such that no two consecutive terms have distance 4. - David Nacin, May 31 2017

Crossrefs

Programs

Formula

a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 6. - Philippe Deléham, Jul 10 2005
From Philippe Deléham, Nov 21 2007: (Start)
a(n) = 8*7^(n-1) for n>=1, a(0)=1 .
G.f.: (1+x)/(1-7x).
The Hankel transform of this sequence is [1,-8,0,0,0,0,0,0,0,0,...]. (End)
a(0)=1, a(1)=8, a(n) = 7*a(n-1). - Vincenzo Librandi, Dec 10 2012
E.g.f.: (8*exp(7*x) - 1)/7. - G. C. Greubel, Sep 24 2019

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A003952 Expansion of g.f.: (1+x)/(1-9*x).

Original entry on oeis.org

1, 10, 90, 810, 7290, 65610, 590490, 5314410, 47829690, 430467210, 3874204890, 34867844010, 313810596090, 2824295364810, 25418658283290, 228767924549610, 2058911320946490, 18530201888518410, 166771816996665690, 1500946352969991210, 13508517176729920890
Offset: 0

Views

Author

N. J. A. Sloane, Mar 15 1996

Keywords

Comments

Coordination sequence for infinite tree with valency 10.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954 m is 2, 3, 4, 5, 6 . - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001
Except 1, all terms are in A033583. - Vincenzo Librandi, May 26 2014
For n>=1, a(n) equals the number of words of length n on alphabet {0,1,...,9} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015 [Corrected by David Nacin, May 31 2017]
a(n) is the number of sequences over the alphabet {0,1,...,9} of length n such that no two consecutive terms have distance 5. - David Nacin, May 31 2017

Crossrefs

Programs

Formula

a(n) = (10*9^n - 0^n)/9. Binomial transform is A000042. - Paul Barry, Jan 29 2004
G.f.: (1+x)/(1-9*x). - Philippe Deléham, Jan 31 2004
a(n) = Sum_{k=0..n} A029653(n, k)*x^k for x = 8. - Philippe Deléham, Jul 10 2005
The Hankel transform of this sequence is: [1,-10,0,0,0,0,0,0,0,...]. - Philippe Deléham, Nov 21 2007
E.g.f.: (10*exp(9*x) - 1)/9. - G. C. Greubel, Sep 24 2019

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A003954 Expansion of g.f.: (1+x)/(1-11*x).

Original entry on oeis.org

1, 12, 132, 1452, 15972, 175692, 1932612, 21258732, 233846052, 2572306572, 28295372292, 311249095212, 3423740047332, 37661140520652, 414272545727172, 4556998002998892, 50126978032987812, 551396758362865932, 6065364341991525252, 66719007761906777772, 733909085380974555492
Offset: 0

Views

Author

Keywords

Comments

Coordination sequence for infinite tree with valency 12.
The n-th term of the coordination sequence of the infinite tree with valency 2m is the same as the number of reduced words of size n in the free group on m generators. In the five sequences A003946, A003948, A003950, A003952, A003954 m is 2, 3, 4, 5, 6 . - Avi Peretz (njk(AT)netvision.net.il), Feb 23 2001 and Ola Veshta (olaveshta(AT)my-deja.com), Mar 30 2001.
For n>=1, a(n) equals the numbers of words of length n-1 on alphabet {0,1,...,11} with no two adjacent letters identical. - Milan Janjic, Jan 31 2015

Crossrefs

Programs

Formula

a(n) = Sum_{k=0..n} A029653(n,k)*x^k for x = 10. - Philippe Deléham, Jul 10 2005
G.f.: (1+x)/(1-11*x). The Hankel transform of this sequence is [1,-12,0,0,0,0,0,0,0,...]. - Philippe Deléham, Nov 21 2007
a(0) = 1; for n>0, a(n) = 12*11^(n-1). - Vincenzo Librandi, Nov 18 2010
a(0) = 1, a(1)=12, a(n) = 11*a(n-1). - Vincenzo Librandi, Dec 10 2012
E.g.f.: (12*exp(11*x) - 1)/11. - Elmo R. Oliveira, Mar 24 2025

Extensions

Edited by N. J. A. Sloane, Dec 04 2009

A036561 Nicomachus triangle read by rows, T(n, k) = 2^(n - k)*3^k, for 0 <= k <= n.

Original entry on oeis.org

1, 2, 3, 4, 6, 9, 8, 12, 18, 27, 16, 24, 36, 54, 81, 32, 48, 72, 108, 162, 243, 64, 96, 144, 216, 324, 486, 729, 128, 192, 288, 432, 648, 972, 1458, 2187, 256, 384, 576, 864, 1296, 1944, 2916, 4374, 6561, 512, 768, 1152, 1728, 2592, 3888, 5832, 8748, 13122, 19683
Offset: 0

Views

Author

Keywords

Comments

The triangle pertaining to this sequence has the property that every row, every column and every diagonal contains a nontrivial geometric progression. More interestingly every line joining any two elements contains a nontrivial geometric progression. - Amarnath Murthy, Jan 02 2002
Kappraff states (pp. 148-149): "I shall refer to this as Nicomachus' table since an identical table of numbers appeared in the Arithmetic of Nicomachus of Gerasa (circa 150 A.D.)" The table was rediscovered during the Italian Renaissance by Leon Battista Alberti, who incorporated the numbers in dimensions of his buildings and in a system of musical proportions. Kappraff states "Therefore a room could exhibit a 4:6 or 6:9 ratio but not 4:9. This ensured that ratios of these lengths would embody musical ratios". - Gary W. Adamson, Aug 18 2003
After Nichomachus and Alberti several Renaissance authors described this table. See for instance Pierre de la Ramée in 1569 (facsimile of a page of his Arithmetic Treatise in Latin in the links section). - Olivier Gérard, Jul 04 2013
The triangle sums, see A180662 for their definitions, link Nicomachus's table with eleven different sequences, see the crossrefs. It is remarkable that these eleven sequences can be described with simple elegant formulas. The mirror of this triangle is A175840. - Johannes W. Meijer, Sep 22 2010
The diagonal sums Sum_{k} T(n - k, k) give A167762(n + 2). - Michael Somos, May 28 2012
Where d(n) is the divisor count function, then d(T(i,j)) = A003991, the rows of which sum to the tetrahedral numbers A000292(n+1). For example, the sum of the divisors of row 4 of this triangle (i = 4), gives d(16) + d(24) + d(36) + d(54) + d(81) = 5 + 8 + 9 + 8 + 5 = 35 = A000292(5). In fact, where p and q are distinct primes, the aforementioned relationship to the divisor function and tetrahedral numbers can be extended to any triangle of numbers in which the i-th row is of form {p^(i-j)*q^j, 0<=j<=i}; i >= 0 (e.g., A003593, A003595). - Raphie Frank, Nov 18 2012, corrected Dec 07 2012
Sequence (or tree) generated by these rules: 1 is in S, and if x is in S, then 2*x and 3*x are in S, and duplicates are deleted as they occur; see A232559. - Clark Kimberling, Nov 28 2013
Partial sums of rows produce Stirling numbers of the 2nd kind: A000392(n+2) = Sum_{m=1..(n^2+n)/2} a(m). - Fred Daniel Kline, Sep 22 2014
A permutation of A003586. - L. Edson Jeffery, Sep 22 2014
Form a word of length i by choosing a (possibly empty) word on alphabet {0,1} then concatenating a word of length j on alphabet {2,3,4}. T(i,j) is the number of such words. - Geoffrey Critzer, Jun 23 2016
Form of Zorach additive triangle (see A035312) where each number is sum of west and northwest numbers, with the additional condition that each number is GCD of the two numbers immediately below it. - Michel Lagneau, Dec 27 2018

Examples

			The start of the sequence as a triangular array read by rows:
   1
   2   3
   4   6   9
   8  12  18  27
  16  24  36  54  81
  32  48  72 108 162 243
  ...
The start of the sequence as a table T(n,k) n, k > 0:
    1    2    4    8   16   32 ...
    3    6   12   24   48   96 ...
    9   18   36   72  144  288 ...
   27   54  108  216  432  864 ...
   81  162  324  648 1296 2592 ...
  243  486  972 1944 3888 7776 ...
  ...
- _Boris Putievskiy_, Jan 08 2013
		

References

  • Jay Kappraff, Beyond Measure, World Scientific, 2002, p. 148.
  • Flora R. Levin, The Manual of Harmonics of Nicomachus the Pythagorean, Phanes Press, 1994, p. 114.

Crossrefs

Cf. A001047 (row sums), A000400 (central terms), A013620, A007318.
Triangle sums (see the comments): A001047 (Row1); A015441 (Row2); A005061 (Kn1, Kn4); A016133 (Kn2, Kn3); A016153 (Fi1, Fi2); A016140 (Ca1, Ca4); A180844 (Ca2, Ca3); A180845 (Gi1, Gi4); A180846 (Gi2, Gi3); A180847 (Ze1, Ze4); A016185 (Ze2, Ze3). - Johannes W. Meijer, Sep 22 2010, Sep 10 2011
Antidiagonal cumulative sum: A000392; square arrays cumulative sum: A160869. Antidiagonal products: 6^A000217; antidiagonal cumulative products: 6^A000292; square arrays products: 6^A005449; square array cumulative products: 6^A006002.

Programs

  • Haskell
    a036561 n k = a036561_tabf !! n !! k
    a036561_row n = a036561_tabf !! n
    a036561_tabf = iterate (\xs@(x:_) -> x * 2 : map (* 3) xs) [1]
    -- Reinhard Zumkeller, Jun 08 2013
    
  • Magma
    /* As triangle: */ [[(2^(i-j)*3^j)/3: j in [1..i]]: i in [1..10]]; // Vincenzo Librandi, Oct 17 2014
  • Maple
    A036561 := proc(n,k): 2^(n-k)*3^k end:
    seq(seq(A036561(n,k),k=0..n),n=0..9);
    T := proc(n,k) option remember: if k=0 then 2^n elif k>=1 then procname(n,k-1) + procname(n-1,k-1) fi: end: seq(seq(T(n,k),k=0..n),n=0..9);
    # Johannes W. Meijer, Sep 22 2010, Sep 10 2011
  • Mathematica
    Flatten[Table[ 2^(i-j) 3^j, {i, 0, 12}, {j, 0, i} ]] (* Flatten added by Harvey P. Dale, Jun 07 2011 *)
  • PARI
    for(i=0,9,for(j=0,i,print1(3^j<<(i-j)", "))) \\ Charles R Greathouse IV, Dec 22 2011
    
  • PARI
    {T(n, k) = if( k<0 || k>n, 0, 2^(n - k) * 3^k)} /* Michael Somos, May 28 2012 */
    

Formula

T(n,k) = A013620(n,k)/A007318(n,k). - Reinhard Zumkeller, May 14 2006
T(n,k) = T(n,k-1) + T(n-1,k-1) for n>=1 and 1<=k<=n with T(n,0) = 2^n for n>=0. - Johannes W. Meijer, Sep 22 2010
T(n,k) = 2^(k-1)*3^(n-1), n, k > 0 read by antidiagonals. - Boris Putievskiy, Jan 08 2013
a(n) = 2^(A004736(n)-1)*3^(A002260(n)-1), n > 0, or a(n) = 2^(j-1)*3^(i-1) n > 0, where i=n-t*(t+1)/2, j=(t*t+3*t+4)/2-n, t=floor[(-1+sqrt(8*n-7))/2]. - Boris Putievskiy, Jan 08 2013
G.f.: 1/((1-2x)(1-3yx)). - Geoffrey Critzer, Jun 23 2016
T(n,k) = (-1)^n * Sum_{q=0..n} (-1)^q * C(k+3*q, q) * C(n+2*q, n-q). - Marko Riedel, Jul 01 2024

A287804 Number of quinary sequences of length n such that no two consecutive terms have distance 1.

Original entry on oeis.org

1, 5, 17, 59, 205, 713, 2481, 8635, 30057, 104629, 364225, 1267923, 4413861, 15365465, 53490097, 186209299, 648230545, 2256616133, 7855718641, 27347281995, 95201200637, 331413874569, 1153716087665, 4016309864843, 13981555011321, 48672509644725
Offset: 0

Views

Author

David Nacin, Jun 01 2017

Keywords

Examples

			For n=2 the a(2)=17=25-8 sequences contain every combination except these eight: 01,10,12,21,23,32,34,43.
		

Crossrefs

Programs

  • Mathematica
    LinearRecurrence[{5, -5, -1}, {1, 5, 17}, 50]
  • Python
    def a(n):
        if n in [0,1,2]:
            return [1,5,17][n]
        return 5*a(n-1)-5*a(n-2)-a(n-3)

Formula

a(n) = 5*a(n-1) - 5a(n-2) - a(n-3), a(0)=1, a(1)=5, a(2)=17.
G.f.: (1 - 3*x^2)/(1 - 5*x + 5*x^2 + x^3).
Showing 1-10 of 126 results. Next