cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-5 of 5 results.

A110294 a(2*n) = A028230(n), a(2*n+1) = -A067900(n+1).

Original entry on oeis.org

1, -8, 15, -112, 209, -1560, 2911, -21728, 40545, -302632, 564719, -4215120, 7865521, -58709048, 109552575, -817711552, 1525870529, -11389252680, 21252634831, -158631825968, 296011017105, -2209456310872, 4122901604639, -30773756526240, 57424611447841
Offset: 0

Views

Author

Creighton Dement, Jul 18 2005

Keywords

Comments

See A110293.

Crossrefs

Programs

  • Magma
    [(3*(-1)^n-1)*Evaluate(ChebyshevSecond(n+1), 2)/2: n in [0..40]]; // G. C. Greubel, Jan 04 2023
    
  • Maple
    seriestolist(series((1-8*x+x^2)/((x^2-4*x+1)*(x^2+4*x+1)), x=0,25));
  • Mathematica
    CoefficientList[Series[(1-8x+x^2)/((1-4x+x^2)(1+4x+x^2)), {x, 0, 24}], x] (* Michael De Vlieger, Nov 01 2016 *)
    LinearRecurrence[{0,14,0,-1},{1,-8,15,-112},30] (* Harvey P. Dale, Dec 16 2024 *)
  • PARI
    Vec((1-8*x+x^2)/((1-4*x+x^2)*(1+4*x+x^2)) + O(x^30)) \\ Colin Barker, Nov 01 2016
    
  • SageMath
    [(3*(-1)^n-1)*chebyshev_U(n,2)/2 for n in range(41)] # G. C. Greubel, Jan 04 2023

Formula

G.f.: (1-8*x+x^2) / ((1-4*x+x^2)*(1+4*x+x^2)).
a(n) = 14*a(n-2) - a(n-4) for n>3. - Colin Barker, Nov 01 2016
a(n) = (3*(-1)^n - 1)*A001353(n+1)/2. - R. J. Mathar, Sep 11 2019

A007655 Standard deviation of A007654.

Original entry on oeis.org

0, 1, 14, 195, 2716, 37829, 526890, 7338631, 102213944, 1423656585, 19828978246, 276182038859, 3846719565780, 53577891882061, 746243766783074, 10393834843080975, 144767444036350576, 2016350381665827089, 28084137899285228670, 391161580208327374291, 5448177985017298011404
Offset: 1

Views

Author

Keywords

Comments

a(n) corresponds also to one-sixth the area of Fleenor-Heronian triangle with middle side A003500(n). - Lekraj Beedassy, Jul 15 2002
a(n) give all (nontrivial, integer) solutions of Pell equation b(n+1)^2 - 48*a(n+1)^2 = +1 with b(n+1)=A011943(n), n>=0.
For n>=3, a(n) equals the permanent of the (n-2) X (n-2) tridiagonal matrix with 14's along the main diagonal, and i's along the superdiagonal and the subdiagonal (i is the imaginary unit). - John M. Campbell, Jul 08 2011
For n>1, a(n) equals the number of 01-avoiding words of length n-1 on alphabet {0,1,...,13}. - Milan Janjic, Jan 25 2015
6*a(n)^2 = 6*S(n-1, 14)^2 is the triangular number Tri((T(n, 7) - 1)/2) with Tri = A000217 and T = A053120. This is instance k = 3 of the general k-identity given in a comment to A001109. - Wolfdieter Lang, Feb 01 2016

Examples

			G.f. = x^2 + 14*x^3 + 195*x^4 + 2716*x^5 + 37829*x^6 + 526890*x^7 + ...
		

References

  • D. A. Benaron, personal communication.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Chebyshev sequence U(n, m): A000027 (m=1), A001353 (m=2), A001109 (m=3), A001090 (m=4), A004189 (m=5), A004191 (m=6), this sequence (m=7), A077412 (m=8), A049660 (m=9), A075843 (m=10), A077421 (m=11), A077423 (m=12), A097309 (m=13), A097311 (m=14), A097313 (m=15), A029548 (m=16), A029547 (m=17), A144128 (m=18), A078987 (m=19), A097316 (m=33).
Cf. A323182.

