A097733 -id:A097733 - cp's OEIS frontend

A097732 Pell equation solutions (7a(n))^2 - 2(5b(n))^2 = -1 with b(n):=A097733(n), n >= 0. Note that D=50=25^2 is not squarefree.

1, 199, 39401, 7801199, 1544598001, 305822602999, 60551330795801, 11988857674965599, 2373733268312392801, 469987198268178808999, 93055091523831091789001, 18424438134520287995413199, 3647945695543493192000024401, 722274823279477131728009418199

Offset: 0

Wolfdieter Lang, Aug 31 2004

Also numbers k such that (7*k+1)^2 + (7*k-1)^2 is a square. - Bruno Berselli, Oct 11 2019

			(x,y) = (7,1), (1393,197), (275807,39005), ... give the positive integer solutions to x^2 - 50*y^2 =-1.

Indranil Ghosh, Table of n, a(n) for n = 0..434
Christian Aebi and Grant Cairns, Lattice equable quadrilaterals III: tangential and extangential cases, Integers (2023) Vol. 23, #A48.
Tanya Khovanova, Recursive Sequences
Giovanni Lucca, Integer Sequences and Circle Chains Inside a Hyperbola, Forum Geometricorum (2019) Vol. 19, 11-16.
H. C. Williams and R. K. Guy, Some fourth-order linear divisibility sequences, Intl. J. Number Theory 7 (5) (2011) 1255-1277.
H. C. Williams and R. K. Guy, Some Monoapparitic Fourth Order Linear Divisibility Sequences Integers, Volume 12A (2012) The John Selfridge Memorial Volume.
Index entries for sequences related to Chebyshev polynomials.
Index entries for linear recurrences with constant coefficients, signature (198,-1).

Cf. A097731 for S(n, 2*99), A100047.

Cf. similar sequences of the type (1/k)*sinh((2*n+1)*arcsinh(k)) listed in A097775.

LinearRecurrence[{198, -1}, {1, 199}, 12] (* Ray Chandler, Aug 11 2015 *)

x='x+O('x^99); Vec((1+x)/(1-2*99*x+x^2)) \\ Altug Alkan, Apr 05 2018

G.f.: (1 + x)/(1 - 2*99*x + x^2).
a(n) = S(n, 2*99) + S(n-1, 2*99) = S(2*n, 10*sqrt(2)), with Chebyshev polynomials of the 2nd kind. See A049310 for the triangle of S(n, x)= U(n, x/2) coefficients. S(-1, x) := 0 =: U(-1, x).
a(n) = ((-1)^n)*T(2*n+1, 7*i)/(7*i) with the imaginary unit i and Chebyshev polynomials of the first kind. See the T-triangle A053120.
a(n) = 198*a(n-1) - a(n-2), n > 1; a(0)=1, a(1)=199. - Philippe Deléham, Nov 18 2008
From Peter Bala, Mar 23 2015: (Start)
a(n) = ( Pell(6*n + 6 - 2*k) + Pell(6*n + 2*k) )/( Pell(6 - 2*k) + Pell(2*k) ), for k an arbitrary integer.
a(n) = ( Pell(6*n + 6 - 2*k - 1) - Pell(6*n + 2*k + 1) )/( Pell(6 - 2*k - 1) - Pell(2*k + 1) ), for k an arbitrary integer, k != 1.
The aerated sequence (b(n))n>=1 = [1, 0, 199, 0, 39401, 0, 7801199, 0, ...] is a fourth-order linear divisibility sequence; that is, if n | m then b(n) | b(m). It is the case P1 = 0, P2 = -196, Q = -1 of the 3-parameter family of divisibility sequences found by Williams and Guy. See A100047 for the connection with Chebyshev polynomials. (End)
a(n) = (1/7)*sinh((2*n + 1)*arcsinh(7)). - Bruno Berselli, Apr 03 2018

