A001541 -id:A001541 - cp's OEIS frontend

A086645 Triangle read by rows: T(n, k) = binomial(2n, 2k).

1, 1, 1, 1, 6, 1, 1, 15, 15, 1, 1, 28, 70, 28, 1, 1, 45, 210, 210, 45, 1, 1, 66, 495, 924, 495, 66, 1, 1, 91, 1001, 3003, 3003, 1001, 91, 1, 1, 120, 1820, 8008, 12870, 8008, 1820, 120, 1, 1, 153, 3060, 18564, 43758, 43758, 18564, 3060, 153, 1, 1, 190, 4845, 38760

Offset: 0

Philippe Deléham, Jul 26 2003

Terms have the same parity as those of Pascal's triangle.
Coefficients of polynomials (1/2)*((1 + x^(1/2))^(2n) + (1 - x^(1/2))^(2n)).
Number of compositions of 2n having k parts greater than 1; example: T(3, 2) = 15 because we have 4+2, 2+4, 3+2+1, 3+1+2, 2+3+1, 2+1+3, 1+3+2, 1+2+3, 2+2+1+1, 2+1+2+1, 2+1+1+2, 1+2+2+1, 1+2+1+2, 1+1+2+2, 3+3. - Philippe Deléham, May 18 2005
Number of binary words of length 2n - 1 having k runs of consecutive 1's; example: T(3,2) = 15 because we have 00101, 01001, 01010, 01011, 01101, 10001, 10010, 10011, 10100, 10110, 10111, 11001, 11010, 11011, 11101. - Philippe Deléham, May 18 2005
Let M_n be the n X n matrix M_n(i, j) = T(i, j-1); then for n > 0, det(M_n) = A000364(n), Euler numbers; example: det([1, 1, 0, 0; 1, 6, 1, 0; 1, 15, 15, 1; 1, 28, 70, 28 ]) = 1385 = A000364(4). - Philippe Deléham, Sep 04 2005
Equals ConvOffsStoT transform of the hexagonal numbers, A000384: (1, 6, 15, 28, 45, ...); e.g., ConvOffs transform of (1, 6, 15, 28) = (1, 28, 70, 28, 1). - Gary W. Adamson, Apr 22 2008
From Peter Bala, Oct 23 2008: (Start)
Let C_n be the root lattice generated as a monoid by {+-2*e_i: 1 <= i <= n; +-e_i +- e_j: 1 <= i not equal to j <= n}. Let P(C_n) be the polytope formed by the convex hull of this generating set. Then the rows of this array are the h-vectors of a unimodular triangulation of P(C_n) [Ardila et al.]. See A127674 for (a signed version of) the corresponding array of f-vectors for these type C_n polytopes. See A008459 for the array of h-vectors for type A_n polytopes and A108558 for the array of h-vectors associated with type D_n polytopes.
The Hilbert transform of this triangle is A142992 (see A145905 for the definition of this term).
(End)
Diagonal sums: A108479. - Philippe Deléham, Sep 08 2009
Coefficients of Product_{k=1..n} (cot(k*Pi/(2n+1))^2 - x) = Sum_{k=0..n} (-1)^k*binomial(2n,2k)*x^k/(2n+1-2k). - David Ingerman (daviddavifree(AT)gmail.com), Mar 30 2010
Generalized Narayana triangle for 4^n (or cosh(2x)). - Paul Barry, Sep 28 2010
Coefficients of the matrix inverse appear to be T^(-1)(n,k) = (-1)^(n+k)*A086646(n,k). - R. J. Mathar, Mar 12 2013
Let E(y) = Sum_{n>=0} y^n/(2*n)! = cosh(sqrt(y)). Then this triangle is the generalized Riordan array (E(y), y) with respect to the sequence (2*n)! as defined in Wang and Wang. Cf. A103327. - Peter Bala, Aug 06 2013
Row 6, (1,66,495,924,495,66,1), plays a role in expansions of powers of the Dedekind eta function. See the Chan link, p. 534, and A034839. - Tom Copeland, Dec 12 2016

			From _Peter Bala_, Oct 23 2008: (Start)
The triangle begins
n\k|..0.....1.....2.....3.....4.....5.....6
===========================================
0..|..1
1..|..1.....1
2..|..1.....6.....1
3..|..1....15....15.....1
4..|..1....28....70....28.....1
5..|..1....45...210...210....45.....1
6..|..1....66...495...924...495....66.....1
...
(End)
From _Peter Bala_, Aug 06 2013: (Start)
Viewed as the generalized Riordan array (cosh(sqrt(y)), y) with respect to the sequence (2*n)! the column generating functions begin
1st col: cosh(sqrt(y)) = 1 + y/2! + y^2/4! + y^3/6! + y^4/8! + ....
2nd col: 1/2!*y*cosh(sqrt(y)) = y/2! + 6*y^2/4! + 15*y^3/6! + 28*y^4/8! + ....
3rd col: 1/4!*y^2*cosh(sqrt(y)) = y^2/4! + 15*y^3/6! + 70*y^4/8! + 210*y^5/10! + .... (End)

A. T. Benjamin and J. J. Quinn, Proofs that really count: the art of combinatorial proof, M.A.A. 2003, id. 224.

Indranil Ghosh, Rows 0.. 120 of triangle, flattened
F. Ardila, M. Beck, S. Hosten, J. Pfeifle and K. Seashore, Root polytopes and growth series of root lattices arXiv:0809.5123 [math.CO], 2008.
H. Chan, S. Cooper and P. Toh, The 26th power of Dedekind's eta function Advances in Mathematics, 207 (2006) 532-543.
T. Copeland, Infinitesimal Generators, the Pascal Pyramid, and the Witt and Virasoro Algebras, 2012.
T. Copeland, Skipping over Dimensions, Juggling Zeros in the Matrix, 2020.
E. Lucas, Théorie des fonctions numériques simplement périodiques, Baltimore, 1878. See page 145 equation (153).
W. Wang and T. Wang, Generalized Riordan array, Discrete Mathematics, Vol. 308, No. 24, 6466-6500.

Cf. A098158, A000384, A081294.
Cf. A008459, A108558, A127674, A142992. - Peter Bala, Oct 23 2008
Cf. A103327 (binomial(2n+1, 2k+1)), A103328 (binomial(2n, 2k+1)), A091042 (binomial(2n+1, 2k)). -Wolfdieter Lang, Jan 06 2013
Cf. A086646 (unsigned matrix inverse), A103327.
Cf. A034839.

/* As triangle: */ [[Binomial(2*n, 2*k): k in [0..n]]: n in [0.. 15]]; // Vincenzo Librandi, Dec 14 2016

A086645:=(n,k)->binomial(2*n,2*k): seq(seq(A086645(n,k),k=0..n),n=0..12);

Table[Binomial[2 n, 2 k], {n, 0, 10}, {k, 0, n}] // Flatten (* Michael De Vlieger, Dec 13 2016 *)

create_list(binomial(2*n,2*k),n,0,12,k,0,n); /* Emanuele Munarini, Mar 11 2011 */

