A123125 -id:A123125 - cp's OEIS frontend

A155585 a(n) = 2^n*E(n, 1) where E(n, x) are the Euler polynomials.

1, 1, 0, -2, 0, 16, 0, -272, 0, 7936, 0, -353792, 0, 22368256, 0, -1903757312, 0, 209865342976, 0, -29088885112832, 0, 4951498053124096, 0, -1015423886506852352, 0, 246921480190207983616, 0, -70251601603943959887872, 0, 23119184187809597841473536, 0

Offset: 0

Paul D. Hanna, Jan 24 2009

Previous name was: a(n) = Sum_{k=0..n-1} (-1)^(k)*C(n-1,k)*a(n-1-k)*a(k) for n>0 with a(0)=1.
Factorials have a similar recurrence: f(n) = Sum_{k=0..n-1} C(n-1,k)*f(n-1-k)*f(k), n > 0.
Related to A102573: letting T(q,r) be the coefficient of n^(r+1) in the polynomial 2^(q-n)/n times Sum_{k=0..n} binomial(n,k)*k^q, then A155585(x) = Sum_{k=0..x-1} T(x,k)*(-1)^k. See Mathematica code below. - John M. Campbell, Nov 16 2011
For the difference table and the relation to the Seidel triangle see A239005. - Paul Curtz, Mar 06 2014
From Tom Copeland, Sep 29 2015: (Start)
Let z(t) = 2/(e^(2t)+1) = 1 + tanh(-t) = e.g.f.(-t) for this sequence = 1 - t + 2*t^3/3! - 16*t^5/5! + ... .
dlog(z(t))/dt = -z(-t), so the raising operators that generate Appell polynomials associated with this sequence, A081733, and its reciprocal, A119468, contain z(-d/dx) = e.g.f.(d/dx) as the differential operator component.
dz(t)/dt = z*(z-2), so the assorted relations to a Ricatti equation, the Eulerian numbers A008292, and the Bernoulli numbers in the Rzadkowski link hold.
From Michael Somos's formula below (drawing on the Edwards link), y(t,1)=1 and x(t,1) = (1-e^(2t))/(1+e^(2t)), giving z(t) = 1 + x(t,1). Compare this to the formulas in my list in A008292 (Sep 14 2014) with a=1 and b=-1,
A) A(t,1,-1) = A(t) = -x(t,1) = (e^(2t)-1)/(1+e^(2t)) = tanh(t) = t + -2*t^3/3! + 16*t^5/5! + -272*t^7/7! + ... = e.g.f.(t) - 1 (see A000182 and A000111)
B) Ainv(t) = log((1+t)/(1-t))/2 = tanh^(-1)(t) = t + t^3/3 + t^5/5 + ..., the compositional inverse of A(t)
C) dA/dt = (1-A^2), relating A(t) to a Weierstrass elliptic function
D) ((1-t^2)d/dt)^n t evaluated at t=0, a generator for the sequence A(t)
F) FGL(x,y)= (x+y)/(1+xy) = A(Ainv(x) + Ainv(y)), a related formal group law corresponding to the Lorentz FGL (Lorentz transformation--addition of parallel velocities in special relativity) and the Atiyah-Singer signature and the elliptic curve (1-t^2)*s = t^3 in Tate coordinates according to the Lenart and Zainoulline link and the Buchstaber and Bunkova link (pp. 35-37) in A008292.
A133437 maps the reciprocal odd natural numbers through the refined faces of associahedra to a(n).
A145271 links the differential relations to the geometry of flow maps, vector fields, and thereby formal group laws. See Mathworld for links of tanh to other geometries and statistics.
Since the a(n) are related to normalized values of the Bernoulli numbers and the Riemann zeta and Dirichlet eta functions, there are links to Witten's work on volumes of manifolds in two-dimensional quantum gauge theories and the Kervaire-Milnor formula for homotopy groups of hyperspheres (see my link below).
See A101343, A111593 and A059419 for this and the related generator (1 + t^2) d/dt and associated polynomials. (End)
With the exception of the first term (1), entries are the alternating sums of the rows of the Eulerian triangle, A008292. - Gregory Gerard Wojnar, Sep 29 2018

			E.g.f.: 1 + x - 2*x^3/3! + 16*x^5/5! - 272*x^7/7! + 7936*x^9/9! -+ ... = exp(x)/cosh(x).
O.g.f.: 1 + x - 2*x^3 + 16*x^5 - 272*x^7 + 7936*x^9 - 353792*x^11 +- ...
O.g.f.: 1 + x/(1+2*x) + 2!*x^2/((1+2*x)*(1+4*x)) + 3!*x^3/((1+2*x)*(1+4*x)*(1+6*x)) + ...

Vincenzo Librandi, Table of n, a(n) for n = 0..200
P. Barry, Riordan Arrays, Orthogonal Polynomials as Moments, and Hankel Transforms, J. Int. Seq. 14 (2011) # 11.2.2, example 28.
Tom Copeland, The Elliptic Lie Triad: Ricatti and KdV Equations, Infinigens, and Elliptic Genera
Tom Copeland, The Kervaire-Milnor formula
H. M. Edwards, A normal form for elliptic curves, Bull. Amer. Math. Soc. (N.S.) 44 (2007), no. 3, 393-422. see section 12, pp. 410-411.
Wenjie Fang, Hsien-Kuei Hwang, and Mihyun Kang, Phase transitions from exp(n^(1/2)) to exp(n^(2/3)) in the asymptotics of banded plane partitions, arXiv:2004.08901 [math.CO], 2020.
G. Rzadkowski, Bernoulli numbers and solitons-revisited, Journal of Nonlinear Math. Physics, 1711, pp. 121-126.

Cf. A001057, A009006, A102573.
Equals row sums of A119879. - Johannes W. Meijer, Apr 20 2011
Cf. A081733, A119468, A008292, A000182, A000111, A101343, A111593, A059419, A133437, A145271.
(-1)^n*a(n) are the alternating row sums of A123125. - Wolfdieter Lang, Jul 12 2017

A155585 := n -> 2^n*euler(n, 1): # Peter Luschny, Jan 26 2009
a := proc(n) option remember; `if`(n::even, 0^n, -(-1)^n - add((-1)^k*binomial(n,k) *a(n-k), k = 1..n-1)) end: # Peter Luschny, Jun 01 2016
# Or via the recurrence of the Fubini polynomials:
F := proc(n) option remember; if n = 0 then return 1 fi;
expand(add(binomial(n, k)*F(n-k)*x, k = 1..n)) end:
a := n -> (-2)^n*subs(x = -1/2, F(n)):
seq(a(n), n = 0..30); # Peter Luschny, May 21 2021

a[m_] := Sum[(-2)^(m - k) k! StirlingS2[m, k], {k, 0, m}] (* Peter Luschny, Apr 29 2009 *)
poly[q_] :=  2^(q-n)/n*FunctionExpand[Sum[Binomial[n, k]*k^q, {k, 0, n}]]; T[q_, r_] :=  First[Take[CoefficientList[poly[q], n], {r+1, r+1}]]; Table[Sum[T[x, k]*(-1)^k, {k, 0, x-1}], {x, 1, 16}] (* John M. Campbell, Nov 16 2011 *)
f[n_] := (-1)^n 2^(n+1) PolyLog[-n, -1]; f[0] = -f[0]; Array[f, 27, 0] (* Robert G. Wilson v, Jun 28 2012 *)

a(n)=if(n==0,1,sum(k=0,n-1,(-1)^(k)*binomial(n-1,k)*a(n-1-k)*a(k)))

a(n)=local(X=x+x*O(x^n));n!*polcoeff(exp(X)/cosh(X),n)

a(n)=polcoeff(sum(m=0,n,m!*x^m/prod(k=1,m,1+2*k*x+x*O(x^n))),n) \\ Paul D. Hanna, Jul 20 2011

{a(n) = local(A); if( n<0, 0, A = x * O(x^n); n! * polcoeff( 1 + sinh(x + A) / cosh(x + A), n))} /* Michael Somos, Jan 16 2012 */

a(n)=local(A=1+x);for(i=1,n,A=sum(k=0,n,intformal(subst(A,x,-x)+x*O(x^n))^k/k!));n!*polcoeff(A,n)
for(n=0,30,print1(a(n),", ")) \\ Paul D. Hanna, Nov 25 2013

from sympy import bernoulli
def A155585(n): return (((2<<(m:=n+1))-2)*bernoulli(m)<>1) if n&1 else (0 if n else 1) # Chai Wah Wu, Apr 14 2023

