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 41 results. Next

A069190 Centered 24-gonal numbers.

Original entry on oeis.org

1, 25, 73, 145, 241, 361, 505, 673, 865, 1081, 1321, 1585, 1873, 2185, 2521, 2881, 3265, 3673, 4105, 4561, 5041, 5545, 6073, 6625, 7201, 7801, 8425, 9073, 9745, 10441, 11161, 11905, 12673, 13465, 14281, 15121, 15985, 16873, 17785, 18721, 19681, 20665, 21673
Offset: 1

Views

Author

Terrel Trotter, Jr., Apr 10 2002

Keywords

Comments

Sequence found by reading the line from 1, in the direction 1, 25, ..., in the square spiral whose vertices are the generalized octagonal numbers A001082. Semi-axis opposite to A135453 in the same spiral. - Omar E. Pol, Sep 16 2011

Examples

			a(5) = 241 because 12*5^2 - 12*5 + 1 = 300 - 60 + 1 = 241.
		

Crossrefs

Programs

Formula

a(n) = 12*n^2 - 12*n + 1.
a(n) = 24*n + a(n-1) - 24 with a(1)=1. - Vincenzo Librandi, Aug 08 2010
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3); a(1)=1, a(2)=25, a(3)=73. - Harvey P. Dale, Jul 17 2011
G.f.: x*(1+22*x+x^2)/(1-x)^3. - Harvey P. Dale, Jul 17 2011
Binomial transform of [1, 24, 24, 0, 0, 0, ...] and Narayana transform (cf. A001263) of [1, 24, 0, 0, 0, ...]. - Gary W. Adamson, Jul 26 2011
From Amiram Eldar, Jun 21 2020: (Start)
Sum_{n>=1} 1/a(n) = Pi*tan(Pi/sqrt(6))/(4*sqrt(6)).
Sum_{n>=1} a(n)/n! = 13*e - 1.
Sum_{n>=1} (-1)^n * a(n)/n! = 13/e - 1. (End)
E.g.f.: exp(x)*(1 + 12*x^2) - 1. - Stefano Spezia, May 31 2022

Extensions

More terms from Harvey P. Dale, Jul 17 2011

A090562 Primes of the form 5k^2 + 5k + 1.

Original entry on oeis.org

11, 31, 61, 101, 151, 211, 281, 661, 911, 1051, 1201, 1361, 1531, 1901, 2311, 2531, 3001, 3251, 3511, 4651, 5281, 6301, 6661, 7411, 9461, 9901, 12251, 13781, 14851, 15401, 18301, 18911, 19531, 20161, 22111, 24151, 24851, 25561, 27011, 27751
Offset: 1

Views

Author

Amarnath Murthy, Dec 11 2003

Keywords

Comments

Or, primes obtained as a concatenation of a triangular number and 1.
Centered decagonal primes. - Paul Muljadi, Oct 04 2005

Crossrefs

