cp's OEIS Frontend

For n>=1, a(n) is equal to the number of functions f: {1,2,...,n+1}->{1,2,3,4} such that Im(f) contains 2 fixed elements. - Aleksandar M. Janjic and Milan Janjic, Feb 27 2007

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. Then a(n+1) = |R|. - Ross La Haye, Mar 19 2009

Number of ordered (n+1)-tuples of positive integers, whose minimum is 0 and maximum 3. - Ovidiu Bagdasar, Sep 19 2014

a(n-2) is the number of possible player-reduced binary games observed by each player in an n X 2 game assuming k < n - 1 players reverse their initially fixed individual strategies and the remaining n - k - 1 players will play as one, either maintaining their status quo strategies or jointly adopting an alternative strategy. - Ambrosio Valencia-Romero, Apr 11 2024

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

Let M = n X n matrix with (i,j)-th entry a(n+1-j, n+1-i), e.g., if n = 3, M = [1 1 1; 3 1 0; 2 0 0]. Given a sequence s = [s(0)..s(n-1)], let b = [b(0)..b(n-1)] be its inverse binomial transform and let c = [c(0)..c(n-1)] = M^(-1)*transpose(b). Then s(k) = Sum_{i=0..n-1} b(i)*binomial(k,i) = Sum_{i=0..n-1} c(i)*k^i, k=0..n-1. - Gary W. Adamson, Nov 11 2001

From Gary W. Adamson, Aug 09 2008: (Start)

Julius Worpitzky's 1883 algorithm generates Bernoulli numbers.

By way of example [Wikipedia]:

B0 = 1;

B1 = 1/1 - 1/2;

B2 = 1/1 - 3/2 + 2/3;

B3 = 1/1 - 7/2 + 12/3 - 6/4;

B4 = 1/1 - 15/2 + 50/3 - 60/4 + 24/5;

B5 = 1/1 - 31/2 + 180/3 - 390/4 + 360/5 - 120/6;

B6 = 1/1 - 63/2 + 602/3 - 2100/4 + 3360/5 - 2520/6 + 720/7;

...

Note that in this algorithm, odd n's for the Bernoulli numbers sum to 0, not 1, and the sum for B1 = 1/2 = (1/1 - 1/2). B3 = 0 = (1 - 7/2 + 13/3 - 6/4) = 0. The summation for B4 = -1/30. (End)

Pursuant to Worpitzky's algorithm and given M = A028246 as an infinite lower triangular matrix, M * [1/1, -1/2, 1/3, ...] (i.e., the Harmonic series with alternate signs) = the Bernoulli numbers starting [1/1, 1/2, 1/6, ...]. - Gary W. Adamson, Mar 22 2012

From Tom Copeland, Oct 23 2008: (Start)

G(x,t) = 1/(1 + (1-exp(x*t))/t) = 1 + 1 x + (2 + t)*x^2/2! + (6 + 6t + t^2)*x^3/3! + ... gives row polynomials for A090582, the f-polynomials for the permutohedra (see A019538).

G(x,t-1) = 1 + 1*x + (1 + t)*x^2 / 2! + (1 + 4t + t^2)*x^3 / 3! + ... gives row polynomials for A008292, the h-polynomials for permutohedra.

G[(t+1)x,-1/(t+1)] = 1 + (1+ t) x + (1 + 3t + 2 t^2) x^2 / 2! + ... gives row polynomials for the present triangle. (End)

The Worpitzky triangle seems to be an apt name for this triangle. - Johannes W. Meijer, Jun 18 2009

If Pascal's triangle is written as a lower triangular matrix and multiplied by A028246 written as an upper triangular matrix, the product is a matrix where the (i,j)-th term is (i+1)^j. For example,

1,0,0,0 1,1,1, 1 1,1, 1, 1

1,1,0,0 * 0,1,3, 7 = 1,2, 4, 8

1,2,1,0 0,0,2,12 1,3, 9,27

1,3,3,1 0,0,0, 6 1,4,16,64

So, numbering all three matrices' rows and columns starting at 0, the (i,j) term of the product is (i+1)^j. - Jack A. Cohen (ProfCohen(AT)comcast.net), Aug 03 2010

