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 51-60 of 79 results. Next

A195318 Centered 44-gonal numbers.

Original entry on oeis.org

1, 45, 133, 265, 441, 661, 925, 1233, 1585, 1981, 2421, 2905, 3433, 4005, 4621, 5281, 5985, 6733, 7525, 8361, 9241, 10165, 11133, 12145, 13201, 14301, 15445, 16633, 17865, 19141, 20461, 21825, 23233, 24685, 26181, 27721, 29305, 30933, 32605, 34321, 36081, 37885, 39733
Offset: 1

Views

Author

Omar E. Pol, Sep 16 2011

Keywords

Comments

Sequence found by reading the line from 1, in the direction 1, 45, ..., in the square spiral whose vertices are the generalized tridecagonal numbers A195313. Semi-axis opposite to A195323 in the same spiral.

Crossrefs

Programs

Formula

a(n) = 22*n^2 - 22*n + 1.
Sum_{n>=1} 1/a(n) = Pi*tan(3*Pi/(2*sqrt(11)))/(6*sqrt(11)). - Amiram Eldar, Feb 11 2022
G.f.: -x*(1+42*x+x^2)/(x-1)^3. - R. J. Mathar, May 07 2024
From Elmo R. Oliveira, Nov 15 2024: (Start)
E.g.f.: exp(x)*(22*x^2 + 1) - 1.
a(n) = 2*A069173(n) - 1.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 3. (End)

A335333 Square array T(n,k), n >= 0, k >= 0, read by antidiagonals, where column k is the expansion of 1/sqrt(1 - 2*(2*k+1)*x + x^2).

Original entry on oeis.org

1, 1, 1, 1, 3, 1, 1, 5, 13, 1, 1, 7, 37, 63, 1, 1, 9, 73, 305, 321, 1, 1, 11, 121, 847, 2641, 1683, 1, 1, 13, 181, 1809, 10321, 23525, 8989, 1, 1, 15, 253, 3311, 28401, 129367, 213445, 48639, 1, 1, 17, 337, 5473, 63601, 458649, 1651609, 1961825, 265729, 1
Offset: 0

Views

Author

Seiichi Manyama, Jun 02 2020

Keywords

Examples

			Square array begins:
  1,    1,     1,      1,      1,       1, ...
  1,    3,     5,      7,      9,      11, ...
  1,   13,    37,     73,    121,     181, ...
  1,   63,   305,    847,   1809,    3311, ...
  1,  321,  2641,  10321,  28401,   63601, ...
  1, 1683, 23525, 129367, 458649, 1256651, ...
		

Crossrefs

Columns k=0..4 give A000012, A001850, A006442, A084768, A084769.
Rows n=0..6 give A000012, A005408, A003154(n+1), A160674, A144124, A335338, A144126.
Main diagonal gives A331656.
T(n,n-1) gives A331657.

Programs

  • Mathematica
    T[n_, k_] := LegendreP[n, 2*k + 1]; Table[T[k, n - k], {n, 0, 9}, {k, 0, n}] // Flatten (* Amiram Eldar, May 03 2021 *)
  • PARI
    T(n, k) = pollegendre(n, 2*k+1);

Formula

T(n,k) is the coefficient of x^n in the expansion of (1 + (2*k+1)*x + k*(k+1)*x^2)^n.
T(n,k) = Sum_{j=0..n} k^j * (k+1)^(n-j) * binomial(n,j)^2.
T(n,k) = Sum_{j=0..n} k^j * binomial(n,j) * binomial(n+j,j).
n * T(n,k) = (2*k+1) * (2*n-1) * T(n-1,k) - (n-1) * T(n-2,k).
T(n,k) = P_n(2*k+1), where P_n is n-th Legendre polynomial.
From Seiichi Manyama, Aug 30 2025: (Start)
T(n,k) = (-1)^n * Sum_{j=0..n} (1/(2*(2*k+1)))^(n-2*j) * binomial(-1/2,j) * binomial(j,n-j).
T(n,k) = Sum_{j=0..floor(n/2)} (k*(k+1))^j * (2*k+1)^(n-2*j) * binomial(n,2*j) * binomial(2*j,j).
E.g.f. of column k: exp((2*k+1)*x) * BesselI(0, 2*sqrt(k*(k+1))*x). (End)