A085478 Triangle read by rows: T(n, k) = binomial(n + k, 2*k).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 6, 5, 1, 1, 10, 15, 7, 1, 1, 15, 35, 28, 9, 1, 1, 21, 70, 84, 45, 11, 1, 1, 28, 126, 210, 165, 66, 13, 1, 1, 36, 210, 462, 495, 286, 91, 15, 1, 1, 45, 330, 924, 1287, 1001, 455, 120, 17, 1, 1, 55, 495, 1716, 3003, 3003, 1820, 680, 153, 19, 1

Offset: 0

Philippe Deléham, Aug 14 2003

Coefficient array for Morgan-Voyce polynomial b(n,x). A053122 (unsigned) is the coefficient array for B(n,x). Reversal of A054142. - Paul Barry, Jan 19 2004
This triangle is formed from even-numbered rows of triangle A011973 read in reverse order. - Philippe Deléham, Feb 16 2004
T(n,k) is the number of nondecreasing Dyck paths of semilength n+1, having k+1 peaks. T(n,k) is the number of nondecreasing Dyck paths of semilength n+1, having k peaks at height >= 2. T(n,k) is the number of directed column-convex polyominoes of area n+1, having k+1 columns. - Emeric Deutsch, May 31 2004
Riordan array (1/(1-x), x/(1-x)^2). - Paul Barry, May 09 2005
The triangular matrix a(n,k) = (-1)^(n+k)*T(n,k) is the matrix inverse of A039599. - Philippe Deléham, May 26 2005
The n-th row gives absolute values of coefficients of reciprocal of g.f. of bottom-line of n-wave sequence. - Floor van Lamoen (fvlamoen(AT)planet.nl), Sep 24 2006
Unsigned version of A129818. - Philippe Deléham, Oct 25 2007
T(n, k) is also the number of idempotent order-preserving full transformations (of an n-chain) of height k >=1 (height(alpha) = |Im(alpha)|) and of waist n (waist(alpha) = max(Im(alpha))). - Abdullahi Umar, Oct 02 2008
A085478 is jointly generated with A078812 as a triangular array of coefficients of polynomials u(n,x): initially, u(1,x) = v(1,x) = 1; for n>1, u(n,x) = u(n-1,x)+x*v(n-1)x and v(n,x) = u(n-1,x)+(x+1)*v(n-1,x). See the Mathematica section. - Clark Kimberling, Feb 25 2012
Per Kimberling's recursion relations, see A102426. - Tom Copeland, Jan 19 2016
Subtriangle of the triangle given by (0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, 1, -1, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Mar 26 2012
T(n,k) is also the number of compositions (ordered partitions) of 2*n+1 into 2*k+1 parts which are all odd. Proof: The o.g.f. of column k, x^k/(1-x)^(2*k+1) for k >= 0, is the o.g.f. of the odd-indexed members of the sequence with o.g.f. (x/(1-x^2))^(2*k+1) (bisection, odd part). Thus T(n,k) is obtained from the sum of the multinomial numbers A048996 for the partitions of 2*n+1 into 2*k+1 parts, all of which are odd. E.g., T(3,1) = 3 + 3 from the numbers for the partitions [1,1,5] and [1,3,3], namely 3!/(2!*1!) and 3!/(1!*2!), respectively. The number triangle with the number of these partitions as entries is A152157. - Wolfdieter Lang, Jul 09 2012
The matrix elements of the inverse are T^(-1)(n,k) = (-1)^(n+k)*A039599(n,k). - R. J. Mathar, Mar 12 2013
T(n,k) = A258993(n+1,k) for k = 0..n-1. - Reinhard Zumkeller, Jun 22 2015
The n-th row polynomial in descending powers of x is the n-th Taylor polynomial of the algebraic function F(x)*G(x)^n about 0, where F(x) = (1 + sqrt(1 + 4*x))/(2*sqrt(1 + 4*x)) and G(x) = ((1 + sqrt(1 + 4*x))/2)^2. For example, for n = 4, (1 + sqrt(1 + 4*x))/(2*sqrt(1 + 4*x)) * ((1 + sqrt(1 + 4*x))/2)^8 = (x^4  + 10*x^3 + 15*x^2 + 7*x + 1) + O(x^5). - Peter Bala, Feb 23 2018
Row n also gives the coefficients of the characteristc polynomial of the tridiagonal n X n matrix M_n given in A332602: Phi(n, x) := Det(M_n - x*1_n) = Sum_{k=0..n} T(n, k)*(-x)^k, for n >= 0, with Phi(0, x) := 1. - Wolfdieter Lang, Mar 25 2020
It appears that the largest root of the n-th degree polynomial is equal to the sum of the distinct diagonals of a (2*n+1)-gon including the edge, 1. The largest root of x^3 - 6*x^2 + 5*x - 1 is 5.048917... = the sum of (1 + 1.80193... + 2.24697...). Alternatively, the largest root of the n-th degree polynomial is equal to the square of sigma(2*n+1). Check: 5.048917... is the square of sigma(7), 2.24697.... Given N = 2*n+1, sigma(N) (N odd) can be defined as 1/(2*sin(Pi/(2*N))). Relating to the 9-gon, the largest root of x^4 - 10*x^3 + 15*x^2 - 7*x + 1 is 8.290859..., = the sum of (1 + 1.879385... + 2.532088... + 2.879385...), and is the square of sigma(9), 2.879385... Refer to A231187 for a further clarification of sigma(7). - Gary W. Adamson, Jun 28 2022
For n >=1, the n-th row is given by the coefficients of the minimal polynomial of -4*sin(Pi/(4*n + 2))^2. - Eric W. Weisstein, Jul 12 2023
Denoting this lower triangular array by L, then L * diag(binomial(2*k,k)^2) * transpose(L) is the LDU factorization of A143007, the square array of crystal ball sequences for the A_n X A_n lattices. - Peter Bala, Feb 06 2024
T(n, k) is the number of occurrences of the periodic substring (01)^k in the periodic string (01)^n (see Proposition 4.7 at page 7 in Fang). - Stefano Spezia, Jun 09 2024

			Triangle begins as:
  1;
  1    1;
  1    3    1;
  1    6    5    1;
  1   10   15    7    1;
  1   15   35   28    9    1;
  1   21   70   84   45   11    1;
  1   28  126  210  165   66   13    1;
  1   36  210  462  495  286   91   15    1;
  1   45  330  924 1287 1001  455  120   17    1;
  1   55  495 1716 3003 3003 1820  680  153   19    1;
...
From _Philippe Deléham_, Mar 26 2012: (Start)
(0, 1, 0, 1, 0, 0, 0, ...) DELTA (1, 0, 1, -1, 0, 0, 0, ...) begins:
  1
  0, 1
  0, 1,  1
  0, 1,  3,   1
  0, 1,  6,   5,   1
  0, 1, 10,  15,   7,   1
  0, 1, 15,  35,  28,   9,  1
  0, 1, 21,  70,  84,  45, 11,  1
  0, 1, 28, 126, 210, 165, 66, 13, 1. (End)

Reinhard Zumkeller, Rows n = 0..125 of triangle, flattened
J.-P. Allouche and M. Mendès France, Stern-Brocot polynomials and power series, arXiv preprint arXiv:1202.0211 [math.NT], 2012. - _N. J. A. Sloane_, May 10 2012
Peter Bala, A 4-parameter family of embedded Riordan arrays
E. Barcucci, A. Del Lungo, S. Fezzi, and R. Pinzani, Nondecreasing Dyck paths and q-Fibonacci numbers, Discrete Math., 170, 1997, 211-217.
Paul Barry, A Catalan Transform and Related Transformations on Integer Sequences, J. Integer Sequ., Vol. 8 (2005), Article 05.4.5.
Paul Barry, Symmetric Third-Order Recurring Sequences, Chebyshev Polynomials, and Riordan Arrays, JIS 12 (2009) 09.8.6.
Paul Barry, Notes on the Hankel transform of linear combinations of consecutive pairs of Catalan numbers, arXiv:2011.10827 [math.CO], 2020.
Paul Barry, The second production matrix of a Riordan array, arXiv:2011.13985 [math.CO], 2020.
Paul Barry and Aoife Hennessy, Notes on a Family of Riordan Arrays and Associated Integer Hankel Transforms , JIS 12 (2009) 09.5.3.
Eduardo H. M. Brietzke, Recurrence Relations for Sums of Binomial Coefficients and Some Generalizations, Journal of Integer Sequences, Vol. 27 (2024), Article 24.3.4.
E. Czabarka et al, Enumerations of peaks and valleys on non-decreasing Dyck paths, Disc. Math. 341 (2018) 2789-2807, Theorem 3.
Emeric Deutsch and H. Prodinger, A bijection between directed column-convex polyominoes and ordered trees of height at most three, Theoretical Comp. Science, 307, 2003, 319-325.
James East and Nicholas Ham, Lattice paths and submonoids of Z^2, arXiv:1811.05735 [math.CO], 2018.
Wenjie Fang, Maximal number of subword occurrences in a word, arXiv:2406.02971 [math.CO], 2024.
A. Laradji and A. Umar, Combinatorial results for semigroups of order-preserving full transformations, Semigroup Forum 72 (2006), 51-62. - _Abdullahi Umar_, Oct 02 2008
Donatella Merlini and Renzo Sprugnoli, Arithmetic into geometric progressions through Riordan arrays, Discrete Mathematics 340.2 (2017): 160-174.
Yidong Sun, Numerical Triangles and Several Classical Sequences, Fib. Quart. 43, no. 4, (2005) 359-370.
Eric Weisstein's World of Mathematics, Morgan-Voyce Polynomials

Row sums: A001519. Signed versions: A123970, A129818.
Cf. A000108, A001006, A007318, A098158, A183160, A078812, A091042.
Cf. A258993, A025581, A051162, A054142, A332602.
Cf. A231187, A143007.

Flat(List([0..12], n-> List([0..n], k-> Binomial(n+k, 2*k) ))); # G. C. Greubel, Aug 01 2019

a085478 n k = a085478_tabl !! n !! k
a085478_row n = a085478_tabl !! n
a085478_tabl = zipWith (zipWith a007318) a051162_tabl a025581_tabl
-- Reinhard Zumkeller, Jun 22 2015

[Binomial(n+k, 2*k): k in [0..n], n in [0..12]]; // G. C. Greubel, Aug 01 2019

T := (n,k) -> binomial(n+k,2*k): seq(seq(T(n,k), k=0..n), n=0..11);

(* First program *)
u[1, x_]:= 1; v[1, x_]:= 1; z = 13;
u[n_, x_]:= u[n-1, x] + x*v[n-1, x];
v[n_, x_]:= u[n-1, x] + (x+1)*v[n-1, x];
Table[Expand[u[n, x]], {n, 1, z/2}]
Table[Expand[v[n, x]], {n, 1, z/2}]
cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
TableForm[cu]
Flatten[%]   (* A085478 *)
Table[Expand[v[n, x]], {n, 1, z}]
cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
TableForm[cv]
Flatten[%]   (* A078812 *) (*Clark Kimberling, Feb 25 2012 *)
(* Second program *)
Table[Binomial[n + k, 2 k], {n, 0, 12}, {k, 0, n}] // Flatten (* G. C. Greubel, Aug 01 2019 *)
CoefficientList[Table[Fibonacci[2 n + 1, Sqrt[x]], {n, 0, 10}], x] // Flatten (* Eric W. Weisstein, Jul 03 2023 *)
Join[{{1}}, CoefficientList[Table[MinimalPolynomial[-4 Sin[Pi/(4 n + 2)]^2, x], {n, 20}], x]] (* Eric W. Weisstein, Jul 12 2023 *)

PARI
```
T(n,k) = binomial(n+k,n-k)
```

[[binomial(n+k,2*k) for k in (0..n)] for n in (0..12)] # G. C. Greubel, Aug 01 2019

T(n, k) = (n+k)!/((n-k)!*(2*k)!).
G.f.: (1-z)/((1-z)^2-tz). - Emeric Deutsch, May 31 2004
Row sums are A001519 (Fibonacci(2n+1)). Diagonal sums are A011782. Binomial transform of A026729 (product of lower triangular matrices). - Paul Barry, Jun 21 2004
T(n, 0) = 1, T(n, k) = 0 if n=0} T(n-1-j, k-1)*(j+1). T(0, 0) = 1, T(0, k) = 0 if k>0; T(n, k) = T(n-1, k-1) + T(n-1, k) + Sum_{j>=0} (-1)^j*T(n-1, k+j)*A000108(j). For the column k, g.f.: Sum_{n>=0} T(n, k)*x^n = (x^k) / (1-x)^(2*k+1). - Philippe Deléham, Feb 15 2004
Sum_{k=0..n} T(n,k)*x^(2*k) = A000012(n), A001519(n+1), A001653(n), A078922(n+1), A007805(n), A097835(n), A097315(n), A097838(n), A078988(n), A097841(n), A097727(n), A097843(n), A097730(n), A098244(n), A097733(n), A098247(n), A097736(n), A098250(n), A097739(n), A098253(n), A097742(n), A098256(n), A097767(n), A098259(n), A097770(n), A098262(n), A097773(n), A098292(n), A097776(n) for x=0,1,2,...,27,28 respectively. - Philippe Deléham, Dec 31 2007
T(2*n,n) = A005809(n). - Philippe Deléham, Sep 17 2009
A183160(n) = Sum_{k=0..n} T(n,k)*T(n,n-k). - Paul D. Hanna, Dec 27 2010
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k). - Philippe Deléham, Feb 06 2012
O.g.f. for column k: x^k/(1-x)^(2*k+1), k >= 0. [See the o.g.f. of the triangle above, and a comment on compositions. - Wolfdieter Lang, Jul 09 2012]
E.g.f.: (2/sqrt(x + 4))*sinh((1/2)*t*sqrt(x + 4))*cosh((1/2)*t*sqrt(x)) = t + (1 + x)*t^3/3! + (1 + 3*x + x^2)*t^5/5! + (1 + 6*x + 5*x^2 + x^3)*t^7/7! + .... Cf. A091042. - Peter Bala, Jul 29 2013
T(n, k) = A065941(n+3*k, 4*k) = A108299(n+3*k, 4*k) = A194005(n+3*k, 4*k). - Johannes W. Meijer, Sep 05 2013
Sum_{k=0..n} (-1)^k*T(n,k)*A000108(k) = A000007(n) for n >= 0. - Werner Schulte, Jul 12 2017
Sum_{k=0..floor(n/2)} T(n-k,k)*A000108(k) = A001006(n) for n >= 0. - Werner Schulte, Jul 12 2017
From  Peter Bala, Jun 26 2025: (Start)
The n-th row polynomial b(n, x) = (-1)^n * U(2*n, (i/2)*sqrt(x)), where U(n,x) is the n-th Chebyshev polynomial of the second kind.
b(n, x) = (-1)^n * Dir(n, -1 - x/2), where Dir(n, x) is the n-th row polynomial of the triangle A244419.
b(n, -1 - x) is the n-th row polynomial of A098493. (End)