Programs

  • GAP
    m:=7;; a:=[0,1];; for n in [3..20] do a[n]:=2*m*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 23 2019
  • Magma
    [n le 2 select n-1 else 14*Self(n-1)-Self(n-2): n in [1..70]]; // Vincenzo Librandi, Feb 02 2016
    
  • Maple
    0,seq(orthopoly[U](n,7),n=0..30); # Robert Israel, Feb 04 2016
  • Mathematica
    Table[GegenbauerC[n, 1, 7], {n,0,20}] (* Vladimir Joseph Stephan Orlovsky, Sep 11 2008 *)
    LinearRecurrence[{14,-1}, {0,1}, 20] (* Vincenzo Librandi, Feb 02 2016 *)
    ChebyshevU[Range[21] -2, 7] (* G. C. Greubel, Dec 23 2019 *)
    Table[Sum[Binomial[n, 2 k - 1]*7^(n - 2 k + 1)*48^(k - 1), {k, 1, n}], {n, 0, 15}] (* Horst H. Manninger, Jan 16 2022 *)
  • PARI
    concat(0, Vec((x^2/(1-14*x+x^2) + O(x^30)))) \\ Michel Marcus, Feb 02 2016
    
  • PARI
    vector(21, n, polchebyshev(n-2, 2, 7) ) \\ G. C. Greubel, Dec 23 2019
    
  • Sage
    [lucas_number1(n,14,1) for n in range(0,20)] # Zerinvary Lajos, Jun 25 2008
    
  • Sage
    [chebyshev_U(n,7) for n in (-1..20)] # G. C. Greubel, Dec 23 2019
    

Formula

a(n) = 14*a(n-1) - a(n-2).
G.f.: x^2/(1-14*x+x^2).
a(n+1) ~ 1/24*sqrt(3)*(2 + sqrt(3))^(2*n). - Joe Keane (jgk(AT)jgk.org), May 15 2002
a(n+1) = S(n-1, 14), n>=0, with S(n, x) := U(n, x/2) Chebyshev's polynomials of the second kind. S(-1, x) := 0. See A049310.
a(n+1) = ( (7+4*sqrt(3))^n - (7-4*sqrt(3))^n )/(8*sqrt(3)).
a(n+1) = sqrt((A011943(n)^2 - 1)/48), n>=0.
Chebyshev's polynomials U(n-2, x) evaluated at x=7.
a(n) = A001353(2n)/4. - Lekraj Beedassy, Jul 15 2002
4*a(n+1) + A046184(n) = A055793(n+2) + A098301(n+1) 4*a(n+1) + A098301(n+1) + A055793(n+2) = A046184(n+1) (4*a(n+1))^2 = A098301(2n+1) (conjectures). - Creighton Dement, Nov 02 2004
(4*a(n))^2 = A103974(n)^2 - A011922(n-1)^2. - Paul D. Hanna, Mar 06 2005
From Mohamed Bouhamida, May 26 2007: (Start)
a(n) = 13*( a(n-1) + a(n-2) ) - a(n-3).
a(n) = 15*( a(n-1) - a(n-2) ) + a(n-3). (End)
a(n) = b such that (-1)^n/4*Integral_{x=-Pi/2..Pi/2} (sin((2*n-2)*x))/(2-sin(x)) dx = c+b*log(3). - Francesco Daddi, Aug 02 2011
a(n+2) = Sum_{k=0..n} A101950(n,k)*13^k. - Philippe Deléham, Feb 10 2012
Product {n >= 1} (1 + 1/a(n)) = 1/3*(3 + 2*sqrt(3)). - Peter Bala, Dec 23 2012
Product {n >= 2} (1 - 1/a(n)) = 1/7*(3 + 2*sqrt(3)). - Peter Bala, Dec 23 2012
a(n) = (A028230(n) - A001570(n))/2. - Richard R. Forberg, Nov 14 2013
E.g.f.: 1 - exp(7*x)*(12*cosh(4*sqrt(3)*x) - 7*sqrt(3)*sinh(4*sqrt(3)*x))/12. - Stefano Spezia, Dec 11 2022

Extensions

Chebyshev comments from Wolfdieter Lang, Nov 08 2002

A045899 Numbers k such that k+1 and 3*k+1 are perfect squares.

Original entry on oeis.org