A006051 Square hex numbers.

Original entry on oeis.org

1, 169, 32761, 6355441, 1232922769, 239180661721, 46399815451081, 9001325016847969, 1746210653453054881, 338755865444875798921, 65716891685652451935769, 12748738231151130799740241, 2473189499951633722697670961, 479786014252385791072548426169
Offset: 1

Views

Author

Keywords

Comments

Numbers n of the form n = y^2 = 3*x^2 - 3*x + 1.

Examples

			G.f. = x + 169*x^2 + 32761*x^3 + 6355441*x^4 + 1232922769*x^5 + ...
		

References

  • M. Gardner, Time Travel and Other Mathematical Bewilderments. Freeman, NY, 1988, p. 19.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A003500.
Intersection of A000290 and A003215.
Values of x are given by A001922, values of y by A001570.

Programs

  • Magma
    [(7*Evaluate(ChebyshevSecond(n),97) - 7*Evaluate(ChebyshevU(n-1), 97) + 1)/8: n in [1..30]]; // G. C. Greubel, Nov 04 2017; Oct 07 2022
    
  • Mathematica
    Rest@ CoefficientList[Series[x(1-26x+x^2)/((1-x)(1-194x+x^2)), {x,0,20}], x] (* Michael De Vlieger, Jan 02 2017 *)
    LinearRecurrence[{195,-195,1},{1,169,32761},20] (* Harvey P. Dale, Nov 03 2017 *)
  • PARI
    {a(n) = sqr( real( (2 + quadgen( 12)) ^ (2*n - 1)) / 2)} /* Michael Somos, Feb 15 2011 */
    
  • SageMath
    def A006051(n): return (7*chebyshev_U(n-1,97) - 7*chebyshev_U(n-2,97) + 1)/8
    [A006051(n) for n in range(1,31)] # G. C. Greubel, Oct 07 2022

Formula

a(n) = A001570(n)^2.
a(1 - n) = a(n).
G.f.: x * (1 - 26*x + x^2) / ((1 - x) * (1 - 194*x + x^2)). - Simon Plouffe in his 1992 dissertation
a(n) = 194*a(n-1) - a(n-2) - 24, a(1)=1, a(2)=169. - James Sellers, Jul 04 2000
a(n+1) = A003215(A001921(n)). - Joerg Arndt, Jan 02 2017
a(n) = (1/8)*(1 + 7*(ChebyshevU(n-1, 97) - ChebyshevU(n-2, 97))). - G. C. Greubel, Oct 07 2022

A211014 Second 14-gonal numbers: n*(6*n+5).

Original entry on oeis.org

0, 11, 34, 69, 116, 175, 246, 329, 424, 531, 650, 781, 924, 1079, 1246, 1425, 1616, 1819, 2034, 2261, 2500, 2751, 3014, 3289, 3576, 3875, 4186, 4509, 4844, 5191, 5550, 5921, 6304, 6699, 7106, 7525, 7956, 8399, 8854, 9321, 9800, 10291, 10794, 11309, 11836, 12375
Offset: 0

Views

Author

Omar E. Pol, Aug 04 2012

Keywords

Comments

Sequence found by reading the line from 0, in the direction 0, 34, ... and the line from 11 in the direction 11, 69, ..., in the square spiral whose vertices are the generalized 14-gonal numbers A195818.

Crossrefs

Bisection of A195818.
Second k-gonal numbers (k=5..14): A005449, A014105, A147875, A045944, A179986, A033954, A062728, A135705, A211013, this sequence.
Cf. A051866.
Cf. A003154.