The Fi1 and Fi2 triangle sums are both given by sequence A000670. For the definition of these triangle sums see A180662. The mirror image of the Worpitzky triangle is A130850. - Johannes W. Meijer, Apr 20 2011

Let S_n(m) = 1^m + 2^m + ... + n^m. Then, for n >= 0, we have the following representation of S_n(m) as a linear combination of the binomial coefficients:

S_n(m) = Sum_{i=1..n+1} a(i+n*(n+1)/2)*C(m,i). E.g., S_2(m) = a(4)*C(m,1) + a(5)*C(m,2) + a(6)*C(m,3) = C(m,1) + 3*C(m,2) + 2*C(m,3). - Vladimir Shevelev, Dec 21 2011

Given the set X = [1..n] and 1 <= k <= n, then a(n,k) is the number of sets T of size k of subset S of X such that S is either empty or else contains 1 and another element of X and such that any two elemements of T are either comparable or disjoint. - Michael Somos, Apr 20 2013

Working with the row and column indexing starting at -1, a(n,k) gives the number of k-dimensional faces in the first barycentric subdivision of the standard n-dimensional simplex (apply Brenti and Welker, Lemma 2.1). For example, the barycentric subdivision of the 2-simplex (a triangle) has 1 empty face, 7 vertices, 12 edges and 6 triangular faces giving row 4 of this triangle as (1,7,12,6). Cf. A053440. - Peter Bala, Jul 14 2014

See A074909 and above g.f.s for associations among this array and the Bernoulli polynomials and their umbral compositional inverses. - Tom Copeland, Nov 14 2014

An e.g.f. G(x,t) = exp[P(.,t)x] = 1/t - 1/[t+(1-t)(1-e^(-xt^2))] = (1-t) * x + (-2t + 3t^2 - t^3) * x^2/2! + (6t^2 - 12t^3 + 7t^4 - t^5) * x^3/3! + ... for the shifted, reverse, signed polynomials with the first element nulled, is generated by the infinitesimal generator g(u,t)d/du = [(1-u*t)(1-(1+u)t)]d/du, i.e., exp[x * g(u,t)d/du] u eval. at u=0 generates the polynomials. See A019538 and the G. Rzadkowski link below for connections to the Bernoulli and Eulerian numbers, a Ricatti differential equation, and a soliton solution to the KdV equation. The inverse in x of this e.g.f. is Ginv(x,t) = (-1/t^2)*log{[1-t(1+x)]/[(1-t)(1-tx)]} = [1/(1-t)]x + [(2t-t^2)/(1-t)^2]x^2/2 + [(3t^2-3t^3+t^4)/(1-t)^3]x^3/3 + [(4t^3-6t^4+4t^5-t^6)/(1-t)^4]x^4/4 + ... . The numerators are signed, shifted A135278 (reversed A074909), and the rational functions are the columns of A074909. Also, dG(x,t)/dx = g(G(x,t),t) (cf. A145271). (Analytic G(x,t) added, and Ginv corrected and expanded on Dec 28 2015.) - Tom Copeland, Nov 21 2014

The operator R = x + (1 + t) + t e^{-D} / [1 + t(1-e^(-D))] = x + (1+t) + t - (t+t^2) D + (t+3t^2+2t^3) D^2/2! - ... contains an e.g.f. of the reverse row polynomials of the present triangle, i.e., A123125 * A007318 (with row and column offset 1 and 1). Umbrally, R^n 1 = q_n(x;t) = (q.(0;t)+x)^n, with q_m(0;t) = (t+1)^(m+1) - t^(m+1), the row polynomials of A074909, and D = d/dx. In other words, R generates the Appell polynomials associated with the base sequence A074909. For example, R 1 = q_1(x;t) = (q.(0;t)+x) = q_1(0;t) + q__0(0;t)x = (1+2t) + x, and R^2 1 = q_2(x;t) = (q.(0;t)+x)^2 = q_2(0:t) + 2q_1(0;t)x + q_0(0;t)x^2 = 1+3t+3t^2 + 2(1+2t)x + x^2. Evaluating the polynomials at x=0 regenerates the base sequence. With a simple sign change in R, R generates the Appell polynomials associated with A248727. - Tom Copeland, Jan 23 2015

For a natural refinement of this array, see A263634. - Tom Copeland, Nov 06 2015

