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.

Previous Showing 11-20 of 27 results. Next

A322699 Square array A(n,k), n >= 0, k >= 0, read by antidiagonals, where A(n,k) is 1/2 * (-1 + Sum_{j=0..k} binomial(2*k,2*j)*(n+1)^(k-j)*n^j).

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 8, 2, 0, 0, 49, 24, 3, 0, 0, 288, 242, 48, 4, 0, 0, 1681, 2400, 675, 80, 5, 0, 0, 9800, 23762, 9408, 1444, 120, 6, 0, 0, 57121, 235224, 131043, 25920, 2645, 168, 7, 0, 0, 332928, 2328482, 1825200, 465124, 58080, 4374, 224, 8, 0
Offset: 0

Views

Author

Seiichi Manyama, Dec 23 2018

Keywords

Examples

			Square array begins:
   0, 0,   0,    0,      0,       0,        0, ...
   0, 1,   8,   49,    288,    1681,     9800, ...
   0, 2,  24,  242,   2400,   23762,   235224, ...
   0, 3,  48,  675,   9408,  131043,  1825200, ...
   0, 4,  80, 1444,  25920,  465124,  8346320, ...
   0, 5, 120, 2645,  58080, 1275125, 27994680, ...
   0, 6, 168, 4374, 113568, 2948406, 76545000, ...
		

Crossrefs

Columns 0-5 give A000004, A001477, A033996, A322675, A322677, A322745.
Main diagonal gives A322746.
Cf. A173175 (A(n,2*n)), A322790.