Programs

Formula

a(n) = -2*Sum_{k=0..n-1} binomial(6*n+5, 6*k+8)*Bernoulli(6*k+8). - Michel Marcus, Jan 11 2016
From G. C. Greubel, Jul 04 2019: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
G.f.: x*(11+x)/(1-x)^3.
E.g.f.: x*(11+6*x)*exp(x). (End)
From Amiram Eldar, Feb 28 2022: (Start)
Sum_{n>=1} 1/a(n) = sqrt(3)*Pi/10 + 6/25 - 3*log(3)/10 - 2*log(2)/5.
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/5 + log(2)/5 - 6/25 - sqrt(3)*log(sqrt(3)+2)/5. (End)
a(n) = A003154(n+1) - n - 1. - Leo Tavares, Jan 29 2023

A257565 Generalized Fubini numbers. Square array read by ascending antidiagonals, A(n,k) = 1 + k*(Sum_{j=1..n-1} C(n,j)*A(j,k)); n>=0 and k>=0.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 13, 5, 1, 1, 1, 75, 37, 7, 1, 1, 1, 541, 365, 73, 9, 1, 1, 1, 4683, 4501, 1015, 121, 11, 1, 1, 1, 47293, 66605, 17641, 2169, 181, 13, 1, 1, 1, 545835, 1149877, 367927, 48601, 3971, 253, 15, 1, 1, 1, 7087261, 22687565, 8952553, 1306809, 108901, 6565, 337, 17, 1, 1
Offset: 0

Views

Author

Peter Luschny, May 08 2015

Keywords

Comments

M. Mureşan defined the generalized Fubini numbers as the enumerators of the k-labeled ordered p partitions of an n-set.

Examples

			      1,       1,       1,       1,        1,         1, ...  A000012
      1,       1,       1,       1,        1,         1, ...  A000012
      1,       3,       5,       7,        9,        11, ...  A005408
      1,      13,      37,      73,      121,       181, ...  A003154
      1,      75,     365,    1015,     2169,      3971, ...  A193252
      1,     541,    4501,   17641,    48601,    108901, ...
      1,    4683,   66605,  367927,  1306809,   3583811, ...
      1,   47293, 1149877, 8952553, 40994521, 137595781, ...
A000012, A000670, A050351, A050352,  A050353,
		

References

  • M. Mureşan, On the generalized Fubini numbers. (Romanian) Stud. Cercet. Mat. 37, 70-76 (1985).

Crossrefs

Programs

  • Maple
    F := proc(n,k) option remember; 1+k*add(binomial(n,j)*F(j,k),j=1..n-1) end:
    seq(print(seq(F(n-k,k),k=0..n)), n=0..7); # triangular form
    egf := k -> 1+1/(1/(exp(z)-1)-k): # egf of column k
    for k from 0 to 4 do seq(j!*coeff(series(egf(k),z,10),z,j),j=0..8) od;
    A := (n,k) -> `if`(n=0,1,add(k^(n-j-1)*(k+1)^j*combinat:-eulerian1(n,j),j=0..n-1)): seq(print(seq(A(n,k),k=0..5)),n=0..7);
  • Mathematica
    A[n_, k_] := A[n, k] = 1 + k Sum[Binomial[n, j] A[j, k], {j, 1, n - 1}]; Table[A[n - k, k], {n, 0, 10}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 30 2016 *)

Formula

E.g.f. of column k: 1+1/(1/(exp(z)-1)-k).
A(n,k) = Sum_{j=0..n-1} k^j*j!*{n,j+1} for n>0, else 1; {n,j} denotes the Stirling subset numbers.
A(n,k) = Sum_{j=0..n-1} k^(n-j-1)*(k+1)^j* for n>0, else 1; denotes the Eulerian numbers.

A146325 Period 3: repeat [1, 4, 1].

Original entry on oeis.org

1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1, 4, 1, 1
Offset: 1

Views

Author

