A072256 -id:A072256 - cp's OEIS frontend

A138288 a(n) = A054320(n) - A001078(n).

1, 9, 89, 881, 8721, 86329, 854569, 8459361, 83739041, 828931049, 8205571449, 81226783441, 804062262961, 7959395846169, 78789896198729, 779939566141121, 7720605765212481, 76426118085983689, 756540575094624409, 7488979632860260401, 74133255753507979601, 733843577902219535609

Offset: 0

Reinhard Zumkeller, Mar 12 2008

Numbers k such that 6*k^2 - 2 is a square. - Bruno Berselli, Feb 10 2014

			1 + 9*x + 89*x^2 + 881*x^3 + 8721*x^4 + 86329*x^5 + ...

H. Brocard, Note #2049, L'Intermédiaire des Mathématiciens, 8 (1901), pp. 212-213. - N. J. A. Sloane, Mar 02 2022

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Jean-Paul Allouche, Jeffrey Shallit, and Manon Stipulanti, Combinatorics on words and generating Dirichlet series of automatic sequences, arXiv:2401.13524 [math.CO], 2024.
Bruno Deschamps, Sur les bonnes valeurs initiales de la suite de Lucas-Lehmer, Journal of Number Theory, Volume 130, Issue 12, December 2010, Pages 2658-2670.
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (I).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (II).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (III).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (IV).
Index entries for linear recurrences with constant coefficients, signature (10, -1).

Cf. A001078, A001079, A072256, A138281.
Cf. similar sequences listed in A238379.
Cf. A041007, A041039, A054320.

CoefficientList[Series[(1 - x)/(1 - 10 x + x^2), {x, 0, 40}], x] (* Vincenzo Librandi, Feb 12 2014 *)
a[c_, n_] := Module[{},
  p := Length[ContinuedFraction[ Sqrt[ c]][[2]]];
  d := Denominator[Convergents[Sqrt[c], n p]];
  t := Table[d[[1 + i]], {i, 0, Length[d] - 1, p}];
  Return[t];
  ] (* Complement of A041007, A041039 *)
a[6, 20] (* Gerry Martens, Jun 07 2015 *)

{a(n) = subst( poltchebi(n+1) + poltchebi(n), x, 5) / 6} /* Michael Somos, Jan 25 2013 */

[lucas_number1(n,10,1)-lucas_number1(n-1,10,1) for n in range(1, 20)] # Zerinvary Lajos, Nov 10 2009

a(n) = A072256(n+1).
a(n) = A001079(n) + 2*A001078(n).
a(n) = 10*a(n-1) - a(n-2). a(-1) = a(0) = 1.
(sqrt(2)+sqrt(3))^(2*n+1) = A054320(n-1)*sqrt(2) + a(n)*sqrt(3).
From Michael Somos, Jan 25 2013: (Start)
G.f.: (1 - x) / (1 - 10*x + x^2).
a(-1-n) = a(n). (End)
a(n) = sqrt(2+(5-2*sqrt(6))^(1+2*n)+(5+2*sqrt(6))^(1+2*n))/(2*sqrt(3)). - Gerry Martens, Jun 04 2015
E.g.f.: exp(5*x)*(3*cosh(2*sqrt(6)*x) + sqrt(6)*sinh(2*sqrt(6)*x))/3. - Stefano Spezia, May 16 2023

A160682 The list of the A values in the common solutions to 13k+1 = A^2 and 17k+1 = B^2.

Original entry on oeis.org

1, 14, 209, 3121, 46606, 695969, 10392929, 155197966, 2317576561, 34608450449, 516809180174, 7717529252161, 115246129602241, 1720974414781454, 25699370092119569, 383769576967012081, 5730844284413061646, 85578894689228912609, 1277952576054020627489

Offset: 1

Paul Weisenhorn, May 23 2009