Programs

  • Magma
    [a: n in [0..100] | IsPrime(a) where a is 5*n^2 + 5*n + 1]; // Vincenzo Librandi, Dec 13 2011
  • Mathematica
    Select[5(#^2 - #) + 1 & /@ Range[75], PrimeQ[ # ] &] (* Robert G. Wilson v, Oct 10 2005 *)

Extensions

Edited and extended by Robert G. Wilson v, Oct 10 2005

A124080 10 times triangular numbers: a(n) = 5*n*(n + 1).

Original entry on oeis.org

0, 10, 30, 60, 100, 150, 210, 280, 360, 450, 550, 660, 780, 910, 1050, 1200, 1360, 1530, 1710, 1900, 2100, 2310, 2530, 2760, 3000, 3250, 3510, 3780, 4060, 4350, 4650, 4960, 5280, 5610, 5950, 6300, 6660, 7030, 7410, 7800, 8200, 8610, 9030, 9460, 9900, 10350
Offset: 0

Views

Author

Zerinvary Lajos, Nov 24 2006

Keywords

Comments

If Y is a 5-subset of an n-set X then, for n >= 5, a(n-4) is equal to the number of 5-subsets of X having exactly three elements in common with Y. Y is a 5-subset of an n-set X then, for n >= 6, a(n-6) is the number of (n-5)-subsets of X having exactly two elements in common with Y. - Milan Janjic, Dec 28 2007
Also sequence found by reading the line from 0, in the direction 0, 10, ... and the same line from 0, in the direction 0, 30, ..., in the square spiral whose vertices are the generalized dodecagonal numbers A195162. Axis perpendicular to A195148 in the same spiral. - Omar E. Pol, Sep 18 2011

Crossrefs

Programs

  • Magma
    [ 5*n*(n+1) : n in [0..50] ]; // Wesley Ivan Hurt, Jun 09 2014
    
  • Maple
    [seq(10*binomial(n,2),n=1..51)];
    seq(n*(n+1)*5, n=0..39); # Zerinvary Lajos, Mar 06 2007
  • Mathematica
    10*Accumulate[Range[0,50]] (* or *) LinearRecurrence[{3,-3,1},{0,10,30},50] (* Harvey P. Dale, Jul 21 2011 *)
  • PARI
    a(n)=5*n*(n+1) \\ Charles R Greathouse IV, Sep 28 2015

Formula

a(n) = 10*C(n,2), n >= 1.
a(n) = A049598(n) - A002378(n). - Zerinvary Lajos, Mar 06 2007
a(n) = 5*n*(n + 1), n >= 0. - Zerinvary Lajos, Mar 06 2007
a(n) = 5*n^2 + 5*n = 10*A000217(n) = 5*A002378(n) = 2*A028895(n). - Omar E. Pol, Dec 12 2008
a(n) = 10*n + a(n-1) (with a(0) = 0). - Vincenzo Librandi, Nov 12 2009
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), a(0) = 0, a(1) = 10, a(2) = 30. - Harvey P. Dale, Jul 21 2011
a(n) = A062786(n+1) - 1. - Omar E. Pol, Oct 03 2011
a(n) = A131242(10*n+9). - Philippe Deléham, Mar 27 2013
From G. C. Greubel, Aug 22 2017: (Start)
G.f.: 10*x/(1 - x)^3.
E.g.f.: 5*x*(x + 2)*exp(x). (End)
From Amiram Eldar, Sep 04 2022: (Start)
Sum_{n>=1} 1/a(n) = 1/5.
Sum_{n>=1} (-1)^(n+1)/a(n) = (2*log(2)-1)/5. (End)
From Amiram Eldar, Feb 21 2023: (Start)
Product_{n>=1} (1 - 1/a(n)) = -(5/Pi)*cos(3*Pi/(2*sqrt(5))).
Product_{n>=1} (1 + 1/a(n)) = (5/Pi)*cos(Pi/(2*sqrt(5))). (End)

A069131 Centered 18-gonal numbers.

Original entry on oeis.org

1, 19, 55, 109, 181, 271, 379, 505, 649, 811, 991, 1189, 1405, 1639, 1891, 2161, 2449, 2755, 3079, 3421, 3781, 4159, 4555, 4969, 5401, 5851, 6319, 6805, 7309, 7831, 8371, 8929, 9505, 10099, 10711, 11341, 11989, 12655, 13339, 14041, 14761, 15499, 16255, 17029, 17821
Offset: 1

Views

Author

Terrel Trotter, Jr., Apr 07 2002

Keywords

Comments

Equals binomial transform of [1, 18, 18, 0, 0, 0, ...]. Example: a(3) = 55 = (1, 2, 1) dot (1, 18, 18) = (1 + 36 + 18). - Gary W. Adamson, Aug 24 2010
Narayana transform (A001263) of [1, 18, 0, 0, 0, ...]. - Gary W. Adamson, Jul 28 2011
From Lamine Ngom, Aug 19 2021: (Start)
Sequence is a spoke of the hexagonal spiral built from the terms of A016777 (see illustration in links section).
a(n) is a bisection of A195042.
a(n) is a trisection of A028387.
a(n) + 1 is promic (A002378).
a(n) + 2 is a trisection of A002061.
a(n) + 9 is the arithmetic mean of its neighbors.
4*a(n) + 5 is a square: A016945(n)^2. (End)

Examples

			a(5) = 181 because 9*5^2 - 9*5 + 1 = 225 - 45 + 1 = 181.
		

Crossrefs

Programs

Formula

a(n) = 9*n^2 - 9*n + 1.
a(n) = 18*n + a(n-1) - 18 (with a(1)=1). - Vincenzo Librandi, Aug 08 2010
G.f.: ( x*(1+16*x+x^2) ) / ( (1-x)^3 ). - R. J. Mathar, Feb 04 2011
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3); a(1)=1, a(2)=19, a(3)=55. - Harvey P. Dale, Jan 20 2014
From Amiram Eldar, Jun 21 2020: (Start)
Sum_{n>=1} 1/a(n) = Pi*tan(sqrt(5)*Pi/6)/(3*sqrt(5)).
Sum_{n>=1} a(n)/n! = 10*e - 1.
Sum_{n>=1} (-1)^n * a(n)/n! = 10/e - 1. (End)
From Lamine Ngom, Aug 19 2021: (Start)
a(n) = 18*A000217(n) + 1 = 9*A002378(n) + 1.
a(n) = 3*A003215(n) - 2.
a(n) = A247792(n) - 9*n.
a(n) = A082040(n) + A304163(n) - a(n-1) = A016778(n) + A016790(n) - a(n-1), n > 0.
a(n) + a(n+1) = 2*A247792(n) = A010008(n), n > 0.
a(n+1) - a(n) = 18*n = A008600(n). (End)
From Leo Tavares, Oct 31 2021: (Start)
a(n)= A000290(n) + A139278(n-1)
a(n) = A069129(n) + A002378(n-1)
a(n) = A062786(n) + 8*A000217(n-1)
a(n) = A062786(n) + A033996(n-1)
a(n) = A060544(n) + 9*A000217(n-1)
a(n) = A060544(n) + A027468(n-1)
a(n) = A016754(n-1) + 10*A000217(n-1)
a(n) = A016754(n-1) + A124080
a(n) = A069099(n) + 11*A000217(n-1)
a(n) = A069099(n) + A152740(n-1)
a(n) = A003215(n-1) + 12*A000217(n-1)
a(n) = A003215(n-1) + A049598(n-1)
a(n) = A005891(n-1) + 13*A000217(n-1)
a(n) = A005891(n-1) + A152741(n-1)
a(n) = A001844(n) + 14*A000217(n-1)
a(n) = A001844(n) + A163756(n-1)
a(n) = A005448(n) + 15*A000217(n-1)
a(n) = A005448(n) + A194715(n-1). (End)
E.g.f.: exp(x)*(1 + 9*x^2) - 1. - Nikolaos Pantelidis, Feb 06 2023