PARI
```
{T(n, k) = binomial(2*n, 2*k)};
```

{T(n, k) = sum( i=0, min(k, n-k), 4^i * binomial(n, 2*i) * binomial(n - 2*i, k-i))}; /* Michael Somos, May 26 2005 */

T(n, k) = (2*n)!/((2*(n-k))!*(2*k)!) row sums = A081294. COLUMNS: A000012, A000384
Sum_{k>=0} T(n, k)*A000364(k) = A000795(n) = (2^n)*A005647(n).
Sum_{k>=0} T(n, k)*2^k = A001541(n). Sum_{k>=0} T(n, k)*3^k = 2^n*A001075(n). Sum_{k>=0} T(n, k)*4^k = A083884(n). - Philippe Deléham, Feb 29 2004
O.g.f.: (1 - z*(1+x))/(x^2*z^2 - 2*x*z*(1+z) + (1-z)^2) = 1 + (1 + x)*z +(1 + 6*x + x^2)*z^2 + ... . - Peter Bala, Oct 23 2008
Sum_{k=0..n} T(n,k)*x^k = A000007(n), A081294(n), A001541(n), A090965(n), A083884(n), A099140(n), A099141(n), A099142(n), A165224(n), A026244(n) for x = 0,1,2,3,4,5,6,7,8,9 respectively. - Philippe Deléham, Sep 08 2009
Product_{k=1..n} (cot(k*Pi/(2n+1))^2 - x) = Sum_{k=0..n} (-1)^k*binomial(2n,2k)*x^k/(2n+1-2k). - David Ingerman (daviddavifree(AT)gmail.com), Mar 30 2010
From Paul Barry, Sep 28 2010: (Start)
G.f.: 1/(1-x-x*y-4*x^2*y/(1-x-x*y)) = (1-x*(1+y))/(1-2*x*(1+y)+x^2*(1-y)^2);
E.g.f.: exp((1+y)*x)*cosh(2*sqrt(y)*x);
T(n,k) = Sum_{j=0..n} C(n,j)*C(n-j,2*(k-j))*4^(k-j). (End)
T(n,k) = 2*T(n-1,k) + 2*T(n-1,k-1) + 2*T(n-2,k-1) - T(n-2,k) - T(n-2,k-2), with T(0,0)=T(1,0)=T(1,1)=1, T(n,k)=0 if k<0 or if k>n. - Philippe Deléham, Nov 26 2013
From Peter Bala, Sep 22 2021: (Start)
n-th row polynomial R(n,x) = (1-x)^n*T(n,(1+x)/(1-x)), where T(n,x) is the n-th Chebyshev polynomial of the first kind. Cf. A008459.
R(n,x) = Sum_{k = 0..n} binomial(n,2*k)*(4*x)^k*(1 + x)^(n-2*k).
R(n,x) = n*Sum_{k = 0..n} (n+k-1)!/((n-k)!*(2*k)!)*(4*x)^k*(1-x)^(n-k) for n >= 1. (End)

A046090 Consider all Pythagorean triples (X,X+1,Z) ordered by increasing Z; sequence gives X+1 values.

Original entry on oeis.org

1, 4, 21, 120, 697, 4060, 23661, 137904, 803761, 4684660, 27304197, 159140520, 927538921, 5406093004, 31509019101, 183648021600, 1070379110497, 6238626641380, 36361380737781, 211929657785304, 1235216565974041, 7199369738058940, 41961001862379597, 244566641436218640

Offset: 0

Eric W. Weisstein

Solution to a*(a-1) = 2b*(b-1) in natural numbers: a = a(n), b = b(n) = A011900(n).
n such that n^2 = (1/2)*(n+floor(sqrt(2)*n*floor(sqrt(2)*n))). - Benoit Cloitre, Apr 15 2003
Place a(n) balls in an urn, of which b(n) = A011900(n) are red; draw 2 balls without replacement; 2*Probability(2 red balls) = Probability(2 balls); this is equivalent to the Pell equation A(n)^2-2*B(n)^2 = -1 with a(n) = (A(n)+1)/2; b(n) = (B(n)+1)/2; and the fundamental solution (7;5) and the solution (3;2) for the unit form. - Paul Weisenhorn, Aug 03 2010
Find base x in which repdigit yy has a square that is repdigit zzzz, corresponding to Diophantine equation zzzz_x = (yy_x)^2; then, solution z = a(n) with x = A002315(n) and y = A001653(n+1) for n >= 1 (see Maurice Protat reference). - Bernard Schott, Dec 21 2022

			For n=4: a(4)=697; b(4)=493; 2*binomial(493,2)=485112=binomial(697,2). - _Paul Weisenhorn_, Aug 03 2010

A. H. Beiler, Recreations in the Theory of Numbers. New York: Dover, pp. 122-125, 1964.
Maurice Protat, Des Olympiades à l'Agrégation, De zzzz_x = (yy_x)^2 à Pell-Fermat, Problème 23, pp. 52-54, Ellipses, Paris, 1997.

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
T. W. Forget and T. A. Larkin, Pythagorean triads of the form X, X+1, Z described by recurrence sequences, Fib. Quart., 6 (No. 3, 1968), 94-104.
L. J. Gerstein, Pythagorean triples and inner products, Math. Mag., 78 (2005), 205-213.
H. J. Hindin, Stars, hexes, triangular numbers and Pythagorean triples, J. Rec. Math., 16 (1983/1984), 191-193. (Annotated scanned copy)
Ron Knott, Pythagorean Triples and Online Calculators
S. Northshield, An Analogue of Stern's Sequence for Z[sqrt(2)], Journal of Integer Sequences, 18 (2015), #15.11.6.
Eric Weisstein's World of Mathematics, Pythagorean Triple
Index entries for two-way infinite sequences
Index entries for linear recurrences with constant coefficients, signature (7,-7,1).

Other 2 sides are A001652 and A001653.
Cf. A001009, A001541, A001542, A011900, A046729, A053141, A084159, A084703, A156035.
See comments in A301383.
Cf. A001109, A001653, A002024, A002315, A029549.

a046090 n = a046090_list !! n
a046090_list = 1 : 4 : map (subtract 2)
   (zipWith (-) (map (* 6) (tail a046090_list)) a046090_list)
-- Reinhard Zumkeller, Jan 10 2012

m:=30; R:=PowerSeriesRing(Integers(), m); Coefficients(R!((1-3*x)/((1-6*x+x^2)*(1-x)))); // G. C. Greubel, Jul 15 2018

Digits:=100: seq(round((1+(7+5*sqrt(2))*(3+2*sqrt(2))^(n-1))/2)/2, n=0..20); # Paul Weisenhorn, Aug 03 2010