Artur Jasinski, Oct 30 2008

Keywords

Comments

Continued fraction of (1 + sqrt(26))/5 = A188659.
Digital roots of the centered triangular numbers A005448. - Ant King, May 08 2012
Also the digital roots of centered 12-gonal numbers A003154. - Peter M. Chema, Dec 20 2023

Crossrefs

Cf. A003154, A005448, A021337, A131534 (square roots), A188659.

Programs

  • Magma
    &cat [[1,4,1]^^40]; // Bruno Berselli, Jun 27 2016
  • Maple
    seq(op([1, 4, 1]), n=1..50); # Wesley Ivan Hurt, Jul 01 2016
  • Mathematica
    Table[Round[N[4 (Cos[(2 n - 1) ArcTan[Sqrt[3]]])^2, 100]], {n, 1, 100}]
    PadLeft[{},111,{1,4,1}] (* Harvey P. Dale, Sep 18 2011 *)
  • PARI
    a(n)=1+3*(n%3==2) \\ Jaume Oliver Lafont, Mar 24 2009
    

Formula

a(n) = 4*(cos((2*n - 1)*Pi/3))^2 = 4 - 4*(sin((2*n - 1)*Pi/3))^2.
a(n+3) = a(n).
a(n) = 2 - cos(2*Pi*n/3) + sqrt(3)*sin(2*Pi*n/3).
O.g.f.: x*(1+4*x+x^2)/(1-x^3). [Richard Choulet, Nov 03 2008]
a(n) = 6 - a(n-1) - a(n-2) for n>2. - Ant King, Jun 12 2012
a(n) = (n mod 3)^(n mod 3). - Bruno Berselli, Jun 27 2016
a(n) = 1 + A021337(n) for n>0. - Wesley Ivan Hurt, Jul 01 2016

A201279 a(n) = 6n^2 + 10n + 5.

Original entry on oeis.org

5, 21, 49, 89, 141, 205, 281, 369, 469, 581, 705, 841, 989, 1149, 1321, 1505, 1701, 1909, 2129, 2361, 2605, 2861, 3129, 3409, 3701, 4005, 4321, 4649, 4989, 5341, 5705, 6081, 6469, 6869, 7281, 7705, 8141, 8589, 9049, 9521, 10005, 10501, 11009, 11529, 12061
Offset: 0

Views

Author

Keywords

Comments

Numbers n where 6n-5 is a square of a number type 6n-1.
Also sequence found by reading the line from 5, in the direction 5, 21,..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. - Omar E. Pol, Jul 18 2012
The spiral mentioned above naturally appears on a "graphene" like lattice (planar net 6^3). The opposite diagonal is A080859. - Yuriy Sibirmovsky, Oct 04 2016
First differences of A048395. - Leo Tavares, Nov 24 2021 [Corrected by Omar E. Pol, Dec 26 2021]

Crossrefs

Programs

  • Magma
    [6*n^2 + 10*n + 5: n in [0..60]]; // Vincenzo Librandi, Dec 01 2011
  • Mathematica
    LinearRecurrence[{3,-3,1},{5,21,49},50] (* Vincenzo Librandi, Dec 01 2011 *)
    Table[6 n^2 + 10 n + 5, {n, 0, 44}] (* or *)
    CoefficientList[Series[(1 + x) (5 + x)/(1 - x)^3, {x, 0, 44}], x] (* Michael De Vlieger, Oct 04 2016 *)
  • PARI
    a(n)=6*n^2+10*n+5 \\ Charles R Greathouse IV, Nov 29 2011
    

Formula

G.f.: (1+x)*(5+x)/(1-x)^3. - Colin Barker, Jan 09 2012
a(n) = 1 + A033579(n+1). - Omar E. Pol, Jul 18 2012
a(n) = (n+1)*A001844(n+1)-n*A001844(n). [Bruno Berselli, Jan 15 2013]
From Leo Tavares, Nov 24 2021: (Start)
a(n) = A003154(n+2) - A022144(n+1). See Diamond Frame Stars illustration.
a(n) = A016754(n) + A046092(n+1). (End)

