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

A126869 a(n) = Sum_{k = 0..n} binomial(n,floor(k/2))*(-1)^(n-k).

Original entry on oeis.org

1, 0, 2, 0, 6, 0, 20, 0, 70, 0, 252, 0, 924, 0, 3432, 0, 12870, 0, 48620, 0, 184756, 0, 705432, 0, 2704156, 0, 10400600, 0, 40116600, 0, 155117520, 0, 601080390, 0, 2333606220, 0, 9075135300, 0, 35345263800, 0, 137846528820, 0, 538257874440, 0, 2104098963720, 0, 8233430727600, 0, 32247603683100, 0, 126410606437752, 0
Offset: 0

Views

Author

Philippe Deléham, Mar 16 2007

Keywords

Comments

Hankel transform is 2^n. Successive binomial transforms are A002426, A000984, A026375, A081671, A098409, A098410.
From Andrew V. Sutherland, Feb 29 2008: (Start)
Counts returning walks of length n on a 1-d integer lattice with step set {-1,+1}.
Moment sequence of the trace of a random matrix in G = SO(2). If X = tr(A) is a random variable (A distributed with Haar measure on G), then a(n) = E[X^n].
Also the moment sequence of the trace of the k-th power of a random matrix in USp(2) = SU(2), for all k > 2.
(End)
From Paul Barry, Aug 10 2009: (Start)
The Hankel transform of 0,1,0,2,0,6,... is 0,-1,0,4,0,-16,0,... with general term I*(-4)^(n/2)(1 - (-1)^n)/4, I = sqrt(-1).
The Hankel transform of 1,1,0,2,0,6,... (which has g.f. 1 + x/sqrt(1 - 4*x^2)) is A164111. (End)
a(n) = A204293(2*n,n): central terms of the triangle in A204293. - Reinhard Zumkeller, Jan 14 2012
a(n) is the total number of closed walks (round trips) of length n on the graph P_N (a line with N nodes and N-1 edges), divided by N, in the limit N -> infinity. See a comment on A198632 and a link under A201198. - Wolfdieter Lang, Oct 10 2012

Examples

			a(4) = 6 {UUDD,UDUD,UDDU,DUUD,DUDU,DDUU}.
		

References

  • Lin Yang and S.-L. Yang, The parametric Pascal rhombus. Fib. Q., 57:4 (2019), 337-346.

Crossrefs

This is A000984 with interspersed zeros. m-th binomial transforms of A000984: A126869 (m = -2), A002426 (m = -1 and m = -3 for signed version), A000984 (m = 0 and m = -4 for signed version), A026375 (m = 1 and m = -5 for signed version), A081671 (m = 2 and m = -6 for signed version), A098409 (m = 3 and m = -7 for signed version), A098410 (m = 4 and m = -8 for signed version), A104454 (m = 5 and m = -9 for signed version).

Programs

  • Haskell
    a126869 n = a204293_row (2*n) !! n  -- Reinhard Zumkeller, Jan 14 2012
    
  • Maple
    seq((-1)^(n/2)*pochhammer(-n,n/2)/(n/2)!, n=0..43); # Peter Luschny, May 17 2013
    seq(n!*coeff(series(hypergeom([],[1],x^2),x,n+1),x,n),n=0..42); # Peter Luschny, Jan 31 2015
  • Mathematica
    Table[(-1)^Floor[n/2] HypergeometricPFQ[{-n,-n},{1},-1],{n,0,30}] (* Peter Luschny, Nov 01 2011 *)
  • Sage
    A126869 = lambda n: (2^(n-1)*((-1)^n+1)*gamma((n+1)/2))/(sqrt(pi)*gamma((n+2)/2))
    [A126869(n) for n in range(44)] # Peter Luschny, Sep 10 2014

Formula