A302329 a(0)=1, a(1)=61; for n>1, a(n) = 62*a(n-1) - a(n-2).

Original entry on oeis.org

1, 61, 3781, 234361, 14526601, 900414901, 55811197261, 3459393815281, 214426605350161, 13290990137894701, 823826961944121301, 51063980650397625961, 3165142973362708688281, 196187800367837541047461, 12160478479832564836254301, 753753477949251182306719201
Offset: 0

Views

Author

Bruno Berselli, Apr 05 2018

Keywords

Comments

Centered hexagonal numbers (A003215) with index in A145607. Example: 35 is a member of A145607, therefore A003215(35) = 3781 is a term of this sequence.
Also, centered 10-gonal numbers (A062786) with index in A182432. Example: 28 is a member of A182432 and A062786(28) = 3781.
a(n) is a solution to the Pell equation (4*a(n))^2 - 15*b(n)^2 = 1. The corresponding b(n) are A258684(n). - Klaus Purath, Jul 19 2025

Crossrefs

Fourth row of the array A188646.
First bisection of A041449, A042859.
Similar sequences of the type cosh((2*n+1)*arccosh(k))/k: A000012 (k=1), A001570 (k=2), A077420 (k=3), this sequence (k=4), A302330 (k=5), A302331 (k=6), A302332 (k=7), A253880 (k=8).

Programs

  • Mathematica
    LinearRecurrence[{62, -1}, {1, 61}, 20]
  • PARI
    x='x+O('x^99); Vec((1-x)/(1-62*x+x^2)) \\ Altug Alkan, Apr 06 2018

Formula

G.f.: (1 - x)/(1 - 62*x + x^2).
a(n) = a(-1-n).
a(n) = cosh((2*n + 1)*arccosh(4))/4.
a(n) = ((4 + sqrt(15))^(2*n + 1) + 1/(4 + sqrt(15))^(2*n + 1))/8.
a(n) = (1/4)*T(2*n+1, 4), where T(n,x) denotes the n-th Chebyshev polynomial of the first kind. - Peter Bala, Jul 08 2022
E.g.f.: exp(31*x)*(4*cosh(8*sqrt(15)*x) + sqrt(15)*sinh(8*sqrt(15)*x))/4. - Stefano Spezia, Jul 25 2025

A075117 Table by antidiagonals of generalized Lucas numbers: T(n,k) = T(n,k-1) + n*T(n,k-2) with T(n,0)=2 and T(n,1)=1.

Original entry on oeis.org

2, 1, 2, 1, 1, 2, 1, 3, 1, 2, 1, 4, 5, 1, 2, 1, 7, 7, 7, 1, 2, 1, 11, 17, 10, 9, 1, 2, 1, 18, 31, 31, 13, 11, 1, 2, 1, 29, 65, 61, 49, 16, 13, 1, 2, 1, 47, 127, 154, 101, 71, 19, 15, 1, 2, 1, 76, 257, 337, 297, 151, 97, 22, 17, 1, 2, 1, 123, 511, 799, 701, 506, 211, 127, 25, 19, 1, 2
Offset: 0