Join[{1},#+1&/@With[{c=3+2Sqrt[2]},NestList[Floor[c #]+3&,3,20]]] (* Harvey P. Dale, Aug 19 2011 *)
LinearRecurrence[{7,-7,1},{1,4,21},25] (* Harvey P. Dale, Apr 13 2012 *)
a[n_] := (2-ChebyshevT[n, 3]+ChebyshevT[n+1, 3])/4; Array[a, 21, 0] (* Jean-François Alcover, Jul 10 2016, adapted from PARI *)

a(n)=(2-subst(poltchebi(abs(n))-poltchebi(abs(n+1)),x,3))/4

x='x+O('x^30); Vec((1-3*x)/((1-6*x+x^2)*(1-x))) \\ G. C. Greubel, Jul 15 2018

a(n) = (-1+sqrt(1+8*b(n)*(b(n)+1)))/2 with b(n) = A011900(n). [corrected by Michel Marcus, Dec 23 2022]
a(n) = 6*a(n-1) - a(n-2) - 2, n >= 2, a(0) = 1, a(1) = 4.
a(n) = (A(n+1) - 3*A(n) + 2)/4 with A(n) = A001653(n).
A001652(n) = -a(-1-n).
From Barry E. Williams, May 03 2000: (Start)
G.f.: (1-3*x)/((1-6*x+x^2)*(1-x)).
a(n) = partial sums of A001541(n). (End)
From Charlie Marion, Jul 01 2003: (Start)
A001652(n)*A001652(n+1) + a(n)*a(n+1) = A001542(n+1)^2 = A084703(n+1).
Let a(n) = A001652(n), b(n) = this sequence and c(n) = A001653(n). Then for k > j, c(i)*(c(k) - c(j)) = a(k+i) + ... + a(i+j+1) + a(k-i-1) + ... + a(j-i) + k - j. For n < 0, a(n) = -b(-n-1). Also a(n)*a(n+2k+1) + b(n)*b(n+2k+1) + c(n)*c(n+2k+1) = (a(n+k+1) - a(n+k))^2; a(n)*a(n+2k) + b(n)*b(n+2k) + c(n)*c(n+2k) = 2*c(n+k)^2. (End)
a(n) = 1/2 + ((1-2^(1/2))/4)*(3 - 2^(3/2))^n + ((1+2^(1/2))/4)*(3 + 2^(3/2))^n. - Antonio Alberto Olivares, Oct 13 2003
2*a(n) = 2*A084159(n) + 1 + (-1)^(n+1) = 2*A046729(n) + 1 - (-1)^(n+1). - Lekraj Beedassy, Jul 16 2004
a(n) = A001109(n+1) - A053141(n). - Manuel Valdivia, Apr 03 2010
From Paul Weisenhorn, Aug 03 2010: (Start)
a(n+1) = round((1+(7+5*sqrt(2))*(3+2*sqrt(2))^n)/2);
b(n+1) = round((2+(10+7*sqrt(2))*(3+2*sqrt(2))^n)/4) = A011900(n+1).
(End)
a(n)*(a(n)-1)/2 = b(n)*b(n+1) and 2*a(n) - 1 = b(n) + b(n+1), where b(n) = A001109. - Kenneth J Ramsey, Apr 24 2011
T(a(n)) = A011900(n)^2 + A001109(n), where T(n) is the n-th triangular number. See also A001653. - Charlie Marion, Apr 25 2011
a(0)=1, a(1)=4, a(2)=21, a(n) = 7*a(n-1) - 7*a(n-2) + a(n-3). - Harvey P. Dale, Apr 13 2012
Limit_{n->oo} a(n+1)/a(n) = 3 + 2*sqrt(2) = A156035. - Ilya Gutkovskiy, Jul 10 2016
a(n) = A001652(n)+1. - Dimitri Papadopoulos, Jul 06 2017
a(n) = (A002315(n) + 1)/2. - Bernard Schott, Dec 21 2022
E.g.f.: (exp(x) + exp(3*x)*(cosh(2*sqrt(2)*x) + sqrt(2)*sinh(2*sqrt(2)*x)))/2. - Stefano Spezia, Mar 16 2024
a(n) = A002024(A029549(n))+1. - Pontus von Brömssen, Sep 11 2024

Additional comments from Wolfdieter Lang

Comment moved to A001653 by Claude Morin, Sep 22 2023

A055792 a(n) and floor(a(n)/2) are both squares; i.e., squares which remain squares when written in base 2 and last digit is removed.

Original entry on oeis.org

0, 1, 9, 289, 9801, 332929, 11309769, 384199201, 13051463049, 443365544449, 15061377048201, 511643454094369, 17380816062160329, 590436102659356801, 20057446674355970889, 681362750825443653409, 23146276081390728245001, 786292024016459316676609

Offset: 0

Henry Bottomley, Jul 14 2000

a(n) > 0 is a square such that a(n) - 1 is a product of powers. - Michel Lagneau, Feb 16 2012

			a(2) = 9 because 9 = 3^2 = 1001_2 and 100_2 = 4 = 2^2.

Charles R Greathouse IV, Table of n, a(n) for n = 0..654
M. F. Hasler, Truncated squares, OEIS wiki, Jan 16 2012
Giovanni Lucca, Integer Sequences and Circle Chains Inside a Circular Segment, Forum Geometricorum, Vol. 18 (2018), 47-55.
Giovanni Lucca, Circle chains inside the arbelos and integer sequences, Int'l J. Geom. (2023) Vol. 12, No. 1, 71-82.
Index to sequences related to truncating digits of squares.
Index entries for linear recurrences with constant coefficients, signature (35,-35,1).

Cf. A023110, A247375.

LinearRecurrence[{35, -35, 1}, {0, 1, 9, 289}, 25] (* Paolo Xausa, Jul 22 2024 *)

concat(0, Vec(-x*(9*x^2-26*x+1)/((x-1)*(x^2-34*x+1)) + O(x^100))) \\ Colin Barker, Sep 15 2014

is(n)=issquare(n) && issquare(n\2) \\ Charles R Greathouse IV, May 07 2015

a(n) = 34*a(n-1) - a(n-2) - 16 = A001541(n-1)^2 = 2*A001542(n-1)^2 + 1 = 8*A001110(n-1) + 1.
From Colin Barker, Sep 15 2014: (Start)
a(n) = 35*a(n-1) - 35*a(n-2) + a(n-3) for n > 3.
G.f.: -x*(9*x^2 - 26*x + 1) / ((x-1)*(x^2 - 34*x + 1)). (End)
a(n) = c*k^n + 1/2 + o(1) with k = 17+sqrt(288) = 33.97... and c = 17/4 - sqrt(18). - Charles R Greathouse IV, May 07 2015
a(n) = (4 + 2*(17 + 12*sqrt(2))^(1-n) + (34 - 24*sqrt(2))*(17 + 12*sqrt(2))^n)/8 for n > 0. - Colin Barker, Mar 02 2016

A023110 Squares which remain squares when the last digit is removed.

Original entry on oeis.org

0, 1, 4, 9, 16, 49, 169, 256, 361, 1444, 3249, 18496, 64009, 237169, 364816, 519841, 2079364, 4678569, 26666896, 92294449, 341991049, 526060096, 749609641, 2998438564, 6746486769, 38453641216, 133088524969, 493150849009, 758578289296, 1080936581761

