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 15 results. Next

A081956 Duplicate of A056182.

Original entry on oeis.org

2, 10, 38, 130, 422, 1330, 4118, 12610, 38342, 116050, 350198, 1054690, 3172262, 9533170, 28632278, 85962370, 258018182, 774316690, 2323474358, 6971471650, 20916512102, 62753730610, 188269580438, 564825518530, 1694510110022
Offset: 1

Views

Author

Amarnath Murthy, Apr 02 2003

Keywords

A001047 a(n) = 3^n - 2^n.

Original entry on oeis.org

0, 1, 5, 19, 65, 211, 665, 2059, 6305, 19171, 58025, 175099, 527345, 1586131, 4766585, 14316139, 42981185, 129009091, 387158345, 1161737179, 3485735825, 10458256051, 31376865305, 94134790219, 282412759265, 847255055011, 2541798719465, 7625463267259, 22876524019505
Offset: 0

Views

Author

Keywords

Comments

a(n+1) is the sum of the elements in the n-th row of triangle pertaining to A036561. - Amarnath Murthy, Jan 02 2002
Number of 2 X n binary arrays with a path of adjacent 1's and no path of adjacent 0's from top row to bottom row. - R. H. Hardin, Mar 21 2002
With offset 1, partial sums of A027649. - Paul Barry, Jun 24 2003
Number of distinct lines through the origin in the n-dimensional lattice of side length 2. A049691 has the values for the 2-dimensional lattice of side length n. - Joshua Zucker, Nov 19 2003
a(n+1)/(n+1)=(3*3^n-2*2^n)/(n+1) is the second binomial transform of the harmonic sequence 1/(n+1). - Paul Barry, Apr 19 2005
a(n+1) is the sum of n-th row of A036561. - Reinhard Zumkeller, May 14 2006
The sequence gives the sum of the lengths of the segments in Cantor's dust generating sequence up to the i-th step. Measurement unit = length of the segment of i-th step. - Giorgio Balzarotti, Nov 18 2006
Let T be a binary relation on the power set P(A) of a set A having n = |A| elements such that for every element x, y of P(A), xTy if x is a proper subset of y. Then a(n) = |T|. - Ross La Haye, Dec 22 2006
From Alexander Adamchuk, Jan 04 2007: (Start)
a(n) is prime for n in A057468.
p divides a(p) - 1 for prime p.
Quotients (3^p - 2^p - 1)/p, where p = prime(n), are listed in A127071.
Numbers k such that k divides 3^k - 2^k - 1 are listed in A127072.
Pseudoprimes in A127072(n) include all powers of primes {2,3,7} and some composite numbers that are listed in A127073, which includes all Carmichael numbers A002997.
Numbers n such that n^2 divides 3^n - 2^n - 1 are listed in A127074.
5 divides a(2n).
5^2 divides a(2*5n).
5^3 divides a(2*5^2n).
5^4 divides a(2*5^3n).
7^2 divides a(6*7n).
13 divides a(4n).
13^2 divides a(4*13n).
19 divides a(3n).
19^2 divides a(3*19n).
23^2 divides a(11n).
23^3 divides a(11*23n).
23^4 divides a(11*23^2n).
29 divides a(7n).
p divides a((p-1)n) for prime p>3.
p divides a((p-1)/2) for prime p in A097934. Also primes p such that 6 is a square mod p, except {2,3}, A038876(n).
p^(k+1) divides a(p^k*(p-1)/2*n) for prime p in A097934.
p^(k+1) divides a(p^k*(p-1)*n) for prime p>3.
Note the exception that for p = 23, p^(k+2) divides a(p^k*(p-1)/2*n).
There are no more such exceptions for primes p up to 600000. (End)
a(n) divides a(q*(n+1)-1), for all q integer. Leonardo Sarasua, Apr 15 2024
Final digits of terms follow sequence 1,5,9,5. - Enoch Haga, Nov 26 2007
This is also the second column sequence of the Sheffer triangle A143494 (2-restricted Stirling2 numbers). See the e.g.f. given below. - Wolfdieter Lang, Oct 08 2011
Partial sums give A000392. - Jon Perry, Apr 05 2014
For n >= 1, this is also row 2 of A281890: when consecutive positive integers are written as a product of primes in nondecreasing order, "3" occurs in n-th position a(n) times out of every 6^n. - Peter Munn, May 17 2017
a(n) is the number of ternary sequences of length n which include the digit 2. For example, a(2)=5 since the sequences are 02,20,12,21,22. - Enrique Navarrete, Apr 05 2021
a(n-1) is the number of ways we can form disjoint unions of two nonempty subsets of [n] such that the union contains n. For example, for n = 3, a(2) = 5 since the disjoint unions are {1}U{3}, {1}U{2,3}, {2}U{3}, {2}U{1,3}, and {1,2}U{3}. Cf. A000392 if we drop the requirement that the union contains n. - Enrique Navarrete, Aug 24 2021
Configures as a composite Koch Snowflake Fractal (see illustration in links) based on the five-fold division of the Cantor Square/Cantor Dust Fractal of (9^n-4^n)/5 see my illustration in (A016153). - John Elias, Oct 13 2021
Number of pairs (A,B) where B is a subset of {1,2,...,n} and A is a proper subset of B. - Jianing Song, Jun 18 2022
From Manfred Boergens, Mar 29 2023: (Start)
With regard to the comments by Ross La Haye and Jianing Song: Omitting "proper" gives A000244.
Number of pairs (A,B) where B is a nonempty subset of {1,2,...,n} and A is a nonempty subset of B. For nonempty proper subsets see a(n+1) in A028243. (End)
a(n) is the number of n-digit numbers whose smallest decimal digit is 7. - Stefano Spezia, Nov 15 2023
a(n-1) is the number of all possible player-reduced binary games observed by each player in an nx2 game assuming the individual strategies of k < n - 1 players are fixed and the remaining n - k - 1 player will play as one, either maintaining their status quo strategies or jointly adopting an alternative strategy. - Ambrosio Valencia-Romero, Apr 11 2024

