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

A016754 Odd squares: a(n) = (2n+1)^2. Also centered octagonal numbers.

Original entry on oeis.org

1, 9, 25, 49, 81, 121, 169, 225, 289, 361, 441, 529, 625, 729, 841, 961, 1089, 1225, 1369, 1521, 1681, 1849, 2025, 2209, 2401, 2601, 2809, 3025, 3249, 3481, 3721, 3969, 4225, 4489, 4761, 5041, 5329, 5625, 5929, 6241, 6561, 6889, 7225, 7569, 7921, 8281, 8649, 9025
Offset: 0

Views

Author

Keywords

Comments

The brown rat (rattus norwegicus) breeds very quickly. It can give birth to other rats 7 times a year, starting at the age of three months. The average number of pups is 8. The present sequence gives the total number of rats, when the intervals are 12/7 of a year and a young rat starts having offspring at 24/7 of a year. - Hans Isdahl, Jan 26 2008
Numbers n such that tau(n) is odd where tau(x) denotes the Ramanujan tau function (A000594). - Benoit Cloitre, May 01 2003
If Y is a fixed 2-subset of a (2n+1)-set X then a(n-1) is the number of 3-subsets of X intersecting Y. - Milan Janjic, Oct 21 2007
Binomial transform of [1, 8, 8, 0, 0, 0, ...]; Narayana transform (A001263) of [1, 8, 0, 0, 0, ...]. - Gary W. Adamson, Dec 29 2007
All terms of this sequence are of the form 8k+1. For numbers 8k+1 which aren't squares see A138393. Numbers 8k+1 are squares iff k is a triangular number from A000217. And squares have form 4n(n+1)+1. - Artur Jasinski, Mar 27 2008
Sequence arises from reading the line from 1, in the direction 1, 25, ... and the line from 9, in the direction 9, 49, ..., in the square spiral whose vertices are the squares A000290. - Omar E. Pol, May 24 2008
Equals the triangular numbers convolved with [1, 6, 1, 0, 0, 0, ...]. - Gary W. Adamson & Alexander R. Povolotsky, May 29 2009
First differences: A008590(n) = a(n) - a(n-1) for n>0. - Reinhard Zumkeller, Nov 08 2009
Central terms of the triangle in A176271; cf. A000466, A053755. - Reinhard Zumkeller, Apr 13 2010
Odd numbers with odd abundance. Odd numbers with even abundance are in A088828. Even numbers with odd abundance are in A088827. Even numbers with even abundance are in A088829. - Jaroslav Krizek, May 07 2011
Appear as numerators in the non-simple continued fraction expansion of Pi-3: Pi-3 = K_{k>=1} (1-2*k)^2/6 = 1/(6+9/(6+25/(6+49/(6+...)))), see also the comment in A007509. - Alexander R. Povolotsky, Oct 12 2011
Ulam's spiral (SE spoke). - Robert G. Wilson v, Oct 31 2011
All terms end in 1, 5 or 9. Modulo 100, all terms are among { 1, 9, 21, 25, 29, 41, 49, 61, 69, 81, 89 }. - M. F. Hasler, Mar 19 2012
Right edge of both triangles A214604 and A214661: a(n) = A214604(n+1,n+1) = A214661(n+1,n+1). - Reinhard Zumkeller, Jul 25 2012
Also: Odd numbers which have an odd sum of divisors (= sigma = A000203). - M. F. Hasler, Feb 23 2013
Consider primitive Pythagorean triangles (a^2 + b^2 = c^2, gcd(a, b) = 1) with hypotenuse c (A020882) and respective even leg b (A231100); sequence gives values c-b, sorted with duplicates removed. - K. G. Stier, Nov 04 2013
For n>1 a(n) is twice the area of the irregular quadrilateral created by the points ((n-2)*(n-1),(n-1)*n/2), ((n-1)*n/2,n*(n+1)/2), ((n+1)*(n+2)/2,n*(n+1)/2), and ((n+2)*(n+3)/2,(n+1)*(n+2)/2). - J. M. Bergot, May 27 2014
Number of pairs (x, y) of Z^2, such that max(abs(x), abs(y)) <= n. - Michel Marcus, Nov 28 2014
Except for a(1)=4, the number of active (ON, black) cells in n-th stage of growth of two-dimensional cellular automaton defined by "Rule 737", based on the 5-celled von Neumann neighborhood. - Robert Price, May 23 2016
a(n) is the sum of 2n+1 consecutive numbers, the first of which is n+1. - Ivan N. Ianakiev, Dec 21 2016
a(n) is the number of 2 X 2 matrices with all elements in {0..n} with determinant = 2*permanent. - Indranil Ghosh, Dec 25 2016
Engel expansion of Pi*StruveL_0(1)/2 where StruveL_0(1) is A197037. - Benedict W. J. Irwin, Jun 21 2018
Consider all Pythagorean triples (X,Y,Z=Y+1) ordered by increasing Z; the segments on the hypotenuse {p = a(n)/A001844(n), q = A060300(n)/A001844(n) = A001844(n) - p} and their ratio p/q = a(n)/A060300(n) are irreducible fractions in Q\Z. X values are A005408, Y values are A046092, Z values are A001844. - Ralf Steiner, Feb 25 2020
a(n) is the number of large or small squares that are used to tile primitive squares of type 2 (A344332). - Bernard Schott, Jun 03 2021
Also, positive odd integers with an odd number of odd divisors (for similar sequence with 'even', see A348005). - Bernard Schott, Nov 21 2021
a(n) is the least odd number k = x + y, with 0 < x < y, such that there are n distinct pairs (x,y) for which x*y/k is an integer; for example, a(2) = 25 and the two corresponding pairs are (5,20) and (10,15). The similar sequence with 'even' is A016742 (see Comment of Jan 26 2018). - Bernard Schott, Feb 24 2023
From Peter Bala, Jan 03 2024: (Start)
The sequence terms are the exponents of q in the series expansions of the following infinite products:
1) q*Product_{n >= 1} (1 - q^(16*n))*(1 + q^(8*n)) = q + q^9 + q^25 + q^49 + q^81 + q^121 + q^169 + ....
2) q*Product_{n >= 1} (1 + q^(16*n))*(1 - q^(8*n)) = q - q^9 - q^25 + q^49 + q^81 - q^121 - q^169 + + - - ....
3) q*Product_{n >= 1} (1 - q^(8*n))^3 = q - 3*q^9 + 5*q^25 - 7*q^49 + 9*q^81 - 11*q^121 + 13*q^169 - + ....
4) q*Product_{n >= 1} ( (1 + q^(8*n))*(1 - q^(16*n))/(1 + q^(16*n)) )^3 = q + 3*q^9 - 5*q^25 - 7*q^49 + 9*q^81 + 11*q^121 - 13*q^169 - 15*q^225 + + - - .... (End)

