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-10 of 13 results. Next

A139597 A139278(n) followed by A139274(n+1).

Original entry on oeis.org

0, 7, 15, 30, 46, 69, 93, 124, 156, 195, 235, 282, 330, 385, 441, 504, 568, 639, 711, 790, 870, 957, 1045, 1140, 1236, 1339, 1443, 1554, 1666, 1785, 1905, 2032, 2160, 2295, 2431, 2574, 2718, 2869, 3021, 3180, 3340, 3507, 3675, 3850
Offset: 0

Views

Author

Omar E. Pol, May 03 2008

Keywords

Comments

Sequence found by reading the line from 0, in the direction 0, 7,... and the line from 15, in the direction 15, 46,..., in the square spiral whose vertices are the triangular numbers A000217.

Examples

			Array begins:
0, 7
15, 30
46, 69
93, 124
		

Crossrefs

Formula

Array read by rows: row n gives 8*n^2 + 7n, 8*(n+1)^2 - (n+1).
a(n) = (3-3*(-1)^n+14*n+8*n^2)/4. a(n) = 2*a(n-1)-2*a(n-3)+a(n-4). G.f.: x*(7+x)/((1-x)^3*(1+x)). [Colin Barker, Jul 22 2012]

A274979 Integers of the form m*(m + 7)/8.

Original entry on oeis.org

0, 1, 15, 18, 46, 51, 93, 100, 156, 165, 235, 246, 330, 343, 441, 456, 568, 585, 711, 730, 870, 891, 1045, 1068, 1236, 1261, 1443, 1470, 1666, 1695, 1905, 1936, 2160, 2193, 2431, 2466, 2718, 2755, 3021, 3060, 3340, 3381, 3675, 3718, 4026, 4071, 4393, 4440, 4776, 4825
Offset: 1

Views

Author

Bruno Berselli, Jul 15 2016

Keywords

Comments

Nonnegative values of m are listed in A047393.
Also, numbers h such that 32*h + 49 is a square.
Equivalently, numbers of the form i*(8*i + 7) with i = 0, -1, 1, -2, 2, -3, 3, ...
Infinitely many squares belong to this sequence.
The first bisection is A139278, and 0 followed by the second bisection gives A051870.
Generalized 18-gonal (or octadecagonal) numbers (see the third comment). - Omar E. Pol, Jun 06 2018
Partial sums of A317314. - Omar E. Pol, Jul 28 2018
Exponents in expansion of Product_{n >= 1} (1 + x^(16*n-15))*(1 + x^(16*n-1))*(1 - x^(16*n)) = 1 + x + x^15 + x^18 + x^46 + .... - Peter Bala, Dec 10 2020
Generalized k-gonal numbers are second k-gonal numbers and positive terms of k-gonal numbers interleaved, k >= 5. They are also the partial sums of the sequence formed by the multiples of (k - 4) and the odd numbers (A005408) interleaved, k >= 5. In this case k = 18. - Omar E. Pol, Apr 25 2021

Examples

			100 is in the sequence because 100 = 25*(25+7)/8 or also 100 = 4*(8*4-7).
From _Omar E. Pol_, Apr 24 2021: (Start)
Illustration of initial terms as vertices of a rectangular spiral:
        46_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _18
         |                                                       |
         |                           0                           |
         |                           |_ _ _ _ _ _ _ _ _ _ _ _ _ _|
         |                           1                           15
         |
        51
More generally, all generalized k-gonal numbers can be represented with this kind of spirals, k >= 5. In this case  k = 18. (End)
		

Crossrefs

Cf. sequences of the form m*(m+k)/(k+1) listed in A274978.
Cf. similar sequences listed in A299645.
Cf. A317314.
Sequences of generalized k-gonal numbers: A001318 (k=5), A000217 (k=6), A085787 (k=7), A001082 (k=8), A118277 (k=9), A074377 (k=10), A195160 (k=11), A195162 (k=12), A195313 (k=13), A195818 (k=14), A277082 (k=15), A274978 (k=16), A303305 (k=17), this sequence (k=18), A303813 (k=19), A218864 (k=20), A303298 (k=21), A303299 (k=22), A303303 (k=23), A303814 (k=24), A303304 (k=25), A316724 (k=26), A316725 (k=27), A303812 (k=28), A303815 (k=29), A316729 (k=30).