References

  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 86-87.
  • 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

a(n) = row sums of A091913, row 2 of A047969, column 1 of A090888 and column 1 of A038719.
Cf. partitions: A241766, A241759.
A diagonal of A262307.

Programs

  • Haskell
    a001047 n = a001047_list !! n
    a001047_list = map fst $ iterate (\(u, v) -> (3 * u + v, 2 * v)) (0, 1)
    -- Reinhard Zumkeller, Jun 09 2013
  • Magma
    [3^n - 2^n: n in [0..30]]; // Vincenzo Librandi, Jul 17 2011
    
  • Maple
    seq(3^n - 2^n, n=0..40); # Giorgio Balzarotti, Nov 18 2006
    A001047:=1/(3*z-1)/(2*z-1); # Simon Plouffe in his 1992 dissertation, dropping the initial zero
  • Mathematica
    Table[ 3^n - 2^n, {n, 0, 25} ]
    LinearRecurrence[{5, -6}, {0, 1}, 25] (* Harvey P. Dale, Aug 18 2011 *)
    Numerator@NestList[(3#+1)/2&,1/2,100] (* Zak Seidov, Oct 03 2011 *)
  • PARI
    {a(n) = 3^n - 2^n};
    
  • Python
    [3**n - 2**n for n in range(25)] # Ross La Haye, Aug 19 2005; corrected by David Radcliffe, Jun 26 2016
    
  • Sage
    [lucas_number1(n, 5, 6) for n in range(26)]  # Zerinvary Lajos, Apr 22 2009
    

Formula

G.f.: x/((1-2*x)*(1-3*x)).
a(n) = 5*a(n-1) - 6*a(n-2).
a(n) = 3*a(n-1) + 2^(n-1). - Jon Perry, Aug 23 2002
Starting 0, 0, 1, 5, 19, ... this is 3^n/3 - 2^n/2 + 0^n/6, the binomial transform of A086218. - Paul Barry, Aug 18 2003
a(n) = A083323(n)-1 = A056182(n)/2 = (A002783(n)-1)/2 = (A003063(n+2)-A003063(n+1))/2. - Ralf Stephan, Jan 12 2004
Binomial transform of A000225. - Ross La Haye, Feb 07 2005
a(n) = Sum_{k=0..n-1} binomial(n, k)*2^k. - Ross La Haye, Aug 20 2005
a(n) = 2^(2n) - A083324(n). - Ross La Haye, Sep 10 2005
a(n) = A112626(n, 1). - Ross La Haye, Jan 11 2006
E.g.f.: exp(3*x) - exp(2*x). - Mohammad K. Azarian, Jan 14 2009
a(n) = A217764(n,1). - Ross La Haye, Mar 27 2013
a(n) = 2*a(n-1) + 3^(n-1). - Toby Gottfried, Mar 28 2013
a(n) = A000244(n) - A000079(n). - Omar E. Pol, Mar 28 2013
a(n) = Sum_{k=0..2} Stirling1(2,k)*(k+1)^n = c_2^{(-n)}, poly-Cauchy numbers. - Takao Komatsu, Mar 28 2013
a(n) = A227048(n,A098294(n)). - Reinhard Zumkeller, Jun 30 2013
a(n+1) = Sum_{k=0..n} 2^k*3^(n-k). - J. M. Bergot, Mar 27 2018
Sum_{n>=1} 1/a(n) = A329064. - Amiram Eldar, Nov 20 2020
a(n) = (1/2)*Sum_{k=0..n} binomial(n, k)*(2^(n-k) + 2^k - 2).
a(n) = A001117(n) + 2*A000918(n) + 1. - Ambrosio Valencia-Romero, Mar 08 2022
a(n) = A000225(n) + A028243(n+1). - Ambrosio Valencia-Romero, Mar 09 2022
From Peter Bala, Jun 27 2025: (Start)
exp(Sum_{n >=1} a(2*n)/a(n)*x^n/n) = Sum_{n >= 0} a(n+1)*x^n.
exp(Sum_{n >=1} a(3*n)/a(n)*x^n/n) = 1 + 19*x + 247*x^2 + ... is the g.f. of A019443.
exp(Sum_{n >=1} a(4*n)/a(n)*x^n/n) = 1 + 65*x + 2743*x^2 + ... is the g.f. of A383754.
The following are all examples of telescoping series:
Sum_{n >= 1} 6^n/(a(n)*a(n+1)) = 2, since 6^n/(a(n)*a(n+1)) = b(n) - b(n+1), where b(n) = 2^n/a(n);
Sum_{n >= 1} 18^n/(a(n)*a(n+1)*a(n+2)) = 22/75, since 18^n/(a(n)*a(n+1)*a(n+2)) = c(n) - c(n+1), where c(n) = (5*6^n - 2*4^n)/(15*a(n)*a(n+1));
Sum_{n >= 1} 54^n/(a(n)*a(n+1)*a(n+2)*a(n+3)) = 634/48735 since 54^n/(a(n)*a(n+1)*a(n+2)*a(n+3)) = d(n) - d(n+1), where d(n) = (57*18^n - 38*12^n + 8*8^n)/(513*a(n)*a(n+1)*a(n+2)).
Sum_{n >= 1} 6^n/(a(n)*a(n+2)) = 14/25; Sum_{n >= 1} (-6)^n/(a(n)*a(n+2)) = -6/25.
Sum_{n >= 1} 6^n/(a(n)*a(n+3)) = 306/1805.
Sum_{n >= 1} 6^n/(a(n)*a(n+4)) = 4282/80275; Sum_{n >= 1} (-6)^n/(a(n)*a(n+4)) = -1698/80275. (End)

Extensions

Edited by Charles R Greathouse IV, Mar 24 2010

A027649 a(n) = 2*(3^n) - 2^n.

Original entry on oeis.org

1, 4, 14, 46, 146, 454, 1394, 4246, 12866, 38854, 117074, 352246, 1058786, 3180454, 9549554, 28665046, 86027906, 258149254, 774578834, 2323998646, 6972520226, 20918609254, 62757924914, 188277969046, 564842295746, 1694543664454, 5083664547794, 15251060752246
Offset: 0

Views

Author

Keywords

Comments

Poly-Bernoulli numbers B_n^(k) with k=-2.
Binomial transform of A007051, if both sequences start at 0. Binomial transform of A000225(n+1). - Paul Barry, Mar 24 2003
Euler expands (1-z)/(1-5z+6z^2) and finds the general term. Section 226 of the Introductio indicates that he could have written down the recursion relation: a(n) = 5 a(n-1)-6 a(n-2). - V. Frederick Rickey (fred-rickey(AT)usma.edu), Feb 10 2006
Let R be a binary relation on the power set P(A) of a set A having n = |A| elements such that for every element x, y of P(A), xRy if x is a subset of y or y is a subset of x. Then a(n) = |R|. - Ross La Haye, Dec 22 2006
With regard to the comment by Ross La Haye: For proper subsets see A056182. - For nonempty subsets see A091344. - For nonempty proper subsets see a(n+1) in A260217. - Manfred Boergens, Aug 02 2023
If x, y are two n-bit binary strings then a(n) gives the number of pairs (x,y) such that XOR(x, y) = ABS(x - y). - Ramasamy Chandramouli, Feb 15 2009
Equals row sums of the triangular version of A038573. - Gary W. Adamson, Jun 04 2009
Inverse binomial transform of A085350. - Paul Curtz, Nov 14 2009
Related to the number of even a's in a nontrivial cycle (should one exist) in the 3x+1 Problem, where a <= floor(log_2(2*(3^n) - 2^n)). The value n correlates to the number of odds in such a nontrivial cycle. See page 1288 of Crandall's paper. Also, this relation gives another proof that the number of odds divided by the number of evens in a nontrivial cycle is bounded by log 2 / log 3 (this observation does not resolve the finite cycles conjecture as the value could be arbitrarily close to this bound). However, the same argument gives that log 2 / log 3 is less than or equal to the number of odds divided by the number of evens in a divergent sequence (should one exist), as log 2 / log 3 is the limit value for a cycle of an arbitrarily large length, where the length is given by the value n. - Jeffrey R. Goodwin, Aug 04 2011
Row sums of Riordan triangle A106516. - Wolfdieter Lang, Jan 09 2015
Number of restricted barred preferential arrangements having 3 bars in which the sections are all restricted sections such that (for fixed sections i and j) section i or section j is empty. - Sithembele Nkonkobe, Oct 12 2015
This is also row 2 of A281891: for n >= 1, when consecutive positive integers are written as a product of primes in nondecreasing order, a factor of 2 or 3 occurs in n-th position a(n) times out of every 6^n. - Peter Munn, May 18 2017
Also row sums of A124929. - Omar E. Pol, Jun 15 2017
This is the sum of A318921(n) for n in the range 2^(k+1) to 2^(k+2)-1. See A318921 for proof. - N. J. A. Sloane, Sep 25 2018
a(n) is also the number of acyclic orientations of the complete bipartite graph K_{2,n}. - Vincent Pilaud, Sep 15 2020
a(n-1) is also the number of n-digit numbers whose largest decimal digit is 2. - Stefano Spezia, Nov 15 2023

References

  • Leonhard Euler, Introductio in analysin infinitorum (1748), section 216.

Crossrefs

Row n = 2 of array A099594.
Also occurs as a row, column, diagonal or as row sums in A038573, A085870, A090888, A106516, A217764, A281891.

Programs

  • Haskell
    a027649 n = a027649_list !! n
    a027649_list = map fst $ iterate (\(u, v) -> (3 * u + v, 2 * v)) (1, 1)
    -- Reinhard Zumkeller, Jun 09 2013
    
  • Magma
    [2*(3^n)-2^n: n in [0..30]]; // Vincenzo Librandi, Jul 17 2011
    
  • Maple
    a(n, k):= (-1)^n*sum( (-1)^'m'*'m'!*Stirling2(n,'m')/('m'+1)^k,'m'=0..n);
    seq(a(n, -2), n=0..30);
  • Mathematica
    Table[2(3^n)-2^n,{n,0,30}] (* or *) LinearRecurrence[ {5,-6},{1,4},31]  (* Harvey P. Dale, Apr 22 2011 *)
  • PARI
    a(n)=2*(3^n)-2^n \\ Charles R Greathouse IV, Jul 16 2011
    
  • PARI
    Vec((1-x)/((1-2*x)*(1-3*x)) + O(x^50)) \\ Altug Alkan, Oct 12 2015
    
  • SageMath
    [2*(3^n - 2^(n-1)) for n in (0..30)] # G. C. Greubel, Aug 01 2022

Formula

G.f.: (1-x)/((1-2*x)*(1-3*x)).
a(n) = 3*a(n-1) + 2^(n-1), with a(0) = 1.
a(n) = Sum_{k=0..n} binomial(n, k)*(2^(k+1) - 1). - Paul Barry, Mar 24 2003
Partial sums of A053581. - Paul Barry, Jun 26 2003
Main diagonal of array (A085870) defined by T(i, 1) = 2^i - 1, T(1, j) = 2^j - 1, T(i, j) = T(i-1, j) + T(i-1, j-1). - Benoit Cloitre, Aug 05 2003
a(n) = A090888(n, 3). - Ross La Haye, Sep 21 2004
a(n) = Sum_{k=0..n} binomial(n+2, k+1)*Sum_{j=0..floor(k/2)} A001045(k-2j). - Paul Barry, Apr 17 2005
a(n) = Sum_{k=0..n} Sum_{j=0..n} binomial(n,j)*binomial(j+1,k+1). - Paul Barry, Sep 18 2006
a(n) = A166060(n+1)/6. - Philippe Deléham, Oct 21 2009
a(n) = 5*a(n-1) - 6*a(n-2), a(0)=1, a(1)=4. - Harvey P. Dale, Apr 22 2011
a(n) = A217764(n,2). - Ross La Haye, Mar 27 2013
For n>0, a(n) = 3 * a(n-1) + 2^(n-1) = 2 * (a(n-1) + 3^(n-1)). - J. Conrad, Oct 29 2015
for n>0, a(n) = 2 * (1 + 2^(n-2) + Sum_{x=1..n-2} Sum_{k=0..x-1} (binomial(x-1,k)*(2^(k+1) + 2^(n-x+k)))). - J. Conrad, Dec 10 2015
E.g.f.: exp(2*x)*(2*exp(x) - 1). - Stefano Spezia, May 18 2024

Extensions

Better formulas from David W. Wilson and Michael Somos
Incorrect formula removed by Charles R Greathouse IV, Mar 18 2010
Duplications (due to corrections to A numbers) removed by Peter Munn, Jun 15 2017

A028243 a(n) = 3^(n-1) - 2^n + 1 (essentially Stirling numbers of second kind).

Original entry on oeis.org

0, 0, 2, 12, 50, 180, 602, 1932, 6050, 18660, 57002, 173052, 523250, 1577940, 4750202, 14283372, 42915650, 128878020, 386896202, 1161212892, 3484687250, 10456158900, 31372671002, 94126401612, 282395982050, 847221500580, 2541731610602, 7625329049532
Offset: 1

Views

Author

N. J. A. Sloane, Doug McKenzie (mckfam4(AT)aol.com)

Keywords

Comments

For n >= 3, a(n) is equal to the number of functions f: {1,2,...,n-1} -> {1,2,3} such that Im(f) contains 2 fixed elements. - Aleksandar M. Janjic and Milan Janjic, Mar 08 2007
Let P(A) be the power set of an n-element set A. Then a(n+1) = the number of pairs of elements {x,y} of P(A) for which x and y are intersecting and for which either x is a proper subset of y or y is a proper subset of x. - Ross La Haye, Jan 02 2008
Let P(A) be the power set of an n-element set A and R be a relation on P(A) such that for all x, y of P(A), xRy if x is not a subset of y and y is not a subset of x and x and y are disjoint. Then a(n+1) = |R|. - Ross La Haye, Mar 19 2009
Let P(A) be the power set of an n-element set A and R be a relation on P(A) such that for all x, y of P(A), xRy if either 0) x is a proper subset of y or y is a proper subset of x, or 1) x is not a subset of y and y is not a subset of x and x and y are disjoint. Then a(n+2) = |R|. - Ross La Haye, Mar 19 2009
In the terdragon curve, a(n) is the number of triple-visited points in expansion level n. The first differences of this sequence (A056182) are the number of enclosed unit triangles since on segment expansion each unit triangle forms a new triple-visited point, and existing triple-visited points are unchanged. - Kevin Ryde, Oct 20 2020
a(n+1) is the number of ternary strings of length n that contain at least one 0 and one 1; for example, for n=3, a(4)=12 since the strings are the 3 permutations of 100, the 3 permutations of 110, and the 6 permutations of 210. - Enrique Navarrete, Aug 13 2021
From Sanjay Ramassamy, Dec 23 2021: (Start)
a(n+1) is the number of topological configurations of n points and n lines where the points lie at the vertices of a convex cyclic n-gon and the lines are the perpendicular bisectors of its sides.
a(n+1) is the number of 2n-tuples composed of n 0's and n 1's which have an interlacing signature. The signature of a 2n-tuple (v_1,...,v_{2n}) is the n-tuple (s_1,...,s_n) defined by s_i=v_i+v_{i+n}. The signature is called interlacing if after deleting the 1's, there are letters remaining and the remaining 0's and 2's are alternating. (End)
a(n+1) is the number of pairs (A,B) where B is a nonempty subset of {1,2,...,n} and A is a nonempty proper subset of B. If either "nonempty" or "proper" is omitted then see A001047. If "nonempty" and "proper" are omitted then see A000244. - Manfred Boergens, Mar 28 2023
a(n) is the number of (n-1) X (n-1) nilpotent Boolean relation matrices with rank equal to 1. a(n) = A060867(n-1) - A005061(n-1) (since every rank 1 matrix is either idempotent or nilpotent). - Geoffrey Critzer, Jul 13 2023
For odd n > 3, a(n) is also the number of minimum vertex colorings in the (n-1)-prism graph. - Eric W. Weisstein, Mar 05 2024