References

  • L. Lorentzen and H. Waadeland, Continued Fractions with Applications, North-Holland 1992, p. 586.

Crossrefs

Cf. A000447 (partial sums).
Cf. A348005, A379481 [= a(A048673(n)-1)].
Partial sums of A022144.
Positions of odd terms in A341528.
Sequences on the four axes of the square spiral: Starting at 0: A001107, A033991, A007742, A033954; starting at 1: A054552, A054556, A054567, A033951.
Sequences on the four diagonals of the square spiral: Starting at 0: A002939 = 2*A000384, A016742 = 4*A000290, A002943 = 2*A014105, A033996 = 8*A000217; starting at 1: A054554, A053755, A054569, A016754.
Sequences obtained by reading alternate terms on the X and Y axes and the two main diagonals of the square spiral: Starting at 0: A035608, A156859, A002378 = 2*A000217, A137932 = 4*A002620; starting at 1: A317186, A267682, A002061, A080335.

Programs

Formula

a(n) = 1 + Sum_{i=1..n} 8*i = 1 + 8*A000217(n). - Xavier Acloque, Jan 21 2003; Zak Seidov, May 07 2006; Robert G. Wilson v, Dec 29 2010
O.g.f.: (1+6*x+x^2)/(1-x)^3. - R. J. Mathar, Jan 11 2008
a(n) = 4*n*(n + 1) + 1 = 4*n^2 + 4*n + 1. - Artur Jasinski, Mar 27 2008
a(n) = A061038(2+4n). - Paul Curtz, Oct 26 2008
Sum_{n>=0} 1/a(n) = Pi^2/8 = A111003. - Jaume Oliver Lafont, Mar 07 2009
a(n) = A000290(A005408(n)). - Reinhard Zumkeller, Nov 08 2009
a(n) = a(n-1) + 8*n with n>0, a(0)=1. - Vincenzo Librandi, Aug 01 2010
a(n) = A033951(n) + n. - Reinhard Zumkeller, May 17 2009
a(n) = A033996(n) + 1. - Omar E. Pol, Oct 03 2011
a(n) = (A005408(n))^2. - Zak Seidov, Nov 29 2011
From George F. Johnson, Sep 05 2012: (Start)
a(n+1) = a(n) + 4 + 4*sqrt(a(n)).
a(n-1) = a(n) + 4 - 4*sqrt(a(n)).
a(n+1) = 2*a(n) - a(n-1) + 8.
a(n+1) = 3*a(n) - 3*a(n-1) + a(n-2).
(a(n+1) - a(n-1))/8 = sqrt(a(n)).
a(n+1)*a(n-1) = (a(n)-4)^2.
a(n) = 2*A046092(n) + 1 = 2*A001844(n) - 1 = A046092(n) + A001844(n).
Limit_{n -> oo} a(n)/a(n-1) = 1. (End)
a(n) = binomial(2*n+2,2) + binomial(2*n+1,2). - John Molokach, Jul 12 2013
E.g.f.: (1 + 8*x + 4*x^2)*exp(x). - Ilya Gutkovskiy, May 23 2016
a(n) = A101321(8,n). - R. J. Mathar, Jul 28 2016
Product_{n>=1} A033996(n)/a(n) = Pi/4. - Daniel Suteu, Dec 25 2016
a(n) = A014105(n) + A000384(n+1). - Bruce J. Nicholson, Nov 11 2017
a(n) = A003215(n) + A002378(n). - Klaus Purath, Jun 09 2020
From Amiram Eldar, Jun 20 2020: (Start)
Sum_{n>=0} a(n)/n! = 13*e.
Sum_{n>=0} (-1)^(n+1)*a(n)/n! = 3/e. (End)
Sum_{n>=0} (-1)^n/a(n) = A006752. - Amiram Eldar, Oct 10 2020
From Amiram Eldar, Jan 28 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = cosh(Pi/2).
Product_{n>=1} (1 - 1/a(n)) = Pi/4 (A003881). (End)
From Leo Tavares, Nov 24 2021: (Start)
a(n) = A014634(n) - A002943(n). See Diamond Triangles illustration.
a(n) = A003154(n+1) - A046092(n). See Diamond Stars illustration. (End)
From Peter Bala, Mar 11 2024: (Start)
Sum_{k = 1..n+1} 1/(k*a(k)*a(k-1)) = 1/(9 - 3/(17 - 60/(33 - 315/(57 - ... - n^2*(4*n^2 - 1)/((2*n + 1)^2 + 2*2^2 ))))).
3/2 - 2*log(2) = Sum_{k >= 1} 1/(k*a(k)*a(k-1)) = 1/(9 - 3/(17 - 60/(33 - 315/(57 - ... - n^2*(4*n^2 - 1)/((2*n + 1)^2 + 2*2^2 - ... ))))).
Row 2 of A142992. (End)
From Peter Bala, Mar 26 2024: (Start)
8*a(n) = (2*n + 1)*(a(n+1) - a(n-1)).
Sum_{n >= 0} (-1)^n/(a(n)*a(n+1)) = 1/2 - Pi/8 = 1/(9 + (1*3)/(8 + (3*5)/(8 + ... + (4*n^2 - 1)/(8 + ... )))). For the continued fraction use Lorentzen and Waadeland, p. 586, equation 4.7.9 with n = 1. Cf. A057813. (End)