Offset: 1

David W. Wilson

This A023110 = A031149^2 is the base 10 version of A001541^2 = A055792 (base 2), A001075^2 = A055793 (base 3), A004275^2 = A055808 (base 4), A204520^2 = A055812 (base 5), A204518^2 = A055851 (base 6), A204516^2 = A055859 (base 7), A204514^2 = A055872 (base 8) and A204502^2 = A204503 (base 9). - M. F. Hasler, Sep 28 2014

For the first 4 terms the square has only one digit. It is understood that deleting this digit yields 0. - Colin Barker, Dec 31 2017

R. K. Guy, Neg and Reg, preprint, Jan 2012.

Jon E. Schoenfield, Table of n, a(n) for n = 1..70 (terms 1..38 from David W. Wilson, terms 39..40 from Robert G. Wilson v, terms 41..67 from Dmitry Petukhov)
M. F. Hasler, Truncated squares, OEIS wiki, Jan 16 2012
Joshua Stucky, Pell's Equation and Truncated Squares, Number Theory Seminar, Kansas State University, Feb 19 2018.
Index to sequences related to truncating digits of squares.

Cf. A023111.
Cf. A031150, A053784, A031149, A055792, A055793, A055808, A055812, A055851, A055859, A055872.
Cf. A001541, A001075, A004275, A204520, A204518, A204516, A204514, A204502, A204503.

count:= 1: A[1]:= 0:
for n from 0 while count < 35 do
  for t in [1,4,6,9] do
    if issqr(10*n^2+t) then
       count:= count+1;
       A[count]:= 10*n^2+t;
    fi
  od
od:
seq(A[i],i=1..count); # Robert Israel, Sep 28 2014

fQ[n_] := IntegerQ@ Sqrt@ Quotient[n^2, 10]; Select[ Range@ 1000000, fQ]^2 (* Robert G. Wilson v, Jan 15 2011 *)

for(n=0,1e7, issquare(n^2\10) & print1(n^2",")) \\  M. F. Hasler, Jan 16 2012

Appears to satisfy a(n)=1444*a(n-7)+a(n-14)-76*sqrt(a(n-7)*a(n-14)) for n >= 16. For n = 15, 14, 13, ... this would require a(1) = 16, a(0) = 49, a(-1) = 169, ... - Henry Bottomley, May 08 2001; edited by Robert Israel, Sep 28 2014
a(n) = A031149(n)^2. - M. F. Hasler, Sep 28 2014
Conjectures from Colin Barker, Dec 31 2017: (Start)
G.f.: x^2*(1 + 4*x + 9*x^2 + 16*x^3 + 49*x^4 + 169*x^5 + 256*x^6 - 1082*x^7 - 4328*x^8 - 9738*x^9 - 4592*x^10 - 6698*x^11 - 6698*x^12 - 4592*x^13 + 361*x^14 + 1444*x^15 + 3249*x^16 + 256*x^17 + 169*x^18 + 49*x^19 + 16*x^20) / ((1 - x)*(1 + x + x^2 + x^3 + x^4 + x^5 + x^6)*(1 - 1442*x^7 + x^14)).
a(n) = 1443*a(n-7) - 1443*a(n-14) + a(n-21) for n>22.
(End)

More terms from M. F. Hasler, Jan 16 2012

A038761 a(n) = 6*a(n-1) - a(n-2), n >= 2, a(0)=1, a(1)=9.

Original entry on oeis.org

1, 9, 53, 309, 1801, 10497, 61181, 356589, 2078353, 12113529, 70602821, 411503397, 2398417561, 13979001969, 81475594253, 474874563549, 2767771787041, 16131756158697, 94022765165141, 548004834832149, 3194006243827753, 18616032628134369, 108502189524978461

Offset: 0

Barry E. Williams, May 02 2000

Bisection of A048654. - Lambert Klasen (lambert.klasen(AT)gmx.de), Nov 24 2004

This gives part of the (increasingly sorted) positive solutions y to the Pell equation x^2 - 2*y^2 = +7. For the x solutions see A038762. For the other part of solutions see A101386 and A253811. - Wolfdieter Lang, Feb 05 2015

			A038762(3)^2 - 2*a(4)^2 = 2547^2 - 2*1801^2 = +7. - _Wolfdieter Lang_, Feb 05 2015

A. H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, pp. 122-125, 194-196.

Vincenzo Librandi, Table of n, a(n) for n = 0..400
M. J. DeLeon, Pell's Equation and Pell Number Triples, Fib. Quart., 14(1976), pp. 456-460.
Tanya Khovanova, Recursive Sequences
Index entries for linear recurrences with constant coefficients, signature (6,-1).
Index entries for sequences related to Chebyshev polynomials.

Cf. A001109, A001541, A001542, A001652, A001653, A038762, A046090, A048654, A055997, A101386, A124124, A253811.

I:=[1, 9]; [n le 2 select I[n] else 6*Self(n-1)-Self(n-2): n in [1..40]]; // Vincenzo Librandi, Nov 16 2011

a[0]:=1: a[1]:=9: for n from 2 to 26 do a[n]:=6*a[n-1]-a[n-2] od: seq(a[n], n=0..19); # Zerinvary Lajos, Jul 26 2006

LinearRecurrence[{6,-1},{1,9},40] (* Vincenzo Librandi, Nov 16 2011 *)

a(n)=([0,1; -1,6]^n*[1;9])[1,1] \\ Charles R Greathouse IV, Apr 08 2016

a(n) = (9*((3+2*sqrt(2))^n -(3-2*sqrt(2))^n)-((3+2*sqrt(2))^(n-1) - (3-2*sqrt(2))^(n-1)))/(4*sqrt(2)).
a(n) = sqrt(2*(A038762(n))^2-14)/2.
For n>1, a(n)-4a(n-1)=A001541(n)-A001542(n-2); e.g. 309-4*53=97=99-2. - Charlie Marion, Nov 12 2003
For n>0, a(n)=A046090(n)+A001653(n)+A001652(n-1)=A055997(n+1)+A001652(n-1); e.g., 309=120+169+20. - Charlie Marion, Oct 11 2006
G.f.: (1+3*x)/(1-6*x+x^2). - Philippe Deléham, Nov 03 2008
a(n) = third binomial transform of 1,6,8,48,64,384. - Al Hakanson (hawkuu(AT)gmail.com), Aug 15 2009
a(n)^2 + 2^2 = A124124(2*n+1)^2 + (A124124(2*n+1)+1)^2. - Hermann Stamm-Wilbrandt, Aug 31 2014
a(n) = irrational part of z(n) = (3 + sqrt(2))*(3 + 2*sqrt(2))^n, n >= 0. z(n) gives only part of the general positive solutions to the Pell equation x^2 - 2*y^2 = 7. See the Nagell reference in A038762 on how to find z(n), and a comment above. - Wolfdieter Lang, Feb 05 2015
a(n) = S(n, 6) + 3*S(n-1, 6), n >= 0, with the Chebyshev S-polynomials evaluated at x=6. See S(n-1, 6) = A001109(n). - Wolfdieter Lang, Mar 30 2015
E.g.f.: exp(3*x)*(2*cosh(2*sqrt(2)*x) + 3*sqrt(2)*sinh(2*sqrt(2)*x))/2. - Stefano Spezia, Mar 16 2024