def A155585(n) :
    if n == 0 : return 1
    return add(add((-1)^(j+1)*binomial(n+1,k-j)*j^n for j in (0..k)) for k in (1..n))
[A155585(n) for n in (0..26)] # Peter Luschny, Jul 23 2012

def A155585_list(n): # Akiyama-Tanigawa algorithm
    A = [0]*(n+1); R = []
    for m in range(n+1) :
        d = divmod(m+3, 4)
        A[m] = 0 if d[1] == 0 else (-1)^d[0]/2^(m//2)
        for j in range(m, 0, -1) :
            A[j - 1] = j * (A[j - 1] - A[j])
        R.append(A[0])
    return R
A155585_list(30) # Peter Luschny, Mar 09 2014

E.g.f.: exp(x)*sech(x) = exp(x)/cosh(x). (See A009006.) - Paul Barry, Mar 15 2006
Sequence of absolute values is A009006 (e.g.f. 1+tan(x)).
O.g.f.: Sum_{n>=0} n! * x^n / Product_{k=1..n} (1 + 2*k*x). - Paul D. Hanna, Jul 20 2011
a(n) = 2^n*E_{n}(1) where E_{n}(x) are the Euler polynomials. - Peter Luschny, Jan 26 2009
a(n) = EL_{n}(-1) where EL_{n}(x) are the Eulerian polynomials. - Peter Luschny, Aug 03 2010
a(n+1) = (4^n-2^n)*B_n(1)/n, where B_{n}(x) are the Bernoulli polynomials (B_n(1) = B_n for n <> 1). - Peter Luschny, Apr 22 2009
G.f.: 1/(1-x+x^2/(1-x+4*x^2/(1-x+9*x^2/(1-x+16*x^2/(1-...))))) (continued fraction). - Paul Barry, Mar 30 2010
G.f.: -log(x/(exp(x)-1))/x = Sum_{n>=0} a(n)*x^n/(2^(n+1)*(2^(n+1)-1)*n!). - Vladimir Kruchinin, Nov 05 2011
E.g.f.: exp(x)/cosh(x) = 2/(1+exp(-2*x)) = 2/(G(0) + 1); G(k) = 1 - 2*x/(2*k + 1 - x*(2*k+1)/(x - (k+1)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Dec 10 2011
E.g.f. is x(t,1) + y(t,1) where x(t,a) and y(t,a) satisfy y(t,a)^2 = (a^2 - x(t,a)^2) / (1 - a^2 * x(t,a)^2) and dx(t,a) / dt = y(t,a) * (1 - a * x(t,a)^2) and are the elliptic functions of Edwards. - Michael Somos, Jan 16 2012
E.g.f.: 1/(1 - x/(1+x/(1 - x/(3+x/(1 - x/(5+x/(1 - x/(7+x/(1 - x/(9+x/(1 +...))))))))))), a continued fraction. - Paul D. Hanna, Feb 11 2012
E.g.f. satisfies: A(x) = Sum_{n>=0} Integral( A(-x) dx )^n / n!. - Paul D. Hanna, Nov 25 2013
a(n) = -2^(n+1)*Li_{-n}(-1). - Peter Luschny, Jun 28 2012
a(n) = Sum_{k=1..n} Sum_{j=0..k} (-1)^(j+1)*binomial(n+1,k-j)*j^n for n > 0. - Peter Luschny, Jul 23 2012
From Sergei N. Gladkovskii, Oct 25 2012 to Dec 16 2013: (Start)
Continued fractions:
G.f.: 1 + x/T(0) where T(k) = 1 + (k+1)*(k+2)*x^2/T(k+1).
E.g.f.: exp(x)/cosh(x) = 1 + x/S(0) where S(k) = (2*k+1) + x^2/S(k+1).
E.g.f.: 1 + x/(U(0)+x) where U(k) = 4*k+1 - x/(1 + x/(4*k+3 - x/(1 + x/U(k+1)))).
E.g.f.: 1 + tanh(x) = 4*x/(G(0)+2*x) where G(k) = 1 - (k+1)/(1 - 2*x/(2*x + (k+1)^2/G(k+1)));
G.f.: 1 + x/G(0) where G(k) = 1 + 2*x^2*(2*k+1)^2 - x^4*(2*k+1)*(2*k+2)^2*(2*k+3)/G(k+1) (due to Stieltjes).
E.g.f.: 1 + x/(G(0) + x) where G(k) = 1 - 2*x/(1 + (k+1)/G(k+1)).
G.f.: 2 - 1/Q(0) where Q(k) = 1 + x*(k+1)/( 1 - x*(k+1)/Q(k+1)).
G.f.: 2 - 1/Q(0) where Q(k) = 1 + x*k^2 + x/(1 - x*(k+1)^2/Q(k+1)).
G.f.: 1/Q(0) where Q(k) = 1 - 2*x + x*(k+1)/(1-x*(k+1)/Q(k+1)).
G.f.: 1/Q(0) where Q(k) = 1 - x*(k+1)/(1 + x*(k+1)/Q(k+1)).
E.g.f.: 1 + x*Q(0) where Q(k) = 1 - x^2/( x^2 + (2*k+1)*(2*k+3)/Q(k+1)).
G.f.: 2 - T(0)/(1+x) where T(k) = 1 - x^2*(k+1)^2/(x^2*(k+1)^2 + (1+x)^2/T(k+1)).
E.g.f.: 1/(x - Q(0)) where Q(k) = 4*k^2 - 1 + 2*x + x^2*(2*k-1)*(2*k+3)/Q(k+1). (End)
G.f.: 1 / (1 - b(1)*x / (1 - b(2)*x / (1 - b(3)*x / ... ))) where b = A001057. - Michael Somos, Jan 03 2013
From Paul Curtz, Mar 06 2014: (Start)
a(2n) = A000007(n).
a(2n+1) = (-1)^n*A000182(n+1).
a(n) is the binomial transform of A122045(n).
a(n) is the row sum of A081658. For fractional Euler numbers see A238800.
a(n) + A122045(n) = 2, 1, -1, -2, 5, 16, ... = -A163982(n).
a(n) - A122045(n) = -A163747(n).
a(n) is the Akiyama-Tanigawa transform applied to 1, 0, -1/2, -1/2, -1/4, 0, ... = A046978(n+3)/A016116(n). (End)
a(n) = 2^(2*n+1)*(zeta(-n,1/2) - zeta(-n, 1)), where zeta(a, z) is the generalized Riemann zeta function. - Peter Luschny, Mar 11 2015
a(n) = 2^(n + 1)*(2^(n + 1) - 1)*Bernoulli(n + 1, 1)/(n + 1). (From Bill Gosper, Oct 28 2015) - N. J. A. Sloane, Oct 28 2015 [See the above comment from Peter Luschny, Apr 22 2009.]
a(n) = -(n mod 2)*((-1)^n + Sum_{k=1..n-1} (-1)^k*C(n,k)*a(n-k)) for n >= 1. - Peter Luschny, Jun 01 2016
a(n) = (-2)^n*F_{n}(-1/2), where F_{n}(x) is the Fubini polynomial. - Peter Luschny, May 21 2021

New name from Peter Luschny, Mar 12 2015

A114088 Triangle read by rows: T(n,k) is number of partitions of n whose tail below its Durfee square has k parts (n >= 1; 0 <= k <= n-1).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 3, 3, 2, 1, 1, 1, 3, 4, 3, 2, 1, 1, 1, 4, 5, 5, 3, 2, 1, 1, 1, 5, 6, 6, 5, 3, 2, 1, 1, 1, 6, 8, 8, 7, 5, 3, 2, 1, 1, 1, 7, 10, 10, 9, 7, 5, 3, 2, 1, 1, 1, 9, 13, 13, 12, 10, 7, 5, 3, 2, 1, 1, 1, 10, 16, 17, 15, 13, 10, 7, 5, 3, 2, 1, 1, 1, 12, 20, 22, 20, 17

Offset: 1

Emeric Deutsch, Feb 12 2006

From Gus Wiseman, May 21 2022: (Start)
Also the number of integer partitions of n with k parts below the diagonal. For example, the partition (3,2,2,1) has two parts (at positions 3 and 4) below the diagonal (1,2,3,4). Row n = 8 counts the following partitions:
  8    71    611    5111    41111    311111   2111111   11111111
  44   332   2222   22211   221111
  53   422   3221   32111
  62   431   3311
       521   4211
Indices of parts below the diagonal are also called strong nonexcedances.
(End)

			T(7,2)=3 because we have [5,1,1], [3,2,1,1] and [2,2,2,1] (the bottom tails are [1,1], [1,1] and [2,1], respectively).
Triangle starts:
  1;
  1, 1;
  1, 1, 1;
  2, 1, 1, 1;
  2, 2, 1, 1, 1;
  3, 3, 2, 1, 1, 1;
  3, 4, 3, 2, 1, 1, 1;

G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976 (pp. 27-28).
G. E. Andrews and K. Eriksson, Integer Partitions, Cambridge Univ. Press, 2004 (pp. 75-78).

John Tyler Rascoe, Rows n = 1..140, flattened
Findstat, St000183: The side length of the Durfee square of an integer partition

Cf. A114087, A114089, A115995, A116365.
Row sums: A000041.
Column k = 0: A003114.
Weak opposite: A115994.
Permutations: A173018, weak A123125.
Ordered: A352521, rank stat A352514, weak A352522.
Opposite ordered: A352524, first col A008930, rank stat A352516.
Weak opposite ordered: A352525, first col A177510, rank stat A352517.
Weak: A353315.
Opposite: A353318.
A000700 counts self-conjugate partitions, ranked by A088902.
A115720 counts partitions by Durfee square, rank stat A257990.
A352490 gives the (strong) nonexcedance set of A122111, counted by A000701.
Cf. A002620, A006918, A219282, A238352, A238874, A325039, A330644, A352487.

g:=sum(z^(k^2)/product((1-z^j)*(1-t*z^j),j=1..k),k=1..20): gserz:=simplify(series(g,z=0,30)): for n from 1 to 14 do P[n]:=coeff(gserz,z^n) od: for n from 1 to 14 do seq(coeff(t*P[n],t^j),j=1..n) od; # yields sequence in triangular form

subdiags[y_]:=Length[Select[Range[Length[y]],#>y[[#]]&]];
Table[Length[Select[IntegerPartitions[n],subdiags[#]==k&]],{n,1,15},{k,0,n-1}] (* Gus Wiseman, May 21 2022 *)

T_qt(max_row) = {my(N=max_row+1, q='q+O('q^N), h = sum(k=1,N, q^(k^2)/prod(j=1,k, (1-q^j)*(1-t*q^j))) ); for(i=1, N-1, print(Vecrev(polcoef(h, i))))}
T_qt(10) \\ John Tyler Rascoe, Oct 24 2024

G.f. = Sum_{k>=1} q^(k^2) / Product_{j=1..k} (1 - q^j)*(1 - t*q^j).

Sum_{k=0..n-1} k*T(n,k) = A114089(n).

A002662 a(n) = 2^n - 1 - n*(n+1)/2.

Original entry on oeis.org

0, 0, 0, 1, 5, 16, 42, 99, 219, 466, 968, 1981, 4017, 8100, 16278, 32647, 65399, 130918, 261972, 524097, 1048365, 2096920, 4194050, 8388331, 16776915, 33554106, 67108512, 134217349, 268435049, 536870476, 1073741358, 2147483151, 4294966767, 8589934030

Offset: 0

N. J. A. Sloane

Number of subsets with at least 3 elements of an n-element set.
For n>4, number of simple rank-(n-1) matroids over S_n.
Number of non-interval subsets of {1,2,3,...,n} (cf. A000124). - Jose Luis Arregui (arregui(AT)unizar.es), Jun 27 2006
The partial sums of the second diagonal of A008292 or third column of A123125. - Tom Copeland, Sep 09 2008
a(n) is the number of binary sequences of length n having at least three 0's. - Geoffrey Critzer, Feb 11 2009
Starting with "1" = eigensequence of a triangle with the tetrahedral numbers (1, 4, 10, 20, ...) as the left border and the rest 1's. - Gary W. Adamson, Jul 24 2010
a(n) is also the number of crossing set partitions of [n+1] with two blocks. - Peter Luschny, Apr 29 2011
The Kn24 sums, see A180662, of triangle A065941 equal the terms (doubled) of this sequence minus the three leading zeros. - Johannes W. Meijer, Aug 14 2011
From L. Edson Jeffery, Dec 28 2011: (Start)
Nonzero terms of this sequence can be found from the row sums of the fourth sub-triangle extracted from Pascal's triangle as indicated below by braces:
                     1;
                  1,    1;
               1,    2,    1;
           {1},   3,    3,    1;
        {1,    4},   6,    4,    1;
     {1,    5,   10},  10,    5,    1;
  {1,    6,   15,   20},  15,    6,    1;
  ... (End)
Partial sums of A000295 (Eulerian Numbers, Column 2).
Second differences equal 2^(n-2) - 1, for n >= 4. - Richard R. Forberg, Jul 11 2013
Starting (0, 0, 1, 5, 16, ...) is the binomial transform of (0, 0, 1, 2, 2, 2, ...). - Gary W. Adamson, Jul 27 2015
a(n - 1) is the rank of the divisor class group of the moduli space of stable rational curves with n marked points, see Keel p. 550. - Harry Richman, Aug 10 2024

			a(4) = 5 is the number of crossing set partitions of {1,2,..,5}, card{13|245, 14|235, 24|135, 25|134, 35|124}. - _Peter Luschny_, Apr 29 2011

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).

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VI: Voronoi Reduction of Three-Dimensional Lattices, Proc. Royal Soc. London, Series A, 436 (1992), 55-68. (See Table 1.)
W. M. B. Dukes, On the number of matroids on a finite set, arXiv:math/0411557 [math.CO], 2004.
J. Eckhoff, Der Satz von Radon in konvexen Produktstrukturen II, Monat. f. Math., 73 (1969), 7-30.
R. K. Guy, Letter to N. J. A. Sloane
Sean Keel, Intersection theory of moduli space of stable n-pointed curves of genus zero, Trans. Amer. Math. Soc. 330 (1992), 545-574.
Pablo Hueso Merino , The first problem from the 55th Spanish Mathematical Olympiad asks to find the value of a(2019) (see comment from Jose Luis Arregui).
Ângela Mestre and José Agapito, Square Matrices Generated by Sequences of Riordan Arrays, J. Int. Seq., Vol. 22 (2019), Article 19.8.4.
Sergey V. Muravyov, Liudmila I. Khudonogova, and Ekaterina Y. Emelyanova, Interval data fusion with preference aggregation, Measurement (2017), see page 5.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Index entries for linear recurrences with constant coefficients, signature (5,-9,7,-2).

a(n) = A055248(n,3).
First differences are A000295.
Cf. A000079, A000124, A000217, A000225, A000247, A002663, A002664, A035038-A035042, A007318, A342352.
Cf. also A000290, A001045.

a002662 n = a002662_list !! n
a002662_list = map (sum . drop 3) a007318_tabl
-- Reinhard Zumkeller, Jun 20 2015

[2^n - 1 - n*(n+1)/2: n in [0..35]]; // Vincenzo Librandi, May 20 2011

A002662 := z**2/(2*z-1)/(z-1)**3; # conjectured by Simon Plouffe in his 1992 dissertation
A002662 := proc(n): 2^n - 1 - n*(n+1)/2 end: seq(A002662(n), n=0..33); # Johannes W. Meijer, Aug 14 2011

With[{nn=40},Join[{0},First[#]-1-Last[#]&/@Thread[{2^Range[nn], Accumulate[ Range[nn]]}]]] (* Harvey P. Dale, May 10 2012 *)
Table[2^n - Binomial[n, 2] - n - 1, {n, 1, 100}] (* Pablo Hueso Merino, Dec 17 2019 *)

a(n)=2^n-1-n*(n+1)/2 \\ Charles R Greathouse IV, Oct 11 2015

def A002662(n): return (1<>1) # Chai Wah Wu, Aug 29 2023

G.f.: x^3/((1-2*x)*(1-x)^3).
a(n) = Sum_{k=0..n} binomial(n,k+3) = Sum_{k=3..n} binomial(n,k). - Paul Barry, Jul 30 2004
a(n+1) = 2*a(n) + binomial(n,2). - Paul Barry, Aug 23 2004
(1, 5, 16, 42, 99, ...) = binomial transform of (1, 4, 7, 8, 8, 8, ...). - Gary W. Adamson, Sep 30 2007
E.g.f.: exp(x)*(exp(x)-x^2/2-x-1). - Geoffrey Critzer, Feb 11 2009
a(n) = n - 2 + 3*a(n-1) - 2*a(n-2), for n >= 2. - Richard R. Forberg, Jul 11 2013
For n>1, a(n) = (1/4)*Sum_{k=1..n-2} 2^k*(n-k-1)*(n-k). For example, (1/4)*(2^1*(4*5) + 2^2*(3*4) + 2^3*(2*3) + 2^4*(1*2)) = 168/4 = 42. - J. M. Bergot, May 27 2014 [edited by Danny Rorabaugh, Apr 19 2015]
Convolution of A001045 and (A000290 shifted by one place). - Oboifeng Dira,  Aug 16 2016
a(n) = Sum_{k=1..n-2} Sum_{i=1..n} (n-k-1) * C(k,i). - Wesley Ivan Hurt, Sep 19 2017
a(n) = 5*a(n-1) - 9*a(n-2) + 7*a(n-3) - 2*a(n-4) for n > 3. - Chai Wah Wu, Apr 03 2021
a(n) = a(n-1) + 1 + A000247(n-1). - Harry Richman, Aug 13 2024

A352523 Number of integer compositions of n with exactly k nonfixed points (parts not on the diagonal).

Original entry on oeis.org

1, 1, 0, 0, 2, 0, 1, 1, 2, 0, 0, 4, 2, 2, 0, 0, 5, 5, 4, 2, 0, 1, 3, 12, 8, 6, 2, 0, 0, 7, 14, 19, 14, 8, 2, 0, 0, 8, 21, 33, 32, 22, 10, 2, 0, 0, 9, 30, 54, 63, 54, 32, 12, 2, 0, 1, 6, 47, 80, 116, 116, 86, 44, 14, 2, 0, 0, 11, 53, 129, 194, 229, 202, 130, 58, 16, 2, 0

Offset: 0

Gus Wiseman, Mar 26 2022

A nonfixed point in a composition c is an index i such that c_i != i.

			Triangle begins:
   1
   1   0
   0   2   0
   1   1   2   0
   0   4   2   2   0
   0   5   5   4   2   0
   1   3  12   8   6   2   0
   0   7  14  19  14   8   2   0
   0   8  21  33  32  22  10   2   0
   0   9  30  54  63  54  32  12   2   0
   1   6  47  80 116 116  86  44  14   2   0
   ...
For example, row n = 6 counts the following compositions (empty column indicated by dot):
  (123)  (6)   (24)    (231)    (2112)   (21111)    .
         (15)  (33)    (312)    (2121)   (111111)
         (42)  (51)    (411)    (3111)
               (114)   (1113)   (11112)
               (132)   (1122)   (11121)
               (141)   (1311)   (11211)
               (213)   (2211)
               (222)   (12111)
               (321)
               (1131)
               (1212)
               (1221)

John Tyler Rascoe, Rows n = 0..130, flattened

Column k = 0 is A010054.
Row sums are A011782.
The version for permutations is A098825.
The corresponding rank statistic is A352513.
Column k = 1 is A352520.
A238349 and A238350 count comps by fixed points, first col A238351, rank stat A352512.
A352486 gives the nonfixed points of A122111, counted by A330644.
A352521 counts comps by strong nonexcedances, first A219282, stat A352514.
A352522 counts comps by weak nonexcedances, first col A238874, stat A352515.
A352524 counts comps by strong excedances, first col A008930, stat A352516.
A352525 counts comps by weak excedances, first col A177510, stat A352517.
Cf. A000700, A064428, A088218, A114088, A115994, A123125, A173018, A238352.

b:= proc(n, i) option remember; expand(`if`(n=0, 1,
      add(`if`(i=j, 1, x)*b(n-j, i+1), j=1..n)))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=0..n))(b(n, 1)):
seq(T(n), n=0..12);  # Alois P. Heinz, Mar 19 2025

pnq[y_]:=Length[Select[Range[Length[y]],#!=y[[#]]&]];
Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],pnq[#]==k&]],{n,0,9},{k,0,n}]

T_xy(max_row) = {my(N=max_row+1, x='x+O('x^N), h= sum(i=0, N, prod(j=1, i, y*(x/(1-x)-x^j)+x^j))); vector(N, n, my(r=Vecrev(polcoeff(h, n-1))); if(n<2, r, concat(r,[0])))}
T_xy(10) \\ John Tyler Rascoe, Mar 21 2025

G.f.: Sum_{i>=0} Product_{j=1..i} y*(x/(1-x) - x^j) + x^j. - John Tyler Rascoe, Mar 19 2025

A009998 Triangle in which j-th entry in i-th row is (j+1)^(i-j).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 3, 1, 1, 8, 9, 4, 1, 1, 16, 27, 16, 5, 1, 1, 32, 81, 64, 25, 6, 1, 1, 64, 243, 256, 125, 36, 7, 1, 1, 128, 729, 1024, 625, 216, 49, 8, 1, 1, 256, 2187, 4096, 3125, 1296, 343, 64, 9, 1, 1, 512, 6561, 16384, 15625, 7776, 2401, 512, 81, 10, 1

Offset: 0

N. J. A. Sloane

Read as a square array this is the Hilbert transform of triangle A123125 (see A145905 for the definition of this term). For example, the fourth row of A123125 is (0,1,4,1) and the expansion (x + 4*x^2 + x^3)/(1-x)^4 = x + 8*x^2 + 27*x^3 + 64*x^4 + ... generates the entries in the fourth row of this array read as a square. - Peter Bala, Oct 28 2008

			Triangle begins:
  1;
  1,  1;
  1,  2,  1;
  1,  4,  3,  1;
  1,  8,  9,  4,  1;
  1, 16, 27, 16,  5,  1;
  1, 32, 81, 64, 25,  6,  1;
  ...
From _Gus Wiseman_, May 01 2021: (Start)
The rows of the triangle are obtained by reading antidiagonals upward in the following table of A(k,n) = n^k, with offset k = 0, n = 1:
         n=1:     n=2:     n=3:     n=4:     n=5:     n=6:
   k=0:   1        1        1        1        1        1
   k=1:   1        2        3        4        5        6
   k=2:   1        4        9       16       25       36
   k=3:   1        8       27       64      125      216
   k=4:   1       16       81      256      625     1296
   k=5:   1       32      243     1024     3125     7776
   k=6:   1       64      729     4096    15625    46656
   k=7:   1      128     2187    16384    78125   279936
   k=8:   1      256     6561    65536   390625  1679616
   k=9:   1      512    19683   262144  1953125 10077696
  k=10:   1     1024    59049  1048576  9765625 60466176
(End)

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 24.

T. D. Noe, Rows n=0..50 of triangle, flattened
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Peter Luschny, Figurate number - a very short introduction

Row sums give A026898.
Cf. A088956, A123125, A179927.
Column n = 2 of the array is A000079.
Column n = 3 of the array is A000244.
Row k = 2 of the array is A000290.
Row k = 3 of the array is A000578.
Diagonal n = k of the array is A000312.
Diagonal n = k + 1 of the array is A000169.
Diagonal n = k + 2 of the array is A000272.
The transpose of the array is A009999.
The numbers of divisors of the entries are A343656 (row sums: A343657).
A007318 counts k-sets of elements of {1..n}.
A059481 counts k-multisets of elements of {1..n}.
Cf. A002109, A062319, A066959, A143773, A176029, A204688, A326358, A327527.

a009998 n k = (k + 1) ^ (n - k)
a009998_row n = a009998_tabl !! n
a009998_tabl = map reverse a009999_tabl
-- Reinhard Zumkeller, Feb 02 2014

E := (n,x) -> `if`(n=0,1,x*(1-x)*diff(E(n-1,x),x)+E(n-1,x)*(1+(n-1)*x));
G := (n,x) -> E(n,x)/(1-x)^(n+1);
A009998 := (n,k) -> coeff(series(G(n-k,x),x,18),x,k);
seq(print(seq(A009998(n,k),k=0..n)),n=0..6);
# Peter Luschny, Aug 02 2010

Flatten[Table[(j+1)^(i-j),{i,0,20},{j,0,i}]] (* Harvey P. Dale, Dec 25 2012 *)

T(i,j)=(j+1)^(i-j) \\ Charles R Greathouse IV, Feb 06 2017

T(n,n) = 1; T(n,k) = (k+1)*T(n-1,k) for k=0..n-1. - Reinhard Zumkeller, Feb 02 2014

T(n,m) = (m+1)*Sum_{k=0..n-m}((n+1)^(k-1)*(n-m)^(n-m-k)*(-1)^(n-m-k)*binomial(n-m-1,k-1)). - Vladimir Kruchinin, Sep 12 2015

a(62) corrected to 512 by T. D. Noe, Dec 20 2007

A046802 T(n, k) = Sum_{j=k..n} binomial(n, j)*E1(j, j-k), where E1 are the Eulerian numbers A173018. Triangle read by rows, T(n, k) for 0 <= k <= n.

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 7, 7, 1, 1, 15, 33, 15, 1, 1, 31, 131, 131, 31, 1, 1, 63, 473, 883, 473, 63, 1, 1, 127, 1611, 5111, 5111, 1611, 127, 1, 1, 255, 5281, 26799, 44929, 26799, 5281, 255, 1, 1, 511, 16867, 131275, 344551, 344551, 131275, 16867, 511, 1, 1, 1023, 52905

Offset: 0

N. J. A. Sloane

T(n,k) is the number of positroid cells of the totally nonnegative Grassmannian G+(k,n) (cf. Postnikov/Williams). It is the triangle of the h-vectors of the stellahedra. - Tom Copeland, Oct 10 2014
See A248727 for a simple transformation of the row polynomials of this entry that produces the umbral compositional inverses of the polynomials of A074909, related to the face polynomials of the simplices. - Tom Copeland, Jan 21 2015
From Tom Copeland, Jan 24 2015: (Start)
The reciprocal of this entry's e.g.f. is [t e^(-xt) - e^(-x)] / (t-1) = 1 - (1+t) x + (1+t+t^2) x^2/2! - (1+t+t^2+t^3) x^3/3! + ... = e^(q.(0;t)x), giving the base sequence (q.(0;t))^n = q_n(0;t) = (-1)^n [1-t^(n+1)] / (1-t) for the umbral compositional inverses (q.(0;t)+z)^n = q_n(z;t) of the Appell polynomials associated with this entry, p_n(z;t) below, i.e., q_n(p.(z;t)) = z^n = p_n(q.(z;t)), in umbral notation. The relations in A133314 then apply between the two sets of base polynomials. (Inserted missing index in a formula - Mar 03 2016.)
The associated o.g.f. for the umbral inverses is Og(x) = x / (1-x q.(0:t)) = x / [(1+x)(1+tx)] = x / [1+(1+t)x+tx^2]. Applying A134264 to h(x) = x / Og(x) = 1 + (1+t) x + t x^2 leads to an o.g.f. for the Narayana polynomials A001263 as the comp. inverse Oginv(x) = [1-(1+t)x-sqrt[1-2(1+t)x+((t-1)x)^2]] / (2xt). Note that Og(x) gives the signed h-polynomials of the simplices and that Oginv(x) gives the h-polynomials of the simplicial duals of the Stasheff polynomials, or type A associahedra. Contrast this with A248727 = A046802 * A007318, which has o.g.f.s related to the corresponding f-polynomials. (End)
The Appell polynomials p_n(x;t) in the formulas below specialize to the Swiss-knife polynomials of A119879 for t = -1, so the Springer numbers A001586 are given by 2^n p_n(1/2;-1). - Tom Copeland, Oct 14 2015
The row polynomials are the h-polynomials associated to the stellahedra, whose f-polynomials are the row polynomials of A248727. Cf. page 60 of the Buchstaber and Panov link. - Tom Copeland, Nov 08 2016
The row polynomials are the h-polynomials of the stellohedra, which enumerate partial permutations according to descents. Cf. Section 10.4 of the Postnikov-Reiner-Williams reference. - Lauren Williams, Jul 05 2022
From p. 60 of the Buchstaber and Panov link, S = P * C / T where S, P, C, and T are the bivariate e.g.f.s of the h vectors of the stellahedra, permutahedra, hypercubes, and (n-1)-simplices, respectively. - Tom Copeland, Jan 09 2017
The number of Le-diagrams of type (k, n) this means the diagram uses the bounding box size k x (n-k), equivalently the number of Grassmann necklaces of type (k, n) and also the number of decorated permutations with k anti-exceedances. - Thomas Scheuerle, Dec 29 2024

			The triangle T(n, k) begins:
n\k 0   1     2      3      4      5      6     7
0:  1
1:  1   1
2:  1   3     1
3:  1   7     7      1
4:  1  15    33     15      1
5:  1  31   131    131     31      1
6:  1  63   473    883    473     63      1
7:  1 127  1611   5111   5111   1611    127     1
... Reformatted. - _Wolfdieter Lang_, Feb 14 2015

L. Comtet, Advanced Combinatorics, Reidel, Holland, 1974, page 245 [From Roger L. Bagula, Nov 21 2009]
D. Singh, The numbers L(m,n) and their relations with prepared Bernoulli and Eulerian numbers, Math. Student, 20 (1952), 66-70.

Michael De Vlieger, Table of n, a(n) for n = 0..11475 (rows 0 <= n <= 150, flattened)
Paul Barry, General Eulerian Polynomials as Moments Using Exponential Riordan Arrays, Journal of Integer Sequences, 16 (2013), #13.9.6.
Paul Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.
Paul Barry, Series reversion with Jacobi and Thron continued fractions, arXiv:2107.14278 [math.NT], 2021.
Carolina Benedetti, Anastasia Chavez, and Daniel Tamayo, Quotients of uniform positroids, arXiv:1912.06873 [math.CO], 2019.
V. Buchstaber and T. Panov Toric Topology, arXiv:1210.2368v3 [math.AT], 2014.
Jan Geuenich, Tilting modules for the Auslander algebra of K[x]/(xn), arXiv:1803.10707 [math.RT], 2018.
Jun Ma, Shi-mei Ma, and Yeong-Nan Yeh, Recurrence relations for binomial-eulerian polynomials, arXiv:1711.09016 [math.CO], 2017, Thm. 2.1.
A. Postnikov, Total positivity, Grassmannians, and networks, arXiv:math/0609764 [math.CO], 2006.
A. Postnikov, V. Reiner, and L. Williams, Faces of generalized permutohedra, arXiv:math/0609184 [math.CO], 2006-2007.
D. Singh, The numbers L(m,n) and their relations with prepared Bernoulli and Eulerian numbers, Math. Student, 20 (1952), 66-70. [Annotated scanned copy]
L. K. Williams, Enumeration of totally positive Grassmann cells, arXiv:math/0307271 [math.CO], 2003-2004.
L. Williams, The Positive Grassmannian (from a mathematician's perspective), 2014

Row sums give A000522.
Cf. A008292, A123125, A248727, A074909, A007318, A000225, A066810, A028246, A001263, A119879, A001586, A019538, A090582, A123125, A130850.
Cf. A122045, A007047, A000522, A026898, A036919.

T := (n, k) -> add(binomial(n, r)*combinat:-eulerian1(r, r-k), r = k .. n):
for n from 0 to 8 do seq(T(n, k), k=0..n) od; # Peter Luschny, Jun 27 2018

t[, 1] = 1; t[n, n_] = 1; t[n_, 2] = 2^(n-1)-1;
t[n_, k_] = Sum[((i-k+1)^i*(k-i)^(n-i-1) - (i-k+2)^i*(k-i-1)^(n-i-1))*Binomial[n-1, i], {i, 0, k-1}];
T[n_, k_] := t[n+1, k+1]; Table[T[n, k], {n, 0, 12}, {k, 0, n}] // Flatten
(* Jean-François Alcover, Jan 22 2015, after Tom Copeland *)
T[ n_, k_] := Coefficient[n! SeriesCoefficient[(1-x) Exp[t] / (1 - x Exp[(1-x) t]), {t, 0, n}] // Simplify, x, k];
Table[T[n, k], {n, 0, 10}, {k, 0, n}] (* Michael Somos, Jan 22 2015 *)

E.g.f.: (y-1)*exp(x*y)/(y-exp((y-1)*x)). - Vladeta Jovovic, Sep 20 2003
p(t,x) = (1 - x)*exp(t)/(1 - x*exp(t*(1 - x))). - Roger L. Bagula, Nov 21 2009
With offset=0, T(n,0)=1 otherwise T(n,k) = sum_{i=0..k-1} C(n,i)((i-k)^i*(k-i+1)^(n-i) - (i-k+1)^i*(k-i)^(n-i)) (cf. Williams). - Tom Copeland, Oct 10 2014
With offset 0, T = A007318 * A123125. Second column is A000225; 3rd, appears to be A066810. - Tom Copeland, Jan 23 2015
A raising operator (with D = d/dx) associated with this entry's row polynomials is R = x + t + (1-t) / [1-t e^{(1-t)D}]  = x + t + 1 + t D + (t+t^2) D^2/2! + (t+4t^2+t^3) D^3/3! + ... , containing the e.g.f. for the Eulerian polynomials of A123125. Then R^n 1 = (p.(0;t)+x)^n = p_n(x;t) are the Appell polynomials with this entry's row polynomials p_n(0;t) as the base sequence. Examples of this formalism are given in A028246 and A248727. - Tom Copeland, Jan 24 2015
With offset 0, T = A007318*(padded A090582)*(inverse of A097805) = A007318*(padded A090582)*(padded A130595) = A007318*A123125 = A007318*(padded A090582)*A007318*A097808*A130595, where padded matrices are of the form of padded A007318, which is A097805. Inverses of padded matrices are just the padded versions of inverses of the unpadded matrices. This relates the face vectors, or f-vectors, and h-vectors of the permutahedra / permutohedra to those of the stellahedra / stellohedra. - Tom Copeland, Nov 13 2016
Umbrally, the row polynomials (offset 0) are r_n(x) = (1 + q.(x))^n, where (q.(x))^k = q_k(x) are the row polynomials of A123125. - Tom Copeland, Nov 16 2016
From the previous umbral statement, OP(x,d/dy) y^n =  (y + q.(x))^n, where OP(x,y) = exp[y * q.(x)] = (1-x)/(1-x*exp((1-x)y)), the e.g.f. of A123125, so OP(x,d/dy) y^n evaluated at y = 1 is  r_n(x), the n-th row polynomial of this entry, with offset 0. - Tom Copeland, Jun 25 2018
Consolidating some formulas in this entry and A248727, in umbral notation for concision, with all offsets 0: Let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of A123125. Then the row polynomials of this entry (A046802, the h-polynomials of the stellahedra) are given by h_n(x) = A_n(x;1); the row polynomials of A248727 (the face polynomials of the stellahedra), by f_n(x) = A_n(1 + x;1); the Swiss-knife polynomials of A119879, by Sw_n(x) = A_n(-1;1 + x); and the row polynomials of the Worpitsky triangle (A130850), by w_n(x) = A(1 + x;0). Other specializations of A_n(x;y) give A090582 (the f-polynomials of the permutohedra, cf. also A019538) and A028246 (another version of the Worpitsky triangle). - Tom Copeland, Jan 24 2020
From Peter Luschny, Apr 30 2021: (Start)
Sum_{k=0..n} (-1)^k*T(n, k) = A122045(n).
Sum_{k=0..n} 2^(n-k)*T(n,k) = A007047(n).
Sum_{k=0..n} T(n, n-k) = A000522(n).
Sum_{k=0..n} T(n-k, k) = Sum_{k=0..n} (n - k)^k = A026898(n-1) for n >= 1.
Sum_{k=0..n} k*T(n, k) = A036919(n) = floor(n*n!*e/2).
(End)

More terms from Vladeta Jovovic, Sep 20 2003
First formula corrected by Wolfdieter Lang, Feb 14 2015
Offset set to 0 and edited by Peter Luschny, Apr 30 2021

A352525 Irregular triangle read by rows where T(n,k) is the number of integer compositions of n with k weak excedances (parts on or above the diagonal), all zeros removed.

Original entry on oeis.org

1, 1, 2, 3, 1, 5, 3, 8, 8, 14, 17, 1, 25, 35, 4, 46, 70, 12, 87, 137, 32, 167, 268, 76, 1, 324, 525, 170, 5, 634, 1030, 367, 17, 1248, 2026, 773, 49, 2466, 3999, 1598, 129, 4887, 7914, 3267, 315, 1, 9706, 15695, 6631, 730, 6, 19308, 31181, 13393, 1631, 23

Offset: 0

Gus Wiseman, Mar 22 2022

			Triangle begins:
     1
     1
     2
     3     1
     5     3
     8     8
    14    17     1
    25    35     4
    46    70    12
    87   137    32
   167   268    76     1
   324   525   170     5
For example, row n = 6 counts the following compositions:
  (6)       (15)     (123)
  (51)      (24)
  (312)     (33)
  (411)     (42)
  (1113)    (114)
  (1122)    (132)
  (2112)    (141)
  (2121)    (213)
  (3111)    (222)
  (11112)   (231)
  (11121)   (321)
  (11211)   (1131)
  (21111)   (1212)
  (111111)  (1221)
            (1311)
            (2211)
            (12111)

Andrew Howroyd, Table of n, a(n) for n = 0..2489 (rows 0..200)
MathOverflow, Why 'excedances' of permutations? [closed].

Row sums are A011782.
The version for partitions is A115994.
The version for permutations is A123125, strong A173018.
Column k = 1 is A177510.
The corresponding rank statistic is A352517.
The strong opposite is A352521, first col A219282, rank statistic A352514.
The opposite version is A352522, first col A238874, rank statistic A352515.
The strong version is A352524, first column A008930, rank statistic A352516.
A008292 is the triangle of Eulerian numbers (version without zeros).
A238349 counts comps by fixed points, first col A238351, rank stat A352512.
A352489 lists the weak excedance set of A122111.
A352523 counts comps by unfixed points, first A352520, rank stat A352513.
Cf. A088218, A098825, A114088, A238352, A319005, A350839, A352488.

pdw[y_]:=Length[Select[Range[Length[y]],#<=y[[#]]&]];
DeleteCases[Table[Length[Select[Join@@ Permutations/@IntegerPartitions[n],pdw[#]==k&]],{n,0,10},{k,0,n}],0,{2}]

T(n)={my(v=vector(n+1, i, i==1), r=v); for(k=1, n, v=vector(#v, j, sum(i=1, j-1, if(k<=i,x,1)*v[j-i])); r+=v); r[1]=x; [Vecrev(p) | p<-r/x]}
{ my(A=T(10)); for(i=1, #A, print(A[i])) } \\ Andrew Howroyd, Jan 19 2023

A090582 T(n, k) = Sum_{j=0..n-k} (-1)^jbinomial(n - k + 1, j)(n - k + 1 - j)^n. Triangle read by rows, T(n, k) for 1 <= k <= n.

Original entry on oeis.org

1, 2, 1, 6, 6, 1, 24, 36, 14, 1, 120, 240, 150, 30, 1, 720, 1800, 1560, 540, 62, 1, 5040, 15120, 16800, 8400, 1806, 126, 1, 40320, 141120, 191520, 126000, 40824, 5796, 254, 1, 362880, 1451520, 2328480, 1905120, 834120, 186480, 18150, 510, 1, 3628800, 16329600, 30240000, 29635200, 16435440, 5103000, 818520, 55980, 1022, 1

Offset: 1

Hugo Pfoertner, Jan 11 2004

Let Q(m, n) = Sum_(k=0..n-1) (-1)^k * binomial(n, k) * (n-k)^m. Then Q(m,n) is the numerator of the probability P(m,n) = Q(m,n)/n^m of seeing each card at least once if m >= n cards are drawn with replacement from a deck of n cards, written in a two-dimensional array read by antidiagonals with Q(m,m) as first row and Q(m,1) as first column.
The sequence is given as a matrix with the first row containing the cases #draws = size_of_deck. The second row contains #draws = 1 + size_of_deck. If "mn" indicates m cards drawn from a deck with n cards then the locations in the matrix are:
  11 22 33 44 55 66 77 ...
  21 32 43 54 65 76 87 ...
  31 42 53 64 75 86 97 ...
  41 52 63 74 85 .. .. ...
read by antidiagonals ->:
  11, 22, 21, 33, 32, 31, 44, 43, 42, 41, 55, 54, 53, 52, ....
The probabilities are given by Q(m,n)/n^m:
.(m,n):.....11..22..21..33..32..31..44..43..42..41...55...54..53..52..51
.....Q:......1...2...1...6...6...1..24..36..14...1..120..240.150..30...1
...n^m:......1...4...1..27...8...1.256..81..16...1.3125.1024.243..32...1
%.Success:.100..50.100..22..75.100...9..44..88.100....4...23..62..94.100
P(n,n) = n!/n^n which can be approximated by sqrt(Pi*(2n+1/3))/e^n (Gosper's approximation to n!).
Let P[n] be the set of all n-permutations. Build a superset Q[n] of P[n] composed of n-permutations in which some (possibly all or none) ascents have been designated. An ascent in a permutation s[1]s[2]...s[n] is a pair of consecutive elements s[i],s[i+1] such that s[i] < s[i+1]. As a triangular array read by rows T(n,k) is the number of elements in Q[n] that have exactly k distinguished ascents, n >= 1, 0 <= k <= n-1. Row sums are A000670. E.g.f. is y/(1+y-exp(y*x)). For example, T(3,1)=6 because there are four 3-permutations with one ascent, with these we would also count 1->2 3, and 1 2->3 where exactly one ascent is designated by "->". (After Flajolet and Sedgewick.) - Geoffrey Critzer, Nov 13 2012
Sum_(k=1..n) Q(n, k)*binomial(x, k) = x^n such that Sum_{k=1..n} Q(n, i)*binomial(x+1,i+1) = Sum_{k=1..x} k^n. - David A. Corneth, Feb 17 2014
A141618(n,n-k+1) = a(n,k) * C(n,k-1) / k. - Tom Copeland, Oct 25 2014
See A074909 and above g.f.s below for associations among this array and the Bernoulli polynomials and their umbral compositional inverses. - Tom Copeland, Nov 14 2014
For connections to toric varieties and Eulerian polynomials (in addition to those noted below), see the Dolgachev and Lunts and the Stembridge links in A019538. - Tom Copeland, Dec 31 2015
See A008279 for a relation between the e.g.f.s enumerating the faces of permutahedra (this entry) and stellahedra. - Tom Copeland, Nov 14 2016
From the Hasan and Franco and Hasan papers: The m-permutohedra for m=1,2,3,4 are the line segment, hexagon, truncated octahedron and omnitruncated 5-cell. The first three are well-known from the study of elliptic models, brane tilings and brane brick models. The m+1 torus can be tiled by a single (m+2)-permutohedron. Relations to toric Calabi-Yau Kahler manifolds are also discussed. - Tom Copeland, May 14 2020

			For m = 4, n = 2, we draw 4 times from a deck of two cards. Call the cards "a" and "b" - of the 16 possible combinations of draws (each of which is equally likely to occur), only two do not contain both a and b: a, a, a, a and b, b, b, b. So the probability of seeing both a and b is 14/16. Therefore Q(m, n) = 14.
Table starts:
  [1] 1;
  [2] 2,      1;
  [3] 6,      6,       1;
  [4] 24,     36,      14,      1;
  [5] 120,    240,     150,     30,      1;
  [6] 720,    1800,    1560,    540,     62,     1;
  [7] 5040,   15120,   16800,   8400,    1806,   126,    1;
  [8] 40320,  141120,  191520,  126000,  40824,  5796,   254,   1;
  [9] 362880, 1451520, 2328480, 1905120, 834120, 186480, 18150, 510, 1.

Reinhard Zumkeller, Rows n = 1..125 of triangle, flattened
Paul Barry, Three Études on a sequence transformation pipeline, arXiv:1803.06408 [math.CO], 2018.
Kevin Buhr, Michael Press and DMB, Collecting a deck of cards on the street. Thread in NG sci.math.
José L. Cereceda, Figurate numbers and sums of powers of integers, arXiv:2001.03208 [math.NT], 2020. See triangle Table 1 p. 5.
Tom Copeland, Generators, Inversion, and Matrix, Binomial, and Integral Transforms
P. Flajolet and R. Sedgewick, Analytic Combinatorics, 2009, page 209.
S. Franco and A. Hasan, Graded Quivers, Generalized Dimer Models and Toric Geometry, arXiv preprint arXiv:1904.07954 [hep-th], 2019.
A. Hasan, Physics and Mathematics of Graded Quivers, dissertation, Graduate Center, City University of New York, 2019.
Hugo Pfoertner, "Hit all" probabilities: Program and results.

Cf. A073593 first m >= n giving at least 50% probability, A085813 ditto for 95%, A055775 n^n/n!, A090583 Gosper's approximation to n!.
Reflected version of A019538.
Cf. A233734 (central terms).
Cf. A074909, A008292, A028246.
Cf. A046802, A119879, A123125, A130850, and A248727.

a090582 n k = a090582_tabl !! (n-1) !! (k-1)
a090582_row n = a090582_tabl !! (n-1)
a090582_tabl = map reverse a019538_tabl
-- Reinhard Zumkeller, Dec 15 2013

T := (n, k) -> add((-1)^j*binomial(n - k + 1, j)*(n - k + 1 - j)^n, j = 0..n-k):
# Or:
T := (n, k) -> (n - k + 1)!*Stirling2(n, n - k + 1):
for n from 1 to 9 do seq( T(n, k), k = 1..n) od; # Peter Luschny, May 21 2021

In[1]:= Table[Table[k! StirlingS2[n, k], {k, n, 1, -1}], {n, 1, 6}] (* Victor Adamchik, Oct 05 2005 *)
nn=6; a=y/(1+y-Exp[y x]); Range[0,nn]! CoefficientList[Series[a, {x,0,nn}], {x,y}]//Grid (* Geoffrey Critzer, Nov 10 2012 *)

a(n)={m=ceil((-1+sqrt(1+8*n))/2);k=m+1+binomial(m,2)-n;k*sum(i=1,k,(-1)^(i+k)*i^(m-1)*binomial(k-1,i-1))} \\ David A. Corneth, Feb 17 2014

T(n, k) = (n - k + 1)!*Stirling2(n, n - k + 1), generated by Stirling numbers of the second kind multiplied by a factorial. - Victor Adamchik, Oct 05 2005
Triangle T(n,k), 1 <= k <= n, read by rows given by [1, 1, 2, 2, 3, 3, 4, 4, 5, 5, ...] DELTA [0, 1, 0, 2, 0, 3, 0, 4, 0, 5, 0, 6, 0, ...] where DELTA is the operator defined in A084938. - Philippe Deléham, Nov 10 2006
From Tom Copeland, Oct 07 2008: (Start)
G(x,t) = 1/ (1 + (1-exp(x*t))/t) = 1 + 1*x + (2 + t)*x^2/2! + (6 + 6*t + t^2)*x^3/3! + ... gives row polynomials of A090582, the f-polynomials for the permutohedra (see A019538).
G(x,t-1) = 1 + 1*x + (1 + t)*x^2/2! + (1 + 4*t + t^2)*x^3/3! + ... gives row polynomials for A008292, the h-polynomials of permutohedra.
G[(t+1)x,-1/(t+1)] = 1 + (1 + t)*x + (1 + 3*t + 2*t^2)*x^2/2! + ... gives row polynomials of A028246. (End)
From Tom Copeland, Oct 11 2011: (Start)
With e.g.f. A(x,t) = G(x,t) - 1, the compositional inverse in x is
B(x,t) = log((t+1)-t/(1+x))/t. Let h(x,t) = 1/(dB/dx) = (1+x)*(1+(1+t)x), then the row polynomial P(n,t) is given by (1/n!)*(h(x,t)*d/dx)^n x, evaluated at x=0, A=exp(x*h(y,t)*d/dy) y, eval. at y=0, and dA/dx = h(A(x,t),t). (End)
k <= 0 or k > n yields Q(n, k) = 0; Q(1,1) = 1; Q(n, k) = k * (Q(n-1, k) + Q(n-1, k-1)). - David A. Corneth, Feb 17 2014
T = A008292*A007318. - Tom Copeland, Nov 13 2016
With all offsets 0 for this entry, let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of A123125 with offsets -1 so that the array becomes A008292; i.e., we ignore the first row and first column of A123125. Then the row polynomials of this entry, A090582, are given by A_n(1 + x;0). Other specializations of A_n(x;y) give A028246, A046802, A119879, A130850, and A248727. - Tom Copeland, Jan 24 2020

A352522 Triangle read by rows where T(n,k) is the number of integer compositions of n with k weak nonexcedances (parts on or below the diagonal).

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 3, 1, 2, 1, 2, 3, 4, 3, 3, 1, 3, 4, 8, 6, 6, 4, 1, 4, 7, 12, 13, 12, 10, 5, 1, 5, 13, 16, 26, 24, 22, 15, 6, 1, 7, 19, 27, 43, 48, 46, 37, 21, 7, 1, 10, 26, 47, 68, 90, 93, 83, 58, 28, 8, 1, 14, 36, 77, 109, 159, 180, 176, 141

Offset: 0

Gus Wiseman, Mar 22 2022

			Triangle begins:
   1
   0   1
   1   0   1
   1   1   1   1
   1   3   1   2   1
   2   3   4   3   3   1
   3   4   8   6   6   4   1
   4   7  12  13  12  10   5   1
   5  13  16  26  24  22  15   6   1
   7  19  27  43  48  46  37  21   7   1
  10  26  47  68  90  93  83  58  28   8   1
For example, row n = 6 counts the following compositions:
  (6)   (15)   (114)  (123)   (1113)   (11112)  (111111)
  (24)  (42)   (132)  (1311)  (1122)   (11121)
  (33)  (51)   (141)  (2112)  (1131)   (11211)
        (231)  (213)  (2121)  (1212)   (12111)
               (222)  (2211)  (1221)
               (312)  (3111)  (21111)
               (321)
               (411)

Andrew Howroyd, Table of n, a(n) for n = 0..1325 (rows 0..50)
MathOverflow, Why 'excedances' of permutations? [closed].

Row sums are A011782.
The strong version for partitions is A114088.
The opposite version for partitions is A115994.
The version for permutations is A123125, strong A173018.
Column k = 0 is A238874.
The corresponding rank statistic is A352515.
The strong version is A352521, first column A219282, rank statistic A352514.
The strong opposite is A352524, first col A008930, rank statistic A352516.
The opposite version is A352525, first col A177510, rank statistic A352517.
A000041 counts integer partitions, strict A000009.
A008292 is the triangle of Eulerian numbers (version without zeros).
A238349 counts comps by fixed points, first col A238351, rank stat A352512.
A352488 lists the weak nonexcedance set of A122111.
A352523 counts comps by unfixed points, first A352520, rank stat A352513.
Cf. A088218, A098825, A238352, A352489.

pw[y_]:=Length[Select[Range[Length[y]],#>=y[[#]]&]];
Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],pw[#]==k&]],{n,0,15},{k,0,n}]

T(n)={my(v=vector(n+1, i, i==1), r=v); for(k=1, n, v=vector(#v, j, sum(i=1, j-1, if(k>=i,x,1)*v[j-i])); r+=v); [Vecrev(p) | p<-r]}
{ my(A=T(10)); for(i=1, #A, print(A[i])) } \\ Andrew Howroyd, Jan 19 2023

A119879 Exponential Riordan array (sech(x),x).

Original entry on oeis.org

1, 0, 1, -1, 0, 1, 0, -3, 0, 1, 5, 0, -6, 0, 1, 0, 25, 0, -10, 0, 1, -61, 0, 75, 0, -15, 0, 1, 0, -427, 0, 175, 0, -21, 0, 1, 1385, 0, -1708, 0, 350, 0, -28, 0, 1, 0, 12465, 0, -5124, 0, 630, 0, -36, 0, 1, -50521, 0, 62325, 0, -12810, 0, 1050, 0, -45, 0, 1

Offset: 0

Paul Barry, May 26 2006

Row sums have e.g.f. exp(x)*sech(x) (signed version of A009006). Inverse of masked Pascal triangle A119467. Transforms the sequence with e.g.f. g(x) to the sequence with e.g.f. g(x)*sech(x).
Coefficients of the Swiss-Knife polynomials for the computation of Euler, tangent and Bernoulli number (triangle read by rows). Another version in A153641. - Philippe Deléham, Oct 26 2013
Relations to Green functions and raising/creation and lowering/annihilation/destruction operators are presented in Hodges and Sukumar and in Copeland's discussion of this sequence and 2020 pdf. - Tom Copeland, Jul 24 2020

			Triangle begins:
     1;
     0,    1;
    -1,    0,     1;
     0,   -3,     0,   1;
     5,    0,    -6,   0,   1;
     0,   25,     0, -10,   0,   1;
   -61,    0,    75,   0, -15,   0,   1;
     0, -427,     0, 175,   0, -21,   0,  1;
  1385,    0, -1708,   0, 350,   0, -28,  0,  1;

G. C. Greubel, Rows n=0..100 of triangle, flattened
Paul Barry, Riordan Arrays, Orthogonal Polynomials as Moments, and Hankel Transforms, J. Int. Seq. 14 (2011) # 11.2.2, example 28.
Tom Copeland, The Elliptic Lie Triad: Ricatti and KdV Equations, Infinigens, and Elliptic Genera, 2015.
Tom Copeland, Discussion of this sequence
Tom Copeland, SKipping over Dimensions, Juggling Zeros in the Matrix, 2020.
A. Hodges and C. V. Sukumar, Bernoulli, Euler, permutations and quantum algebras, Proc. R. Soc. A Oct. 2007 vol 463 no. 463 2086 2401-2414.
Peter Luschny, Additive decompositions of classical numbers via the Swiss-Knife polynomials [_Tom Copeland_, Oct 20 2015]
Miguel Méndez and Rafael Sánchez, On the combinatorics of Riordan arrays and Sheffer polynomials: monoids, operads and monops, arXiv:1707.00336 [math.CO], 2017, Section 4.3, Example 4.
Miguel A. Méndez and Rafael Sánchez Lamoneda, Monops, Monoids and Operads: The Combinatorics of Sheffer Polynomials, The Electronic Journal of Combinatorics 25(3) (2018), #P3.25.

Row sums are A155585. - Johannes W. Meijer, Apr 20 2011
Rows reversed: A081658.
Cf. A109449, A153641, A162660. - Philippe Deléham, Oct 26 2013
Cf. A000182, A046802, A119467, A133314.
Cf. A019538, A028246, A090582, A123125, A130850, A248727.

T := (n,k) -> binomial(n,k)*2^(n-k)*euler(n-k,1/2): # Peter Luschny, Jan 25 2009

T[n_, k_] := Binomial[n, k] 2^(n-k) EulerE[n-k, 1/2];
Table[T[n, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 20 2018, after Peter Luschny *)

{T(n,k) = binomial(n,k)*2^(n-k)*(2/(n-k+1))*(subst(bernpol(n-k+1, x), x, 1/2) - 2^(n-k+1)*subst(bernpol(n-k+1, x), x, 1/4))};
for(n=0,5, for(k=0,n, print1(T(n,k), ", "))) \\ G. C. Greubel, Feb 25 2019

@CachedFunction
def A119879_poly(n, x) :
    return 1 if n == 0  else add(A119879_poly(k, 0)*binomial(n, k)*(x^(n-k)-1+n%2) for k in range(n)[::2])
def A119879_row(n) :
    R = PolynomialRing(ZZ, 'x')
    return R(A119879_poly(n,x)).coeffs()  # Peter Luschny, Jul 16 2012
# Alternatively:

# uses[riordan_array from A256893]
riordan_array(sech(x), x, 9, exp=true) # Peter Luschny, Apr 19 2015

Number triangle whose k-th column has e.g.f. sech(x)*x^k/k!.
T(n,k) = C(n,k)*2^(n-k)*E_{n-k}(1/2) where C(n,k) is the binomial coefficient and E_{m}(x) are the Euler polynomials. - Peter Luschny, Jan 25 2009
The coefficients in ascending order of x^i of the polynomials p{0}(x) = 1 and p{n}(x) = Sum_{k=0..n-1; k even} binomial(n,k)*p{k}(0)*((n mod 2) - 1 + x^(n-k)). - Peter Luschny, Jul 16 2012
E.g.f.: exp(x*z)/cosh(x). - Peter Luschny, Aug 01 2012
Sum_{k=0..n} T(n,k)*x^k = A122045(n), A155585(n), A119880(n), A119881(n) for x = 0, 1, 2, 3 respectively. - Philippe Deléham, Oct 27 2013
With all offsets 0, let A_n(x;y) = (y + E.(x))^n, an Appell sequence in y where E.(x)^k = E_k(x) are the Eulerian polynomials of A123125. Then the row polynomials of A046802 (the h-polynomials of the stellahedra) are given by h_n(x) = A_n(x;1); the row polynomials of A248727 (the face polynomials of the stellahedra), by f_n(x) = A_n(1 + x;1); the Swiss-knife polynomials of this entry, A119879, by Sw_n(x) = A_n(-1;1 + x); and the row polynomials of the Worpitsky triangle (A130850), by w_n(x) = A(1 + x;0). Other specializations of A_n(x;y) give A090582 (the f-polynomials of the permutohedra, cf. also A019538) and A028246 (another version of the Worpitsky triangle). - Tom Copeland, Jan 24 2020
Triangle equals P*((I + P^2)/2)^(-1), where P denotes Pascal's triangle A007318. - Peter Bala, Mar 07 2024

cp's OEIS Frontend

A155585 a(n) = 2^n*E(n, 1) where E(n, x) are the Euler polynomials.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

Python

Sage

Sage

Formula

Extensions

A114088 Triangle read by rows: T(n,k) is number of partitions of n whose tail below its Durfee square has k parts (n >= 1; 0 <= k <= n-1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A002662 a(n) = 2^n - 1 - n*(n+1)/2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Formula

A352523 Number of integer compositions of n with exactly k nonfixed points (parts not on the diagonal).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A009998 Triangle in which j-th entry in i-th row is (j+1)^(i-j).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

A090582 T(n, k) = Sum_{j=0..n-k} (-1)^jbinomial(n - k + 1, j)(n - k + 1 - j)^n. Triangle read by rows, T(n, k) for 1 <= k <= n.