Programs

  • Mathematica
    Unprotect[Power]; 0^0 := 1; Protect[Power]; Table[(-1 + Sum[Binomial[2 k, 2 j] (# + 1)^(k - j)*#^j, {j, 0, k}])/2 &[n - k], {n, 0, 9}, {k, n, 0, -1}] // Flatten (* Michael De Vlieger, Jan 01 2019 *)
    nmax = 9; row[n_] := LinearRecurrence[{4n+3, -4n-3, 1}, {0, n, 4n(n+1)}, nmax+1]; T = Array[row, nmax+1, 0]; A[n_, k_] := T[[n+1, k+1]];
    Table[A[n-k, k], {n, 0, nmax}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Jan 06 2019 *)
  • Ruby
    def ncr(n, r)
      return 1 if r == 0
      (n - r + 1..n).inject(:*) / (1..r).inject(:*)
    end
    def A(k, n)
      (0..n).map{|i| (0..k).inject(-1){|s, j| s + ncr(2 * k, 2 * j) * (i + 1) ** (k - j) * i ** j} / 2}
    end
    def A322699(n)
      a = []
      (0..n).each{|i| a << A(i, n - i)}
      ary = []
      (0..n).each{|i|
        (0..i).each{|j|
          ary << a[i - j][j]
        }
      }
      ary
    end
    p A322699(10)

Formula

sqrt(A(n,k)+1) + sqrt(A(n,k)) = (sqrt(n+1) + sqrt(n))^k.
sqrt(A(n,k)+1) - sqrt(A(n,k)) = (sqrt(n+1) - sqrt(n))^k.
A(n,0) = 0, A(n,1) = n and A(n,k) = (4*n+2) * A(n,k-1) - A(n,k-2) + 2*n for k > 1.
A(n,k) = (T_{k}(2*n+1) - 1)/2 where T_{k}(x) is a Chebyshev polynomial of the first kind.
T_1(x) = x. So A(n,1) = (2*n+1-1)/2 = n.

A173130 a(n) = Cosh[(2 n - 1) ArcCosh[n]].

Original entry on oeis.org

0, 1, 26, 3363, 937444, 456335045, 343904160606, 371198523608647, 543466014742175624, 1036834190110356583689, 2499384905955651114739810, 7429238104512325157021090411, 26695718139185294187938997247212
Offset: 0

Views

Author

Artur Jasinski, Feb 10 2010

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Round[Cosh[(2 n - 1) ArcCosh[n]]], {n, 0, 20}] (* Artur Jasinski *)

Formula

a(n) ~ 2^(2*n-2) * n^(2*n-1). - Vaclav Kotesovec, Apr 05 2016

A173131 a(n) = (Cosh[(2n-1)ArcSinh[n]])^2.

Original entry on oeis.org

1, 2, 1445, 19740250, 1361599599377, 298514762397852026, 160545187370375075046277, 179656719395983409634002348450, 373368546362937441101158606899394625
Offset: 0

Views

Author

Artur Jasinski, Feb 10 2010

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Round[Cosh[(2 n - 1) ArcSinh[n]]^2], {n, 0, 10}] (* Artur Jasinski *)

Formula

a(n) ~ 2^(4*n-4) * n^(4*n-2). - Vaclav Kotesovec, Apr 05 2016

A173148 a(n) = cos(2*n*arccos(sqrt(n))).

Original entry on oeis.org

1, 1, 17, 485, 18817, 930249, 55989361, 3974443213, 325142092801, 30122754096401, 3117419602578001, 356452534779818421, 44627167107085622401, 6071840759403431812825, 892064955046043465408177, 140751338790698080509966749, 23737154316161495960243527681
Offset: 0

Views

Author

Artur Jasinski, Feb 11 2010

Keywords

Comments

The Chebyshev polynomial T_n is defined by cos(nx) = T_n(cos(x)). So T_2n(cos(x)) = cos(2nx) = cos^2(nx) - 1 = (T_n(x))^2 - 1 consists of only even powers of x. As a result, a(n) = T_2n(sqrt(n)) is an integer. - Michael B. Porter, Jan 01 2019

Crossrefs

Programs

  • GAP
    a:=List([0..20],n->Sum([0..n],k->Binomial(2*n,2*k)*(n-1)^(n-k)*n^k));; Print(a); # Muniru A Asiru, Jan 03 2019
    
  • Magma
    [&+[Binomial(2*n,2*k)*(n-1)^(n-k)*n^k: k in [0..n]]: n in [0..20]]; // Vincenzo Librandi, Jan 03 2019
  • Mathematica
    Table[Round[Cos[2 n ArcCos[Sqrt[n]]]], {n, 0, 30}] (* Artur Jasinski, Feb 11 2010 *)
  • PARI
    {a(n) = sum(k=0, n, binomial(2*n, 2*k)*(n-1)^(n-k)*n^k)} \\ Seiichi Manyama, Dec 27 2018
    
  • PARI
    {a(n) = round(cosh(2*n*acosh(sqrt(n))))} \\ Seiichi Manyama, Dec 27 2018
    
  • PARI
    {a(n) = polchebyshev(n, 1, 2*n-1)} \\ Seiichi Manyama, Dec 29 2018
    

Formula

a(n) ~ exp(-1/2) * 2^(2*n-1) * n^n. - Vaclav Kotesovec, Apr 05 2016
a(n) = Sum_{k=0..n} binomial(2*n,2*k)*(n-1)^(n-k)*n^k. - Seiichi Manyama, Dec 27 2018
a(n) = cosh(2*n*arccosh(sqrt(n))). - Seiichi Manyama, Dec 27 2018
a(n) = T_{2*n}(sqrt(n)) = T_{n}(2*n-1) where T_{n}(x) is a Chebyshev polynomial of the first kind. - Seiichi Manyama, Dec 29 2018
a(n) = A322790(n-1, n) for n > 0. - Seiichi Manyama, Dec 29 2018

Extensions

Minor edits by Vaclav Kotesovec, Apr 05 2016

A173133 a(n) = Sinh[(2n-1) ArcSinh[n]].

Original entry on oeis.org

0, 1, 38, 4443, 1166876, 546365045, 400680904674, 423859315570607, 611038907405197432, 1151555487914640463209, 2748476184146759127540190, 8102732939160371170806346243, 28915133156938367486730067779348
Offset: 0

Views

Author

Artur Jasinski, Feb 10 2010

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Round[Sinh[(2 n - 1) ArcSinh[n]]], {n, 0, 20}] (* Artur Jasinski *)
    Round[Table[1/2 (n - Sqrt[1 + n^2])^(2 n - 1) + 1/2 (n + Sqrt[1 + n^2])^(2 n - 1), {n, 0, 10}]] (* Artur Jasinski, Feb 14 2010 *)

Formula

a(n) = 1/2 (n - sqrt(1 + n^2))^(2 n - 1) + 1/2 (n + sqrt(1 + n^2))^(2 n - 1). - Artur Jasinski, Feb 14 2010

Extensions

Minor edits by Vaclav Kotesovec, Apr 05 2016

A173134 a(n) = Sinh[(2n-1)ArcCosh[n]]^2.

Original entry on oeis.org

-1, 0, 675, 11309768, 878801253135, 208241673295152024, 118270071682117442287235, 137788343929239264227213170608, 295355309179742652677310128859789375
Offset: 0

Views

Author

Artur Jasinski, Feb 10 2010

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Round[Sinh[(2 n - 1) ArcCosh[n]]^2], {n, 0, 20}]