Crossrefs

Cf. A000392, A008277, A163626, A056182 (first differences), A000244, A001047.

Programs

  • Magma
    [3^(n-1) - 2*2^(n-1) + 1: n in [1..30]]; // G. C. Greubel, Nov 19 2017
    
  • Mathematica
    Table[2 StirlingS2[n, 3], {n, 24}] (* or *)
    Table[3^(n - 1) - 2*2^(n - 1) + 1, {n, 24}] (* or *)
    Rest@ CoefficientList[Series[-2 x^3/(-1 + x)/(-1 + 3 x)/(-1 + 2 x), {x, 0, 24}], x] (* Michael De Vlieger, Sep 24 2016 *)
  • PARI
    a(n) = 3^(n-1) - 2*2^(n-1) + 1 \\ G. C. Greubel, Nov 19 2017
  • Sage
    [stirling_number2(i,3)*2 for i in range(1,30)] # Zerinvary Lajos, Jun 26 2008
    

Formula

a(n) = 2*S(n, 3) = 2*A000392(n). - Emeric Deutsch, May 02 2004
G.f.: -2*x^3/(-1+x)/(-1+3*x)/(-1+2*x) = -1/3 - (1/3)/(-1+3*x) + 1/(-1+2*x) - 1/(-1+x). - R. J. Mathar, Nov 22 2007
E.g.f.: (exp(3*x) - 3*exp(2*x) + 3*exp(x) - 1)/3. - Wolfdieter Lang, May 03 2017
E.g.f. with offset 0: exp(x)*(exp(x)-1)^2. - Enrique Navarrete, Aug 13 2021
a(n) = Sum_{k = 1..n-2} binomial(n-1, k) * (2^(n-k-1)-1). - Ocean Wong, Jan 03 2025