From Wolfdieter Lang, Mar 13 2017: (Start)

The e.g.f. E(n, x) for {S(n, m)}{m>=0} with S(n, m) = Sum{k=1..m} k^n, n >= 0, (with undefined sum put to 0) is exp(x)*R(n+1, x) with the exponential row polynomials R(n, x) = Sum_{k=1..n} a(n, k)*x^k/k!. E.g., e.g.f. for n = 2, A000330: exp(x)*(1*x/1!+3*x^2/2!+2*x^3/3!).

The o.g.f. G(n, x) for {S(n, m)}{m >=0} is then found by Laplace transform to be G(n, 1/p) = p*Sum{k=1..n} a(n+1, k)/(p-1)^(2+k).

Hence G(n, x) = x/(1 - x)^(n+2)*Sum_{k=1..n} A008292(n,k)*x^(k-1).

E.g., n=2: G(2, 1/p) = p*(1/(p-1)^2 + 3/(p-1)^3 + 2/(p-1)^4) = p^2*(1+p)/(p-1)^4; hence G(2, x) = x*(1+x)/(1-x)^4.

This works also backwards: from the o.g.f. to the e.g.f. of {S(n, m)}_{m>=0}. (End)

a(n,k) is the number of k-tuples of pairwise disjoint and nonempty subsets of a set of size n. - Dorian Guyot, May 21 2019

From Rajesh Kumar Mohapatra, Mar 16 2020: (Start)

a(n-1,k) is the number of chains of length k in a partially ordered set formed from subsets of an n-element set ordered by inclusion such that the first term of the chains is either the empty set or an n-element set.

Also, a(n-1,k) is the number of distinct k-level rooted fuzzy subsets of an n-set ordered by set inclusion. (End)

The relations on p. 34 of Hasan (also p. 17 of Franco and Hasan) agree with the relation between A019538 and this entry given in the formula section. - Tom Copeland, May 14 2020

T(n,k) is the size of the Green's L-classes in the D-classes of rank (k-1) in the semigroup of partial transformations on an (n-1)-set. - Geoffrey Critzer, Jan 09 2023

T(n,k) is the number of strongly connected binary relations on [n] that have period k (A367948) and index 1. See Theorem 5.4.25(6) in Ki Hang Kim reference. - Geoffrey Critzer, Dec 07 2023

Stirling transform of A052849(n-1) = [1,2,4,12,48,...] is a(n-1) =[1,3,11,51,299,...]. - Michael Somos, Mar 04 2004

It is interesting to note that a chain in the power set of a set X can be thought of as a fuzzy subset of X and conversely. Chains originating with empty set are fuzzy subsets with empty core and those chains not ending with the whole set are with support strictly contained in X. - Venkat Murali (v.murali(AT)ru.ac.za), May 18 2005

Equals the binomial transform of A000629: (1, 2, 6, 26, 150, 1082, ...) and the double binomial transform of A000670: (1, 1, 3, 13, 75, 541, ...). - Gary W. Adamson, Aug 04 2009

Row sums of A038719. - Peter Bala, Jul 09 2014

Also the number of restricted barred preferential arrangements of an n-set having two bars, where one fixed section is a free section and the other two sections are restricted sections. - Sithembele Nkonkobe, Jun 16 2015

The number of all predictable outcomes of a race between a given number registered competitors, where clean finishes, dead heats (ties), disqualifications, cancellations and their combinations are all counted. (11 outcomes for two competitors, 51 for three, 299 for four, etc.. Example for two competitors shown below.) - Gergely Földvári, Jul 28 2024

Next-to-last diagonal of A038719.

a(n-1) is the sum of the n-th entries in all cycles of all permutations of [n]. a(2) = 5 because the sum of the third entries in all cycles of all permutations of [3] ((123), (132), (12)(3), (13)(2), (1)(23), (1)(2)(3)) is 3+2+0+0+0+0 = 5. - Alois P. Heinz, May 03 2017

Previous name was: A simple grammar.

Stirling transform of A005359(n-1)=[0,0,2,0,24,0,...] is a(n-1)=[0,0,2,12,74,...]. - Michael Somos, Mar 04 2004