This summarizes the case C=13 of common solutions to C*k+1=A^2, (C+4)*k+1=B^2.
The 2 equations are equivalent to the Pell equation x^2-C*(C+4)*y^2=1,
with x=(C*(C+4)*k+C+2)/2; y=A*B/2 and with smallest values x(1) = (C+2)/2, y(1)=1/2.
Generic recurrences are:
A(j+2)=(C+2)*A(j+1)-A(j) with A(1)=1; A(2)=C+1.
B(j+2)=(C+2)*B(j+1)-B(j) with B(1)=1; B(2)=C+3.
k(j+3)=(C+1)*(C+3)*( k(j+2)-k(j+1) )+k(j) with k(1)=0; k(2)=C+2; k(3)=(C+1)*(C+2)*(C+3).
x(j+2)=(C^2+4*C+2)*x(j+1)-x(j) with x(1)=(C+2)/2; x(2)=(C^2+4*C+1)*(C+2)/2;
Binet-type of solutions of these 2nd order recurrences are:
R=C^2+4*C; S=C*sqrt(R); T=(C+2); U=sqrt(R); V=(C+4)*sqrt(R);
A(j)=((R+S)*(T+U)^(j-1)+(R-S)*(T-U)^(j-1))/(R*2^j);
B(j)=((R+V)*(T+U)^(j-1)+(R-V)*(T-U)^(j-1))/(R*2^j);
x(j)+sqrt(R)*y(j)=((T+U)*(C^2*4*C+2+(C+2)*sqrt(R))^(j-1))/2^j;
k(j)=(((T+U)*(R+2+T*U)^(j-1)+(T-U)*(R+2-T*U)^(j-1))/2^j-T)/R. [Paul Weisenhorn, May 24 2009]
.C -A----- -B----- -k-----
01 A001519 A002878 A058038
02 A001653 A002315 A045899/2
03 A004253 A030221 A160695
04 A001653 A002315 A078522/4
05 A049685 A033890 A161582
06 A070997 A057080 A159683/2
07 A070998 A057081 A161585
08 A072256 A054320 A045502/4
09 A078922 A097783 A161586
10 A077417 A077416 A159681/2
11 A085260 A126816 A161588
12 A001570 A028230 A059989/4
13 A160682 A161591 A161584
14 A157456 A159678 A159679/2
15 A161595 A161599 A161583
16 A007805 A049629 A157459/4
For n>=2, a(n) equals the permanent of the (2n-2)X(2n-2) tridiagonal matrix with sqrt(13)'s along the main diagonal, and 1's along the superdiagonal and the subdiagonal. [John M. Campbell, Jul 08 2011]
Positive values of x (or y) satisfying x^2 - 15xy + y^2 + 13 = 0. - Colin Barker, Feb 11 2014

Vincenzo Librandi, Table of n, a(n) for n = 1..200
J.-C. Novelli, J.-Y. Thibon, Hopf Algebras of m-permutations,(m+1)-ary trees, and m-parking functions, arXiv preprint arXiv:1403.5962 [math.CO], 2014.
Index entries for the Pell equation
Index entries for linear recurrences with constant coefficients, signature (15,-1).

Cf. similar sequences listed in A238379.