Views

Author

Henry Bottomley, Sep 02 2002

Keywords

Examples

			Array starts as:
  2, 1,  1,  1,  1,   1, ...;
  2, 1,  3,  4,  7,  11, ...;
  2, 1,  5,  7, 17,  31, ...;
  2, 1,  7, 10, 31,  61, ...;
  2, 1,  9, 13, 49, 101, ...;
  2, 1, 11, 16, 71, 151, ...; etc.
		

Crossrefs

Cf. A060959.
Columns include: A007395, A000012, A005408, A016777, A056220, A062786.

Programs

  • Magma
    [2^(1+k-n)*(&+[Binomial(n-k,2*j)*(1+4*k)^j: j in [0..Floor((n-k)/2)]]): k in [0..n], n in [0..13]]; // G. C. Greubel, Jan 27 2020
    
  • Maple
    seq(seq( 2^(1+k-n)*add( binomial(n-k, 2*j)*(1+4*k)^j, j=0..floor((n-k)/2)), k=0..n), n=0..13); # G. C. Greubel, Jan 27 2020
  • Mathematica
    T[n_, k_]:= ((1 + Sqrt[1+4n])/2)^k + ((1 - Sqrt[1+4n])/2)^k; Table[If[n==0 && k==0, 2, T[k, n-k]]//Simplify, {n,0,13}, {k,0,n}]//Flatten (* G. C. Greubel, Jan 27 2020 *)
  • Sage
    def T(n, k): return 2^(1-k)*sum( binomial(k, 2*j)*(1+4*n)^j for j in (0..floor(k/2)) )
    [[T(k,n-k) for k in (0..n)] for n in (0..13)] # G. C. Greubel, Jan 27 2020

Formula

T(n, k) = ((1+sqrt(4*n+1))/2)^k + ((1-sqrt(4*n+1))/2)^k = 2*A060959(n, k+1) - A060959(n, k).
T(n, k) = 2^(1-k)*Sum_{j=0..floor(k/2)} binomial(k, 2*j)*(1+4*n)^j. - G. C. Greubel, Jan 27 2020

A273366 a(n) = 10*n^2 + 10*n + 2.

Original entry on oeis.org

2, 22, 62, 122, 202, 302, 422, 562, 722, 902, 1102, 1322, 1562, 1822, 2102, 2402, 2722, 3062, 3422, 3802, 4202, 4622, 5062, 5522, 6002, 6502, 7022, 7562, 8122, 8702, 9302, 9922, 10562, 11222, 11902, 12602, 13322, 14062, 14822, 15602
Offset: 0

Views

Author

Keywords

Comments

These are the numbers k such that 10*k+5 is a perfect square.

Crossrefs

Cf. A033583 (perfect squares ending in 0 in base 10 with final 0 removed).

Programs

Formula

G.f.: 2*(x^2+8x+1)/(1-x)^3.
From G. C. Greubel, May 20 2016: (Start)
E.g.f.: 2*(1 + 10*x + 5*x^2)*exp(x).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). (End)
a(n) = 2*A062786(n+1). - R. J. Mathar, Jun 03 2016
Sum_{n>=0} 1/a(n) = Pi/(2*sqrt(5)) * tan(Pi/(2*sqrt(5))) (A350760). - Amiram Eldar, Jan 20 2022

A247643 a(n) = ( 10*n*(n+1)+(2*n+1)*(-1)^n+7 )/8.

Original entry on oeis.org

1, 3, 9, 15, 27, 37, 55, 69, 93, 111, 141, 163, 199, 225, 267, 297, 345, 379, 433, 471, 531, 573, 639, 685, 757, 807, 885, 939, 1023, 1081, 1171, 1233, 1329, 1395, 1497, 1567, 1675, 1749, 1863, 1941, 2061, 2143, 2269, 2355, 2487, 2577, 2715, 2809, 2953, 3051
Offset: 0

Views

Author

N. J. A. Sloane, Sep 23 2014

Keywords

Comments

From Paul Curtz, Jan 01 2020: (Start)
In the following pentagonal spiral of odd numbers
101
99 61 63
97 59 31 33 65
95 57 29 11 13 35 67
93 55 27 9 1 3 15 37 69
91 53 25 7 5 17 39 71
89 51 23 21 19 41 73
87 49 47 45 43 75
85 83 81 79 77
the terms of this sequence appear on the x axis. A062786 and A172043 are in the spiral as well. (End)