Programs

  • Magma
    [t: m in [0..200] | IsIntegral(t) where t is m*(m+7)/8];
    
  • Mathematica
    Select[m = Range[0, 200]; m (m + 7)/8, IntegerQ] (* Jean-François Alcover, Jul 21 2016 *)
    Select[Table[(m(m+7))/8,{m,0,200}],IntegerQ] (* or *) LinearRecurrence[ {1,2,-2,-1,1},{0,1,15,18,46},50] (* Harvey P. Dale, May 07 2019 *)
  • Python
    def A274979(n): return (n>>1)*((n<<2)+(3 if n&1 else -7)) # Chai Wah Wu, Mar 11 2025
  • Sage
    def A274979_list(len):
        h = lambda m: m*(m+7)/8
        return [h(m) for m in (0..len) if h(m) in ZZ]
    print(A274979_list(199)) # Peter Luschny, Jul 18 2016
    

Formula

O.g.f.: x^2*(1 + 14*x + x^2)/((1 + x)^2*(1 - x)^3).
E.g.f.: (3*(2*x + 1)*exp(-x) + (8*x^2 - 3)*exp(x))/4.
a(n) = (8*(n-1)*n - 3*(2*n-1)*(-1)^n - 3)/4.
a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5) for n >= 6. - Wesley Ivan Hurt, Dec 18 2020
From Amiram Eldar, Feb 28 2022: (Start)
Sum_{n>=2} 1/a(n) = (8 + 7*(sqrt(2)+1)*Pi)/49.
Sum_{n>=2} (-1)^n/a(n) = 8*log(2)/7 + 2*sqrt(2)*log(sqrt(2)+1)/7 - 8/49. (End)
a(n) = (n-1)*(4*n+3)/2 if n is odd and a(n) = n*(4*n-7)/2 if n is even. - Chai Wah Wu, Mar 11 2025

A051870 18-gonal (or octadecagonal) numbers: a(n) = n*(8*n-7).

Original entry on oeis.org

0, 1, 18, 51, 100, 165, 246, 343, 456, 585, 730, 891, 1068, 1261, 1470, 1695, 1936, 2193, 2466, 2755, 3060, 3381, 3718, 4071, 4440, 4825, 5226, 5643, 6076, 6525, 6990, 7471, 7968, 8481, 9010, 9555, 10116, 10693, 11286, 11895, 12520
Offset: 0

Views

Author

N. J. A. Sloane, Dec 15 1999

Keywords

Comments

Also, sequence found by reading the segment (0, 1) together with the line from 1, in the direction 1, 18, ..., in the square spiral whose vertices are the triangular numbers A000217. - Omar E. Pol, Apr 26 2008
This sequence does not contain any triangular numbers other than 0 and 1. See A188892. - T. D. Noe, Apr 13 2011
Also sequence found by reading the line from 0, in the direction 0, 18, ... and the parallel line from 1, in the direction 1, 51, ..., in the square spiral whose vertices are the generalized 18-gonal numbers. - Omar E. Pol, Jul 18 2012
Partial sums of 16n + 1 (with offset 0), compare A005570. - Jeremy Gardiner, Aug 04 2012
All x values for Diophantine equation 32*x + 49 = y^2 are given by this sequence and A139278. - Bruno Berselli, Nov 11 2014
This is also a star enneagonal number: a(n) = A001106(n) + 9*A000217(n-1). - Luciano Ancora, Mar 30 2015

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, p. 189.
  • Elena Deza and Michel Marie Deza, Figurate numbers, World Scientific Publishing, 2012, page 6.

Crossrefs

Programs

Formula

G.f.: x*(1+15*x)/(1-x)^3. - Bruno Berselli, Feb 04 2011
a(n) = 16*n + a(n-1) - 15, with n > 0, a(0) = 0. - Vincenzo Librandi, Aug 06 2010
a(16*a(n)+121*n+1) = a(16*a(n)+121*n) + a(16*n+1). - Vladimir Shevelev, Jan 24 2014
E.g.f.: (8*x^2 + x)*exp(x). - G. C. Greubel, Jul 18 2017
Sum_{n>=1} 1/a(n) = ((1+sqrt(2))*Pi + 2*sqrt(2)*arccoth(sqrt(2)) + 8*log(2))/14. - Amiram Eldar, Oct 20 2020
Product_{n>=2} (1 - 1/a(n)) = 8/9. - Amiram Eldar, Jan 22 2021

A139275 a(n) = n*(8*n+1).