0, 8, 120, 1680, 23408, 326040, 4541160, 63250208, 880961760, 12270214440, 170902040408, 2380358351280, 33154114877520, 461777249934008, 6431727384198600, 89582406128846400, 1247721958419651008, 17378525011746267720, 242051628206028097080, 3371344269872647091408
Offset: 1

Views

Author

Andrej Dujella (duje(AT)math.hr)

Keywords

Comments

Essentially the same as A051047.
It appears that a(n) = A046175(n)-A046174(n), that is, the triangular index of the n-th pentagonal triangular number minus its pentagonal index. - Jonathan Vos Post, Feb 28 2011
Sequence lists the nonnegative x solutions when (x + 1)*(3*x + 1) is a square. Positive x solutions when (x - 1)*(3*x - 1) is a square are in A011922. - Bruno Berselli, Feb 20 2018

Crossrefs

Programs

  • Mathematica
    f[n_] := FullSimplify[((Sqrt[3] + 2)*(7 + 4*Sqrt[3])^n - (Sqrt[3] - 2) (7 - 4 Sqrt[3])^n - 4)/6]; Array[f, 18, 0] (* Joseph Biberstine (jrbibers(AT)indiana.edu), Apr 23 2006 *)
    Rest[CoefficientList[Series[-8*x^2/((x - 1)*(x^2 - 14*x + 1)), {x,0,50}], x]] (* G. C. Greubel, Jun 07 2017 *)
    LinearRecurrence[{15,-15,1},{0,8,120},20] (* Harvey P. Dale, Jul 14 2024 *)
  • PARI
    x='x+O('x^50); concat([0], Vec(-8*x^2/((x - 1)*(x^2 - 14*x + 1)))) \\ G. C. Greubel, Jun 07 2017

Formula

a(n) = A046184(n+1) - 1.
a(n) = 14*a(n-1) - a(n-2) + 8.
a(n) = ((2 + sqrt(3))*(7 + 4*sqrt(3))^n + (2 - sqrt(3))*(7 - 4*sqrt(3))^n - 4)/6. - Joseph Biberstine (jrbibers(AT)indiana.edu), Apr 23 2006
a(n) = 8*A076139(n-1) = 4*A217855(n-1) = 2*A123480(n-1) = 8/3*A076140(n-1). - Peter Bala, Dec 31 2012
From Colin Barker, Jul 30 2013: (Start)
G.f.: -8*x^2 / ((x - 1)*(x^2 - 14*x + 1)).
a(n) = 15*a(n-1) - 15*a(n-2) + a(n-3). (End)
E.g.f.: (-4*exp(x) + (2 + sqrt(3))*exp((7-4*sqrt(3))*x) + (2 - sqrt(3))*exp((7+4*sqrt(3))*x))/6. - Ilya Gutkovskiy, Apr 28 2016

A067902 a(n) = 14*a(n-1) - a(n-2); a(0) = 2, a(1) = 14.

Original entry on oeis.org

2, 14, 194, 2702, 37634, 524174, 7300802, 101687054, 1416317954, 19726764302, 274758382274, 3826890587534, 53301709843202, 742397047217294, 10340256951198914, 144021200269567502, 2005956546822746114, 27939370455248878094, 389145229826661547202, 5420093847118012782734
Offset: 0

Views

Author

Lekraj Beedassy, May 13 2003

Keywords

Comments

Solves for x in x^2 - 3*y^2 = 4. [Complete nonnegative solutions are in A003500 and A052530. - Wolfdieter Lang, Sep 05 2021]
For n>0, a(n)+2 is the number of dimer tilings of a 4 X 2n Klein bottle (cf. A103999).
This is the Lucas sequence V(14,1). In addition to the comment above: If x = a(n) then y(n) = (a(n+1) - a(n-1))/24, n >= 1. - Klaus Purath, Aug 17 2021

Examples

			G.f. = 2 + 14*x + 194*x^2 + 2702*x^3 + 37634*x^4 + 524174*x^5 + ...
		

Crossrefs

Row 2 * 2 of array A188644.