A188647 Array read by antidiagonals of a(n) = a(n-1)*k-((k-1)/(k^n)) where a(0)=1 and k=(sqrt(x^2+1)+x)^2 for integers x>=1.

Original entry on oeis.org

1, 5, 1, 29, 17, 1, 169, 305, 37, 1, 985, 5473, 1405, 65, 1, 5741, 98209, 53353, 4289, 101, 1, 33461, 1762289, 2026009, 283009, 10301, 145, 1, 195025, 31622993, 76934989, 18674305, 1050601, 21169, 197, 1, 1136689, 567451585, 2921503573, 1232221121, 107151001, 3090529, 39005, 257, 1

Offset: 0

Charles L. Hohn, Apr 06 2011

Conjecture: Given function f(x, y)=(sqrt(x^2+y)+x)^2; constant k=f(x, y); and initial term a(0)=1; then for all integers x>=1 and y=[+-]1, k may be irrational, but sequence a(n)=a(n-1)*k-((k-1)/(k^n)) always produces integer sequences; y=1 results shown here; y=-1 results are A188646.

Also square array A(n,k), n >= 1, k >= 0, read by antidiagonals, where A(n,k) is (1/sqrt(n^2+1)) * T_{2*k+1}(sqrt(n^2+1)), with T the Chebyshev polynomial of the first kind. - Seiichi Manyama, Jan 02 2019

			Square array begins:
     | 0    1       2          3             4
