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-2 of 2 results.

A219865 First differences of A219795.

Original entry on oeis.org

0, 0, 0, 1, 0, 2, 2, 3, 1, 5, 7, 10, 11, 14, 23, 33, 44, 54, 81, 114, 158, 212, 287, 407, 565, 777, 1050, 1471, 2036, 2813, 3863, 5310, 7370, 10183, 14046, 19296, 26726, 36909, 50955, 70251, 96867, 133886, 184841, 255092, 351673
Offset: 0

Views

Author

Paul Curtz, Nov 30 2012

Keywords

Comments

2,2,3, 5,7,10,11, 23,33,44, 81,114,158,212, 407,565,777, 1471,2036,2813,3863, 7370,10183,14046, 26726,36909,50955,70251, is a subsequence of A219795(n).
Consider A219795(n) and its successive differences:
2, 2, 2, 2, 3, 3, 5, 7, 10, 11, 16, 23, 33, 44; A219795(n)
0, 0, 0, 1, 0, 2, 2, 3, 1, 5, 7, 10, 11, 14; a(n)
0, 0, 1, -1, 2, 0, 1, -2, 4, 2, 3, 1, 3, 9;
0, 1, -2, 3, -2, 1, -3, 6, -2, 1, -2, 2, 6, 1;

Programs

Formula

a(n) = A219795(n+1) - A219795(n).

A135929 Triangle read by rows: row n gives coefficients of polynomial P_n(x)= U_{n}(x,1) + 3 * U_{n-2}(x,1) where U is the Chebyshev polynomial of the second kind, in order of decreasing exponents.

Original entry on oeis.org

1, 1, 0, 1, 0, 2, 1, 0, 1, 0, 1, 0, 0, 0, -2, 1, 0, -1, 0, -3, 0, 1, 0, -2, 0, -3, 0, 2, 1, 0, -3, 0, -2, 0, 5, 0, 1, 0, -4, 0, 0, 0, 8, 0, -2, 1, 0, -5, 0, 3, 0, 10, 0, -7, 0, 1, 0, -6, 0, 7, 0, 10, 0, -15, 0, 2, 1, 0, -7, 0, 12, 0, 7, 0, -25, 0, 9, 0, 1, 0, -8, 0, 18, 0, 0, 0, -35, 0, 24, 0, -2
Offset: 0

Views

Author

N. J. A. Sloane, Mar 09 2008

Keywords

Comments

Take a(0)=-2 instead of 1. The recurrence begins immediately (at the third instead of the fourth polynomial). Companion: A192011(n). - Paul Curtz, Sep 20 2011

Examples

			The coefficients and polynomials are
  1;                                 1
  1, 0;                              x
  1, 0,  2;                          x^2 + 2
  1, 0,  1, 0;                       x^3 +   x
  1, 0,  0, 0, -2;                   x^4 - 2
  1, 0, -1, 0, -3, 0;                x^5 -   x^3 - 3*x
  1, 0, -2, 0, -3, 0,  2;            x^6 - 2*x^4 - 3*x^2 + 2
  1, 0, -3, 0, -2, 0,  5, 0;         x^7 - 3*x^5 - 2*x^3 + 5*x
  1, 0, -4, 0,  0, 0,  8, 0, -2;     x^8 - 4*x^6 + 8*x^2 - 2
  1, 0, -5, 0,  3, 0, 10, 0, -7, 0;  x^9 - 5*x^7 + 3*x^5 + 10*x^3 - 7*x
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972; see Chapter 22.

Crossrefs

Programs

  • Magma
    A053119:= func< n,k | (1/2)*(-1)^Floor(3*k/2)*(1+(-1)^k)*Binomial(n - Floor(k/2), n-k) >;
    A135929:= func< n,k | n eq 0 select 1 else A053119(n, k) + 3*A053119(n-2, k-2) >;
    [A135929(n,k): k in [0..n], n in [0..16]]; // G. C. Greubel, Apr 24 2023
    
  • Maple
    A135929 := proc(n, m) coeftayl( coeftayl( (1+3*t^2)/(1-x*t+t^2), t=0, n), x=0, n-m) ; end proc: seq(seq(A135929(n,m), m=0..n),n=0..14) ; # R. J. Mathar, Nov 03 2009
  • Mathematica
    p[0, ]= 1; p[1, x]:= x; p[2, x_]:= x^2+2; p[n_, x_]:= p[n, x] = x*p[n-1, x] - p[n-2, x]; row[n_]:= CoefficientList[p[n, x], x]; Table[row[n]//Reverse, {n, 0, 13}]//Flatten (* Jean-François Alcover, Nov 26 2012, after Paul Curtz's formula *)
    (* Second program *)
    p=1; q=2; t[, 0]=p; t[2, 2]=q; t[, ?OddQ]=0; t[n, k_] /; k > n = 0; t[n_ /; n >= 0, k_ /; k >= 0]:= t[n, k] = t[n-1, k] - t[n-2, k-2]; Table[t[n, k], {n, 0, 13}, {k, 0, n}]//Flatten (* Jean-François Alcover, Nov 27 2012, from recurrence *)
  • SageMath
    def A053119(n,k): return (-1)^(3*k/2)*((k+1)%2)*binomial(n-k/2, n-k)
    def A135929(n,k): return 1 if (n==0) else A053119(n, k) + 3*A053119(n-2, k-2)
    flatten([[A135929(n,k) for k in range(n+1)] for n in range(16)]) # G. C. Greubel, Apr 24 2023

Formula

G.f.: (1+3*t^2)/(1-x*t+t^2).
P_n(x) = U_{n}(x,1) + 3 * U_{n-2}(x,1) for n>=2. - Max Alekseyev, Dec 04 2009
P_n(x) = S_{n}(x) + 3*S_{n-2}(x), with Chebyshev Polynomials S_n(x) defined in A049310 and A053119. - R. J. Mathar, Dec 07 2009
P_0(x)=1, P_1(x)=x, P_2(x)=x^2+2, and P_n(x)= x*P_{n-1}(x) - P_{n-2}(x) for n>=3. - Paul Curtz, Aug 14 2011
From G. C. Greubel, Apr 24 2023: (Start)
T(n, k) = A053119(n, k) + 3*A053119(n-2, k-2), with T(0,0) = 1.
Sum_{k=0..n} T(n, k) = A138034(n). (End)

Extensions

Extended by R. J. Mathar, Nov 03 2009
Showing 1-2 of 2 results.