%I #65 Jan 01 2019 12:01:14
%S 1,1,1,1,7,1,1,97,17,1,1,1351,577,31,1,1,18817,19601,1921,49,1,1,
%T 262087,665857,119071,4801,71,1,1,3650401,22619537,7380481,470449,
%U 10081,97,1,1,50843527,768398401,457470751,46099201,1431431,18817,127,1
%N Array of (k^n + k^(-n))/2 where k = (sqrt(x^2-1) + x)^2 for integers x >= 1.
%C Conjecture: Given the function f(x,y) = (sqrt(x^2+y) + x)^2 and constant k=f(x,y), then for all integers x >= 1 and y=[+-]1, k may be irrational, but (k^n + k^(-n))/2 always produces integer sequences; y=-1 results shown here; y=1 results are A188645.
%C Also square array A(n,k), n >= 1, k >= 0, read by antidiagonals, where A(n,k) is Chebyshev polynomial of the first kind T_{2*k}(x), evaluated at x=n. - _Seiichi Manyama_, Dec 30 2018
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Chebyshev_polynomials">Chebyshev polynomials</a>.
%H <a href="/index/Ch#Cheby">Index entries for sequences related to Chebyshev polynomials.</a>
%F A(n,k) = (A188646(n,k-1) + A188646(n,k))/2.
%F A(n,k) = Sum_{j=0..k} binomial(2*k,2*j)*(n^2-1)^(k-j)*n^(2*j). - _Seiichi Manyama_, Jan 01 2019
%e Row 2 gives {( (2+sqrt(3))^(2*n) + (2-sqrt(3))^(2*n) )/2}.
%e Square array begins:
%e | 0 1 2 3 4
%e -----+---------------------------------------------
%e 1 | 1, 1, 1, 1, 1, ...
%e 2 | 1, 7, 97, 1351, 18817, ...
%e 3 | 1, 17, 577, 19601, 665857, ...
%e 4 | 1, 31, 1921, 119071, 7380481, ...
%e 5 | 1, 49, 4801, 470449, 46099201, ...
%e 6 | 1, 71, 10081, 1431431, 203253121, ...
%e 7 | 1, 97, 18817, 3650401, 708158977, ...
%e 8 | 1, 127, 32257, 8193151, 2081028097, ...
%e 9 | 1, 161, 51841, 16692641, 5374978561, ...
%e 10 | 1, 199, 79201, 31521799, 12545596801, ...
%e 11 | 1, 241, 116161, 55989361, 26986755841, ...
%e 12 | 1, 287, 164737, 94558751, 54276558337, ...
%e 13 | 1, 337, 227137, 153090001, 103182433537, ...
%e 14 | 1, 391, 305761, 239104711, 186979578241, ...
%e 15 | 1, 449, 403201, 362074049, 325142092801, ...
%e ...
%t max = 9; y = -1; t = Table[k = ((x^2 + y)^(1/2) + x)^2; ((k^n) + (k^(-n)))/2 // FullSimplify, {n, 0, max - 1}, {x, 1, max}]; Table[ t[[n - k + 1, k]], {n, 1, max}, {k, 1, n}] // Flatten (* _Jean-François Alcover_, Jul 17 2013 *)
%Y Row 2 is A011943, row 3 is A056771, row 8 is A175633, (row 2)*2 is A067902, (row 9)*2 is A089775.
%Y Column 0-5 give A000012, A056220, A144130, A243132, A243134, A243136.
%Y (column 1)*2 is A060626.
%Y Cf. A188645 (f(x, y) as above with y=1).
%Y Diagonals give A173129, A322899.
%Y Cf. A188646, A322836.
%K nonn,tabl
%O 0,5
%A _Charles L. Hohn_, Apr 06 2011
%E Edited by _Seiichi Manyama_, Dec 30 2018
%E More terms from _Seiichi Manyama_, Jan 01 2019