Original entry on oeis.org

0, 9, 34, 75, 132, 205, 294, 399, 520, 657, 810, 979, 1164, 1365, 1582, 1815, 2064, 2329, 2610, 2907, 3220, 3549, 3894, 4255, 4632, 5025, 5434, 5859, 6300, 6757, 7230, 7719, 8224, 8745, 9282, 9835, 10404, 10989, 11590, 12207, 12840
Offset: 0

Views

Author

Omar E. Pol, Apr 26 2008

Keywords

Comments

Sequence found by reading the line from 0, in the direction 0, 9,..., in the square spiral whose vertices are the triangular numbers A000217.

Crossrefs

Programs

  • Mathematica
    Table[n (8 n + 1), {n, 0, 40}] (* Bruno Berselli, Sep 21 2016 *)
    LinearRecurrence[{3,-3,1},{0,9,34},50] (* Harvey P. Dale, Apr 21 2020 *)
  • PARI
    a(n) = n*(8*n+1); \\ Altug Alkan, Sep 21 2016

Formula

a(n) = 8*n^2 + n.
Sequences of the form a(n) = 8*n^2+c*n have generating functions x{c+8+(8-c)x} / (1-x)^3 and recurrence a(n) = 3*a(n-1)-3*a(n-2)+a(n-3). The inverse binomial transform is 0, c+8, 16, 0, 0, ... (0 continued). This applies to A139271-A139278, positive or negative c. - R. J. Mathar, May 12 2008
a(n) = 16*n + a(n-1) - 7 with n>0, a(0)=0. - Vincenzo Librandi, Aug 03 2010
a(n) = A000217(5*n) - A000217(3*n). - Bruno Berselli, Sep 21 2016
Sum_{n>=1} 1/a(n) = 8 - (1+sqrt(2))*Pi/2 - 4*log(2) - sqrt(2) * log(1+sqrt(2)) = 0.1887230016056779928... . - Vaclav Kotesovec, Sep 21 2016
From G. C. Greubel, Jul 18 2017: (Start)
G.f.: x*(7*x + 9)/(1-x)^3.
E.g.f.: (8*x^2 + 9*x)*exp(x). (End)

A139271 a(n) = 2*n*(4*n-3).

Original entry on oeis.org

0, 2, 20, 54, 104, 170, 252, 350, 464, 594, 740, 902, 1080, 1274, 1484, 1710, 1952, 2210, 2484, 2774, 3080, 3402, 3740, 4094, 4464, 4850, 5252, 5670, 6104, 6554, 7020, 7502, 8000, 8514, 9044, 9590, 10152, 10730, 11324, 11934, 12560, 13202, 13860, 14534, 15224
Offset: 0

Views

Author

Omar E. Pol, Apr 26 2008

Keywords

Comments

Sequence found by reading the line from 0, in the direction 0, 2, ..., in the square spiral whose vertices are the triangular numbers A000217. Opposite numbers to the members of A033585 in the same spiral.
Twice decagonal numbers (or twice 10-gonal numbers). - Omar E. Pol, May 15 2008
a(n) is the number of walks in a cubic lattice of n dimensions that reach the point of origin for the first time after 4 steps. - Shel Kaphan, Mar 20 2023

Crossrefs

Cf. A001107.
Cf. numbers of the form n*(n*k-k+4)/2 listed in A226488 (this sequence is the case k=16). - Bruno Berselli, Jun 10 2013
Row n=2 of A361397.

Programs

Formula

a(n) = 8*n^2 - 6*n.
Sequences of the form a(n) = 8*n^2 + c*n have generating functions x{c+8+(8-c)x} / (1-x)^3 and recurrence a(n) = 3a(n-1) - 3a(n-2) + a(n-3). The inverse binomial transform is 0, c+8, 16, 0, 0, ... (0 continued). This applies to A139271-A139278, positive or negative c. - R. J. Mathar, May 12 2008
a(n) = A001107(n)*2. - Omar E. Pol, May 15 2008
a(n) = 16*n + a(n-1) - 14 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
From G. C. Greubel, Jul 18 2017: (Start)
G.f.: (2*x)*(7*x+1)/(1-x)^3.
E.g.f.: (8*x^2 + 2*x)*exp(x). (End)
Sum_{n>=1} 1/a(n) = Pi/12 + log(2)/2. - Amiram Eldar, Mar 28 2023

Extensions

Corrected by Harvey P. Dale, Sep 26 2016