Extensions

Additional description from Terrel Trotter, Jr., Apr 06 2002

A064038 Numerator of average number of swaps needed to bubble sort a string of n distinct letters.

Original entry on oeis.org

0, 1, 3, 3, 5, 15, 21, 14, 18, 45, 55, 33, 39, 91, 105, 60, 68, 153, 171, 95, 105, 231, 253, 138, 150, 325, 351, 189, 203, 435, 465, 248, 264, 561, 595, 315, 333, 703, 741, 390, 410, 861, 903, 473, 495, 1035, 1081, 564, 588, 1225, 1275, 663, 689, 1431, 1485, 770
Offset: 1

Views

Author

Antti Karttunen, Aug 23 2001

Keywords

Comments

Denominators are given by the simple periodic sequence [1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, ...] (= A014695) thus we get an average of 1/2, 3/2, 3, 5, 15/2, 21/2, 14, 18, etc. swappings required to bubble sort a string of 2, 3, 4, 5, 6, ... letters.

References

  • E. Reingold, J. Nievergelt and N. Deo, Combinatorial Algorithms, Prentice-Hall, 1977, section 7.1, p. 287.

Crossrefs

Programs

  • Magma
    [Numerator(n*(n-1)/4): n in [1..100]]; // G. C. Greubel, Sep 21 2018
  • Maple
    [seq(numer((n*(n-1))/4), n=1..120)];
  • Mathematica
    f[n_] := Numerator[n (n - 1)/4]; Array[f, 56]
    f[n_] := n/GCD[n, 4]; Array[f[#] f[# - 1] &, 56]
    LinearRecurrence[{3,-6,10,-12,12,-10,6,-3,1},{0,1,3,3,5,15,21,14,18},80] (* Harvey P. Dale, Jan 23 2023 *)
  • PARI
    vector(100, n, numerator(n*(n-1)/4)) \\ G. C. Greubel, Sep 21 2018
    

Formula

a(n) = numerator(A001809(n)/(n!)).
a(4n) = A033991(n).
a(4n+1) = A007742(n).
a(4n+2) = A014634(n).
a(4n+3) = A033567(n+1).
a(n+1) = A061041(8*n-4). - Paul Curtz, Jan 03 2011
G.f.: -x^2*(1+4*x^3+x^6) / ( (x-1)^3*(1+x^2)^3 ). - R. J. Mathar, Jan 03 2011
a(n+1) = A060819(n)*A060819(n+1).
a(n+1) = A000217(n)/(period 4:repeat 2,1,1,2=A014695(n+2)=A130658(n+3)).
a(n) = 3*a(n-4) -3*a(n-8) +a(n-12). - Paul Curtz, Mar 04 2011
a(n) = +3*a(n-1) -6*a(n-2) +10*a(n-3) -12*a(n-4) +12*a(n-5) -10*a(n-6) +6*a(n-7) -3*a(n-8) +1*a(n-9). - Joerg Arndt, Mar 04 2011
a(n+1) = A026741(A000217(n)). - Paul Curtz, Apr 04 2011
a(n) = numerator(Sum_{k=0..n-1} k/2). - Arkadiusz Wesolowski, Aug 09 2012
a(n) = n*(n-1)*(3-i^(n*(n-1)))/8, where i=sqrt(-1). - Bruno Berselli, Oct 01 2012, corrected by Vaclav Kotesovec, Aug 09 2022
Sum_{n>=2} 1/a(n) = 4 - Pi/2. - Amiram Eldar, Aug 09 2022
E.g.f.: x^2*(3*exp(x) + cos(x) + sin(x))/8. - Stefano Spezia, Aug 23 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

A188358 T(n,k)=Number of arrangements of n+1 numbers x(i) in -k..k with the sum of x(i)*x(i+1) equal to zero.

Original entry on oeis.org

5, 9, 15, 13, 45, 29, 17, 91, 117, 91, 21, 153, 257, 565, 197, 25, 231, 497, 1775, 1953, 627, 29, 325, 749, 4201, 7369, 9153, 1477, 33, 435, 1205, 8051, 22229, 49747, 37269, 4671, 37, 561, 1569, 14477, 45801, 181449, 262925, 173409, 11693, 41, 703, 2161
Offset: 1

Views

Author

R. H. Hardin Mar 28 2011

Keywords

Comments

Table starts
.....5.......9.......13........17.........21..........25..........29
....15......45.......91.......153........231.........325.........435
....29.....117......257.......497........749........1205........1569
....91.....565.....1775......4201.......8051.......14477.......22583
...197....1953.....7369.....22229......45801......100885......163745
...627....9153....49747....181449.....480363.....1171345.....2295051
..1477...37269...262925...1264985....3882729....11535493....25217201
..4671..173409..1763151..10396729...40629247...136061137...355442823
.11693..766869.10536621..81850801..383541193..1540834589..4602176797
.36487.3589357.70819767.684867741.4033546763.18592832393.65556183495

Examples

			Some solutions for n=6 k=4
.-4...-1...-4....0....3....0...-1....0...-4...-1...-1...-1...-2...-1....1...-3
..0...-2...-3...-2...-4...-4...-3...-4....0...-1...-2...-3...-2....0...-2...-1
.-4....2....3....1....3....4...-3...-1....4....0....1....1...-2....4....1....1
..4...-1...-1....1....3....4....0...-2....0....0...-3...-4....1...-4...-4....2
..4....2....2....1....1....0....1....3...-2...-3...-1...-1...-2...-4....3...-3
..3....1...-2....0...-4....1...-3...-2....0...-1....0....2....1...-4....4...-1
.-4....4...-3...-1...-4....0....3...-3...-4....4....3....1...-2....4....2....1
		

Crossrefs

Row 1 is A016813
Row 2 is A014634

A033567 a(n) = (2*n-1)*(4*n-1).

Original entry on oeis.org

1, 3, 21, 55, 105, 171, 253, 351, 465, 595, 741, 903, 1081, 1275, 1485, 1711, 1953, 2211, 2485, 2775, 3081, 3403, 3741, 4095, 4465, 4851, 5253, 5671, 6105, 6555, 7021, 7503, 8001, 8515, 9045, 9591, 10153, 10731, 11325, 11935, 12561, 13203, 13861, 14535, 15225
Offset: 0

Views

Author

Keywords

Comments

a(n+1) = A005563(1), A061037(3), A061039(5), A061041(7), A061043(9), A061045(11), A061047(13), A061049(15). Lyman, Balmer, Paschen, Brackett, Pfund, Humphreys, Hansen-Strong, ... spectra of hydrogen. - Paul Curtz, Oct 08 2008
Sequence found by reading the segment [1, 3] together with the line from 3, in the direction 3, 21, ..., in the square spiral whose vertices are the triangular numbers A000217. - Omar E. Pol, Sep 03 2011

Crossrefs

Programs

  • Magma
    [(2*n-1)*(4*n-1): n in [0..50]]; // G. C. Greubel, Sep 19 2018
  • Mathematica
    Table[(2*n - 1)*(4*n - 1), {n, 0, 50}] (* G. C. Greubel, Jul 06 2017 *)
    LinearRecurrence[{3,-3,1},{1,3,21},50] (* Harvey P. Dale, Aug 25 2019 *)
  • PARI
    vector(60, n, n--; (2*n-1)*(4*n-1)) \\ Michel Marcus, Apr 12 2015
    

Formula

a(n) = a(n-1) + 16*n - 14 (with a(0)=1). - Vincenzo Librandi, Nov 17 2010
From G. C. Greubel, Jul 06 2017: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-2).
E.g.f.: (1 + 2*x + 8*x^2)*exp(x).
G.f.: (1 + 15*x^2)/(1 - x)^3. (End)
From Amiram Eldar, Jan 03 2022: (Start)
Sum_{n>=0} 1/a(n) = 1 + Pi/4 - log(2)/2.
Sum_{n>=0} (-1)^n/a(n) = 1 + (sqrt(2)-1)*Pi/4 + log(sqrt(2)-1)/sqrt(2). (End)

Extensions

More terms from Michel Marcus, Apr 12 2015

A139278 a(n) = n*(8*n+7).

Original entry on oeis.org

0, 15, 46, 93, 156, 235, 330, 441, 568, 711, 870, 1045, 1236, 1443, 1666, 1905, 2160, 2431, 2718, 3021, 3340, 3675, 4026, 4393, 4776, 5175, 5590, 6021, 6468, 6931, 7410, 7905, 8416, 8943, 9486, 10045, 10620, 11211, 11818, 12441, 13080
Offset: 0

Views

Author

Omar E. Pol, Apr 26 2008

Keywords

Comments

Sequence found by reading the segment (0, 15) together with the line from 15, in the direction 15, 46, ..., in the square spiral whose vertices are the triangular numbers A000217.

Crossrefs

Programs

  • Mathematica
    Table[n (8 n + 7), {n, 0, 40}] (* or *) LinearRecurrence[{3, -3, 1}, {0, 15, 46}, 50] (* Harvey P. Dale, Oct 07 2015 *)
  • PARI
    a(n)=n*(8*n+7) \\ Charles R Greathouse IV, Jun 17 2017

Formula

a(n) = 8*n^2 + 7*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)-1 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
From G. C. Greubel, Jul 18 2017: (Start)
G.f.: x*(x+15)/(1-x)^3.
E.g.f.: (8*x^2 + 15*x)*exp(x). (End)
Sum_{n>=1} 1/a(n) = 8/49 + (sqrt(2)+1)*Pi/14 - 4*log(2)/7 - sqrt(2)*log(sqrt(2)+1)/7. - Amiram Eldar, Mar 17 2022

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
Showing 1-10 of 24 results. Next