Stirling transform of -(-1)^n*A052566(n-1)=[1,-1,4,-6,48,...] is a(n-1)=[1,0,2,12,74,...]. - Michael Somos, Mar 04 2004

Stirling transform of A000142(n)=[0,2,6,24,120,...] is a(n)=[0,2,12,74,...]. - Michael Somos, Mar 04 2004

Stirling transform of A007680(n)=[2,10,42,216,...] is a(n+1)=[2,12,74,...]. - Michael Somos, Mar 04 2004

a(n) is the number of chains in the power set of {1,2,...,n} that do not contain the empty set and do not contain {1,2,...,n}. Equivalently, a(n) is the number of ordered set partitions of {1,2,...,n} into at least 2 classes. - Geoffrey Critzer, Sep 01 2014

For n >= 1, a(n) is the number of chains of n X n (0, 1) matrices.

a(n) is also the number of chains in the power set of n^2 elements.

a(n) is the n^2-th term of A007047.

A chain of binary (crisp or Boolean or logical) matrices of order n can be thought of as a fuzzy matrix of order n.

a(n) is the number of distinct n X n fuzzy matrices.

a(n) is the sum of the n^2-th row of triangle A038719.

For n >= 1, a(n) is the number of chains of binary reflexive symmetric matrices of order n.

The number of chains of strictly upper triangular or strictly lower triangular matrices.

Also, number of chains in power set of (n^2-n)/2 elements.

a(n) is the number of distinct reflexive symmetric fuzzy matrices of order n.

Let T(n,k) = sum_{v=0..k} (-1)^v*v*binomial(k,v)*(v+1)^(n-1) for n >= 1, k >= 1 and additionally T(0,0) = 1. Then a(n) = sum_{k=0..n} abs(T(n,k)).

a(n) = A073146(n,n-1) for n >= 1.

a(n) appears to be the total number of subsets over all chains of the poset on the powerset of {1,2,...,n-1} ordered by set inclusion. That is, a(n) = Sum_{k=0..n} A038719(n,k)*(k+1). For example a(2)=4 because there are three chains: {}; {1}; {},{1}; and there are 4 total subsets. - Geoffrey Critzer, Nov 28 2014

T(n,k) = sum_{v=0..k} (-1)^v*v*binomial(k,v)*(v+1)^(n-1)

for n >= 1, k >= 1; by convention T(0,0) = 1.

Gives a representation of the Bernoulli numbers B_{n} = B_{n}(1) (with B_1 = 1/2)

B_{n} = sum_{j=0..n} sum_{k=0..j} T(j,k)/(k+1)

T(n,1) = -2^(n-1) (n>=1)

T(n,n) = (-1)^n*n! (n>=1)

sum_{k=0..n} T(n,k) = -A000007(n-1) = -1,0,0,0,0,... (n>=1)

sum_{k=0..n} abs(T(n,k)) = A162509(n) = A073146(n,n-1) (n>=1)

sum_{k=0..n} T(n,k)/(k+1) = Bernoulli(n,1)-Bernoulli(n-1,1) (n>=1)

numer(sum(T(n,k)/(k+1),k=0..n)) = A051716(n) (n>=0)

denom(sum(T(n,k)/(k+1),k=0..n)) = A051717(n) (n>=0)

Contribution from Peter Luschny, Jul 08 2009: (Start)

More generally, define the polynomials (assume p[0,0](x)=1 and 0^0=1)

p[n,k](x) = sum_{v=0..k} (-1)^v*v*binomial(k,v)*(v+1+x)^(n-1)

[1], [0, -1], [0, -2-x, 2], [0, -4-4x-x^2, 10+4x, -6], ...

then T(n,k)=p[n,k](0) and (-1)^k*k!*Stirling2(n,k)=p[n,k](-1) (cf. A019538).

Assume now k >= 1 and read by rows. Then

p[n,k](1) = -1,-3,2,-9,14,-6,-27,74,-72,24,-81,350,-582,432,-120,...

(-1)^n*(-2)^(n-k)*p[n,k](-1/2))=1,3,2,9,16,6,27,98,90,24,81,544,924,576,120,..

(-1)^n*(-2)^(n-k)*p[n,k](-3/2))=1,1,2,1,8,6,1,26,54,24,1,80,348,384,120,... (End)

Variant of A199400.