A202804 a(n) = n*(6*n+4).

Original entry on oeis.org

0, 10, 32, 66, 112, 170, 240, 322, 416, 522, 640, 770, 912, 1066, 1232, 1410, 1600, 1802, 2016, 2242, 2480, 2730, 2992, 3266, 3552, 3850, 4160, 4482, 4816, 5162, 5520, 5890, 6272, 6666, 7072, 7490, 7920, 8362, 8816, 9282, 9760, 10250, 10752, 11266, 11792, 12330
Offset: 0

Views

Author

Jeremy Gardiner, Dec 24 2011

Keywords

Comments

Sequence found by reading the line from 0, in the direction 0, 10, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. Opposite numbers to the members of A033579 in the same spiral. - Omar E. Pol, Jul 17 2012
Partial sums give A163815. - Leo Tavares, Feb 25 2022

Crossrefs

Programs

  • Maple
    A202804:=n->n*(6*n+4): seq(A202804(n), n=0..100); # Wesley Ivan Hurt, Apr 09 2017
  • Mathematica
    Table[n(6n+4),{n,0,50}] (* or *) LinearRecurrence[{3,-3,1},{0,10,32},50] (* Harvey P. Dale, Dec 28 2015 *)
  • PARI
    x='x + O('x^50); concat([0], Vec(-2*x*(5 + x)/(x - 1)^3)) \\ Indranil Ghosh, Apr 10 2017

Formula

a(n) = 2*n(3*n+2) = 6*n^2 + 4*n = 2*A045944(n).
a(n) = A080859(n) - 1. - Omar E. Pol, Jul 18 2012
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, Dec 28 2015
G.f.: 2*x*(5 + x)/(1 - x)^3. - Indranil Ghosh, Apr 10 2017
a(n) = A003154(n+1) - A005408(n). - Leo Tavares, Feb 25 2022
From Amiram Eldar, Mar 01 2022: (Start)
Sum_{n>=1} 1/a(n) = (Pi/sqrt(3) - 3*log(3) + 3)/8.
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi/(4*sqrt(3)) - 3/8. (End)
E.g.f.: 2*exp(x)*x*(5 + 3*x). - Elmo R. Oliveira, Dec 12 2024

A144535 Numerators of continued fraction convergents to sqrt(3)/2.

Original entry on oeis.org

0, 1, 6, 13, 84, 181, 1170, 2521, 16296, 35113, 226974, 489061, 3161340, 6811741, 44031786, 94875313, 613283664, 1321442641, 8541939510, 18405321661, 118973869476, 256353060613, 1657092233154, 3570537526921, 23080317394680, 49731172316281, 321467351292366
Offset: 0

Views

Author

N. J. A. Sloane, Dec 29 2008

Keywords

Examples

			0, 1, 6/7, 13/15, 84/97, 181/209, 1170/1351, 2521/2911, 16296/18817, 35113/40545, ...
		

Crossrefs

Bisections give A001570, A011945.

