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

A052150 Partial sums of A000340, second partial sums of A003462.

Original entry on oeis.org

1, 6, 24, 82, 261, 804, 2440, 7356, 22113, 66394, 199248, 597822, 1793557, 5380776, 16142448, 48427480, 145282593, 435847950, 1307544040, 3922632330, 11767897221, 35303691916, 105911076024, 317733228372, 953199685441
Offset: 0

Views

Author

Barry E. Williams, Jan 23 2000

Keywords

References

  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, pp. 189, 194-196.
  • P. Ribenhoim, The Little Book of Big Primes, Springer-Verlag, N.Y., 1991, p. 53.

Crossrefs

Programs

  • Mathematica
    LinearRecurrence[{6,-12,10,-3},{1,6,24,82},40] (* Harvey P. Dale, Sep 05 2013 *)

Formula

a(n) = ((3^(n+3)) - (2*(n^2) + 12n + 19))/8.
a(n) = 3a(n-1)+C(n+2,2); a(0)=1.
a(n) = sum{k=0..n, binomial(n+3, k+3)2^k}. - Paul Barry, Aug 20 2004
From Colin Barker, Dec 18 2012: (Start)
a(n) = 6*a(n-1) - 12*a(n-2) + 10*a(n-3) - 3*a(n-4).
G.f.: 1/((x-1)^3*(3*x-1)). (End)

A014824 a(0) = 0; for n>0, a(n) = 10*a(n-1) + n.

Original entry on oeis.org

0, 1, 12, 123, 1234, 12345, 123456, 1234567, 12345678, 123456789, 1234567900, 12345679011, 123456790122, 1234567901233, 12345679012344, 123456790123455, 1234567901234566, 12345679012345677, 123456790123456788, 1234567901234567899, 12345679012345679010, 123456790123456790121
Offset: 0

Views

Author

Keywords

Comments

The square roots of these numbers have some remarkable properties - see the link to Schizophrenic numbers.
Partial sums of A002275. - Jonathan Vos Post, Apr 25 2010
This sequence is the particular case of a(0) = 0, a(n) = r*a(n-1) + n, when r = 10. If now the first N terms are computed for (r > N) then the resulting set of numbers is readable as the smallest k-digits permutations (1 <= k <= N): Those built from the concatenation of the first k digits in base-r (see links). - R. J. Cano, Jan 09 2013
There is also an interesting structure to the decimal expansion of 1/sqrt(a(2*n+1)), which has long strings of 0's (that gradually shorten in length until they disappear) interspersed with strings of what, at first sight, appear to be random digits. However, if we factorize these blocks of 'random' digits we find they are related to each other. An example illustrating this is given below. - Peter Bala, Sep 13 2015
From Peter Bala, Sep 15 2015: (Start)
Extending the previous empirical observation, it appears that the decimal expansions of the numbers 1/ (a(4*n+1))^(1/2), 1/a((4*n+3))^(1/4), 1/(10*a(4*n))^(1/2), 1/(10*a(4*n+2))^(1/4), and their powers, have the same pattern of beginning with long strings of 0's (that gradually shorten in length until they disappear) interlaced with strings of digits which, when read as integers and factorized, turn out to be related to each other.
The following result, which is a consequence of Bottomley's explicit formula for a(n), should be helpful in explaining these observations: the decimal expansion of 1/a(2*n-1) for n >= 5 begins with long strings of 0's interlaced successively with the digits of the numbers 81*(18*n + 1)^k for k going from 0 up to approximately n/log_10(18*n). For example, for n = 7 we have 1/a(13) = 0.00000000000081000000000102870000000130644900000 165919023..., with 10287 = 81*127, 1306449 = 81*127^2 and 165919023 = 81*127^3. A similar result holds for the decimal expansion of the number 1/a(2*n).
It appears that a(4*n+3)^(1/4), a(4*n+3)^(3/4), (10*a(4*n))^(1/2), (10*a(4*n+2))^(1/4) and (10*a(4*n+2))^(3/4) are further examples of Brown's schizophrenic numbers.
A theorem of Kuzmin in the measure theory of continued fractions says that for a random real number alpha, the probability that some given partial quotient of alpha is equal to a positive integer k is given by 1/log(2)*( log(1 + 1/k) - log(1 + 1/(k+1)) ). Thus large partial quotients are the exception in continued fraction expansions. Empirically, we observe the presence of unexpectedly large partial quotients early in the continued fraction expansions of the numbers (a(4*n+1))^(1/2), (a(4*n+3))^(1/4), (10*a(4*n))^(1/2), (10*a(4*n+2))^(1/4) and their powers. An example is given below. (End)
From Ya-Ping Lu, Dec 21 2024: (Start)
To get a(n), concatenate the first n digits in the cyclic string '123456790' and subtract the number of occurrences of '9' from the concatenated number. For example, a(8) = 12345679 - 1 = 12345678.
There are 2 prime terms for n <= 20000: a(2497) and a(3301). (End)

Examples

			From _Peter Bala_, Sep 13 2015: (Start)