I:=[1,14]; [n le 2 select I[n] else 15*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Feb 12 2014

LinearRecurrence[{15,-1},{1,14},20] (* Harvey P. Dale, Oct 08 2012 *)
CoefficientList[Series[(1 - x)/(1 - 15 x + x^2), {x, 0, 40}], x] (* Vincenzo Librandi, Feb 12 2014 *)

a(n) = round((2^(-1-n)*((15-sqrt(221))^n*(13+sqrt(221))+(-13+sqrt(221))*(15+sqrt(221))^n))/sqrt(221)) \\ Colin Barker, Jul 25 2016

a(n) = 15*a(n-1)-a(n-2).
G.f.: (1-x)*x/(1-15*x+x^2).
a(n) = (2^(-1-n)*((15-sqrt(221))^n*(13+sqrt(221))+(-13+sqrt(221))*(15+sqrt(221))^n))/sqrt(221). - Colin Barker, Jul 25 2016

Edited, extended by R. J. Mathar, Sep 02 2009

First formula corrected by Harvey P. Dale, Oct 08 2012

A077416 Chebyshev S-sequence with Diophantine property.

Original entry on oeis.org

1, 13, 155, 1847, 22009, 262261, 3125123, 37239215, 443745457, 5287706269, 63008729771, 750817050983, 8946795882025, 106610733533317, 1270382006517779, 15137973344680031, 180385298129642593

Offset: 0

Wolfdieter Lang, Nov 29 2002

7*b(n)^2 - 5*a(n)^2 = 2 with companion sequence b(n) = A077417(n), n>=0.
a(n) = L(n,-12)*(-1)^n, where L is defined as in A108299; see also A077417 for L(n,+12). - Reinhard Zumkeller, Jun 01 2005
The aerated sequence (b(n))n>=1 = [1, 0, 13, 0, 155, 0, 1857, 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 = -10, Q = -1 of the 3-parameter family of divisibility sequences found by Williams and Guy. See A100047. - Peter Bala, May 12 2025

Ivan Panchenko, Table of n, a(n) for n = 0..200
K. Andersen, L. Carbone, and D. Penta, Kac-Moody Fibonacci sequences, hyperbolic golden ratios, and real quadratic fields, Journal of Number Theory and Combinatorics, Vol 2, No. 3 pp 245-278, 2011. See Section 9.
Alex Fink, Richard K. Guy, and Mark Krusemeyer, Partitions with parts occurring at most thrice, Contributions to Discrete Mathematics, Vol 3, No 2 (2008), pp. 76-114. See Section 13.
Tanya Khovanova, Recursive Sequences
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 (12,-1).

Cf. A054320(n-1) with companion A072256(n), n>=1.

I:=[1, 13]; [n le 2 select I[n] else 12*Self(n-1) - Self(n-2): n in [1..30]]; // G. C. Greubel, Jan 18 2018

LinearRecurrence[{12,-1},{1,13},30] (* Harvey P. Dale, Apr 03 2013 *)

x='x+O('x^30); Vec((1+x)/(1-12*x+x^2)) \\ G. C. Greubel, Jan 18 2018

[(lucas_number2(n,12,1)-lucas_number2(n-1,12,1))/10 for n in range(1, 18)] # Zerinvary Lajos, Nov 10 2009

a(n) = 12*a(n-1) - a(n-2), a(-1)=-1, a(0)=1.
a(n) = S(n, 12) + S(n-1, 12) = S(2*n, sqrt(14)) with S(n, x) := U(n, x/2) Chebyshev's polynomials of the second kind. See A049310. S(-1, x)=0, S(n, 12) = A004191(n).
G.f.: (1+x)/(1-12*x+x^2).
a(n) = (ap^(2*n+1) - am^(2*n+1))/(ap - am) with ap := (sqrt(7)+sqrt(5))/sqrt(2) and am := (sqrt(7)-sqrt(5))/sqrt(2).
a(n) = Sum_{k=0..n} (-1)^k * binomial(2*n-k,k) * 14^(n-k).
a(n) = sqrt((7*A077417(n)^2 - 2)/5).
From Peter Bala, May 09 2025: (Start)
a(n) = Dir(n, 6), where Dir(n, x) denotes the n-th row polynomial of the triangle A244419.
a(n)^2 - 12*a(n)*a(n+1) + a(n+1)^2 = 14.
More generally, for real x, a(n+x)^2 - 12*a(n+x)*a(n+x+1) + a(n+x+1)^2 = 14, where a(n) := (ap^(2*n+1) - am^(2*n+1))/(ap - am), ap := sqrt(7/2) + sqrt(5/2) and am := sqrt(7/2) - sqrt(5/2), as given above.
Sum_{n >= 1} (-1)^(n+1)/(a(n) - 1/a(n)) = 1/14 (telescoping series).
Product_{n >= 1} (a(n) + 1)/(a(n) - 1) = sqrt(7/5) (telescoping product). (End)

A133607 Triangle read by rows: T(n, k) = qStirling2(n, k, q) for q = -1, with 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 1, -1, 0, 1, -1, -1, 0, 1, -1, -2, 1, 0, 1, -1, -3, 2, 1, 0, 1, -1, -4, 3, 3, -1, 0, 1, -1, -5, 4, 6, -3, -1, 0, 1, -1, -6, 5, 10, -6, -4, 1, 0, 1, -1, -7, 6, 15, -10, -10, 4, 1, 0, 1, -1, -8, 7, 21, -15, -20, 10, 5, -1, 0, 1, -1, -9, 8, 28, -21, -35, 20, 15, -5, -1

Offset: 0

Philippe Deléham, Dec 27 2007

Previous name: Triangle T(n,k), 0<=k<=n, read by rows given by [0, 1, 0, 0, 0, 0, 0, 0, 0, ...] DELTA [1, -2, 1, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938.

			Triangle begins:
  1;
  0, 1;
  0, 1, -1;
  0, 1, -1, -1;
  0, 1, -1, -2, 1;
  0, 1, -1, -3, 2, 1;
  0, 1, -1, -4, 3, 3, -1;
  0, 1, -1, -5, 4, 6, -3, -1;
  0, 1, -1, -6, 5, 10, -6, -4, 1;
  0, 1, -1, -7, 6, 15, -10, -10, 4, 1;
  0, 1, -1, -8, 7, 21, -15, -20, 10, 5, -1;
  0, 1, -1, -9, 8, 28, -21, -35, 20, 15, -5, -1;
  0, 1, -1, -10, 9, 36, -28, -56, 35, 35, -15, -6, 1;
  ...
Triangle A103631 begins:
  1;
  0, 1;
  0, 1, 1;
  0, 1, 1, 1;
  0, 1, 1, 2, 1;
  0, 1, 1, 3, 2, 1;
  0, 1, 1, 4, 3, 3, 1;
  0, 1, 1, 5, 4, 6, 3, 1;
  0, 1, 1, 6, 5, 10, 6, 4, 1;
  0, 1, 1, 7, 6, 15, 10, 10, 4, 1;
  0, 1, 1, 8, 7, 21, 15, 20, 10, 5, 1;
  0, 1, 1, 9, 8, 28, 21, 35, 20, 15, 5, 1;
  0, 1, 1, 10, 9, 36, 28, 56, 35, 35, 15, 6, 1;
  ...
Triangle A108299 begins:
  1;
  1, -1;
  1, -1, -1;
  1, -1, -2, 1;
  1, -1, -3, 2, 1;
  1, -1, -4, 3, 3, -1;
  1, -1, -5, 4, 6, -3, -1;
  1, -1, -6, 5, 10, -6, -4, 1;
  1, -1, -7, 6, 15, -10, -10, 4, 1;
  1, -1, -8, 7, 21, -15, -20, 10, 5, -1;
  1, -1, -9, 8, 28, -21, -35, 20, 15, -5, -1;
  1, -1, -10, 9, 36, -28, -56, 35, 35, -15, -6, 1;
  ...

Another version is A108299.

Unsigned version is A103631 (T(n,k) = A103631(n,k)*A057077(k)).

m = 13
(* DELTA is defined in A084938 *)
DELTA[Join[{0, 1}, Table[0, {m}]], Join[{1, -2, 1}, Table[0, {m}]], m] // Flatten (* Jean-François Alcover, Feb 19 2020 *)
qStirling2[n_, k_, q_] /; 1 <= k <= n := q^(k-1) qStirling2[n-1, k-1, q] + Sum[q^j, {j, 0, k-1}] qStirling2[n-1, k, q];
qStirling2[n_, 0, _] := KroneckerDelta[n, 0];
qStirling2[0, k_, _] := KroneckerDelta[0, k];
qStirling2[, , _] = 0;
Table[qStirling2[n, k, -1], {n, 0, 11}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 10 2020 *)

from sage.combinat.q_analogues import q_stirling_number2
for n in (0..9):
    print([q_stirling_number2(n,k).substitute(q=-1) for k in [0..n]])
# Peter Luschny, Mar 09 2020

Sum_{k, 0<=k<=n}T(n,k)*x^(n-k)= A057077(n), A010892(n), A000012(n), A001519(n), A001835(n), A004253(n), A001653(n), A049685(n-1), A070997(n-1), A070998(n-1), A072256(n), A078922(n), A077417(n-1), A085260(n), A001570(n-1) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14 respectively .
Sum_{k, 0<=k<=n}T(n,k)*x^k = A000007(n), A010892(n), A133631(n), A133665(n), A133666(n), A133667(n), A133668(n), A133669(n), A133671(n), A133672(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively .
G.f.: (1-x+y*x)/(1-x+y^2*x^2). - Philippe Deléham, Mar 14 2012
T(n,k) = T(n-1,k) - T(n-2,k-2), T(0,0) = T(1,1) = T(2,1) = 1, T(1,0) = T(2,0) = 0, T(2,2) = -1 and T(n,k) = 0 if k<0 or if k>n. - Philippe Deléham, Mar 14 2012

New name from Peter Luschny, Mar 09 2020

A244419 Coefficient triangle of polynomials related to the Dirichlet kernel. Rising powers. Riordan triangle ((1+z)/(1+z^2), 2*z/(1+z^2)).

Original entry on oeis.org

1, 1, 2, -1, 2, 4, -1, -4, 4, 8, 1, -4, -12, 8, 16, 1, 6, -12, -32, 16, 32, -1, 6, 24, -32, -80, 32, 64, -1, -8, 24, 80, -80, -192, 64, 128, 1, -8, -40, 80, 240, -192, -448, 128, 256, 1, 10, -40, -160, 240, 672, -448, -1024, 256, 512, -1, 10, 60, -160, -560, 672, 1792, -1024, -2304, 512, 1024

Offset: 0

Wolfdieter Lang, Jul 29 2014

This is the row reversed version of A180870. See also A157751 and A228565.
The Dirichlet kernel is D(n,x) = Sum_{k=-n..n} exp(i*k*x) = 1 + 2*Sum_{k=1..n} T(n,x) = S(n, 2*y) + S(n-1, 2*y) = S(2*n, sqrt(2*(1+y))) with y = cos(x), n >= 0, with the Chebyshev polynomials T (A053120) and S (A049310). This triangle T(n, k) gives in row n the coefficients of the polynomial Dir(n,y) = D(n,x=arccos(y)) = Sum_{m=0..n} T(n,m)*y^m. See A180870, especially the Peter Bala comments and formulas.
This is the Riordan triangle ((1+z)/(1+z^2), 2*z/(1+z^2)) due to the o.g.f. for Dir(n,y) given by (1+z)/(1 - 2*y*z + z^2) = G(z)/(1 - y*F(z)) with G(z) = (1+z)/(1+z^2) and F(z) = 2*z/(1+z^2) (see the Peter Bala formula under A180870). For Riordan triangles and references see the W. Lang link 'Sheffer a- and z- sequences' under A006232.
The A- and Z- sequences of this Riordan triangle are (see the mentioned W. Lang link in the preceding comment also for the references): The A-sequence has o.g.f. 1+sqrt(1-x^2) and is given by A(2*k+1) = 0 and A(2*k) [2, -1/2, -1/8, -1/16, -5/128, -7/256, -21/1024, -33/2048, -429/32768, -715/65536, ...], k >= 0. (See A098597 and A046161.)
  The Z-sequence has o.g.f. sqrt((1-x)/(1+x)) and is given by
  [1, -1, 1/2, -1/2, 3/8, -3/8, 5/16, -5/16, 35/128, -35/128, ...]. (See A001790 and A046161.)
The column sequences are A057077, 2*(A004526 with even numbers signed), 4*A008805 (signed), 8*A058187 (signed), 16*A189976 (signed), 32*A189980 (signed) for m = 0, 1, ..., 5.
The row sums give A005408 (from the o.g.f. due to the Riordan property), and the alternating row sums give A033999.
The row polynomials Dir(n, x), n >= 0, give solutions to the diophantine equation (a + 1)*X^2 - (a - 1)*Y^2 = 2 by virtue of the identity (a + 1)*Dir(n, -a)^2 - (a - 1)*Dir(n, a)^2 = 2, which is easily proved inductively using the recurrence Dir(n, a) = (1 + a)*(-1)^(n-1)*Dir(n-1, -a) + a*Dir(n-1, a) given below by Wolfdieter Lang. - Peter Bala, May 08 2025

			The triangle T(n,m) begins:
  n\m  0   1   2    3    4    5    6     7     8    9    10 ...
  0:   1
  1:   1   2
  2:  -1   2   4
  3:  -1  -4   4    8
  4:   1  -4 -12    8   16
  5:   1   6 -12  -32   16   32
  6:  -1   6  24  -32  -80   32   64
  7:  -1  -8  24   80  -80 -192   64   128
  8:   1  -8 -40   80  240 -192 -448   128   256
  9:   1  10 -40 -160  240  672 -448 -1024   256  512
  10: -1  10  60 -160 -560  672 1792 -1024 -2304  512  1024
  ...
Example for A-sequence recurrence: T(3,1) = Sum_{j=0..2} A(j)*T(2,j) = 2*(-1) + 0*2 + (-1/2)*4 = -4. Example for Z-sequence recurrence: T(4,0) = Sum_{j=0..3} Z(j)*T(3,j) = 1*(-1) + (-1)*(-4) + (1/2)*4 + (-1/2)*8 = +1. (For the A- and Z-sequences see a comment above.)
Example for the alternate recurrence: T(4,2) = 2*T(3,1) - T(3,2) = 2*(-4) - 4 = -12. T(4,3) = 0*T(3,2) + T(3,3) = T(3,3) = 8. - _Wolfdieter Lang_, Jul 30 2014

Wikipedia, Dirichlet kernel.

Dir(n, x) : A005408 (x = 1), A002878 (x = 3/2), A001834 (x = 2), A030221 (x = 5/2), A002315 (x = 3), A033890 (x = 7/2), A057080 (x = 4), A057081 (x = 9/2), A054320 (x = 5), A077416 (x = 6), A028230 (x = 7), A159678 (x = 8), A049629 (x = 9), A083043 (x = 10),
(-1)^n * Dir(n, x): A122367 (x = -3/2); A079935 (x = -2), A004253 (x = -5/2), A001653 (x = -3), A049685 (x = -7/2), A070997 (x = -4), A070998 (x = -9/2), A072256(n+1) (x = -5).
Cf. A180870 (row reversed), A157751, A228565, A049310, A053120, A057077, A004526, A008805, A058187, A189976, A189980, A005408, A033999.

T[n_, k_] := T[n, k] = Which[k == 0, (-1)^Quotient[n, 2], (0 <= n && n < k) || (n == -1 && k == 1), 0, True, 2 T[n-1, k-1] - T[n-2, k]];
Table[T[n, k], {n, 0, 11}, {k, 0, n}] // Flatten (* Jean-François Alcover, Jun 28 2019, from Sage *)

def T(n, k):
    if k == 0: return (-1)^(n//2)
    if (0 <= n and n < k) or (n == -1 and k == 1): return 0
    return 2*T(n-1, k-1) - T(n-2, k)
for n in range(11): [T(n,k) for k in (0..n)] # Peter Luschny, Jul 29 2014

T(n, m) = [y^m] Dir(n,y) for n >= m >= 0 and 0 otherwise, with the polynomials Dir(y) defined in a comment above.
T(n, m) = 2^m*(S(n,m) + S(n-1,m)) with the entries S(n,m) of A049310 given there explicitly.
O.g.f. for polynomials Dir(y) see a comment above (Riordan triangle ((1+z)/(1+z^2), 2*z/(1+z^2))).
O.g.f. for column m: ((1 + x)/(1 + x^2))*(2*x/(1 + x^2))^m, m >= 0, (Riordan property).
Recurrence for the polynomials: Dir(n, y) = 2*y*Dir(n-1, y) - Dir(n-2, y), n >= 1, with input D(-1, y) = -1 and D(0, y) = 1.
Triangle three-term recurrence: T(n,m) = 2*T(n-1,m-1) - T(n-2,m) for n >= m >= 1 with T(n,m) = 0 if 0 <= n < m, T(0,0) = 1, T(-1,1) = 0 and T(n,0) = A057077(n) = (-1)^(floor(n/2)).
From Wolfdieter Lang, Jul 30 2014: (Start)
In analogy to A157751 one can derive a recurrence for the row polynomials Dir(n, y) = Sum_{m=0..n} T(n,m)*y^m also using a negative argument but only one recursive step: Dir(n,y) = (1+y)*(-1)^(n-1)*Dir(n-1,-y) + y*Dir(n-1,y), n >= 1, Dir(0,y) = 1 (Dir(-1,y) = -1). See also A180870 from where this formula can be derived by row reversion.
This entails another triangle recurrence T(n,m) = (1 + (-1)^(n-m))*T(n-1,m-1) - (-1)^(n-m)*T(n-1,m), for n >= m >= 1 with T(n,m) = 0 if n < m and T(n,0) = (-1)^floor(n/2). (End)
From Peter Bala, Aug 14 2022: (Start)
The row polynomials Dir(n,x), n >= 0, are related to the Chebyshev polynomials of the first kind T(n,x) by the binomial transform as follows:
(2^n)*(x - 1)^(n+1)*Dir(n,x) = (-1) * Sum_{k = 0..2*n+1} binomial(2*n+1,k)*T(k,-x).
Note that Sum_{k = 0..2*n} binomial(2*n,k)*T(k,x) = (2^n)*(1 + x)^n*T(n,x). (End)
From Peter Bala, May 04 2025: (Start)
For n >= 1, the n-th row polynomial Dir(n, x) = (-1)^n * (U(n, -x) - U(n-1, -x)) = U(2*n, sqrt((1+x)/2)), where U(n, x) denotes the n-th Chebyshev polynomial of the second kind.
For n >= 1 and x < 1, Dir(n, x) = (-1)^n * sqrt(2/(1 - x )) * T(2*n+1, sqrt((1 - x)/2)), where T(n, x) denotes the n-th Chebyshev polynomial of the first kind.
Dir(n, x)^2 - 2*x*Dir(n, x)*Dir(n+1, x) + Dir(n+1, x)^2 = 2*(1 + x).
Dir(n, x) = (-1)^n * R(n, -2*(x+1)), where R(n, x) is the n-th row polynomial of the triangle A085478.
Dir(n, x) = Sum_{k = 0..n} (-1)^(n+k) * binomial(n+k, 2*k) * (2*x + 2)^k. (End)

A041007 Denominators of continued fraction convergents to sqrt(6).

Original entry on oeis.org

1, 2, 9, 20, 89, 198, 881, 1960, 8721, 19402, 86329, 192060, 854569, 1901198, 8459361, 18819920, 83739041, 186298002, 828931049, 1844160100, 8205571449, 18255302998, 81226783441, 180708869880

Offset: 0

N. J. A. Sloane

sqrt(6) = 4/2 + 4/9 + 4/(9*89) + 4/(89*881) + 4/(881*8721), ...; where sqrt(6) = 2.4494897427... and the sum of the first 5 terms of this series = 2.449489737... - Gary W. Adamson, Dec 21 2007
sqrt(6) = 2 + continued fraction [2, 4, 2, 4, 2, 4, ...] = 4/2 + 4/9 + 4/(9*89) + 4/(89*881) + 4/(881*8721) + ... - Gary W. Adamson, Dec 21 2007
Interspersion of 2 sequences, A072256 and 2*A004189. - Gerry Martens, Jun 10 2015
For n > 0, a(n) equals the permanent of the n X n tridiagonal matrix with the main diagonal alternating sequence [2, 4, 2, 4, ...] and 1's along the superdiagonal and the subdiagonal. - Rogério Serôdio, Apr 01 2018

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (0,10,0,-1).

Cf. A010464, A041006.

Cf. A072256, A004189.

Table[Denominator[FromContinuedFraction[ContinuedFraction[Sqrt[6],n]]],{n,1,50}] (* Vladimir Joseph Stephan Orlovsky, Mar 16 2011 *)

G.f.: (1+2*x-x^2)/(1-10*x^2+x^4). - Colin Barker, Dec 31 2011
From Rogério Serôdio, Apr 01 2018: (Start)
Recurrence formula: a(n) = (3 + (-1)^n)*a(n-1) + a(n-2), a(0) = 1, a(1) = 2.
Some properties:
(1) a(n)^2 - a(n-2)^2 = (3+(-1)^n)*a(2*n-1), for n > 1;
(2) a(2*n+1) = a(n)*(a(n+1) + a(n-1)), for n > 0;
(3) a(2*n) = A142239(2*n), for n >= 0;
(4) a(2*n+1) = A041007(2*n+1)/2, for n >= 0;
(5) a(2*n-1)*A142239(2*n+1) = a(n)^2 - 1, for n > 0;
(6) a(2*n) = a(n)*A142239(n) + a(n-1)*A142239(n-1), for n > 0;
(7) Sum_{k=0..n} a(2*k+1)*(A142239(2*k) + A142239(2*(k+1))) = Sum_{k=0..n} a(3+4*k);
(8) Sum_{k=0..n} (a(2*k-1) + a(2*k+1))*A142239(2*k) = Sum_{k=0..n} A142239(3+4*k). (End)
a(n) = ((2 + sqrt(6))^(n+1) - (2 - sqrt(6))^(n+1))/(sqrt(6) * 2^(ceiling(n/2) + 1)). - Robert FERREOL, Oct 14 2024

A105038 Nonnegative n such that 6n^2 + 6n + 1 is a square.

Original entry on oeis.org

0, 4, 44, 440, 4360, 43164, 427284, 4229680, 41869520, 414465524, 4102785724, 40613391720, 402031131480, 3979697923084, 39394948099364, 389969783070560, 3860302882606240, 38213059042991844, 378270287547312204, 3744489816430130200, 37066627876753989800

Offset: 0

Gerald McGarvey, Apr 03 2005

Harvey P. Dale, Table of n, a(n) for n = 0..1000
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (I).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (II).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (III).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (IV).
Index entries for linear recurrences with constant coefficients, signature (11,-11,1).

Cf. A001652, A001570, A049629, A105040, A104240, A077288, A105036, A103200, A105037.

CoefficientList[Series[4x/(1-11x+11x^2-x^3),{x,0,30}],x] (* or *) LinearRecurrence[{11,-11,1},{0,4,44},30] (* Harvey P. Dale, Sep 29 2013 *)

for(n=0,427284,if(issquare(6*n*(n+1)+1),print1(n,",")))

Vec(4*x/(1-11*x+11*x^2-x^3)+O(x^99)) \\ Charles R Greathouse IV, Nov 13 2012

G.f.: 4*x/(1-11*x+11*x^2-x^3).
a(0)=0, a(1)=4, a(2)=44, a(n)=11*a(n-1)-11*a(n-2)+a(n-3). - Harvey P. Dale, Sep 29 2013
a(n) = (-6-(5-2*sqrt(6))^n*(-3+sqrt(6))+(3+sqrt(6))*(5+2*sqrt(6))^n)/12. - Colin Barker, Mar 05 2016
a(n) = (A072256(n+1) - 1)/2.

More terms from Vladeta Jovovic, Apr 05 2005

Incorrect conjectures deleted by N. J. A. Sloane, Nov 24 2010

A046172 Indices of pentagonal numbers (A000326) that are also squares (A000290).

Original entry on oeis.org

1, 81, 7921, 776161, 76055841, 7452696241, 730288175761, 71560788528321, 7012226987599681, 687126683996240401, 67331402804643959601, 6597790348171111800481, 646516122717964312487521, 63351982236012331511976561, 6207847743006490523861215441

Offset: 1

Eric W. Weisstein

if P_x = y^2 is a pentagonal number that is also a square, the least both pentagonal and square number that is greater as P_x, is P_(49*x + 40*y - 8) = (60*x + 49*y - 10)^2 (in fact, P_(49*x + 40*y - 8) - (60*x + 49*y - 10)^2 = (3/2)*x^2 - (1/2)*x - y^2). - Richard Choulet, Apr 28 2009
a(n)*(3*a(n)-1)/2 = m^2 is equivalent to the Pell equation (6*a(n)-1)^2 - 6*(2*m)^2 = 1 or x(n)^2 - 6*y(n)^2 = 1. - Paul Weisenhorn, May 15 2009
As n increases, this sequence is approximately geometric with common ratio r = lim_{n -> oo} a(n)/a(n-1) = (sqrt(2) + sqrt(3))^4 = 49 + 20*sqrt(6). - Ant King, Nov 07 2011
Numbers k such that the k-th pentagonal number is equal to the sum of two consecutive triangular numbers. - Colin Barker, Dec 11 2014
Indices of pentagonal numbers (A000326) that are also centered octagonal numbers (A016754). - Colin Barker, Jan 11 2015

Muniru A. Asiru, All square chiliagonal numbers, International Journal of Mathematical Education in Science and Technology, Volume 47, 2016 - Issue 7; http://dx.doi.org/10.1080/0020739X.2016.1164346

Colin Barker, Table of n, a(n) for n = 1..503
Leonhard Euler, De solutione problematum diophanteorum per numeros integros, section 21.
Giovanni Lucca, Integer Sequences and Circle Chains Inside a Circular Segment, Forum Geometricorum, Vol. 18 (2018), 47-55.
W. Sierpiński, Sur les nombres pentagonaux, Bull. Soc. Roy. Sci. Liege 33 (1964) 513-517.
Eric Weisstein's World of Mathematics, Pentagonal Square Number.
Index entries for linear recurrences with constant coefficients, signature (99,-99,1).

Cf. A036353, A046173.

Cf. A000217, A000290, A000326, A251914, A248205.

LinearRecurrence[{99, -99, 1}, {1, 81, 7921}, 13] (* Ant King, Nov 07 2011 *)
Table[Round[(1 + x^(2*n+1))^2 / (12*x^(2*n+1)) /. x->5+2*Sqrt@6],{n,0,99}] (* Federico Provvedi, Apr 24 2023 *)

a(n) = 98*a(n-1) - a(n-2) - 16; g.f.: x*(1 - 18*x + x^2)/((1-x)*(1 - 98*x + x^2)). - Warut Roonguthai Jan 05 2001 - Corrected by Colin Barker, Jan 11 2015
a(n+1) = 49*a(n) - 8 + 10*sqrt(8*(3*a(n)^2 - a(n))) with a(1) = 1. - Richard Choulet, Apr 28 2009
a(n) = 1/6+((5 + 2*sqrt(6))^(2*n+1)/12) + ((5 - 2*sqrt(6))^(2*n+1)/12) for n >= 0. - Richard Choulet, Apr 29 2009
From Paul Weisenhorn, May 15 2009: (Start)
x(n+2) = 98*x(n+1) - x(n) with x(1)=5, x(2)=485;
y(n+2) = 98*y(n+1) - y(n) with y(n)=A046173(n)*2;
m(n+2) = 98*m(n+1) - m(n) with m(n)=A046173(n);
a(n) = A072256(n)^2.
(End)
a(n) = b(n)*b(n), b(n) = 10*b(n-1)- b(n-2), b(1)=1, b(2)=9, b(n)=((5 + sqrt(24))^n - (5 - sqrt(24))^n)/(2*sqrt(24)). - Sture Sjöstedt, Sep 21 2009
From Ant King, Nov 07 2011: (Start)
a(n) = 99*a(n-1) - 99*a(n-2) + a(n-3).
a(n) = ceiling((1/12)*(sqrt(3) + sqrt(2))^(4*n-2)).
(End)
a(n) = (1 + x^(2n+1))^2 / (12*x^(2*n+1)), with x = 5 + 2*sqrt(6). - Federico Provvedi, Apr 24 2023

A200993 Triangular numbers, T(m), that are two-thirds of another triangular number; T(m) such that 3T(m) = 2T(k) for some k.

Original entry on oeis.org

0, 10, 990, 97020, 9506980, 931587030, 91286021970, 8945098566040, 876528373449960, 85890835499530050, 8416425350580494950, 824723793521388975060, 80814515339745539060940, 7918997779501541438997070, 775980967875811315482651930, 76038215854050007375860892080

Offset: 0

Charlie Marion, Dec 15 2011

For n>1, a(n) = 98*a(n-1) - a (n-2) + 10. In general, for m>0, let b(n) be those triangular numbers such that for some triangular number c(n), (m+1)*b(n) = m*c(n). Then b(0) = 0, b(1)= A014105(m) and for n>1, b(n) = 2*A069129(m+1)*b(n-1) - b(n-2) + A014105(m).

Further, c(0) = 0, c(1) = A000384(m+1) and for n>1, c(n) = 2*A069129(m+1)*c(n-1) - c(n-2) + A000384(m+1).

			3*0 = 2*0.
3*10 = 2*15.
3*990 = 2*1485.
3*97020 = 2*145530.

Colin Barker, Table of n, a(n) for n = 0..500
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (I).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (II).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (III).
Editors, L'Intermédiaire des Mathématiciens, Query 4500: The equation x(x+1)/2 = y*(y+1)/3, L'Intermédiaire des Mathématiciens, 22 (1915), 255-260 (IV).
Index entries for linear recurrences with constant coefficients, signature (99,-99,1).

Cf. A001652, A029549, A053141, A075528, A200994-A201008.

m:=30; R:=PowerSeriesRing(Integers(), m); [0] cat Coefficients(R!(10*x/((1-x)*(x^2-98*x+1)))); // G. C. Greubel, Jul 15 2018

LinearRecurrence[{99,-99,1},{0,10,990},20] (* Harvey P. Dale, Feb 25 2018 *)

concat(0, Vec(10*x/((1-x)*(1-98*x+x^2)) + O(x^40))) \\ Colin Barker, Mar 02 2016

G.f. 10*x / ((1-x)*(x^2-98*x+1)). - R. J. Mathar, Dec 20 2011
a(n) = 99*a(n-1)-99*a(n-2)+a(n-3) for n>2. - Colin Barker, Mar 02 2016
a(n) = (-10+(5-2*sqrt(6))*(49+20*sqrt(6))^(-n)+(5+2*sqrt(6))*(49+20*sqrt(6))^n)/96. - Colin Barker, Mar 07 2016

A237610 Positive integers k such that x^2 - 10xy + y^2 + k = 0 has integer solutions.

Original entry on oeis.org

8, 15, 20, 23, 24, 32, 47, 60, 71, 72, 80, 87, 92, 95, 96, 116, 128, 135, 152, 159, 167, 180, 188, 191, 200, 207, 212, 215, 216, 239, 240, 263, 276, 284, 288, 303, 311, 320, 335, 344, 348, 359, 368, 375, 380, 383, 384, 392, 404, 423, 431, 447, 456, 464, 479

Offset: 1

Colin Barker, Feb 10 2014

nonn

			15 is in the sequence because x^2 - 10xy + y^2 + 15 = 0 has integer solutions, for example (x, y) = (2, 19).

Cf. A072256 (k = 8), A129445 (k = 15), A080806 (k = 20), A074061 (k = 23), A001079 (k = 24).

Cf. A031363, A084917, A237351, A237599, A237606, A237609.

is(n)=m=bnfisintnorm(bnfinit(x^2-10*x+1),-n);#m>0&&denominator(polcoeff(m[1],1))==1 \\ Ralf Stephan, Feb 11 2014

cp's OEIS Frontend

A138288 a(n) = A054320(n) - A001078(n).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Sage

Formula

A160682 The list of the A values in the common solutions to 13*k+1 = A^2 and 17*k+1 = B^2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

Extensions

A077416 Chebyshev S-sequence with Diophantine property.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Sage

Formula

A133607 Triangle read by rows: T(n, k) = qStirling2(n, k, q) for q = -1, with 0 <= k <= n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Sage

Formula

Extensions

A244419 Coefficient triangle of polynomials related to the Dirichlet kernel. Rising powers. Riordan triangle ((1+z)/(1+z^2), 2*z/(1+z^2)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Sage

Formula

A041007 Denominators of continued fraction convergents to sqrt(6).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

Formula

A105038 Nonnegative n such that 6*n^2 + 6*n + 1 is a square.

Views

Author

Keywords

Links

A160682 The list of the A values in the common solutions to 13k+1 = A^2 and 17k+1 = B^2.

A105038 Nonnegative n such that 6n^2 + 6n + 1 is a square.

A200993 Triangular numbers, T(m), that are two-thirds of another triangular number; T(m) such that 3T(m) = 2T(k) for some k.