Programs

  • Magma
    I:=[0, 1, 6, 13]; [n le 4 select I[n] else 14*Self(n-2)-Self(n-4): n in [1..30]]; // Vincenzo Librandi, Dec 10 2013
    
  • Maple
    with(numtheory); Digits:=200: cf:=convert(evalf(sqrt(3)/2,confrac); [seq(nthconver(cf,i), i=0..100)];
  • Mathematica
    CoefficientList[Series[x (1 + 6 x - x^2)/((1 - 4 x + x^2) (1 + 4 x + x^2)), {x, 0, 40}], x] (* Vincenzo Librandi, Dec 10 2013 *)
    Numerator[Convergents[Sqrt[3]/2,30]] (* or *) LinearRecurrence[{0,14,0,-1},{0,1,6,13},30] (* Harvey P. Dale, Feb 10 2014 *)
  • PARI
    Vec(x*(1+6*x-x^2)/((1-4*x+x^2)*(1+4*x+x^2)) + O(x^30)) \\ Colin Barker, Mar 27 2016

Formula

From Colin Barker, Apr 14 2012: (Start)
a(n) = 14*a(n-2) - a(n-4).
G.f.: x*(1 + 6*x - x^2)/((1 - 4*x + x^2)*(1 + 4*x + x^2)). (End)
a(n) = ((-(-2-sqrt(3))^n*(-3+sqrt(3)) + (2-sqrt(3))^n*(-3+sqrt(3)) - (3+sqrt(3))*((-2+sqrt(3))^n - (2+sqrt(3))^n)))/(8*sqrt(3)). - Colin Barker, Mar 27 2016
a(2*n) = 6*a(2*n-1) + a(2*n-2). a(2*n+1) = A003154(A101265(n+1)). - John Elias, Dec 10 2021

A319384 a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5), a(0)=1, a(1)=5, a(2)=9, a(3)=21, a(4)=29.

Original entry on oeis.org

1, 5, 9, 21, 29, 49, 61, 89, 105, 141, 161, 205, 229, 281, 309, 369, 401, 469, 505, 581, 621, 705, 749, 841, 889, 989, 1041, 1149, 1205, 1321, 1381, 1505, 1569, 1701, 1769, 1909, 1981, 2129, 2205, 2361, 2441, 2605, 2689, 2861, 2949, 3129, 3221, 3409, 3505, 3701, 3801, 4005, 4109, 4321, 4429, 4649, 4761, 4989, 5105, 5341, 5461
Offset: 0

Views

Author

Paul Curtz, Sep 18 2018

Keywords

Comments

The two bisections A136392(n+1)=1,9,29,61, ... and A201279(n)=5,21,49, ... are in the hexagonal spiral based on 2*n+1:
.
67--65--63--61
/ \
69 33--31--29 59
/ / \ \
71 35 11---9 27 57
/ / / \ \ \
73 37 13 1 7 25 55
/ / / / / /
39 15 3---5 23 53
\ \ / /
41 17--19--21 51
\ /
43--45--47--49
.
A201279(n) - A136892(n) = 20*n.

Crossrefs

In the spiral: A003154(n+1), A080859, A126587, A136392, A201279, A227776.
Partial sums of A382154.

Programs

  • Magma
    [(6*n^2 + 6*n + 5 - (2*n + 1)*(-1)^n)/4 : n in [0..50]]; // Wesley Ivan Hurt, Jan 19 2021
    
  • Mathematica
    Table[(6 n^2 + 6 n + 5 - (2 n + 1)*(-1)^n)/4, {n, 0, 80}] (* Wesley Ivan Hurt, Jan 07 2021 *)
  • PARI
    Vec((1 + x^2)*(1 + 4*x + x^2) / ((1 - x)^3*(1 + x)^2) + O(x^50)) \\ Colin Barker, Jun 05 2019
    
  • Python
    def A319384(n): return (n*(3*n+4)+3 if n&1 else n*(3*n+2)+2)>>1 # Chai Wah Wu, Mar 25 2025

Formula

a(2*n) = A136392(n+1), a(2*n+1) = A201279(n).
a(-n) = a(n).
a(2*n) + a(2*n+1) = 6*A001844(n).
a(n) = (6*n^2 + 6*n + 5 - (2*n + 1)*(-1)^n)/4. - Wesley Ivan Hurt, Oct 04 2018
G.f.: (1 + x^2)*(1 + 4*x + x^2) / ((1 - x)^3*(1 + x)^2). - Colin Barker, Jun 05 2019
a(n) = A104585(n) + A032766(n+1). - Alex W. Nowak, Jan 08 2021

Extensions

More terms from N. J. A. Sloane, Mar 23 2025
Previous Showing 51-60 of 79 results. Next