A139273 a(n) = n*(8*n - 3).

Original entry on oeis.org

0, 5, 26, 63, 116, 185, 270, 371, 488, 621, 770, 935, 1116, 1313, 1526, 1755, 2000, 2261, 2538, 2831, 3140, 3465, 3806, 4163, 4536, 4925, 5330, 5751, 6188, 6641, 7110, 7595, 8096, 8613, 9146, 9695, 10260, 10841, 11438, 12051, 12680
Offset: 0

Views

Author

Omar E. Pol, Apr 26 2008

Keywords

Comments

Sequence found by reading the line from 0, in the direction 0, 5, ..., in the square spiral whose vertices are the triangular numbers A000217. Opposite numbers to the members of A139277 in the same spiral.
Also, sequence of numbers of the form d*A000217(n-1) + 5*n with generating functions x*(5+(d-5)*x)/(1-x)^3; the inverse binomial transform is 0,5,d,0,0,.. (0 continued). See Crossrefs. - Bruno Berselli, Feb 11 2011
Even decagonal numbers divided by 2. - Omar E. Pol, Aug 19 2011

Crossrefs

Programs

  • Magma
    [ n*(8*n-3) : n in [0..40] ];  // Bruno Berselli, Feb 11 2011
    
  • Mathematica
    Table[n (8 n - 3), {n, 0, 40}] (* or *) LinearRecurrence[{3, -3, 1}, {0, 5, 26}, 40] (* Harvey P. Dale, Feb 02 2012 *)
  • PARI
    a(n)=n*(8*n-3) \\ Charles R Greathouse IV, Sep 24 2015

Formula

a(n) = 8*n^2 - 3*n.
Sequences of the form a(n) = 8*n^2 + c*n have generating functions x{c+8+(8-c)x} / (1-x)^3 and recurrence a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). The inverse binomial transform is 0, c+8, 16, 0, 0, ... (0 continued). This applies to A139271-A139278, positive or negative c. - R. J. Mathar, May 12 2008
a(n) = 16*n + a(n-1) - 11 for n>0, a(0)=0. - Vincenzo Librandi, Aug 03 2010
From Bruno Berselli, Feb 11 2011: (Start)
G.f.: x*(5 + 11*x)/(1 - x)^3.
a(n) = 4*A000217(n) + A051866(n). (End)
a(n) = A028994(n)/2. - Omar E. Pol, Aug 19 2011
a(0)=0, a(1)=5, a(2)=26; for n>2, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, Feb 02 2012
E.g.f.: (8*x^2 + 5*x)*exp(x). - G. C. Greubel, Jul 18 2017
Sum_{n>=1} 1/a(n) = 4*log(2)/3 - (sqrt(2)-1)*Pi/6 - sqrt(2)*arccoth(sqrt(2))/3. - Amiram Eldar, Jul 03 2020

A139272 a(n) = n*(8*n-5).

Original entry on oeis.org

0, 3, 22, 57, 108, 175, 258, 357, 472, 603, 750, 913, 1092, 1287, 1498, 1725, 1968, 2227, 2502, 2793, 3100, 3423, 3762, 4117, 4488, 4875, 5278, 5697, 6132, 6583, 7050, 7533, 8032, 8547, 9078, 9625, 10188, 10767, 11362, 11973, 12600
Offset: 0

Views

Author

Omar E. Pol, Apr 26 2008

Keywords

Comments

Sequence found by reading the line from 0, in the direction 0, 3, ..., in the square spiral whose vertices are the triangular numbers A000217. Opposite numbers to the members of A139276 in the same spiral.

Crossrefs

Cf. numbers of the form n*(n*k-k+6)/2, this sequence is the case k=16: see Comments lines of A226492.

Programs

Formula

a(n) = 8*n^2 - 5*n.
Sequences of the form a(n) = 8*n^2 + c*n have generating functions x{c+8+(8-c)x} / (1-x)^3 and recurrence a(n) = 3a(n-1) - 3a(n-2) + a(n-3). The inverse binomial transform is 0, c+8, 16, 0, 0, ... (0 continued). This applies to A139271-A139278, positive or negative c. - R. J. Mathar, May 12 2008
a(n) = 16*n + a(n-1) - 13 with n>0, a(0)=0. - Vincenzo Librandi, Aug 03 2010
From G. C. Greubel, Jul 18 2017: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
G.f.: x*(13*x + 3)/(1-x)^3.
E.g.f.: (8*x^2 + 3*x)*exp(x). (End)
Sum_{n>=1} 1/a(n) = ((sqrt(2)-1)*Pi + 8*log(2) - 2*sqrt(2)*log(sqrt(2)+1))/10. - Amiram Eldar, Mar 17 2022