-----+---------------------------------------------
   1 | 1,   5,     29,       169,          985, ...
   2 | 1,  17,    305,      5473,        98209, ...
   3 | 1,  37,   1405,     53353,      2026009, ...
   4 | 1,  65,   4289,    283009,     18674305, ...
   5 | 1, 101,  10301,   1050601,    107151001, ...
   6 | 1, 145,  21169,   3090529,    451196065, ...
   7 | 1, 197,  39005,   7722793,   1529074009, ...
   8 | 1, 257,  66305,  17106433,   4413393409, ...
   9 | 1, 325, 105949,  34539049,  11259624025, ...
  10 | 1, 401, 161201,  64802401,  26050404001, ...
  11 | 1, 485, 235709, 114554089,  55673051545, ...
  12 | 1, 577, 333505, 192765313, 111418017409, ...
  13 | 1, 677, 459005, 311204713, 210996336409, ...
  14 | 1, 785, 617009, 484968289, 381184458145, ...
  15 | 1, 901, 812701, 733055401, 661215159001, ...
  ...

Row 1 is A001653, row 2 is A007805, row 3 is A097315, row 4 is A078988, row 5 is A097727, row 6 is A097730, row 7 is A097733, row 8 is A097736, row 9 is A097739, row 10 is A097742, row 11 is A097767, row 12 is A097770, row 13 is A097773.
Column 1 is A053755.
A(n,n) gives A323012.
Cf. A188645, A188646 (f(x, y) as above with y=-1).

A(n,k) = 2 * A188645(n,k) - A(n,k-1).

A(n,k) = Sum_{j=0..k} binomial(2*k+1,2*j)*(n^2+1)^(k-j)*n^(2*j). - Seiichi Manyama, Jan 02 2019

Edited and extended by Seiichi Manyama, Jan 02 2019

cp's OEIS Frontend

A097732 Pell equation solutions (7*a(n))^2 - 2*(5*b(n))^2 = -1 with b(n):=A097733(n), n >= 0. Note that D=50=2*5^2 is not squarefree.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A085478 Triangle read by rows: T(n, k) = binomial(n + k, 2*k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

Sage

Formula

A188647 Array read by antidiagonals of a(n) = a(n-1)*k-((k-1)/(k^n)) where a(0)=1 and k=(sqrt(x^2+1)+x)^2 for integers x>=1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

Extensions

A097731 Chebyshev U(n,x) polynomial evaluated at x=99 gives 2*7^2+1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A278476 a(n) = floor((1 + sqrt(2))^3*a(n-1)) for n>0, a(0) = 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A097732 Pell equation solutions (7a(n))^2 - 2(5b(n))^2 = -1 with b(n):=A097733(n), n >= 0. Note that D=50=25^2 is not squarefree.