A003063 a(n) = 3^(n-1) - 2^n.

Original entry on oeis.org

-1, -1, 1, 11, 49, 179, 601, 1931, 6049, 18659, 57001, 173051, 523249, 1577939, 4750201, 14283371, 42915649, 128878019, 386896201, 1161212891, 3484687249, 10456158899, 31372671001, 94126401611, 282395982049, 847221500579, 2541731610601, 7625329049531, 22876255584049
Offset: 1

Views

Author

Henrik Johansson (Henrik.Johansson(AT)Nexus.SE)

Keywords

Comments

Binomial transform of A000918: (-1, 0, 2, 6, 14, 30, ...). - Gary W. Adamson, Mar 23 2012
This sequence demonstrates 2^n as a loose lower bound for g(n) in Waring's problem. Since 3^n > 2(2^n) for all n > 2, the number 2^(n + 1) - 1 requires 2^n n-th powers for its representation since 3^n is not available for use in the sum: the gulf between the relevant powers of 2 and 3 widens considerably as n gets progressively larger. - Alonso del Arte, Feb 01 2013

Examples

			a(3) = 1 because 3^2 - 2^3 = 9 - 8 = 1.
a(4) = 11 because 3^3 - 2^4 = 27 - 16 = 11.
a(5) = 49 because 3^4 - 2^5 = 81 - 32 = 49.
		