Edited: Replaced the unspecific Pell comment. Moved a formula from the comment section to the formula section. - Wolfdieter Lang, Feb 05 2015

A132592 X-values of solutions to the equation X(X + 1) - 8Y^2 = 0.

Original entry on oeis.org

0, 8, 288, 9800, 332928, 11309768, 384199200, 13051463048, 443365544448, 15061377048200, 511643454094368, 17380816062160328, 590436102659356800, 20057446674355970888, 681362750825443653408, 23146276081390728245000, 786292024016459316676608, 26710782540478226038759688

Offset: 0

Mohamed Bouhamida, Nov 14 2007

Equivalently, numbers k such that both k/2 and k+1 are squares. - Karl-Heinz Hofmann, Sep 20 2022

Equivalently, numbers k such that the k-dimensional volume and total (k-1)-dimensional volume are equal, with side length being a positive integer, for all regular polyhedra constructible in k dimensions. - Matt Moir, Jul 09 2024

Seiichi Manyama, Table of n, a(n) for n = 0..500
Index entries for linear recurrences with constant coefficients, signature (35,-35,1).

Intersection between A132411 and A001105.
Cf. A007654.
Cf. A001541, A058331, A001079, A037270, A055792, A071253, A108741, A132592, A146311, A146312, A146313, A173115, A173116, A173121.

I:=[0,8,288]; [n le 3 select I[n] else 35*Self(n-1)-35*Self(n-2)+ Self(n-3): n in [1..30]]; // Vincenzo Librandi, Dec 24 2018

Table[Round[N[Sinh[2 n ArcCosh[Sqrt[2]]]^2, 100]], {n, 0, 20}] (* Artur Jasinski, Feb 10 2010 *)
LinearRecurrence[{35, -35, 1}, {0, 8, 288}, 30] (* Vincenzo Librandi, Dec 24 2018 *)

A132592 = [0, 8]
for n in range(2, 18): A132592.append(34 * A132592[-1] - A132592[-2] + 16)
print(A132592) # Karl-Heinz Hofmann, Sep 20 2022

a(0)=0, a(1)=8 and a(n) = 34*a(n-1) - a(n-2) + 16.
a(n) = (A056771(n) - 1)/2. - Max Alekseyev, Nov 13 2009
a(n) = sinh(2*n*arccosh(sqrt(2))^2) (n=0,1,2,3,...). - Artur Jasinski, Feb 10 2010
G.f.: -8*x*(x+1)/((x-1)*(x^2-34*x+1)). - Colin Barker, Oct 24 2012
a(n) = A055792(n+1)-1 = A001541(n)^2 - 1. - Antti Karttunen, Oct 03 2016

More terms from Max Alekseyev, Nov 13 2009

A079291 Squares of Pell numbers.

Original entry on oeis.org

0, 1, 4, 25, 144, 841, 4900, 28561, 166464, 970225, 5654884, 32959081, 192099600, 1119638521, 6525731524, 38034750625, 221682772224, 1292061882721, 7530688524100, 43892069261881, 255821727047184, 1491038293021225

Offset: 0

Ralf Stephan, Feb 08 2003

(-1)^(n+1)*a(n) is the r=-4 member of the r-" of sequences S_r(n), n>=1, defined in A092184 where more information can be found.
Binomial transform of A086346. - Johannes W. Meijer, Aug 01 2010
In general, squaring the terms of a Horadam sequence with signature (c,d) will result in a third-order recurrence with signature (c^2+d, c^2*d+d^2, -d^3). - Gary Detlefs, Nov 11 2021
(Conjectured) For any primitive Pythagorean triple of the form (X, Y, Z=Y+1), it appears that Y or Z will always be (and only be) a square Pell number if X = A001333(n), for n > 1. If n is even, Y is always a square Pell number, and if n is odd, then Z is always a square Pell number. For example: (3, 4, 5), (7, 24, 25), (17, 144, 145), (41, 840, 841), (99, 4900, 4901). - Jules Beauchamp, Feb 02 2022
a(n+1) is the number of tilings of an n-board (a board with dimensions n X 1) using (1/2,1/2)-fences, black half-squares (1/2 X 1 pieces, always placed so that the shorter sides are horizontal), and white half-squares. A (w,g)-fence is a tile composed of two w X 1 pieces separated by a gap of width g. a(n+1) also equals the number of tilings of an n-board using black (1/4,1/4)-fences, white (1/4,1/4)-fences, and (1/4,3/4)-fences. - Michael A. Allen, Dec 29 2022

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Michael A. Allen and Kenneth Edwards, Fence tiling derived identities involving the metallonacci numbers squared or cubed, Fib. Q. 60:5 (2022) 5-17.
Joerg Arndt, Matters Computational (The Fxtbook), sect. 32.1.5, pp. 626-627
Sergio Falcon, Some series of reciprocal k-Fibonacci numbers, Asian Journal of Mathematics and Computer Research, Vol. 11, No. 3 (2016), pp. 184-191; ResearchGate link.
Toufik Mansour, A note on sum of k-th power of Horadam's sequence, arXiv:math/0302015 [math.CO], 2003.
Toufik Mansour, Squaring the terms of an ell-th order linear recurrence, arXiv:math/0303138 [math.CO], 2003.
Pantelimon Stanica, Generating functions, weighted and non-weighted sums for powers of second-order recurrence sequences, arXiv:math/0010149 [math.CO], 2000.
Index entries for linear recurrences with constant coefficients, signature (5,5,-1).
Index entries for sequences related to Chebyshev polynomials.

Cf. A000129, A001254, A001109, A001541, A001653, A174968.

Cf. A002203, A007598, A084158 (partial sums), A086346, A092184.

I:=[0,1,4]; [n le 3 select I[n] else 5*Self(n-1)+ 5*Self(n-2) - Self(n-3): n in [1..31]]; // Vincenzo Librandi, May 17 2013

with(combinat):seq(fibonacci(i,2)^2, i=0..31); # Zerinvary Lajos, Mar 20 2008

CoefficientList[Series[x(1-x)/((1+x)*(1-6x+x^2)), {x, 0, 30}], x] (* Vincenzo Librandi, May 17 2013 *)
LinearRecurrence[{5,5,-1},{0,1,4},40] (* Harvey P. Dale, Dec 20 2015 *)
Fibonacci[Range[0, 30], 2]^2 (* G. C. Greubel, Sep 17 2021 *)

[lucas_number1(n, 2, -1)^2 for n in (0..30)] # G. C. Greubel, Sep 17 2021