Crossrefs

A diagonal of triangle in A247646.

Programs

  • Maple
    f:=n->(10*n*(n+1)+(2*n+1)*(-1)^n+7)/8;
  • Mathematica
    Table[(10 n (n + 1) + (2 n + 1) (-1)^n + 7)/8, {n, 0, 60}] (* Vincenzo Librandi, Sep 26 2014 *)
  • PARI
    Vec(-(x^4+2*x^3+4*x^2+2*x+1) / ((x-1)^3*(x+1)^2) + O(x^100)) \\ Colin Barker, Sep 25 2014

Formula

a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5). - Colin Barker, Sep 25 2014
G.f.: -(x^4+2*x^3+4*x^2+2*x+1) / ((x-1)^3*(x+1)^2). - Colin Barker, Sep 25 2014
From Paul Curtz, Jan 01 2020: (Start)
a(n) = 1 + 2*A085787(n).
a(n+1) = a(n-1) + A090772(n+1). (End)
E.g.f.: (1/4)*((1 + x)*(4 + 5*x)*cosh(x) + (3 + x*(11 + 5*x))*sinh(x)). - Stefano Spezia, Jan 01 2020

Extensions

More terms from Colin Barker, Sep 25 2014

A128922 Numbers simultaneously 10-gonal and centered 10-gonal.

Original entry on oeis.org

1, 451, 145351, 46802701, 15070324501, 4852597686751, 1562521384809451, 503127033310956601, 162005342204743216201, 52165217062894004660251, 16797037888909664757384751
Offset: 0

Views

Author

Steven Schlicker, Apr 24 2007

Keywords

Examples

			a(1) = 451 because 451 is the tenth centered 10-gonal number and the eleventh 10-gonal number.
		

Crossrefs

Programs

  • Maple
    CP := n -> 1+1/2*10*(n^2-n): N:=10: u:=9: v:=1: x:=10: y:=1: k_pcp:=[1]: for i from 1 to N do tempx:=x; tempy:=y; x:=tempx*u+80*tempy*v: y:=tempx*v+tempy*u: s:=(y+1)/2: k_pcp:=[op(k_pcp),CP(s)]: end do: k_pcp;

Formula

Let x(n) + y(n)*sqrt(80) =: (10+sqrt(80))*(9+sqrt(80))^n, s(n) = (y(n)+1)/2; then a(n) = (1/2)*(2+10*(s(n)^2-s(n))).
From Richard Choulet, Oct 01 2007: (Start)
a(n+2) = 322*a(n+1)-a(n)+130.
a(n+1) = 161*a(n)+65+9*(320*a(n)^2+260*a(n)+45)^0.5.
G.f.: z*(1+128*z+z^2)/((1-z)*(1-322*z+z^2)). (End)

A133273 Indices of centered decagonal numbers which are also decagonal numbers.

Original entry on oeis.org

1, 10, 171, 3060, 54901, 985150, 17677791, 317215080, 5692193641, 102142270450, 1832868674451, 32889493869660, 590178020979421, 10590314883759910, 190035489886698951, 3410048503076821200, 61190837565496082641, 1098025027675852666330, 19703259660599851911291
Offset: 1

Views

Author

Richard Choulet, Oct 16 2007

Keywords

Comments

Numbers k such that 80*k^2 - 80*k + 25 is a square.
Also the indices of centered square numbers which are also centered pentagonal numbers. - Colin Barker, Jan 01 2015

Crossrefs

Programs

  • Mathematica
    LinearRecurrence[{19,-19,1},{1,10,171},20] (* Harvey P. Dale, Oct 09 2020 *)
  • PARI
    Vec(x*(-1+9*x)/((-1+x)*(1-18*x+x^2)) + O(x^100)) \\ Colin Barker, Jan 01 2015

Formula

a(n+2) = 18*a(n+1) - a(n) - 8.
a(n+1) = 9*a(n) - 4 + sqrt(80*a(n)^2 - 80*a(n) + 25).
G.f.: x*(-1+9*x)/(-1+x)/(1 - 18*x + x^2). - R. J. Mathar, Nov 14 2007
a(n) = 19*a(n-1) - 19*a(n-2) + a(n-3). - Colin Barker, Jan 01 2015
Product_{n>=2} (1 - 1/a(n)) = 2/sqrt(5) (= A010532 / 10). - Amiram Eldar, Dec 02 2024

Extensions

More terms from Paolo P. Lava, Nov 25 2008
Previous Showing 11-20 of 41 results. Next