The decimal expansion of 1/sqrt(a(51)) begins 9.0...0211050...07423683750...02901423065625000... x 10^(-26). The long strings of 0's gradually shorten in length and are interspersed with eleven blocks of digits  [9, 21105, 742368375, 2901423065625, 1190671490555859375, 5025824361636282421875, 216068565680679841787109375, 940978603539360710982861328125, 4137365297437126626102768402099609375, 18326229731370116994398540261077880859375, 816525165681195562685426961332324981689453125]. Read as ordinary integers these numbers factorize as [3^2, (3^2)*5*7*67, (3^3)*(5^3)*(7^2)*(67^2), (3^2)*(5^5)*(7^3)*(67^3), (3^2)*(5^8)*(7^5)*(67^4), (3^4)*(5^8)*(7^6)*(67^5), (3^3)*(5^10)*(7^7)*11*(67^6), (3^3)*(5^11)*(7^7)*11*13*(67^7), (3^4)*(5^16)*(7^8)*11*13*(67^8), (3^2)*(5^17)*(7^9)*11*13*17*(67^9), (3^2)*(5^18)*(7^10)*11*13*17*19*(67^10)]. (End)
From _Peter Bala_, Sep 15 2015: (Start)
The continued fraction expansion of 1/sqrt(a(51)) begins [0; 11111111111111111111111111, 9, 47382136934375740345889, 2, 21, 3, 1, 7, 2, 1, 101028010521057015662, 5, 14, 9, 1, 1, 2, 2, 8, 5, 1, 1, 1, 1, 215411536292232442, 5, 1, 5, 1, 1, 2, 1, 1, 8, 1, 4, 3, 1, 4, 2, 1, 8, 1, 1, 3, 10, 459299650942926, 4, 1, 1, 4, 1, 20, 64, 5, 9, 2, 2, 1, 2, 1, 1, 1, 1, 30, 1, 11, 3, 979316952969, 1, 2, 93, 1, 5, 1, 1, 11, 1, 1, 1, 1, 5, 1, 29, 1, 29, 1, 1, 1, 2, 4, 1, 37, 1, 1, 2, 8, 2, 2088095848, 12, 1, 3, 1, 3, 2, 2, 3, 1, 5, 6, 1, 3, 1, 4, 2, 2, 1, 2, 2, 14, 4, 1, 2, 1, 50, 2, 6, 1, 11, 135, 4452229, 1, ...] and has several unexpectedly large partial quotients early on. (End)
For n=5, a(5) = 1*15 + 9*20 + 9^2*15 + 9^3*6 + 9^4*1 + 9^5*0 = 12345. - _Bruno Berselli_, Nov 13 2015
		

Crossrefs

Cf. A060011.
Cf. A002275. - Jonathan Vos Post, Apr 25 2010
Similar sequences in other bases are: (base-2) A000295, (base-3) A000340, (base-4) A014825, (base-5) A014827, (base-6) A014829. - R. J. Cano, Jan 11 2013
Differs from A007908, A035239, A057137, A060555, A138957 from n=10 on. - M. F. Hasler, Jan 17 2013
Cf. A030512.

Programs

  • Magma
    [(10^n-1)*(10/81)-n/9: n in [0..20]]; // Vincenzo Librandi, Aug 23 2011
    
  • Maple
    a:=n->sum((10^(n-j)-1^(n-j))/9,j=0..n): seq(a(n), n=0..17); # Zerinvary Lajos, Jan 15 2007
    a:=n->sum(10^(n-j)*j,j=0..n): seq(a(n), n=0..16); # Zerinvary Lajos, Jun 05 2008
  • Mathematica
    Table[Sum[10^i - 1, {i, n}]/9, {n, 18}] (* Robert G. Wilson v, Nov 20 2004 *)
    CoefficientList[Series[x/(1 - 12*x + 21*x^2 - 10*x^3), {x, 0, 20}], x] (* Wesley Ivan Hurt, Sep 15 2015 *)
  • PARI
    linrec01(p,u,base)={my(r=!p,A=1);for(j=2,u,A=A*base+r+p*j); A};
    a(n)=(n!=0)*linrec01(1, n, 10); \\ R. J. Cano, Jan 09 2011; With (0, n, 10) it generates repunit numbers.
    
  • PARI
    A014824(n)=(10^(n+1)\9-n)\9  \\ M. F. Hasler, Jan 17 2013
    
  • Python
    def A014824(n): s = ''.join('123456790'[i%9] for i in range(n)); q, r = divmod(n, 9); return int(s) - q - r//8 # Ya-Ping Lu, Dec 21 2024

Formula

a(n) = (10^n-1)*(10/81) - n/9. - Henry Bottomley, Jul 04 2000
a(n)/10^n converges to 10/81 = 0.123456790123456790...
Let b(n) = if(n = 0, 1, if(n = 1, 10, 10*9^(n-2))). Then a(n) = Sum_{k=0..n} C(n, k)*b(k) (Binomial transform). - Paul Barry, Jan 29 2004
G.f.: x/(1-12*x+21*x^2-10*x^3). - Colin Barker, Jan 08 2012
a(n) = 12*a(n-1) - 21*a(n-2) + 10*a(n-3), n>2. - Wesley Ivan Hurt, Sep 15 2015
a(n) = Sum_{i=0..n} 9^i*binomial(n+1,n-1-i). - Bruno Berselli, Nov 13 2015
a(n) = Sum_{i=0..n} 10^(n-i)*i. - Ya-Ping Lu, Dec 21 2024
E.g.f.: exp(x)*(10*exp(9*x) - 9*x - 10)/81. - Elmo R. Oliveira, Mar 29 2025

A122803 Powers of -2: a(n) = (-2)^n.

Original entry on oeis.org