Crossrefs

Cf. A000918, A056182 (first differences), A064686, A083313, A214091, A369490.
Cf. A363024 (prime terms).
From the third term onward the first differences of A005173.
Difference between two leftmost columns of A090888.
A diagonal in A254027.
Right edge of irregular triangle A252750.

Programs

Formula

Let b(n) = 2*(3/2)^n - 1. Then a(n) = -b(1-n)*3^(n-1) for n > 0. A083313(n) = A064686(n) = b(n)*2^(n-1) for n > 0. - Michael Somos, Aug 06 2006
From Colin Barker, May 27 2013: (Start)
a(n) = 5*a(n-1) - 6*a(n-2).
G.f.: -x*(1-4*x) / ((1-2*x)*(1-3*x)). (End)
E.g.f.: (1/3)*(2 - 3*exp(2*x) + exp(3*x)). - G. C. Greubel, Nov 03 2022

Extensions

A few more terms from Alonso del Arte, Feb 01 2013

A293181 Irregular triangle read by rows: T(n,k) is the number of k-partitions of {1..2n} that are invariant under a permutation consisting of n 2-cycles (1 <= k <= 2n).

Original entry on oeis.org

1, 1, 1, 3, 2, 1, 1, 7, 10, 9, 3, 1, 1, 15, 38, 53, 34, 18, 4, 1, 1, 31, 130, 265, 261, 195, 80, 30, 5, 1, 1, 63, 422, 1221, 1700, 1696, 1016, 515, 155, 45, 6, 1, 1, 127, 1330, 5369, 10143, 13097, 10508, 6832, 2926, 1120, 266, 63, 7, 1
Offset: 1