G.f.: x*(1-x)/((1+x)*(1-6*x+x^2)).
a(n) = (r^n + (1/r)^n - 2*(-1)^n)/8, with r = 3 + sqrt(8).
a(n+3) = 5*a(n+2) + 5*a(n+1) - a(n).
L.g.f.: (1/8)*log((1+2*x+x^2)/(1-6*x+x^2)) = Sum_{n>=0} (a(n)/n)*x^n, see p. 627 of the Fxtbook link; special case of the following: let v(0)=0, v(1)=1, and v(n) = u*v(n-1) + v(n-2), then (1/A)*log((1+2*x+x^2)/(1-(2-A)*x+x^2)) = Sum_{n>=0} v(n)^2/n*x^n where A = u^2 + 4. - Joerg Arndt, Apr 08 2011
a(n+1) = Sum_{k=0..n} ( (-1)^(n-k)*A001653(k) ); e.g., 144 = -1 + 5 - 29 + 169; 25 = 1 - 5 + 29. - Charlie Marion, Jul 16 2003
a(n) = A000129(n)^2.
a(n) = (T(n, 3) - (-1)^n)/4 with Chebyshev's polynomials of the first kind evaluated at x=3: T(n, 3) = A001541(n) = ((3 + 2*sqrt(2))^n + (3 - 2*sqrt(2))^n )/2. - Wolfdieter Lang, Oct 18 2004
a(n) is the rightmost term of M^n * [1 0 0] where M is the 3 X 3 matrix [4 4 1 / 2 1 0 / 1 0 0]. a(n+1) = leftmost term. E.g., a(6) = 4900, a(5) = 841 since M^5 * [1 0 0] = [4900 2030 841]. - Gary W. Adamson, Oct 31 2004
a(n) = ( (-1)^(n+1) + A001109(n+1) - 3*A001109(n) )/4. - R. J. Mathar, Nov 16 2007
a(n) = ( (((1 - sqrt(2))^n + (1 + sqrt(2))^n) /2 )^2 + (-1)^(n+1) )/2. - Antonio Pane (apane1(AT)spc.edu), Dec 15 2007
Lim_{k -> infinity} ( a(n+k)/a(k) ) = A001541(n) + 2*A001109(n)*sqrt(2). - Johannes W. Meijer, Aug 01 2010
For n>0, a(2*n) = 6*a(2*n-1) - a(2*n-2) - 2, a(2*n+1) = 6*a(2*n) - a(2*n-1) + 2. - Charlie Marion, Sep 24 2011
a(n) = (1/8)*(A002203(2*n) - 2*(-1)^n). - G. C. Greubel, Sep 17 2021
Conjectured formula for (X, Y, Z) for primitive Pythagorean triple of the form (X, Y, Z=Y+1) is (A001333(n)^2, A079291(n)^2, A079291(n)^2-1) or (A001333(n)^2, A079291(n)^2-1, A079291(n)^2). As a closed formula (X, Y, Z) = (((1-sqrt(2))^n + (1+sqrt(2))^n)/2, (((1-sqrt(2))^n + (1+sqrt(2))^n)^2 - 4)/8, (((1-sqrt(2))^n + (1+sqrt(2))^n)^2 + 4)/8). - Jules Beauchamp, Feb 02 2022
From Michael A. Allen, Dec 29 2022: (Start)
a(n+1) = 6*a(n) - a(n-1) + 2*(-1)^n.
a(n+1) = (1 + (-1)^n)/2 + 4*Sum_{k=1..n} ( k*a(n+1-k) ). (End)
Product_{n>=2} (1 + (-1)^n/a(n)) = (1 + sqrt(2))/2 (A174968) (Falcon, 2016, p. 189, eq. (3.1)). - Amiram Eldar, Dec 03 2024

A008310 Triangle of coefficients of Chebyshev polynomials T_n(x).

Original entry on oeis.org

1, 1, -1, 2, -3, 4, 1, -8, 8, 5, -20, 16, -1, 18, -48, 32, -7, 56, -112, 64, 1, -32, 160, -256, 128, 9, -120, 432, -576, 256, -1, 50, -400, 1120, -1280, 512, -11, 220, -1232, 2816, -2816, 1024, 1, -72, 840, -3584, 6912, -6144, 2048, 13, -364, 2912, -9984, 16640, -13312, 4096

Offset: 0

N. J. A. Sloane

The row length sequence of this irregular array is A008619(n), n >= 0. Even or odd powers appear in increasing order starting with 1 or x for even or odd row numbers n, respectively. This is the standard triangle A053120 with 0 deleted. - Wolfdieter Lang, Aug 02 2014

Let T* denote the triangle obtained by replacing each number in this triangle by its absolute value. Then T* gives the coefficients for cos(nx) as a polynomial in cos x. - Clark Kimberling, Aug 04 2024

			Rows are: (1), (1), (-1,2), (-3,4), (1,-8,8), (5,-20,16) etc., since if c = cos(x): cos(0x) = 1, cos(1x) = 1c; cos(2x) = -1+2c^2; cos(3x) = -3c+4c^3, cos(4x) = 1-8c^2+8c^4, cos(5x) = 5c-20c^3+16c^5, etc.
From _Wolfdieter Lang_, Aug 02 2014: (Start)
This irregular triangle a(n,k) begins:
  n\k   0    1     2      3      4      5      6      7 ...
  0:    1
  1:    1
  2:   -1    2
  3:   -3    4
  4:    1   -8     8
  5:    5  -20    16
  6:   -1   18   -48     32
  7:   -7   56  -112     64
  8:    1  -32   160   -256    128
  9:    9 -120   432   -576    256
 10:   -1   50  -400   1120  -1280    512
 11:  -11  220 -1232   2816  -2816   1024
 12:    1  -72   840  -3584   6912  -6144   2048
 13:   13 -364  2912  -9984  16640 -13312   4096
 14:   -1   98 -1568   9408 -26880  39424 -28672   8192
 15:  -15  560 -6048  28800 -70400  92160 -61440  16384
  ...
T(4,x) = 1 - 8*x^2 + 8*x^4, T(5,x) = 5*x - 20*x^3 +16*x^5.
(End)

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 795.
E. A. Guilleman, Synthesis of Passive Networks, Wiley, 1957, p. 593.
Yaroslav Zolotaryuk, J. Chris Eilbeck, "Analytical approach to the Davydov-Scott theory with on-site potential", Physical Review B 63, p543402, Jan. 2001. The authors write, "Since the algebra of these is 'hyperbolic', contrary to the usual Chebyshev polynomials defined on the interval 0 <= x <= 1, we call the set of functions (21) the hyperbolic Chebyshev polynomials." (This refers to the triangle T* described in Comments.)

R. J. Mathar, Table of n, a(n) for n = 0..2600
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
Renato Ferreira Pinto Jr. and Nathaniel Harms, Testing Support Size More Efficiently Than Learning Histograms, arXiv:2410.18915 [cs.DS], 2024. See p. 40.
D. Foata and G.-N. Han, Nombres de Fibonacci et polynomes orthogonaux
C. Lanczos, Applied Analysis (Annotated scans of selected pages)
I. Rivin, Growth in free groups (and other stories), arXiv:math/9911076 [math.CO], 1999.
Eric Weisstein's World of Mathematics, Chebyshev Polynomial of the First Kind.
Wikipedia, Chebyshev polynomials.
Index entries for sequences related to Chebyshev polynomials.