1, -2, 4, -8, 16, -32, 64, -128, 256, -512, 1024, -2048, 4096, -8192, 16384, -32768, 65536, -131072, 262144, -524288, 1048576, -2097152, 4194304, -8388608, 16777216, -33554432, 67108864, -134217728, 268435456, -536870912, 1073741824, -2147483648, 4294967296, -8589934592, 17179869184
Offset: 0

Views

Author

Keywords

Comments

The number -2 can be used as a base of numeration (see the Weisstein link). - Alonso del Arte, Mar 30 2014
Contribution from M. F. Hasler, Oct 21 2014: (Start)
This is the inverse binomial transform of A033999 = n->(-1)^n, and the binomial transform of A033999*A000244 = n->(-3)^n, see also A141413.
Prefixed with one 0, i.e., (0,1,-2,4,...) = -A033999*A131577, it is the binomial transform of (0, 1, -4, 13, -40, 121,...) = -A033999*A003462, and inverse binomial transform of (0,1,0,1,0,1,...) = A000035.
Prefixed with two 0's, i.e., (0,0,1,-2,4,-8,...), it is the binomial transform of (0,0,1,-5,18,-58,179,-543,...) (cf. A000340) and inverse binomial transform of (0,0,1,1,2,2,3,3,...) = A004526. (End)
Prefixed with three 0's, this is the inverse binomial difference of (0, 0, 0, 1, 2, 4, 6, 9, 12, 16,...) = concat(0, A002620), which has as successive differences (0, 0, 1, 1, 2, 2,...) = A004526, then (0, 1, 0, 1,...) = A000035, then (1, -1, 1, -1,...) = A033999, and then (-2)^k*A033999 with k=1,2,3,... - Paul Curtz, Oct 16 2014, edited by M. F. Hasler, Oct 21 2014
Stirling-Bernoulli transform of triangular numbers: 1, 3, 6, 10, 15, 21, 28, ... - Philippe Deléham, May 25 2015

Crossrefs

Programs

Formula

a(n) = (-2)^n = (-1)^n * 2^n.
a(n) = -2*a(n-1), n > 0; a(0) = 1. G.f.: 1/(1+2x). - Philippe Deléham, Nov 19 2008
Sum_{n >= 0} 1/a(n) = 2/3. - Jaume Oliver Lafont, Mar 01 2009
E.g.f.: 1/exp(2*x). - Arkadiusz Wesolowski, Aug 13 2012
a(n) = Sum_{k = 0..n} (-2)^(n-k)*binomial(n, k)*A030195(n+1). - R. J. Mathar, Oct 15 2012
G.f.: 1/(1+2x). A122803 = A033999 * A000079. - M. F. Hasler, Oct 21 2014
a(n) = Sum_{k = 0..n} A163626(n,k)*A000217(k+1). - Philippe Deléham, May 25 2015

A156919 Table of coefficients of polynomials related to the Dirichlet eta function.

Original entry on oeis.org

1, 2, 1, 4, 10, 1, 8, 60, 36, 1, 16, 296, 516, 116, 1, 32, 1328, 5168, 3508, 358, 1, 64, 5664, 42960, 64240, 21120, 1086, 1, 128, 23488, 320064, 900560, 660880, 118632, 3272, 1, 256, 95872, 2225728
Offset: 0

Views

Author

Johannes W. Meijer, Feb 20 2009, Jun 24 2009

Keywords

Comments

Essentially the same as A185411. Row reverse of A185410. - Peter Bala, Jul 24 2012
The SF(z; n) formulas, see below, were discovered while studying certain properties of the Dirichlet eta function.
From Peter Bala, Apr 03 2011: (Start)
Let D be the differential operator 2*x*d/dx. The row polynomials of this table come from repeated application of the operator D to the function g(x) = 1/sqrt(1 - x). For example,
D(g) = x*g^3
D^2(g) = x*(2 + x)*g^5
D^3(g) = x*(4 + 10*x + x^2)*g^7
D^4(g) = x*(8 + 60*x + 36*x^2 + x^3)*g^9.
Thus this triangle is analogous to the triangle of Eulerian numbers A008292, whose row polynomials come from the repeated application of the operator x*d/dx to the function 1/(1 - x). (End)

Examples

			The first few rows of the triangle are:
  [1]
  [2, 1]
  [4, 10, 1]
  [8, 60, 36, 1]
  [16, 296, 516, 116, 1]
The first few P(z;n) are:
  P(z; n=0) = 1
  P(z; n=1) = 2 + z
  P(z; n=2) = 4 + 10*z + z^2
  P(z; n=3) = 8 + 60*z + 36*z^2 + z^3
The first few SF(z;n) are:
  SF(z; n=0) = (1/2)*(1)/(1-z)^(3/2);
  SF(z; n=1) = (1/4)*(2+z)/(1-z)^(5/2);
  SF(z; n=2) = (1/8)*(4+10*z+z^2)/(1-z)^(7/2);
  SF(z; n=3) = (1/16)*(8+60*z+36*z^2+z^3)/(1-z)^(9/2);
In the Savage-Viswanathan paper, the coefficients appear as
  1;
  1,    2;
  1,   10,     4;
  1,   36,    60,     8;
  1,  116,   516,   296,    16;
  1,  358,  3508,  5168,  1328,   32;
  1, 1086, 21120, 64240, 42960, 5664, 64;
  ...
		

Crossrefs