Views

Author

Andrew Howroyd, Oct 01 2017

Keywords

Comments

See A002872 for detailed description.
T(m,k) is the number of achiral color patterns in a row or loop of length 2m using exactly k different colors. Two color patterns are equivalent if we permute the colors. - Robert A. Russell, Apr 24 2018
T(n,k) = coefficient of x^k for A(2,n)(x) in Gilbert and Riordan's article. - Robert A. Russell, Jun 14 2018

Examples

			Triangle begins:
  1,   1;
  1,   3,    2,    1;
  1,   7,   10,    9,     3,     1;
  1,  15,   38,   53,    34,    18,     4,    1;
  1,  31,  130,  265,   261,   195,    80,   30,    5,    1;
  1,  63,  422, 1221,  1700,  1696,  1016,  515,  155,   45,   6,  1;
  1, 127, 1330, 5369, 10143, 13097, 10508, 6832, 2926, 1120, 266, 63, 7, 1;
  ...
For T(2,2)=3, the row patterns are AABB, ABAB, and ABBA.  The loop patterns are AAAB, AABB, and ABAB. - _Robert A. Russell_, Apr 24 2018
		

Crossrefs

Row sums are A002872.
Maximum row values are A002873.
Number of achiral color patterns of length odd n in A140735.
Column k=3 gives A056182.