From Andrew V. Sutherland, Feb 29 2008: (Start)
a(2*n) = binomial(2*n,n) = A000984(n); a(2*n+1) = 0.
a(n) = Sum_{k = 0..n} A107430(n,k)*(-1)^(n-k).
a(n) = Sum_{k = 0..n} A061554(n,k)*(-1)^k.
a(n) = (1/Pi)*Integral_{t = 0..Pi} cos^n(t) dt. (End)
E.g.f.: I_0(2*x) where I_n(x) is the modified Bessel function as a function of x. - Benjamin Phillabaum, Mar 10 2011
G.f.: A(x) = 1/sqrt(1 - 4*x^2). - Vladimir Kruchinin, Apr 16 2011
a(n) = (1/Pi)*Integral{x = -2..2} x^n/sqrt((2 - x)*(2 + x)). - Peter Luschny, Sep 12 2011
a(n) = (-1)^floor(n/2) * Hypergeometric([-n,-n],[1], -1). - Peter Luschny, Nov 01 2011
E.g.f.: E(0)/(1 - x) where E(k) = 1 - x/(1 - x/(x - (k+1)^2/E(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Apr 05 2013
E.g.f.: 1 + x^2/(Q(0) - x^2), where Q(k) = x^2 + (k+1)^2 - x^2*(k+1)^2/Q(k+1); (continued fraction). - Sergei N. Gladkovskii, Apr 28 2013
G.f.: 1/(1 - 2*x^2*Q(0)), where Q(k) = 1 + (4*k+1)*x^2/(k+1 - x^2*(2*k+2)*(4*k+3)/(2*x^2*(4*k+3) + (2*k+3)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 15 2013
G.f.: G(0)/2, where G(k) = 1 + 1/(1 - 2*x/(2*x + (k+1)/(x*(2*k+1))/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, May 24 2013
G.f.: G(0)/(1+x), where G(k) = 1 + x*(2+5*x)*(4*k+1)/((4*k+2)*(1+x)^2 - 2*(2*k+1)*(4*k+3)*x*(2+5*x)*(1+x)^2/((4*k+3)*x*(2+5*x) + 4*(k+1)*(1+x)^2/G(k+1) )); (continued fraction). - Sergei N. Gladkovskii, Jan 19 2014
a(n) = 2^n*JacobiP(n,0,-1/2-n,-3). - Peter Luschny, Aug 02 2014
a(n) = (2^(n-1)*((-1)^n+1)*Gamma((n+1)/2))/(sqrt(Pi)*Gamma((n+2)/2)). - Peter Luschny, Sep 10 2014
a(n) = n!*[x^n]hypergeom([],[1],x^2). - Peter Luschny, Jan 31 2015
a(n) = 2^n*hypergeom([1/2,-n],[1],2). - Peter Luschny, Feb 03 2015
From Peter Bala, Jul 25 2016: (Start)
a(n) = (-1)^floor(n/2)*Sum_{k = 0..n} (-1)^k*binomial(n,k)^2.
D-finite with recurrence: a(n) = 4*(n - 1)/n * a(n-2) with a(0) = 1, a(1) = 0. (End)
From Ilya Gutkovskiy, Jul 25 2016: (Start)
Inverse binomial transform of A002426.
a(n) = Sum_{k=0..n} (-1)^k*A128014(k).
a(n) ~ 2^n*((-1)^n + 1)/sqrt(2*Pi*n). (End)

A198632 Triangle version of the array of the number of closed paths of even length on the graph P_n (n vertices, n-1 edges).

Original entry on oeis.org

1, 0, 2, 0, 2, 3, 0, 2, 4, 4, 0, 2, 8, 6, 5, 0, 2, 16, 14, 8, 6, 0, 2, 32, 36, 20, 10, 7, 0, 2, 64, 94, 56, 26, 12, 8, 0, 2, 128, 246, 164, 76, 32, 14, 9, 0, 2, 256, 644, 488, 234, 96, 38, 16, 10, 0, 2, 512, 1686, 1460, 740, 304, 116, 44, 18, 11, 0, 2, 1024, 4414, 4376, 2372, 992, 374, 136, 50, 20, 12, 0, 2, 2048, 11556, 13124, 7654, 3296, 1244, 444, 156, 56, 22, 13
Offset: 0

Views

Author

Wolfdieter Lang, Nov 02 2011

Keywords

Comments

This array is an example of counting walks on a graph whose adjacency matrix is given by a special symmetric tridiagonal matrix with nonnegative integer entries, appropriate for orthogonal polynomials. These are quadratic Jacobi matrices J_n with nonnegative entries. The corresponding graphs could be called Jacobi graphs. Here Chebyshev S-polynomials (coefficients A049310) are considered, which belong to the Jacobi class of the classical orthogonal polynomial systems. The corresponding graph has adjacency matrix [[0,1,0,...],[1,0,1,...],[0,1,0,1,...]...[0,...0,1,0]] (n rows and n columns), with characteristic polynomial S(n,x) (see also a comment by Michael Somos on A049310).
w(n,l;p_k->p_m) = ((J_n)^l)(k,m) is the number of walks of length l from vertex p_k to vertex p_m on such a Jacobi graph. w(n,0; p_k->p_m) = delta(k,m), with the Kronecker symbol delta. The total number of closed walks of length l is w(n,l):=Sum_{i=1..n} w(n,l; p_i->p_i) = trace(J_n^l), which is the l-th power sum of the eigenvalues of J_n, i.e., the zeros of the characteristic polynomial for J_n. There are theorems for the o.g.f. of the normalized power sums of these zeros. See, e.g., the given W. Lang reference, p. 244. This results for the S-polynomial in the o.g.f. G(n,x) = Sum_{l=0..infinity} w(n,l)*x^l = y*(d/dy)S(n,y)/S(y) with y=1/x. This can be rewritten in the form given in the formula section (this results from eq. (3.8b) of the W. Lang reference, and in eq. (3.8d) it should be coth, not tanh).
From Wolfdieter Lang, Oct 10 2012: (Start)
For an accompanying paper on path counting on Jacobi graphs see the W. Lang link under A201198.
The total number of round trips of length L on the graph P_n, taken per site, becomes for n -> infinity A126869(L). See the just mentioned link, p. 8. This limit is derived from the limit of G(n,x)/n with G(n,x) given in the formula section.
Thanks go to Clyde P. Kruskal for asking a question which led to this comment.
(End)

Examples

			The array w(n,2*k) is
n\k  0  1   2   3   4    5    6     7     8      9 ...
1:   1  0   0   0   0    0    0     0     0      0
2:   2  2   2   2   2    2    2     2     2      2
3:   3  4   8  16  32   64  128   256   512   1024
4:   4  6  14  36  94  246  644  1686  4414  11556
5:   5  8  20  56 164  488 1460  4376 13124  39368
6:   6 10  26  76 234  740 2372  7654 24778  80338
7:   7 12  32  96 304  992 3296 11072 37440 127104
8:   8 14  38 116 374 1244 4220 14504 50294 175454
9:   9 16  44 136 444 1496 5144 17936 63164 224056
...
The triangle is
k\n 1  2    3    4    5    6   7    8   9 10 11 12 ...
0:  1
1:  0  2
2:  0  2    3
3:  0  2    4    4
4:  0  2    8    6    5
5:  0  2   16   14    8    6
6:  0  2   32   36   20   10   7
7:  0  2   64   94   56   26  12    8
8:  0  2  128  246  164   76  32   14   9
9:  0  2  256  644  488  234  96   38  16 10
10: 0  2  512 1686 1460  740 304  116  44 18 11
11: 0  2 1024 4414 4376 2372 992  374 136 50 20 12
...
n=3, l=2*k = 4: graph P_3 as 1-2-3, with eight walks of length 4, namely 12121, 12321, 21212, 23232, 21232, 23212, 32323 and 32123.
		

Crossrefs

Column sequences: A000007, 2*A000012, A198633, 2*A005248, A198635, ...

Formula

a(k,n)=w(n,2*(k-n+2)), the total number of closed walks (paths) of length 2*(k-n+2) on the graph P_n, which looks like o-o-o...-o, with n vertices (nodes) and n-1 edges (lines), k+1>=n>=1.
O.g.f. G(n,x) for w(n,l), which vanishes for odd l, is
((n+1)*coth((n+1)*log((2*x)/(1-sqrt(1-(2*x)^2)))) - 1/sqrt(1-(2*x)^2))/sqrt(1-(2*x)^2). See the comment above for a version with Chebyshev S-polynomials.
Conjecture: For the array w(n,2*k) in the example below, w(2*q,2*k)/2 = A185095(q,k), q >= 1, k >= 0. - L. Edson Jeffery, Nov 23 2013

A201199 Triangle version of the array w(N,L) of the total number of round trips of length L on closed Laguerre graphs Lc_N.

Original entry on oeis.org

1, 1, 2, 1, 4, 3, 1, 18, 9, 4, 1, 76, 53, 16, 5, 1, 322, 357, 120, 25, 6, 1, 1364, 2489, 1024, 233, 36, 7, 1, 5778, 17509, 9424, 2545, 404, 49, 8, 1, 24476, 123449, 89536, 29985, 5400, 645, 64, 9, 1, 103682, 870893, 862560, 367505, 78392, 10213, 968, 81, 10
Offset: 0

Views

Author

Wolfdieter Lang, Nov 30 2011

Keywords

Comments

For Laguerre graphs (open and closed ones) see the W. Lang link on Jacobi graphs under A201198. There one also finds a sketch of the closed Laguerre graph Lc_4 as Fig.4.
The total number of round trips on the closed Laguerre graph Lc_N, for N>=3, with N vertices N^2 loops, binomial(N,2) lines between neighboring vertices and two lines between the first and the last vertex (in total (3*N-1)*N/2+2 = (3*N^2-N+4)/2 lines) is w(N,L) = sum(w(N,L;p_n->p_n),n=1..N) = Trace((L_N)^L) = sum((x_n^{(N)})^L,n=1..N), with the N x N symmetric adjacency matrix, also called Lc_N, having non-vanishing elements (Lc_N)[n,n] = 2*n-1, n=1..N, (Lc_N)[n,n+1] = (Lc_N)[n+1,n] = n, n=1..N-1, and (Lc_N)[1,N]= 2=(Lc_N)[N,1]. The eigenvalues of Lc_N are x_n^{(N)}. They are the zeros of the characteristic polynomial Lac_N(x):=Det(x*1_N -Lc_N) with the N x N unit matrix 1_N. These are the polynomials Lac_N(x) = La(N,x) - 4*La1(N-2,x) - 4*(N-1)!, with the ordinary monic Laguerre polynomials La(N,x) with coefficient array given by A021009(n,m)*(-1)^n and the first associated monic Laguerre polynomials La1(N-2,x) with coefficient array given by A199577(n,m). For N=1 one has Lc_1=L_1 (Laguerre graph with one vertex and one loop) with L_1(x)=x-1, and for N=2 one has a graph where one vertex has one loop, the other three, and there are two lines joining these vertices, hence Lc_2(x)= x^2-4*x-1.

Examples

			The array w(N,L) starts:
N\L 0   1    2     3      4        5         6  ...
1:  1   1    1     1      1        1         1
2:  2   4   12    40    136      464      1584
3:  3   9   53   357   2489    17509    123449
4:  4  16  120  1024   9424    89536    862560
5:  5  25  233  2545  29985   367505   4599521
6:  6  36  404  5400  78392  1188336  18460016
7:  7  49  645 10213 176473  3195829  59473593
8:  8  64  968 17728 355536  7493504 162671840
9:  9  81 1385 28809 657953 15826041 392792273
...The triangle a(K,N) = w(N,K-N+1) starts:
K\N 1      2       3      4      5     6     7   8  9..
0:  1
1:  1      2
2:  1      4       3
3:  1     18       9      4
4:  1     76      53     16      5
5:  1    322     357    120     25     6
6:  1   1364    2489   1024    233    36     7
7:  1   5778   17509   9424   2545   404    49   8
8:  1  24476  123449  89536  29985  5400   645  64  9
...
For the graph Lc_4, shown in the W. Lang link as Figure 4, the counting for round trips of length L=2 for each of the four vertices V_i, i=1..4, read from left to right, is as follows.
V_1: 1+1+(1+1+2*1), V_2: 3+2*binomial(3,2)+1+(1+1+2*1),
V_3: 5+2*binomial(5,2)+(1+1+2*1)+(3+2*binomial(3,2)),
V_4: 7+2*binomial(7,2)+(3+2*binomial(3,2))+(1+1+2*1),
this sums to the total number  w(4,2)= 120  =  a(5,4).
Compared to the open L_4 graph (see the corresponding A201198 entry 4*28 = 112) one has to add 2*(1+1+2*1)=8 from the new two lines joining V_1 and V_4.
		

Crossrefs

Cf. A201198 (open Laguerre graphs).

Formula

a(K,N) = w(N,K-N+1),K>=0, N=1,...,K+1, with w(N,L) the total number of round trips of length L on the closed Laguerre graph Lc_N described above in the comment section.
The o.g.f. of w(N,L) is: G(N,x)=y*(d/dx)Lac_N(x)/Lac_N(x) with y=1/x.
The characteristic polynomial Lac_N(x) has also been given in the comment section above.

A199579 Average number of round trips of length n on the Laguerre graph L_4.

Original entry on oeis.org

1, 4, 28, 232, 2056, 18784, 174112, 1625152, 15220288, 142777600, 1340416768, 12588825088, 118252556800, 1110898849792, 10436554713088, 98050271875072, 921180638875648, 8654518327066624, 81309636020912128
Offset: 0

Views

Author

Wolfdieter Lang, Dec 02 2011

Keywords

Comments

See the general array and triangle for the average number of round trips of length L on (open) Laguerre graphs L_N given in A201198. Here a(n) = w(4,L=n), n>=0, the fourth row in this array. In the corresponding triangle this is the column no. N=4 without leading zeros: a(n) = A201198(n+3,4), n>=0.
For a sketch of this Laguerre graph L_4 see Figure 3 of the W. Lang link. The o.g.f. is also given there.
By definition the number of zero length round trips of length 0 for a vertex is put to 1 in order to count vertices.

Examples

			n=0: a(0)=1 because the average number of vertices is 4/4=1.
a(1)= (1+3+5+7)/4 = 4, from the sum of the self-loops of L_4 divided by the number of vertices 4.
The counting for n=2, a(2)= 112/4 = 28, has been given as an example to A201198.
		

Crossrefs

Cf. A201198, A201199 (closed Laguerre graphs), A201200 (closed L_4 graph).

Programs

  • Magma
    I:=[1, 4, 28, 232]; [n le 4 select I[n] else 16*Self(n-1) - 72*Self(n-2) + 96*Self(n-3) -24*Self(n-4): n in [1..30]]; // G. C. Greubel, May 14 2018
  • Mathematica
    LinearRecurrence[{16, -72, 96, -24}, {1, 4, 28, 232}, 50] (* G. C. Greubel, May 14 2018 *)
  • PARI
    x='x+O('x^30); Vec((1-12*x+36*x^2-24*x^3)/(1-16*x+72*x^2- 96*x^3 +24*x^4)) \\ G. C. Greubel, May 14 2018
    

Formula

a(n) = A201198(n+3,4), n>=0.
O.g.f.: (1-12*x+36*x^2-24*x^3)/(1-16*x+72*x^2-96*x^3+24*x^4).

A201200 Total number of round trips of length n on the closed Laguerre graph Lc_4 divided by 4.

Original entry on oeis.org

1, 4, 30, 256, 2356, 22384, 215640, 2090176, 20315536, 197702464, 1925042400, 18749072896, 182629124416, 1779030655744, 17330352562560, 168824779580416, 1644626142474496, 16021353180980224, 156074394613317120, 1520422660926324736
Offset: 0

Views

Author

Wolfdieter Lang, Dec 02 2011

Keywords

Comments

For the general array and triangle for the total number of round trips of length L on closed Laguerre graphs Lc_N see A201199. Here a(n)=w(4,L=n)/4, n>=0, the fourth row of this array divided by 4. In the corresponding triangle a(n) = A201199(n+3,4)/4, n>=0.
For a sketch of the closed Laguerre graph Lc_4 see Figure 4 of the given W. Lang link. The o.g.f. is also found there.
By definition the number of length 0 round trips for a vertex is put to 1 in order to count vertices.
The average number of round trips of length n on a closed Laguerre graph Lc_N is in general a fraction. Therefore A201199 tabulates the total number of round trips.

Crossrefs

Cf. A201199, A201198 (open Laguerre graphs). A199579 (open L_4 graph).

Programs

  • Magma
    I:=[1, 4, 30, 256]; [n le 4 select I[n] else 16*Self(n-1) - 68*Self(n-2) + 64*Self(n-3) + 44*Self(n-4): n in [1..30]]; // G. C. Greubel, May 13 2018
  • Mathematica
    LinearRecurrence[{16,-68,64,44}, {1, 4, 30, 256}, 30] (* G. C. Greubel, May 13 2018 *)
  • PARI
    Vec((1-8*x)*(1-4*x+2*x^2)/((1-4*x-2*x^2)*(1-12*x+22*x^2)) + O(x^50)) \\ Colin Barker, Apr 27 2016
    

Formula

a(n) = A201199(n+3,4)/4, n>=0.
O.g.f.: (8*x-1)*(2*x^2-4*x+1) / ( (22*x^2-12*x+1)*(2*x^2+4*x-1) ).
From Colin Barker, Apr 27 2016: (Start)
a(n) = 16*a(n-1)-68*a(n-2)+64*a(n-3)+44*a(n-4) for n>3.
a(n) = ((2-sqrt(6))^n+(2+sqrt(6))^n+(6-sqrt(14))^n+(6+sqrt(14))^n)/4.
(End)
E.g.f.: (exp((2-sqrt(6))*x) + exp((2+sqrt(6))*x) + exp((6-sqrt(14))*x) + exp((6+sqrt(14))*x))/4. - Ilya Gutkovskiy, Apr 27 2016

Extensions

Typo in formula fixed by Colin Barker, Apr 27 2016
Showing 1-5 of 5 results.