A139274 a(n) = n*(8*n-1).

Original entry on oeis.org

0, 7, 30, 69, 124, 195, 282, 385, 504, 639, 790, 957, 1140, 1339, 1554, 1785, 2032, 2295, 2574, 2869, 3180, 3507, 3850, 4209, 4584, 4975, 5382, 5805, 6244, 6699, 7170, 7657, 8160, 8679, 9214, 9765, 10332, 10915, 11514, 12129, 12760, 13407, 14070, 14749
Offset: 0

Views

Author

Omar E. Pol, Apr 26 2008

Keywords

Comments

Sequence found by reading the line from 0, in the direction 0, 7, ..., in the square spiral whose vertices are the triangular numbers A000217.
Polygonal number connection: 2*P_n + 5*S_n where P_n is the n-th pentagonal number and S_n is the n-th square. - William A. Tedeschi, Sep 12 2010

Examples

			a(1) = 16*1 + 0 - 9 = 7; a(2) = 16*2 + 7 - 9 = 30; a(3) = 16*3 + 30 - 9 = 69. - _Vincenzo Librandi_, Aug 03 2010
		

Crossrefs

Programs

Formula

Sequences of the form a(n) = 8*n^2 + c*n have generating functions x{c+8+(8-c)x} / (1-x)^3 and recurrence a(n) = 3a(n-1) - 3a(n-2) + a(n-3). The inverse binomial transform is 0, c+8, 16, 0, 0, ... (0 continued). This applies to A139271-A139278, positive or negative c. - R. J. Mathar, May 12 2008
a(n) = 16*n + a(n-1) - 9 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
a(n) = (1/3) * Sum_{i=n..(7*n-1)} i. - Wesley Ivan Hurt, Dec 04 2016
From G. C. Greubel, Jul 18 2017: (Start)
G.f.: x*(9*x+7)/(1-x)^3.
E.g.f.: (8*x^2 + 7*x)*exp(x). (End)
Sum_{n>=1} 1/a(n) = 4*log(2) + sqrt(2)*log(sqrt(2)+1) - (sqrt(2)+1)*Pi/2. - Amiram Eldar, Mar 18 2022

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

A139276 a(n) = n*(8*n+3).

Original entry on oeis.org

0, 11, 38, 81, 140, 215, 306, 413, 536, 675, 830, 1001, 1188, 1391, 1610, 1845, 2096, 2363, 2646, 2945, 3260, 3591, 3938, 4301, 4680, 5075, 5486, 5913, 6356, 6815, 7290, 7781, 8288, 8811, 9350, 9905, 10476, 11063, 11666, 12285, 12920
Offset: 0

Views

Author

Omar E. Pol, Apr 26 2008

Keywords

Comments

Sequence found by reading the line from 0, in the direction 0, 11,..., in the square spiral whose vertices are the triangular numbers A000217. Opposite numbers to the members of A139272 in the same spiral.

Examples

			a(1)=16*1+0-5=11; a(2)=16*2+11-5=38; a(3)=16*3+38-5=81. - _Vincenzo Librandi_, Aug 03 2010
		

Crossrefs

Programs

Formula

a(n) = 8*n^2 + 3*n.
Sequences of the form a(n)=8*n^2+c*n have generating functions x{c+8+(8-c)x} / (1-x)^3 and recurrence a(n)= 3a(n-1)-3a(n-2)+a(n-3). The inverse binomial transform is 0, c+8, 16, 0, 0, ... (0 continued). This applies to A139271-A139278, positive or negative c. - R. J. Mathar, May 12 2008
a(n) = 16*n+a(n-1)-5 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
From G. C. Greubel, Jul 18 2017: (Start)
G.f.: x*(5*x + 11)/(1-x)^3.
E.g.f.: (8*x^2 + 11*x)*exp(x). (End)
Sum_{n>=1} 1/a(n) = 8/9 - (sqrt(2)-1)*Pi/6 - 4*log(2)/3 + sqrt(2)*log(sqrt(2)+1)/3. - Amiram Eldar, Mar 17 2022
Showing 1-10 of 13 results. Next