Programs

  • GAP
    m:=7;; a:=[2,14];; for n in [3..20] do a[n]:=2*m*a[n-1]-a[n-2]; od; a; # G. C. Greubel, Dec 23 2019
  • Magma
    [Floor((2+Sqrt(3))^(2*n)+(1+Sqrt(3))^(-n)): n in [0..19]]; // Vincenzo Librandi, Mar 31 2011
    
  • Maple
    a := proc(n) option remember: if n=0 then RETURN(2) fi: if n=1 then RETURN(14) fi: 14*a(n-1)-a(n-2): end: for n from 0 to 30 do printf(`%d,`,a(n)) od:
    seq( simplify(2*ChebyshevT(n, 7)), n=0..20); # G. C. Greubel, Dec 23 2019
  • Mathematica
    a[0]=2; a[1]=14; a[n_]:= 14a[n-1] -a[n-2]; Table[a[n], {n,0,20}] (* Robert G. Wilson v, Jan 30 2004 *)
    2*ChebyshevT[Range[21] -1, 7] (* G. C. Greubel, Dec 23 2019 *)
  • PARI
    vector( 21, n, 2*polchebyshev(n-1, 1, 7) ) \\ G. C. Greubel, Dec 23 2019
    
  • Sage
    [lucas_number2(n,14,1) for n in range(0,20)] # Zerinvary Lajos, Jun 26 2008
    
  • Sage
    [2*chebyshev_T(n,7) for n in (0..20)] # G. C. Greubel, Dec 23 2019
    

Formula

G.f.: 2*(1-7*x)/(1-14*x+x^2). - N. J. A. Sloane, Nov 22 2006
a(n) = p^n + q^n, where p = 7 + 4*sqrt(3) and q = 7 - 4*sqrt(3). - Tanya Khovanova, Feb 06 2007
a(n) = 2*A011943(n+1). - R. J. Mathar, Sep 27 2014
From Peter Bala, Oct 16 2019: (Start)
Let F(x) = Product_{n >= 0} (1 + x^(4*n+1))/(1 + x^(4*n+3)). Let alpha = 7 - 4*sqrt(3). This sequence gives the partial denominators in the simple continued fraction expansion of 1 + F(alpha) = 2.07140228197873197080... = 2 + 1/(14 + 1/(194 + 1/(2702 + ...))). Cf. A005248.
12*Sum_{n >= 1} 1/(a(n) - 16/a(n)) = 1.
16*Sum_{n >= 1} (-1)^(n+1)/(a(n) + 12/a(n)) = 1.
Series acceleration formula for sum of reciprocals:
Sum_{n >= 1} 1/a(n) = 1/12 - 16*Sum_{n >= 1} 1/(a(n)*(a(n)^2 - 16)).
Sum_{n >= 1} 1/a(n) = ( (theta_3(7-4*sqrt(3)))^2 - 1 )/4, where theta_3(x) = 1 + 2*Sum_{n >= 1} x^(n^2) (see A000122). Cf. A153415 and A003499.
(End)
From Klaus Purath, Aug 17 2021: (Start)
a(n) = (a(n-1)*a(n-2) + 2688)/a(n-3), n >= 3.
a(n) = (a(n-1)^2 + 192)/a(n-2), n >= 2.
a(2*n) = A302332(n-1) + A302332(n), n >= 1.
a(2*n+1) = 14*A302332(n). (End)
a(n) = A003500(2*n) = S(2*n,4) - S(2*n-2, 4) = 2*T(2*n,2), for n >= 0, with Chebyshev S and T. S(n, 4) = A001353(n+1) and T(n, 2) = A001075(n). - Wolfdieter Lang, Sep 06 2021

A232771 Values of x satisfying x^2 = floor(y^2/3 + y).

Original entry on oeis.org

0, 1, 3, 6, 17, 43, 84, 237, 599, 1170, 3301, 8343, 16296, 45977, 116203, 226974, 640377, 1618499, 3161340, 8919301, 22542783, 44031786, 124229837, 313980463, 613283664, 1730298417, 4373183699, 8541939510, 24099948001
Offset: 1

Views

Author

Richard R. Forberg, Nov 30 2013

Keywords

Comments

Corresponding values of y are given by A232765.
a(n) are also the values of x satisfying x^2 = floor(y^2/3 - y).
Let b(n) equal the second differences of a(n), where b(1) = 1, then b(3n-2) = b(3n-1) = A028230(n) and b(3n) = A067900(n) for n>0.

Crossrefs

Formula

Empirical g.f.: x^2*(x^4+3*x^3+6*x^2+3*x+1) / (x^6-14*x^3+1). - Colin Barker, Dec 30 2014
Showing 1-5 of 5 results.