A039991 is a row reversed version, but has zeros which enable the triangle to be seen. Columns/diagonals are A011782, A001792, A001793, A001794, A006974, A006975, A006976 etc.
Reflection of A028297. Cf. A008312, A053112.
Row sums are one. Polynomial evaluations include A001075 (x=2), A001541 (x=3), A001091, A001079, A023038, A011943, A001081, A023039, A001085, A077422, A077424, A097308, A097310, A068203.
Cf. A053120.

A008310 := proc(n,m) local x ; coeftayl(simplify(ChebyshevT(n,x),'ChebyshevT'),x=0,m) ; end: i := 0 : for n from 0 to 100 do for m from n mod 2 to n by 2 do printf("%d %d ",i,A008310(n,m)) ; i := i+1 ; od ; od ; # R. J. Mathar, Apr 20 2007
# second Maple program:
b:= proc(n) b(n):= `if`(n<2, 1, expand(2*b(n-1)-x*b(n-2))) end:
T:= n-> (p-> (d-> seq(coeff(p, x, d-i), i=0..d))(degree(p)))(b(n)):
seq(T(n), n=0..15);  # Alois P. Heinz, Sep 04 2019

Flatten[{1, Table[CoefficientList[ChebyshevT[n, x], x], {n, 1, 13}]}]//DeleteCases[#, 0, Infinity]& (* or *) Flatten[{1, Table[Table[((-1)^k*2^(n-2 k-1)*n*Binomial[n-k, k])/(n-k), {k, Floor[n/2], 0, -1}], {n, 1, 13}]}] (* Eugeniy Sokol, Sep 04 2019 *)

a(n,m) = 2^(m-1) * n * (-1)^((n-m)/2) * ((n+m)/2-1)! / (((n-m)/2)! * m!) if n>0. - R. J. Mathar, Apr 20 2007
From Paul Weisenhorn, Oct 02 2019: (Start)
T_n(x) = 2*x*T_(n-1)(x) - T_(n-2)(x), T_0(x) = 1, T_1(x) = x.
T_n(x) = ((x+sqrt(x^2-1))^n + (x-sqrt(x^2-1))^n)/2. (End)
From Peter Bala, Aug 15 2022: (Start)
T(n,x) = [z^n] ( z*x + sqrt(1 + z^2*(x^2 - 1)) )^n.
Sum_{k = 0..2*n} binomial(2*n,k)*T(k,x) = (2^n)*(1 + x)^n*T(n,x).
exp( Sum_{n >= 1} T(n,x)*t^n/n ) = Sum_{n >= 0} P(n,x)*t^n, where P(n,x) denotes the n-th Legendre polynomial. (End)

Additional comments and more terms from Henry Bottomley, Dec 13 2000

Edited: Corrected Cf. A039991 statement. Cf. A053120 added. - Wolfdieter Lang, Aug 06 2014

A001601 a(n) = 2*a(n-1)^2 - 1, if n>1. a(0)=1, a(1)=3.

Original entry on oeis.org

1, 3, 17, 577, 665857, 886731088897, 1572584048032918633353217, 4946041176255201878775086487573351061418968498177

Offset: 0

N. J. A. Sloane

Reduced numerators of Newton's iteration for sqrt(2). - Eric W. Weisstein
An infinite coprime sequence defined by recursion. - Michael Somos, Mar 14 2004
Evaluation of the 2^n - 1 degree interpolating polynomial of 1/x at Chebyshev nodes in the interval (1,2): v = 1.0; for(i = 1, n, v *= 4*(a(i) - x*v)); v *= 2/a(n+1). - Jose Hortal, Apr 07 2012
Smallest positive integer x satisfying the Pell equation x^2 - 2^(2*n+1) * y^2 = 1, for n > 0. - A.H.M. Smeets, Sep 29 2017

L. E. Dickson, History of the Theory of Numbers. Carnegie Institute Public. 256, Washington, DC, Vol. 1, 1919; Vol. 2, 1920; Vol. 3, 1923, see vol. 1, p. 376.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Dennis Martin, Table of n, a(n) for n = 0..11
Dov Jarden, Recurring Sequences, Riveon Lematematika, Jerusalem, 1966. [Annotated scanned copy]
M. Mendes France and A. J. van der Poorten, From geometry to Euler identities, Theoret. Comput. Sci., 65 (1989), 213-220.
Jeffrey Shallit, Rational numbers with non-terminating, non-periodic modified Engel-type expansions, Fib. Quart., 31 (1993), 37-40.
Eric Weisstein's World of Mathematics, Newtons Iteration.
Eric Weisstein's World of Mathematics, Square Root.
Eric Weisstein's World of Mathematics, Pythagoras's Constant.
H. S. Wilf, D. C. B. Marsh and J. V. Whittaker, Problem E1093, Amer. Math. Monthly, 61 (1954), 424-425.
Index entries for sequences related to Engel expansions

Cf. A051009, A003423.

a(n) = A001333(2^n).