A142963 and this sequence can be mapped onto the A156920 triangle.
FP1 sequences A000340, A156922, A156923, A156924.
FP2 sequences A050488, A142965, A142966, A142968.
Appears in A162005, A000182, A162006 and A162007.
Cf. A185410 (row reverse), A185411.

Programs

  • Maple
    A156919 := proc(n,m) if n=m then 1; elif m=0 then 2^n ; elif m<0 or m>n then 0; else 2*(m+1)*procname(n-1,m)+(2*n-2*m+1)*procname(n-1,m-1) ; end if; end proc: seq(seq(A156919(n,m), m=0..n), n=0..7); # R. J. Mathar, Feb 03 2011
  • Mathematica
    g[0] = 1/Sqrt[1-x]; g[n_] := g[n] = 2x*D[g[n-1], x]; p[n_] := g[n] / g[0]^(2n+1) // Cancel; row[n_] := CoefficientList[p[n], x] // Rest; Table[row[n], {n, 0, 9}] // Flatten (* Jean-François Alcover, Aug 09 2012, after Peter Bala *)
    Flatten[Table[Rest[CoefficientList[Nest[2 x D[#, x] &, (1 - x)^(-1/2), k] (1 - x)^(k + 1/2), x]], {k, 10}]] (* Jan Mangaldan, Mar 15 2013 *)

Formula

SF(z; n) = Sum_{m >= 1} m^(n-1)*4^(-m)*z^(m-1)*Gamma(2*m+1)/(Gamma(m)^2) = P(z;n) / (2^(n+1)*(1-z)^((2*n+3)/2)) for n >= 0. The polynomials P(z;n) = Sum_{k = 0..n} a(k)*z^k generate the a(n) sequence.
If we write the sequence as a triangle the following relation holds: T(n,m) = (2*m+2)*T(n-1,m) + (2*n-2*m+1)*T(n-1,m-1) with T(n,m=0) = 2^n and T(n,n) = 1, n >= 0 and 0 <= m <= n.
G.f.: 1/(1-xy-2x/(1-3xy/(1-4x/(1-5xy/(1-6x/(1-7xy/(1-8x/(1-... (continued fraction). - Paul Barry, Jan 26 2011
From Peter Bala, Apr 03 2011 (Start)
E.g.f.: exp(z*(x + 2)) * (1 - x)/(exp(2*x*z) - x*exp(2*z))^(3/2) = Sum_{n >= 0} P(x,n)*z^n/n! = 1 + (2 + x)*z + (4 + 10*x + x^2)*z^2/2! + (8 + 60*x + 36*x^2 + x^3)*z^3/3! + ... .
Explicit formula for the row polynomials:
P(x,n-1) = Sum_{k = 1..n} 2^(n-2*k)*binomial(2k,k)*k!*Stirling2(n,k)*x^(k-1)*(1 - x)^(n-k).
The polynomials x*(1 + x)^n * P(x/(x + 1),n) are the row polynomials of A187075.
The polynomials x^(n+1) * P((x + 1)/x,n) are the row polynomials of A186695.
Row sums are A001147(n+1). (End)
Sum_{k = 0..n} (-1)^k*T(n,k) = (-1)^binomial(n,2)*A012259(n+1). - Johannes W. Meijer, Sep 27 2011

Extensions

Minor edits from Johannes W. Meijer, Sep 27 2011

A126885 T(n,k) = n*T(n,k-1) + k, with T(n,1) = 1, square array read by ascending antidiagonals (n >= 0, k >= 1).

Original entry on oeis.org

1, 1, 2, 1, 3, 3, 1, 4, 6, 4, 1, 5, 11, 10, 5, 1, 6, 18, 26, 15, 6, 1, 7, 27, 58, 57, 21, 7, 1, 8, 38, 112, 179, 120, 28, 8, 1, 9, 51, 194, 453, 543, 247, 36, 9, 1, 10, 66, 310, 975, 1818, 1636, 502, 45, 10, 1, 11, 83, 466, 1865, 4881, 7279, 4916, 1013, 55, 11
Offset: 0

Views

Author

Gary W. Adamson, Dec 30 2006

Keywords

Examples

			Square array begins:
  n\k | 1   2   3   4    5     6      7       8 ...
  -------------------------------------------------
    0 | 1   2   3   4    5     6      7       8 ... A000027
    1 | 1   3   6  10   15    21     28      36 ... A000217
    2 | 1   4  11  26   57   120    247     502 ... A000295
    3 | 1   5  18  58  179   543   1636    4916 ... A000340
    4 | 1   6  27 112  453  1818   7279   29124 ... A014825
    5 | 1   7  38 194  975  4881  24412  122068 ... A014827
    6 | 1   8  51 310 1865 11196  67183  403106 ... A014829
    7 | 1   9  66 466 3267 22875 160132 1120932 ... A014830
    8 | 1  10  83 668 5349 42798 342391 2739136 ... A014831
    ...
		

Crossrefs

Antidiagonal sums are A134195.
Main diagonal gives A062805.

Programs

  • Maxima
    T(n, k) := if k = 1 then 1 else n*T(n, k - 1) + k$
    create_list(T(n - k + 1, k), n, 0, 20, k, 1, n + 1);
    /* Franck Maminirina Ramaharo, Jan 26 2019 */

Formula

T(1,k) = k*(k + 1)/2, and T(n,k) = (k - (k + 1)*n + n^(k + 1))/(n^2 - 2*n + 1) elsewhere.
T(n,k) = third entry in the vector M^k * (1, 0, 0), where M is the following 3 X 3 matrix:
1, 0, 0
1, 1, 0
1, 1, n.

Extensions

Edited and name clarified by Franck Maminirina Ramaharo, Jan 26 2019

A156920 Triangle of the normalized A142963 and A156919 sequences.

Original entry on oeis.org

1, 1, 1, 1, 5, 1, 1, 15, 18, 1, 1, 37, 129, 58, 1, 1, 83, 646, 877, 179, 1, 1, 177, 2685, 8030, 5280, 543, 1, 1, 367, 10002, 56285, 82610, 29658, 1636, 1, 1, 749, 34777, 335162, 919615, 756218, 159742, 4916, 1
Offset: 0

Views

Author

Johannes W. Meijer, Feb 20 2009

Keywords

Comments

The originator sequences are A142963 and A156919.
The Flower Triangle seems to be an appropriate name for the triangular array of this sequence. The zero patterns of the Flower Polynomials of the first, see A156921, the second, see A156925, the third, see A156927, and the fourth kind, see A156933, look like flowers.
The first Maple program generates the Flower Triangle sequence.
The second program generates the Right Hand Columns sequences and the third one generates the Left Hand Column sequences. For an explanation of these two algorithms see A142963.

Examples

			The first few rows of the triangle are:
  [1]
  [1, 1]
  [1, 5, 1]
  [1, 15, 18, 1]
  [1, 37, 129, 58, 1]
  [1, 83, 646, 877, 179, 1]
		

Crossrefs

Originator sequences A142963, A156919.
Related sequences A156921, A156925, A156927, A156933.
Left hand column sequences A050488, A142965, A142966, A142968.
Right hand column sequences A000340, A156922, A156923, A156924.
Row sums A014307(n+1).

Programs

  • Maple
    A156920 := proc(n,m): if n=m then 1; elif m=0 then 1 ; elif m<0 or m>n then 0; else (m+1)*procname(n-1, m)+(2*n-2*m+1)*procname(n-1, m-1) ; end if; end proc: seq(seq(A156920(n, m), m=0..n), n=0..8);
    RHCnr:=5; RHCmax:=10; RHCend:=RHCnr+RHCmax: for k from RHCnr to RHCend do for n from 0 to k do S2[k,n]:=sum((-1)^(n+i)*binomial(n,i)*i^k/n!,i=0..n) end do: G(k,x):= sum(S2[k,p]*((2*p)!/p!) *x^p/(1-4*x)^(p+1),p=0..k)/(((-1)^(k+1)*2*x)/(-1+4*x)^(k+1)): fx:=simplify(G(k,x)): nmax:=degree(fx); RHC[k-RHCnr+1]:= coeff(fx,x,k-RHCnr)/2^(k-RHCnr) end do: a:=n-> RHC[n]: seq(a(n), n=1..RHCend-RHCnr);
    LHCnr:=5; LHCmax:=10: LHCend:=LHCnr+LHCmax: for k from LHCnr to LHCend do for n from 0 to k do S2[k,n]:=sum((-1)^(n+i)*binomial(n,i)*i^k/n!,i=0..n) end do: G(k,x):= sum(S2[k,p]*((2*p)!/p!)*x^p/(1-4*x)^(p+1),p=0..k)/ (((-1)^(k+1)*2*x)/(-1+4*x)^(k+1)): fx:=simplify(G(k,x)): nmax:=degree(fx); for n from 0 to nmax do d[n]:= coeff(fx,x,n)/2^n end do: LHC[n]:=d[LHCnr-1] end do: a:=n-> LHC[n]: seq(a(n), n=LHCnr..LHCend-1);
  • Mathematica
    T[, 0] = 1; T[n, n_] = 1; T[n_, m_] := T[n, m] = (m + 1)*T[n - 1, m] + (2*n - 2*m + 1)*T[n - 1, m - 1];
    Table[T[n, m], {n, 0, 8}, {m, 0, n}] // Flatten (* Jean-François Alcover, Nov 14 2017 *)

Formula

T(n,m) = (m+1)*T(n-1,m) + (2*n-2*m+1)*T(n-1,m-1) with T(n,m=0) = 1 and T(n,n) = 1, n>=0 and 0 <= m <= n.
From Peter Bala, Jul 22 2012: (Start)
T(n,k) = 1/(2^(n-k))*A156919(n,k).
E.g.f.: 1 + t*x + (t+t^2)*x^2/2! + (t+5*t^2+t^3)*x^3/3! + ... = sqrt(E(x,2*t)), where E(x,t) = (1-t)*exp(x*t)/(exp(x*t)-t*exp(x)) = 1 + t*x + (t+t^2)*x^2/2! + (t+4*t^2+t^3)*x^3/3! + ... is the e.g.f. for the Eulerian numbers A008292.
The row polynomials R(n,x) satisfy 1/sqrt(1-2*x)*(x*d/dx)^n(1/sqrt(1-2*x)) = R(n,x)/(1-2*x)^(n+1). (End)

Extensions

Minor edits by Johannes W. Meijer, Sep 28 2011

A142963 Triangle read by rows, coefficients of the polynomials P(k, x) = (1/2) Sum_{p=0..k-1} Stirling2(k, p+1)*x^p*(1-4*x)^(k-1-p)*(2*p+2)!/(p+1)!.

Original entry on oeis.org

1, 1, 2, 1, 10, 4, 1, 30, 72, 8, 1, 74, 516, 464, 16, 1, 166, 2584, 7016, 2864, 32, 1, 354, 10740, 64240, 84480, 17376, 64, 1, 734, 40008, 450280, 1321760, 949056, 104704, 128, 1, 1498, 139108, 2681296, 14713840, 24198976, 10223488, 629248, 256, 1, 3030, 462264, 14341992
Offset: 1

Views

Author

Wolfdieter Lang, Sep 15 2008

Keywords

Comments

Previous name: Table of coefficients of row polynomials of certain o.g.f.s.
The o.g.f.s G(k, x) for the k-family of sequences S(k, n):= Sum_{p=0..n} p^k*binomial(2*p, p)*binomial(2*(n-p), n-p), k=0,1,... (convolution of two sequences involving the central binomial coefficients) are 1/(1-4*x) for k=0 and 2*x*P(k, x)/(1-4*x)^(k+1) for k=1,2,..., with the row polynomials P(k, x) = Sum_{m=0..k-1} a(n,m)*x^m).
The author was led to compute the sums S(k, n) by a question asked by M. Greiter, Jun 27 2008.
In order to keep the index k>=1 of Sigma(k, n) also for the polynomials P(k, x), their degree is then k-1.

Examples

			Triangle starts:
[1]
[1,   2]
[1,  10,     4]
[1,  30,    72,      8]
[1,  74,   516,    464,      16]
[1, 166,  2584,   7016,    2864,     32]
[1, 354, 10740,  64240,   84480,  17376,     64]
[1, 734, 40008, 450280, 1321760, 949056, 104704, 128]
...
P(3,x) = 1+10*x+4*x^2.
G(3,x) = 2*x*(1+10*x+4*x^2)/(1-4*x)^4.
		

Crossrefs

Left hand column sequences 2*A142964, 4*A142965, 8*A142966, 16*A142968.
Row sums A142967.
From Johannes W. Meijer, Feb 20 2009: (Start)
A156919 and this sequence can be mapped onto A156920.
Right hand column sequences 2^n*A000340, 2^n*A156922, 2^n*A156923, 2^n*A156924. (End)

Programs

  • Maple
    A142963 := proc(n,m): if n=m+1 then 2^(n-1); elif m=0 then 1 ; elif m<0 or m>n-1 then 0; else (m+1)*procname(n-1, m)+(4*n-4*m-2)*procname(n-1, m-1); end if; end proc: seq(seq(A142963(n,m), m=0..n-1), n=1..9); # Johannes W. Meijer, Sep 28 2011
    # Alternatively (assumes offset 0):
    p := (n,x) -> (1/2)*add(Stirling2(n+1,k+1)*x^k*(1-4*x)^(n-k)*(2*k+2)!/(k+1)!, k=0..n): for n from 0 to 7 do [n], PolynomialTools:-CoefficientList(p(n,x), x) od;
    # Peter Luschny, Jun 18 2017
  • Mathematica
    t[, 0] = 1; t[n, m_] /; m == n-1 := 2^m; t[n_, m_] := (m+1)*t[n-1, m] + (4*n-4*m-2)*t[n-1, m-1]; Table[t[n, m], {n, 1, 10}, {m, 0, n-1}] // Flatten (* Jean-François Alcover, Jun 21 2013, after Johannes W. Meijer *)

Formula

G(k, x) = Sum_{p=0..k} S2(k, p)*((2*p)!/p!)*x^p/(1-4*x)^(p+1), k >= 0 (here k >= 1), with the Stirling2 triangle S2(k, p):=A048993(k, p). (Proof from the product of the o.g.f.s of the two convoluted sequences and the normal ordering (x^d_x)^k = Sum_{p=0..k} S2(k, p)*x^p*d_x^p, with the derivative operator d_x.)
a(k,m) = [x^m]P(k, x) = [x^m] ((1-4*x)^(k+1))*G(k,x)/(2*x), k>=1, m=0,1,...,k-1.
For the triangle coefficients the following relation holds: T(n,m) = (m+1)*T(n-1,m) + (4*n-4*m-2)*T(n-1,m-1) with T(n,m=0) = 1 and T(n,m=n-1) = 2^(n-1), n >= 1 and 0 <= m <= n-1. - Johannes W. Meijer, Feb 20 2009
From Peter Bala, Jan 18 2018: (Start)
(x*d/dx)^n (1/(sqrt(1 - 4*x)) = 2*x*P(n,x)/sqrt(1 - 4*x)^(n+1/2) for n >= 1.
x*P(n,x)/(1 - 4*x)^(n+1/2) = (1/2)*Sum_{k >= 1} binomial(2*k,k)* k^n*x^k for n >= 1.
P(n+1,x) = ((4*n - 2)*x + 1)*P(n,x) - x*(4*x - 1)*d/dx(P(n,x)).
Hence the polynomial P(n,x) has all real zeros by Liu et al., Theorem 1.1, Corollary 1.2. (End)

Extensions

Minor edits by Johannes W. Meijer, Sep 28 2011
A more precise name by Peter Luschny, Jun 18 2017
Name reformulated with offset corrected, edited by Wolfdieter Lang, Aug 23 2019

A156921 FP1 polynomials related to the generating functions of the right hand columns of the A156920 triangle.

Original entry on oeis.org

1, 1, 1, 1, -6, 1, 7, -79, 119, 126, -270, 1, 28, -515, 1654, 8689, -65864, 142371, -82242, -99090, 113400, 1, 86, -2255, 5784, 300930, -3904584, 20663714, -41517272, -80232259, 657717054
Offset: 0

Views

Author

Johannes W. Meijer, Feb 20 2009

Keywords

Comments

The FP1 polynomials appear in the numerators of the GF1 o.g.f.s. of the right hand columns of A156920. The FP1 can be calculated with the formula for the RHC sequence, see A156920, and the formula for the general structure of the generating function GF1, see below.
An appropriate name for the FP1 polynomials seems to be the flower polynomials of the first kind because the zero patterns of these polynomials look like flowers. The zero patterns of the FP2, see A156925, and the FP1 resemble each other closely.
A Maple program that generates for a right hand column with a certain RHCnr its GF1 and FP1 can be found below. RHCnr stands for right hand column number and starts from 1.

Examples

			The first few rows of the "triangle" of the coefficients of the FP1 polynomials.
In the columns the coefficients of the powers of z^m, m=0,1,2,... , appear.
  [1]
  [1]
  [1, 1, -6]
  [1, 7, -79, 119, 126, -270]
  [1, 28, -515, 1654, 8689, -65864, 142371, -82242, -99090, 113400]
Matrix of the coefficients of the FP1 polynomials. The coefficients in the columns of this matrix are the powers of z^m, m=0,1,2,.. .
  [1, 0 ,0, 0, 0, 0, 0, 0, 0, 0]
  [1, 0 ,0, 0, 0, 0, 0, 0, 0, 0]
  [1, 1, -6, 0 ,0, 0, 0, 0, 0, 0]
  [1, 7, -79, 119, 126, -270, 0, 0, 0, 0]
  [1, 28, -515, 1654, 8689, -65864, 142371, -82242, -99090, 113400]
The first few FP1 polynomials are:
  FP1(z; RHCnr=1) = 1
  FP1(z; RHCnr=2) = 1
  FP1(z; RHCnr =3) = 1+z-6*z^2
Some GF1(z;RHCnr) are:
  GF1(z;RHCnr= 3) = (1+z-6*z^2)/((1-5*z)*(1-3*z)^2*(1-z)^3)
  GF1(z;RHCnr= 4) = (1+7*z-79*z^2+119*z^3+126*z^4-270*z^5)/((1-7*z)*(1-5*z)^2*(1-3*z)^3*(1-z)^4)
		

Crossrefs

For the first few GF1's see A000340, A156922, A156923, A156924.
The number of FP1 terms follow the triangular numbers A000217, with quite surprisingly one exception here a(0)=1.
Abs(Row sums (n)) = A098695(n).
For the polynomials in the denominators of the GF1(z;RHCnr) see A157702.

Programs

  • Maple
    RHCnr:=4: if RHCnr=1 then RHCmax :=1; else RHCmax:=(RHCnr-1)*(RHCnr)/2 end if: RHCend:=RHCnr+RHCmax: for k from RHCnr to RHCend do for n from 0 to k do S2[k,n]:=sum((-1)^(n+i)*binomial(n,i)*i^k/n!,i=0..n) end do: G(k,x):= sum(S2[k,p]*((2*p)!/p!) *x^p/(1-4*x)^(p+1),p=0..k)/(((-1)^(k+1)*2*x)/(-1+4*x)^(k+1)): fx:=simplify(G(k,x)): nmax:=degree(fx); RHC[k-RHCnr+1]:= coeff(fx,x,k-RHCnr)/2^(k-RHCnr) end do: a:=n-> RHC[n]: seq(a(n), n=1..RHCend-RHCnr+1); for nx from 0 to RHCmax do num:=sort(sum(A[t]*z^t, t=0..RHCmax)); nom:=Product((1-(2*u-1)*z)^(RHCnr-u+1),u=1..RHCnr): RHCa:= series(num/nom,z,nx+1); y:=coeff(RHCa,z,nx)-A[nx]; x:=RHC[nx+1]; A[nx]:=x-y; end do: FP1[RHCnr]:=sort(num,z, ascending); GenFun[RHCnr] :=FP1[RHCnr]/product((1-(2*m-1)*z)^(RHCnr-m+1),m=1..RHCnr);

Formula

G.f.: GF1(z;RHCnr) := FP1(z;RHCnr)/product((1-(2*m-1)*z)^(RHCnr+1-m),m=1..RHCnr)
Row sums (n) = (-1)^(1+(n+1)*(n+2)/2)*A098695(n).

A106516 A Pascal-like triangle based on 3^n.

Original entry on oeis.org

1, 3, 1, 9, 4, 1, 27, 13, 5, 1, 81, 40, 18, 6, 1, 243, 121, 58, 24, 7, 1, 729, 364, 179, 82, 31, 8, 1, 2187, 1093, 543, 261, 113, 39, 9, 1, 6561, 3280, 1636, 804, 374, 152, 48, 10, 1, 19683, 9841, 4916, 2440, 1178, 526, 200, 58, 11, 1, 59049, 29524, 14757, 7356, 3618, 1704, 726, 258, 69, 12, 1
Offset: 0

Views

Author

Paul Barry, May 05 2005

Keywords

Comments

Row sums are A027649. Antidiagonal sums are A106517.
From Wolfdieter Lang, Jan 09 2015: (Start)
Alternating row sums give A025192. The A-sequence of this Riordan lower triangular matrix is [1, 1, repeat(0, )] (leading to the Pascal recurrence for T(n,k) for n >= k >= 1. The Z-sequence is [3, repeat(0, )] (leading to the recurrence T(n,0) = 3*T(n-1,0), n >= 1. For A- and Z-sequences see the W. Lang link under A006232.
The inverse of this Riordan matrix is Tinv = ((1 - 2*x)/(1 + x), x/(1 + x)) given as a signed version of A093560: Tinv(n,m) = (-1)^(n-m)*A093560(n,m). (End)

Examples

			The triangle T(n,k) begins:
n\k     0     1     2    3    4    5   6   7  8  9 10 ...
0:      1
1:      3     1
2:      9     4     1
3:     27    13     5    1
4:     81    40    18    6    1
5:    243   121    58   24    7    1
6:    729   364   179   82   31    8   1
7:   2187  1093   543  261  113   39   9   1
8:   6561  3280  1636  804  374  152  48  10  1
9:  19683  9841  4916 2440 1178  526 200  58 11  1
10: 59049 29524 14757 7356 3618 1704 726 258 69 12  1
... reformatted and extended. - _Wolfdieter Lang_, Jan 06 2015
----------------------------------------------------------
With the array M(k) as defined in the Formula section, the infinite product M(0)*M(1)*M(2)*... begins
/ 1        \/1           \/1        \       /1         \
| 3  1     ||0  1        ||0 1      |      | 3  1      |
| 9  4 1   ||0  3  1     ||0 0 1    |... = | 9  7  1   |
|27 13 5 1 ||0  9  4 1   ||0 0 3 1  |      |27 37 12 1 |
|...       ||0 27 13 5 1 ||0 0 9 4 1|      |...        |
|...       ||...         ||...      |      |...        |
= A143495. - _Peter Bala_, Dec 23 2014
		

Crossrefs

Columns 1, 2, 3, 4, 5: A003462, A000340, A052150, A097786, A097787.

Programs

  • Mathematica
    a106516[n_] := Block[{a, k},
    a[x_] := Flatten@ Last@ Reap[For[k = -1, k < x, Sow[Binomial[x, k] +
    2 Sum[3^(i - 1)*Binomial[x - i, k], {i, 1, x}]], k++]]; Flatten@Array[a, n, 0]]; a106516[11] (* Michael De Vlieger, Dec 23 2014 *)

Formula

Riordan array (1/(1-3x), x/(1-x)); Number triangle T(n, 0)=A000244(n), T(n, k)=T(n-1, k-1)+T(n-1, k); T(n, k)=sum{j=0..n, binomial(n, k+j)2^j}.
From Peter Bala, Jul 16 2013: (Start)
T(n,k) = binomial(n,k) + 2*sum {i = 1..n} 3^(i-1)*binomial(n-i,k).
O.g.f.: (1 - t)/( (1 - 3*t)*(1 - (1 + x)*t) ) = 1 + (3 + x)*t + (9 + 4*x + x^2)*t^2 + ....
The n-th row polynomial R(n,x) = 1/(x - 2)*( x*(x + 1)^n - 2*3^n ). (End)
Closed-form formula for arbitrary left and right borders of Pascal-like triangle see A228196. - Boris Putievskiy, Aug 19 2013
T(n,k) = 4*T(n-1,k) + T(n-1,k-1) - 3*T(n-2,k) - 3*T(n-2,k-1), T(0,0)=1, T(1,0)=3, T(1,1)=1, T(n,k)=0 if k<0 or if k>n. - Philippe Deléham, Dec 26 2013
From Peter Bala, Dec 23 2014: (Start)
exp(x) * e.g.f. for row n = e.g.f. for diagonal n. For example, for n = 3 we have exp(x)*(27 + 13*x + 5*x^2/2! + x^3/3!) = 27 + 40*x + 58*x^2/2! + 82*x^3/3! + 113*x^4/4! + .... The same property holds more generally for Riordan arrays of the form ( f(x), x/(1 - x) ).
Let M denote the present triangle. For k = 0,1,2,... define M(k) to be the lower unit triangular block array
/I_k 0\
\ 0 M/ having the k X k identity matrix I_k as the upper left block; in particular, M(0) = M. The infinite product M(0)*M(1)*M(2)*..., which is clearly well-defined, is equal to A143495 (but with a different offset). See the Example section. Cf. A055248. (End)
n-th row polynomial R(n, x) = (2*3^n - x*(1 + x)^n)/(2 - x). - Peter Bala, Mar 05 2025

A156922 Third right hand column (n-m=2) of the A156920 triangle.

Original entry on oeis.org

1, 15, 129, 877, 5280, 29658, 159742, 838038, 4323003, 22057825, 111750519, 563535483, 2833221970, 14216330916, 71243079660, 356731958812, 1785306330981, 8931761831331, 44675371382365
Offset: 0

Views

Author

Johannes W. Meijer, Feb 20 2009

Keywords

Crossrefs

Other columns A000340, A156923, A156924.
Equals A156920 third right hand column.
Equals A156919 third right hand column divided by 4.
Equals A142963 third right hand column divided by 2^n.

Formula

a(n) = 14*a(n-1)-75*a(n-2)+196*a(n-3)-263*a(n-4)+174*a(n-5)-45*a(n-6)
a(n) = (4*n^2-108*n*3^n+24*n-378*3^n+375*5^n+35)/32
G.f.: GF1(z;RHCnr=3) = (1+z-6*z^2)/((1-5*z)*(1-3*z)^2*(1-z)^3)
Showing 1-10 of 28 results. Next