%I #53 Dec 18 2017 22:45:40
%S 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,
%T 94,56,26,12,8,0,2,128,246,164,76,32,14,9,0,2,256,644,488,234,96,38,
%U 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
%N Triangle version of the array of the number of closed paths of even length on the graph P_n (n vertices, n-1 edges).
%C 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).
%C 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).
%C From _Wolfdieter Lang_, Oct 10 2012: (Start)
%C For an accompanying paper on path counting on Jacobi graphs see the W. Lang link under A201198.
%C 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.
%C Thanks go to Clyde P. Kruskal for asking a question which led to this comment.
%C (End)
%H S. Barbero, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL17/Barbero/barbero15.html">Dickson Polynomials, Chebyshev Polynomials, and Some Conjectures of Jeffery</a>, Journal of Integer Sequences, 17 (2014), #14.3.8.
%H S. Barbero, U. Cerruti, N. Murru, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL16/Barbero/barbero11.html">Identities Involving Zeros of Ramanujan and Shanks Cubic Polynomials</a>, J. Integer Seq., Vol. 16 (2013), Article 13.8.1.
%H Carlos M. da Fonseca, M. Lawrence Glasser, Victor Kowalenko, <a href="http://arxiv.org/abs/1601.07839">Basic trigonometric power sums with applications</a>, arXiv:1601.07839 [math.NT], 2016. See Theorem 5.1.
%H W. Lang, <a href="http://dx.doi.org/10.1016/S0377-0427(97)00240-9">On sums of powers of zeros of polynomials</a>, J. Comp. Appl. Math., 89 (1998) 237-256.
%F 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.
%F O.g.f. G(n,x) for w(n,l), which vanishes for odd l, is
%F ((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.
%F 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
%e The array w(n,2*k) is
%e n\k 0 1 2 3 4 5 6 7 8 9 ...
%e 1: 1 0 0 0 0 0 0 0 0 0
%e 2: 2 2 2 2 2 2 2 2 2 2
%e 3: 3 4 8 16 32 64 128 256 512 1024
%e 4: 4 6 14 36 94 246 644 1686 4414 11556
%e 5: 5 8 20 56 164 488 1460 4376 13124 39368
%e 6: 6 10 26 76 234 740 2372 7654 24778 80338
%e 7: 7 12 32 96 304 992 3296 11072 37440 127104
%e 8: 8 14 38 116 374 1244 4220 14504 50294 175454
%e 9: 9 16 44 136 444 1496 5144 17936 63164 224056
%e ...
%e The triangle is
%e k\n 1 2 3 4 5 6 7 8 9 10 11 12 ...
%e 0: 1
%e 1: 0 2
%e 2: 0 2 3
%e 3: 0 2 4 4
%e 4: 0 2 8 6 5
%e 5: 0 2 16 14 8 6
%e 6: 0 2 32 36 20 10 7
%e 7: 0 2 64 94 56 26 12 8
%e 8: 0 2 128 246 164 76 32 14 9
%e 9: 0 2 256 644 488 234 96 38 16 10
%e 10: 0 2 512 1686 1460 740 304 116 44 18 11
%e 11: 0 2 1024 4414 4376 2372 992 374 136 50 20 12
%e ...
%e 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.
%Y Column sequences: A000007, 2*A000012, A198633, 2*A005248, A198635, ...
%K nonn,easy,tabl
%O 0,3
%A _Wolfdieter Lang_, Nov 02 2011