Programs

  • Mathematica
    (* Ach[n, k] is the number of achiral color patterns for a row or loop of n
      colors containing k different colors *)
    Ach[n_, k_] := Ach[n, k] = Which[0==k, Boole[0==n], 1==k, Boole[n>0],
      OddQ[n], Sum[Binomial[(n-1)/2, i] Ach[n-1-2i, k-1], {i, 0, (n-1)/2}],
      True, Sum[Binomial[n/2-1, i] (Ach[n-2-2i, k-1]
      + 2^i Ach[n-2-2i, k-2]), {i, 0, n/2-1}]]
    Table[Ach[n, k], {n, 2, 14, 2}, {k, 1, n}] // Flatten
    (* Robert A. Russell, Feb 06 2018 *)
    Table[Drop[MatrixPower[Table[Switch[j-i, 0, i-1, 1, 1, 2, 1, _, 0],
      {i, 1, 2n+1}, {j, 1, 2n+1}], n][[1]], 1], {n, 1, 10}] // Flatten
    (* Robert A. Russell, Apr 14 2018 *)
    Aeven[m_, k_] := Aeven[m, k] = If[m>0, k Aeven[m-1, k] + Aeven[m-1, k-1]
      + Aeven[m-1, k-2], Boole[m == 0 && k == 0]]
    Table[Aeven[m, k], {m, 1, 10}, {k, 1, 2m}] // Flatten (* Robert A. Russell, Apr 24 2018 *)
  • PARI
    \\ see A056391 for Polya enumeration functions
    T(n,k) = 2*NonequivalentStructsExactly(CylinderPerms(2,n),k) - stirling(2*n,k,2);
    
  • PARI
    seq(n)={Vec(serlaplace(exp(y*(exp(x + O(x*x^n))-1)+(1/2)*y^2*(exp(2*x + O(x*x^n))-1))) - 1)}
    {my(T=seq(10)); for(n=1, #T, for(k=1, 2*n, print1(polcoeff(T[n], k), ", ")); print)} \\ Andrew Howroyd, Jan 31 2018

Formula

T(n,k) = coefficient of t^k x^n/n! in exp(t*(exp(x)-1)+(1/2)*t^2*(exp(2*x)-1)). - Ira M. Gessel, Jan 30 2018
T(m,k) = [m>0]*(k*T(m-1,k)+T(m-1,k-1)+T(m-1,k-2)) + [m==0]*[k==0]. - Robert A. Russell, Apr 24 2018
Conjecture: T(n,k) = R(n,k)-R(n,k-1), with R(n,k) = Sum_{m=0..k} m^n*A000085(m)*A038205(k-m)/(m!*(k-m)!). - Mikhail Kurkov, Jun 26 2018

A091344 a(n) = 2*3^n - 3*2^n + 1.

Original entry on oeis.org

0, 1, 7, 31, 115, 391, 1267, 3991, 12355, 37831, 115027, 348151, 1050595, 3164071, 9516787, 28599511, 85896835, 257887111, 774054547, 2322950071, 6970423075, 20914414951, 62749536307, 188261191831, 564808741315, 1694476555591
Offset: 0

Views

Author

Mario Catalani (mario.catalani(AT)unito.it), Jan 01 2004

Keywords

Comments

Starting with offset 1 = binomial transform of A068293: (1, 6, 18, 42, 90, ...) and double binomial transform of (1, 5, 7, 5, 7, 5, ...). - Gary W. Adamson, Jan 13 2009
Number of pairs (A,B) where A and B are nonempty subsets of {1,2,...,n} and one of these subsets is a subset of the other. - For the case that one of these subsets is a proper subset of the other see a(n+1) in A260217. - If empty subsets are included, see A027649 (all subsets) and A056182 (proper subsets). - Manfred Boergens, Aug 02 2023

Crossrefs

Programs

  • Maple
    a:=n->sum((3^(n-j-1)-2^(n-2-j))*12, j=0..n): seq(a(n), n=-1..24); # Zerinvary Lajos, Feb 11 2007
    with (combinat):a:=n->stirling2(n,3)+stirling2(n+1,3): seq(a(n), n=1..26); # Zerinvary Lajos, Oct 07 2007
  • Mathematica
    Table[Sum[i!i^2 StirlingS2[n, i](-1)^(n - i), {i, 1, n}], {n, 0, 30}]
    Table[2*3^n-3*2^n+1,{n,0,30}] (* or *) LinearRecurrence[{6,-11,6},{0,1,7},30] (* Harvey P. Dale, Dec 31 2013 *)

Formula

a(n) = Sum_{i=1..n} i!*i^2*Stirling2(n,i)*(-1)^(n-i).
From Christian Ballot via R. K. Guy, Jan 13 2009: (Start)
a(n) = 6*a(n-1) - 11*a(n-2) + 6*a(n-3);
G.f.: x*(1+x)/((1-x)*(2-x)*(3-x)). (End)
a(n) = 5*a(n-1) - 6*a(n-2) + 2, a(0)=0, a(1)=1. - Vincenzo Librandi, Nov 25 2010
E.g.f.: exp(x)*(1 - 3*exp(x) + 2*exp(2*x)). - Stefano Spezia, May 18 2024

Extensions

Edited by N. J. A. Sloane, Jan 13 2009 at the suggestion of R. K. Guy; the concise definition was provided by Vladeta Jovovic, Jan 01 2004

A217764 Array defined by a(n,k) = floor((k+2)/2)*3^n - floor((k+1)/2)*2^n, read by antidiagonals.

Original entry on oeis.org

1, 3, 0, 9, 1, 1, 27, 5, 4, 0, 81, 19, 14, 2, 1, 243, 65, 46, 10, 5, 0, 729, 211, 146, 38, 19, 3, 1, 2187, 665, 454, 130, 65, 15, 6, 0, 6561, 2059, 1394, 422, 211, 57, 24, 4, 1, 19683, 6305, 4246, 1330, 665, 195, 84, 20, 7, 0, 59049, 19171, 12866, 4118, 2059, 633, 276, 76, 29, 5, 1
Offset: 0

Views

Author

Ross La Haye, Mar 23 2013

Keywords

Comments

Columns 0,1,2,3 respectively correspond to relations R_3, R_4, R_0, R_1 defined in La Haye paper listed below.

Examples

			a(4,4) = 211 because floor((4+2)/2)*3^4 - floor((4+1)/2)*2^4 = 3*3^4 - 2*2^4 = 243 - 32 = 211.
		

Crossrefs

Cf. a(1,k) = A084964(k+2); a(n,0) = A000244(n); a(n,1) = A001047(n); a(n,2) = A027649(n); a(n,3) = A056182(n); a(n,4) = A001047(n+1); a(n,5) = A210448(n); a(n,6) = A166060(n); a(n,7) = A145563(n); a(n,8) = A102485(n).

Formula

a(n,k) = floor((k+2)/2)*3^n - floor((k+1)/2)*2^n. a(n,k) = 5*a(n-1,k) - 6*a(n-2,k); a(0,k) = floor((k+2)/2) - floor((k+1)/2), a(1,k) = floor((k+2)/2)*3 - floor((k+1)/2)*2.

A056151 Distribution of maximum inversion table entry.

Original entry on oeis.org

1, 1, 1, 1, 3, 2, 1, 7, 10, 6, 1, 15, 38, 42, 24, 1, 31, 130, 222, 216, 120, 1, 63, 422, 1050, 1464, 1320, 720, 1, 127, 1330, 4686, 8856, 10920, 9360, 5040, 1, 255, 4118, 20202, 50424, 80520, 91440, 75600, 40320, 1, 511, 12610, 85182, 276696, 558120, 795600, 851760, 685440, 362880
Offset: 1

Views

Author

Wouter Meeussen, Aug 05 2000

Keywords

Comments

T(n,k) is the number of permutations p of [n] such that max(p(i) - i) = k. Example: T(3,0) = 1 because for p = 123 we have max(p(i) - i) = 0; T(3,1) = 3 because for p = 132, 213, 231 we have max(p(i) - i) = 1; T(3,2) = 2 because for p = 312, 321 we have max(p(i) - i) = 2. - Emeric Deutsch, Nov 12 2004
T(n,k) is the number of permutations of [n] for which the first subcedance occurs at position n + 2 - k. A subcedance of pi occurs at position i if i>pi(i). For example, with n = 3 and k = 2, T(n,k) = 3 counts 132, 231, 321 in each of which the first subcedance occurs at position n+2-k = 3. - David Callan, Dec 14 2021

Examples

			Triangle starts:
  1;
  1,   1;
  1,   3,    2;
  1,   7,   10,     6;
  1,  15,   38,    42,    24;
  1,  31,  130,   222,   216,   120;
  1,  63,  422,  1050,  1464,  1320,   720;
  1, 127, 1330,  4686,  8856, 10920,  9360,  5040;
  1, 255, 4118, 20202, 50424, 80520, 91440, 75600, 40320;
		

References

  • R. Sedgewick and Ph. Flajolet, "An Introduction to the Analysis of Algorithms", Addison-Wesley, 1996, ISBN 0-201-40009-X, table 6.10 (page 356)

Crossrefs

Columns and diagonals give A000225, A018927, A056182, A000142, A056197.
Cf. A343237.

Programs

  • Maple
    T:=proc(n,k) if k>0 and k<=n then k!*(k+1)^(n-k)-(k-1)!*k^(n-k+1) elif k=0 then 1 else 0 fi end: TT:=(n,k)->T(n,k-1): matrix(10,10,TT);
    # Alternative, assuming offset 0:
    egf := exp(exp(x)*y + x)*(exp(x)*y - y + 1): ser := series(egf, x, 12):
    cx := n -> series(coeff(ser, x, n), y, 12):
    T := (n, k) -> k!^2 * (n-k)! * coeff(cx(n - k), y, k):
    for n from 0 to 6 do seq(T(n, k), k=0..n) od; # Peter Luschny, Dec 14 2021
  • Mathematica
    T[, 0] = 1; T[n, k_] := k! (k + 1)^(n - k) - (k - 1)! k^(n - k + 1);
    Table[T[n, k], {n, 1, 10}, {k, 0, n - 1}] // Flatten (* Jean-François Alcover, May 03 2017 *)

Formula

Table[ -((-1 + k)^(1 - k + n)*(-1 + k)!) + k^(-k + n)*k!, {n, 1, 9}, {k, 1, n} ]
T(n, k) = k!(k+1)^(n-k) - (k-1)!k^(n-k+1) if 0Emeric Deutsch, Nov 12 2004
From Peter Luschny, Dec 14 2021: (Start)
We assume T with offset 0.
T(n, k) = k!^2 * (n-k)! * [y^k] [x^(n-k)] (exp(exp(x)*y + x)*(exp(x)*y - y + 1)).
T(n, k) = k!*A343237(n, k). (End)

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Oct 03 2000

A210204 Triangle of coefficients of polynomials v(n,x) jointly generated with A210203; see the Formula section.

Original entry on oeis.org

1, 3, 2, 7, 8, 2, 15, 24, 12, 2, 31, 64, 48, 16, 2, 63, 160, 160, 80, 20, 2, 127, 384, 480, 320, 120, 24, 2, 255, 896, 1344, 1120, 560, 168, 28, 2, 511, 2048, 3584, 3584, 2240, 896, 224, 32, 2, 1023, 4608, 9216, 10752, 8064, 4032, 1344, 288, 36, 2, 2047
Offset: 1

Views

Author

Clark Kimberling, Mar 18 2012

Keywords

Comments

Column 1: -1+2^n.
Row sums: A048473.
Alternating row sums: 1,1,1,1,1,1,1,1,1,...
For a discussion and guide to related arrays, see A208510.
Row sums without first column give A056182. - Alois P. Heinz, Jan 14 2022

Examples

			First five rows:
1
3....2
7....8....2
15...24...12...2
31...64...48...16...2
First three polynomials v(n,x): 1, 3 + 2x , 7 + 8x + 2x^2.
		

Crossrefs

Programs

  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := u[n - 1, x] + v[n - 1, x] + 1;
    v[n_, x_] := (x + 1)*u[n - 1, x] + (x + 1)*v[n - 1, x] + 1;
    Table[Expand[u[n, x]], {n, 1, z/2}]
    Table[Expand[v[n, x]], {n, 1, z/2}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%]   (* A210203 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]   (* A210204 *)

Formula

u(n,x)=u(n-1,x)+v(n-1,x)+1,
v(n,x)=(x+1)*u(n-1,x)+(x+1)*v(n-1,x)+1,
where u(1,x)=1, v(1,x)=1.
Showing 1-10 of 15 results. Next