Table[Simplify[Expand[(1/2) ((1 + Sqrt[2])^(2^n) + (1 - Sqrt[2])^(2^n))]], {n, 0, 7}]  (* Artur Jasinski, Oct 10 2008 *)
Join[{1},NestList[2#^2-1&,3,7]]  (* Harvey P. Dale, Mar 24 2011 *)

PARI
```
a(n)=if(n<1,n==0,2*a(n-1)^2-1)
```

From Mario Catalani (mario.catalani(AT)unito.it), May 27 2003, May 30 2003: (Start)
a(n) = a(n-1)^2 + 2*A051009(n)^2 for n > 0.
a(n)^2 = 2*A051009(n+1)^2 + 1.
a(n) = Sum_{r=0..2^(n-1)} binomial(2^n, 2*r)*2^r. (End)
Expansion of 1/sqrt(2) as an infinite product: 1/sqrt(2) = Product_{k>=1} (1 - 1/(a(n)+1)). a(1)=3; a(n) = floor(1/(1-1/(sqrt(2)*Product_{k=1..n-1} 1-1/(a(k)+1)))). - Thomas Baruchel, Nov 06 2003
2*a(n+1) = A003423(n).
a(n) = (1/2)*((1 + sqrt(2))^(2^n) + (1 - sqrt(2))^(2^n)). - Artur Jasinski, Oct 10 2008
For n > 1: a(n) - 1 = 4^n * Product_{i=1..n-2} a(i)^2. - Jose Hortal, Apr 13 2012
From Peter Bala, Nov 11 2012: (Start)
4*sqrt(2)/7 = Product_{n>=1} (1 - 1/(2*a(n))).
sqrt(2) = Product_{n>=1} (1 + 1/a(n)).
a(n) = (1/2)*A003423(n-1). (End)
a(n) = cos(2^(n-1) * arccos(3)) = cosh(2^(n-1) * log(3 + 2*sqrt(2))) for n >= 1. - Daniel Suteu, Jul 28 2016
a(n+1) = T(2^n,3), where T(n,x) denotes the n-th Chebyshev polynomial of the first kind. - Peter Bala, Feb 01 2017
a(n) = A001541(2^(n-1)). - A.H.M. Smeets, May 28 2017

A005319 a(n) = 6*a(n-1) - a(n-2).

Original entry on oeis.org

0, 4, 24, 140, 816, 4756, 27720, 161564, 941664, 5488420, 31988856, 186444716, 1086679440, 6333631924, 36915112104, 215157040700, 1254027132096, 7309005751876, 42600007379160, 248291038523084, 1447146223759344, 8434586304032980

Offset: 0

N. J. A. Sloane

Solutions y of the equation 2x^2-y^2=2; the corresponding x values are given by A001541. - N-E. Fahssi, Feb 25 2008
The lower intermediate convergents to 2^(1/2) beginning with 4/3, 24/17, 140/99, 816/577, form a strictly increasing sequence; essentially, numerators=A005319 and denominators=A001541. - Clark Kimberling, Aug 26 2008
Numbers n such that (ceiling(sqrt(n*n/2)))^2 = 1 + n*n/2. - Ctibor O. Zizka, Nov 09 2009
All nonnegative solutions of the indefinite binary quadratic form X^2 + 4*X*Y -4*Y^2 of discriminant 32, representing -4 are (X(n), Y(n)) = (a(n), A001653(n+1)), for n >= 0. - Wolfdieter Lang, Jun 13 2018
Also the number of edge covers in the n-triangular snake graph. - Eric W. Weisstein, Jun 08 2019
All of the positive integer solutions of a*b+1=x^2, a*c+1=y^2, b*c+1=z^2, x+z=2*y, 0A001542(n), b=A005319(n), c=A001542(n+1), x=A001541(n), y=A001653(n+1), z=A002315(n) with 0Michael Somos, Jun 26 2022
a(n) is the sum of 4*n consecutive powers of the silver ratio 1+sqrt(2), starting at (1+sqrt(2))^(-2*n) and ending at (1+sqrt(2))^(2*n-1). - Greg Dresden and Ruxin Sheng, Jul 25 2024

			G.f. = 4*x + 24*x^2 + 140*x^3 + 816*x^4 + 4756*x^5 + ... - _Michael Somos_, Jun 26 2022

P. de la Harpe, Topics in Geometric Group Theory, Univ. Chicago Press, 2000, p. 160, middle display.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Marius A. Burtea, Table of n, a(n) for n = 0..100
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.
Hacène Belbachir, Soumeya Merwa Tebtoub, and László Németh, Ellipse Chains and Associated Sequences, J. Int. Seq., Vol. 23 (2020), Article 20.8.5.
John M. Campbell, An Integral Representation of Kekulé Numbers, and Double Integrals Related to Smarandache Sequences, arXiv preprint arXiv:1105.3399 [math.GM], 2011.
Tanya Khovanova, Recursive Sequences
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992
Soumeya M. Tebtoub, Hacène Belbachir, and László Németh, Integer sequences and ellipse chains inside a hyperbola, Proceedings of the 1st International Conference on Algebras, Graphs and Ordered Sets (ALGOS 2020), hal-02918958 [math.cs], 17-18.
Eric Weisstein's World of Mathematics, Edge Cover
Eric Weisstein's World of Mathematics, Triangular Snake Graph
Index entries for linear recurrences with constant coefficients, signature (6,-1).

Cf. A000129, A001109, A001541, A001542, A001652, A001653, A046090, A046729, A090390.

a:=[0,4]; [n le 2 select a[n] else 6*Self(n-1) - Self(n-2):n in [1..22]]; // Marius A. Burtea, Sep 19 2019

LinearRecurrence[{6, -1}, {0, 4}, 22] (* Jean-François Alcover, Sep 26 2017 *)
Table[((3 + 2 Sqrt[2])^n - (3 - 2 Sqrt[2])^n)/Sqrt[2], {n, 20}] // Expand (* Eric W. Weisstein, Jun 08 2019 *)
CoefficientList[Series[(4 x)/(1 - 6 x + x^2), {x, 0, 20}], x] (* Eric W. Weisstein, Jun 08 2019 *)
a[ n_] := 4*ChebyshevU[n-1, 3]; (* Michael Somos, Jun 26 2022 *)

{a(n) = 4*polchebyshev(n-1, 2, 3)}; /* Michael Somos, Jun 26 2022 */

G.f.: 4*x / ( 1-6*x+x^2 ). - Simon Plouffe in his 1992 dissertation.
G.f. for signed version beginning with 1: (1+2*x+x^2)/(1+6*x+x^2).
For any term n of the sequence, 2*n^2 + 4 is a perfect square. Limit_{n->infinity} a(n)/a(n-1) = 3 + 2*sqrt(2). - Gregory V. Richardson, Oct 06 2002
a(n) = ((3+2*sqrt(2))^n - (3-2*sqrt(2))^n) / sqrt(2). - Gregory V. Richardson, Oct 06 2002
(-1)^(n+1) = A090390(n+1) + A001542(n+1) + A046729(n) - a(n) (conjectured). - Creighton Dement, Nov 17 2004
For n > 0, a(n) = A000129(n+1)^2 - A000129(n-1)^2; a(n) = A046090(n-1) + A001652(n); e.g., 816 = 120 + 696; a(n) = A001653(n) - A001653(n-1); e.g., 816 = 985 - 169. - Charlie Marion Jul 22 2005
a(n) = 4*A001109(n). - M. F. Hasler, Mar 2009
For n > 1, a(n) is the denominator of continued fraction [1,4,1,4,...,1,4] with (n-1) repetitions of 1,4. For the numerators, see A001653. - Greg Dresden, Sep 10 2019
1/a(n) - 1/a(n+1) = 1/(Pell(2*n+1) - 1/Pell(2*n+1)) for n >= 1, where Pell(n) = A000129(n). - Peter Bala, Aug 21 2022
E.g.f.: sqrt(2)*exp(3*x)*sinh(2*sqrt(2)*x). - Stefano Spezia, Nov 25 2022
a(n) = 2*A000129(2*n). - Tanya Khovanova and MIT PRIMES STEP senior group, Apr 17 2024

cp's OEIS Frontend

A086645 Triangle read by rows: T(n, k) = binomial(2n, 2k).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

Maxima

PARI

PARI

Formula

A046090 Consider all Pythagorean triples (X,X+1,Z) ordered by increasing Z; sequence gives X+1 values.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Formula

Extensions

A055792 a(n) and floor(a(n)/2) are both squares; i.e., squares which remain squares when written in base 2 and last digit is removed.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

A023110 Squares which remain squares when the last digit is removed.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A038761 a(n) = 6*a(n-1) - a(n-2), n >= 2, a(0)=1, a(1)=9.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

Extensions

A132592 X-values of solutions to the equation X*(X + 1) - 8*Y^2 = 0.

A132592 X-values of solutions to the equation X(X + 1) - 8Y^2 = 0.