Formula

a(n) ~ 2^(4*n-4) * n^(4*n-2). - Vaclav Kotesovec, Apr 05 2016

A097726 Pell equation solutions (5*a(n))^2 - 26*b(n)^2 = -1 with b(n):=A097727(n), n >= 0.

Original entry on oeis.org

1, 103, 10505, 1071407, 109273009, 11144775511, 1136657829113, 115927953794015, 11823514629160417, 1205882564220568519, 122988198035868828521, 12543590317094399940623, 1279323224145592925115025, 130478425272533383961791927, 13307520054574259571177661529
Offset: 0

Views

Author

Wolfdieter Lang, Aug 31 2004

Keywords

Comments

a(-1) = -1. - Artur Jasinski, Feb 10 2010
5*a(n) gives the x-values in the solution to the Pell equation x^2 - 26*y^2 = -1. - Colin Barker, Aug 24 2013

Examples

			(x,y) = (5,1), (515,101), (52525,10301), ... give the positive integer solutions to x^2 - 26*y^2 = -1.
		

Crossrefs

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

Programs

  • Mathematica
    Table[(1/5) Round[N[Sinh[(2 n - 1) ArcSinh[5]], 100]], {n, 1, 50}] (* Artur Jasinski, Feb 10 2010 *)
    CoefficientList[Series[(1 + x)/(1 - 102 x + x^2), {x, 0, 40}], x] (* Vincenzo Librandi, Apr 13 2014 *)
    LinearRecurrence[{102,-1},{1,103},20] (* Harvey P. Dale, Aug 20 2017 *)
  • PARI
    x='x+O('x^99); Vec((1+x)/(1-102*x+x^2)) \\ Altug Alkan, Apr 05 2018

Formula

G.f.: (1 + x)/(1 - 102*x + x^2).
a(n) = S(n, 2*51) + S(n-1, 2*51) = S(2*n, 2*sqrt(26)), 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, 5*i)/(5*i) with the imaginary unit i and Chebyshev polynomials of the first kind. See the T-triangle A053120.
a(n) = 102*a(n-1) - a(n-2) for n > 1; a(0)=1, a(1)=103. - Philippe Deléham, Nov 18 2008
a(n) = (1/5)*sinh((2*n-1)*arcsinh(5)), n >= 1. - Artur Jasinski, Feb 10 2010

Extensions

More terms from Harvey P. Dale, Aug 20 2017

A276575 After a(0)=0, the first differences of A276573.

Original entry on oeis.org

0, 3, 3, 2, 3, 4, 1, 2, 3, 3, 3, 3, 2, 3, 3, 2, 3, 2, 3, 3, 2, 3, 3, 4, 1, 3, 3, 2, 3, 3, 2, 3, 2, 3, 2, 3, 3, 3, 3, 3, 3, 4, 3, 2, 3, 3, 3, 2, 3, 3, 2, 3, 4, 1, 3, 2, 3, 3, 3, 2, 2, 3, 3, 3, 2, 3, 3, 4, 3, 3, 3, 3, 3, 2, 3, 4, 3, 3, 3, 3, 2, 3, 3, 4, 2, 3, 4, 3, 2, 3, 3, 4, 1, 3, 2, 3, 3, 3, 2, 3, 2, 3, 3, 3, 2, 3, 3, 2, 2, 3, 3, 2, 3, 3, 3, 3, 2, 3, 3, 2, 3
Offset: 0

Views

Author

Antti Karttunen, Sep 07 2016

Keywords

Crossrefs

Formula

a(n) = A002828(A276573(n)).
a(0) = 0; for n >= 1, a(n) = A276573(n) - A276573(n-1).
Other identities.
For all n >= 1, a(A260731(A132592(n))) = a(A260733(A001541(n))) = 2. [This is implied by the fact observed in the Polster video. Of course 2's occur at other points too.]

A065113 Sum of the squares of the a(n)-th and the (a(n)+1)st triangular numbers (A000217) is a perfect square.

Original entry on oeis.org

6, 40, 238, 1392, 8118, 47320, 275806, 1607520, 9369318, 54608392, 318281038, 1855077840, 10812186006, 63018038200, 367296043198, 2140758220992, 12477253282758, 72722761475560, 423859315570606, 2470433131948080, 14398739476117878, 83922003724759192
Offset: 1

Views

Author

Robert G. Wilson v, Nov 12 2001

Keywords

Comments

The sequence of square roots of the sum of the squares of the n-th and the (n+1)st triangular numbers is A046176.

Examples

			T6 = 21 and T7 = 28, 21^2 + 28^2 = 441 + 784 = 1225 = 35^2.
		

Crossrefs

Cf. A001652, A002315, A003499 (first differences), A065651.

Programs

  • Mathematica
    CoefficientList[ Series[2*(x - 3)/(-1 + 7x - 7x^2 + x^3), {x, 0, 24} ], x]
    LinearRecurrence[{7,-7,1},{6,40,238},41] (* Harvey P. Dale, Dec 27 2011 *)
  • PARI
    a(n)=-1+subst(poltchebi(abs(n+1))-poltchebi(abs(n)),x,3)/2
    
  • PARI
    Vec(2*x*(3-x)/((1-6*x+x^2)*(1-x)) + O(x^40)) \\ Colin Barker, Mar 05 2016

Formula

a(n) = 2*A001652(n) = -1 + A002315(n).
a(n) - a(n-1) = A003499(n).
From Michael Somos, Apr 07 2003: (Start)
G.f.: 2*x*(3-x)/((1-6*x+x^2)*(1-x)).
a(n) = 6*a(n-1) - a(n-2) + 4.
a(-1-n) = -a(n) - 2. (End)
a(1)=6, a(2)=40, a(3)=238, a(n) = 7*a(n-1)-7*a(n-2)+a(n-3). - Harvey P. Dale, Dec 27 2011
a(n)^2 + (a(n)+2)^2 = A075870(n+1)^2 = A165518(n+1). - Joerg Arndt, Feb 15 2012
a(n) = (-2-(3-2*sqrt(2))^n*(-1+sqrt(2))+(1+sqrt(2))*(3+2*sqrt(2))^n)/2. - Colin Barker, Mar 05 2016
From Klaus Purath, Sep 05 2021: (Start)
(a(n+1) - a(n) - a(n-1) + a(n-2))/8 = A005319(n), for n >= 3.
((a(n) - a(n-1))^2)/2 - 2 = A005319(n)^2 = 2*A132592(n), for n>= 2.
a(n) = A265278(2*n+1).
a(n) = A293004(2*n+1).
a(n) = A213667(2*n).
a(n) = Sum_{k=1..n} A003499(k). (End)

A173170 a(n) = sin^2((2n-1)*arcsin(sqrt n)) = 1 - sin^2( (2n-1)*arccos(sqrt n)).

Original entry on oeis.org

0, 1, 50, 23763, 25421764, 48225038405, 142786923879606, 608447515452613207, 3527836867501829594888, 26710782540478226038759689, 255922222218837615280903143610, 3026917140685147530327256796600411
Offset: 0

Views

Author

Artur Jasinski, Feb 11 2010

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Round[Sin[(2 n - 1) ArcSin[Sqrt[n]]]^2], {n, 0, 20}] (* Artur Jasinski, Feb 11 2010 *)

Formula

a(n) ~ exp(-1) * 2^(4*n-4) * n^(2*n-1). - Vaclav Kotesovec, Apr 05 2016

Extensions

Minor edits by Vaclav Kotesovec, Apr 05 2016
Previous Showing 11-20 of 27 results. Next