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 10 results.

A001318 Generalized pentagonal numbers: m*(3*m - 1)/2, m = 0, +-1, +-2, +-3, ....

Original entry on oeis.org

0, 1, 2, 5, 7, 12, 15, 22, 26, 35, 40, 51, 57, 70, 77, 92, 100, 117, 126, 145, 155, 176, 187, 210, 222, 247, 260, 287, 301, 330, 345, 376, 392, 425, 442, 477, 495, 532, 551, 590, 610, 651, 672, 715, 737, 782, 805, 852, 876, 925, 950, 1001, 1027, 1080, 1107, 1162, 1190, 1247, 1276, 1335
Offset: 0

Views

Author

Keywords

Comments

Partial sums of A026741. - Jud McCranie; corrected by Omar E. Pol, Jul 05 2012
From R. K. Guy, Dec 28 2005: (Start)
"Conway's relation twixt the triangular and pentagonal numbers: Divide the triangular numbers by 3 (when you can exactly):
0 1 3 6 10 15 21 28 36 45 55 66 78 91 105 120 136 153 ...
0 - 1 2 .- .5 .7 .- 12 15 .- 22 26 .- .35 .40 .- ..51 ...
.....-.-.....+..+.....-..-.....+..+......-...-.......+....
"and you get the pentagonal numbers in pairs, one of positive rank and the other negative.
"Append signs according as the pair have the same (+) or opposite (-) parity.
"Then Euler's pentagonal number theorem is easy to remember:
"p(n-0) - p(n-1) - p(n-2) + p(n-5) + p(n-7) - p(n-12) - p(n-15) ++-- = 0^n
where p(n) is the partition function, the left side terminates before the argument becomes negative and 0^n = 1 if n = 0 and = 0 if n > 0.
"E.g. p(0) = 1, p(7) = p(7-1) + p(7-2) - p(7-5) - p(7-7) + 0^7 = 11 + 7 - 2 - 1 + 0 = 15."
(End)
The sequence may be used in order to compute sigma(n), as described in Euler's article. - Thomas Baruchel, Nov 19 2003
Number of levels in the partitions of n + 1 with parts in {1,2}.
a(n) is the number of 3 X 3 matrices (symmetrical about each diagonal) M = {{a, b, c}, {b, d, b}, {c, b, a}} such that a + b + c = b + d + b = n + 2, a,b,c,d natural numbers; example: a(3) = 5 because (a,b,c,d) = (2,2,1,1), (1,2,2,1), (1,1,3,3), (3,1,1,3), (2,1,2,3). - Philippe Deléham, Apr 11 2007
Also numbers a(n) such that 24*a(n) + 1 = (6*m - 1)^2 are odd squares: 1, 25, 49, 121, 169, 289, 361, ..., m = 0, +-1, +-2, ... . - Zak Seidov, Mar 08 2008
From Matthew Vandermast, Oct 28 2008: (Start)
Numbers n for which A000326(n) is a member of A000332. Cf. A145920.
This sequence contains all members of A000332 and all nonnegative members of A145919. For values of n such that n*(3*n - 1)/2 belongs to A000332, see A145919. (End)
Starting with offset 1 = row sums of triangle A168258. - Gary W. Adamson, Nov 21 2009
Starting with offset 1 = Triangle A101688 * [1, 2, 3, ...]. - Gary W. Adamson, Nov 27 2009
Starting with offset 1 can be considered the first in an infinite set generated from A026741. Refer to the array in A175005. - Gary W. Adamson, Apr 03 2010
Vertex number of a square spiral whose edges have length A026741. The two axes of the spiral forming an "X" are A000326 and A005449. The four semi-axes forming an "X" are A049452, A049453, A033570 and the numbers >= 2 of A033568. - Omar E. Pol, Sep 08 2011
A general formula for the generalized k-gonal numbers is given by n*((k - 2)*n - k + 4)/2, n=0, +-1, +-2, ..., k >= 5. - Omar E. Pol, Sep 15 2011
a(n) is the number of 3-tuples (w,x,y) having all terms in {0,...,n} and 2*w = 2*x + y. - Clark Kimberling, Jun 04 2012
Generalized k-gonal numbers are second k-gonal numbers and positive terms of k-gonal numbers interleaved, k >= 5. - Omar E. Pol, Aug 04 2012
a(n) is the sum of the largest parts of the partitions of n+1 into exactly 2 parts. - Wesley Ivan Hurt, Jan 26 2013
Conway's relation mentioned by R. K. Guy is a relation between triangular numbers and generalized pentagonal numbers, two sequences from different families, but as triangular numbers are also generalized hexagonal numbers in this case we have a relation between two sequences from the same family. - Omar E. Pol, Feb 01 2013
Start with the sequence of all 0's. Add n to each value of a(n) and the next n - 1 terms. The result is the generalized pentagonal numbers. - Wesley Ivan Hurt, Nov 03 2014
(6k + 1) | a(4k). (3k + 1) | a(4k+1). (3k + 2) | a(4k+2). (6k + 5) | a(4k+3). - Jon Perry, Nov 04 2014
Enge, Hart and Johansson proved: "Every generalised pentagonal number c >= 5 is the sum of a smaller one and twice a smaller one, that is, there are generalised pentagonal numbers a, b < c such that c = 2a + b." (see link theorem 5). - Peter Luschny, Aug 26 2016
The Enge, et al. result for c >= 5 also holds for c >= 2 if 0 is included as a generalized pentagonal number. That is, 2 = 2*1 + 0. - Michael Somos, Jun 02 2018
Suggestion for title, where n actually matches the list and b-file: "Generalized pentagonal numbers: k(n)*(3*k(n) - 1)/2, where k(n) = A001057(n) = [0, 1, -1, 2, -2, 3, -3, ...], n >= 0" - Daniel Forgues, Jun 09 2018 & Jun 12 2018
Generalized k-gonal numbers are the partial sums of the sequence formed by the multiples of (k - 4) and the odd numbers (A005408) interleaved, with k >= 5. - Omar E. Pol, Jul 25 2018
The last digits form a symmetric cycle of length 40 [0, 1, 2, 5, ..., 5, 2, 1, 0], i.e., a(n) == a(n + 40) (mod 10) and a(n) == a(40*k - n - 1) (mod 10), 40*k > n. - Alejandro J. Becerra Jr., Aug 14 2018
Only 2, 5, and 7 are prime. All terms are of the form k*(k+1)/6, where 3 | k or 3 | k+1. For k > 6, the value divisible by 3 must have another factor d > 2, which will remain after the division by 6. - Eric Snyder, Jun 03 2022
8*a(n) is the product of two even numbers one of which is n + n mod 2. - Peter Luschny, Jul 15 2022
a(n) is the dot product of [1, 2, 3, ..., n] and repeat[1, 1/2]. a(5) = 12 = [1, 2, 3, 4, 5] dot [1, 1/2, 1, 1/2, 1] = [1 + 1 + 3 + 2 + 5]. - Gary W. Adamson, Dec 10 2022
Every nonnegative number is the sum of four terms of this sequence [S. Realis]. - N. J. A. Sloane, May 07 2023
From Peter Bala, Jan 06 2025: (Start)
The sequence terms are the exponents in the expansions of the following infinite products:
1) Product_{n >= 1} (1 - s(n)*q^n) = 1 + q + q^2 + q^5 + q^7 + q^12 + q^15 + ..., where s(n) = (-1)^(1 + mod(n+1,3)).
2) Product_{n >= 1} (1 - q^(2*n))*(1 - q^(3*n))^2/((1 - q^n)*(1 - q^(6*n))) = 1 + q + q^2 + q^5 + q^7 + q^12 + q^15 + ....
3) Product_{n >= 1} (1 - q^n)*(1 - q^(4*n))*(1 - q^(6*n))^5/((1 - q^(2*n))*(1 - q^(3*n))*(1 - q^(12*n)))^2 = 1 - q + q^2 - q^5 - q^7 + q^12 - q^15 + q^22 + q^26 - q^35 + ....
4) Product_{n >= 1} (1 - q^(2*n))^13/((1 - (-1)^n*q^n)*(1 - q^(4*n)))^5 = 1 - 5*q + 7*q^2 - 11*q^5 + 13*q^7 - 17*q^12 + 19*q^15 - + .... See Oliver, Theorem 1.1. (End)

Examples

			G.f. = x + 2*x^2 + 5*x^3 + 7*x^4 + 12*x^5 + 15*x^6 + 22*x^7 + 26*x^8 + 35*x^9 + ...
		

References

  • Enoch Haga, A strange sequence and a brilliant discovery, chapter 5 of Exploring prime numbers on your PC and the Internet, first revised ed., 2007 (and earlier ed.), pp. 53-70.
  • Ross Honsberger, Ingenuity in Mathematics, Random House, 1970, p. 117.
  • Donald E. Knuth, The Art of Computer Programming, vol. 4A, Combinatorial Algorithms, (to appear), section 7.2.1.4, equation (18).
  • Ivan Niven and Herbert S. Zuckerman, An Introduction to the Theory of Numbers, 2nd ed., Wiley, NY, 1966, p. 231.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A080995 (characteristic function), A026741 (first differences), A034828 (partial sums), A165211 (mod 2).
Cf. A000326 (pentagonal numbers), A005449 (second pentagonal numbers), A000217 (triangular numbers).
Indices of nonzero terms of A010815, i.e., the (zero-based) indices of 1-bits of the infinite binary word to which the terms of A068052 converge.
Union of A036498 and A036499.
Sequences of generalized k-gonal numbers: this sequence (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), A274979 (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).
Column 1 of A195152.
Squares in APs: A221671, A221672.
Quadrisection: A049453(k), A033570(k), A033568(k+1), A049452(k+1), k >= 0.
Cf. A002620.

Programs

  • GAP
    a:=[0,1,2,5];; for n in [5..60] do a[n]:=2*a[n-2]-a[n-4]+3; od; a; # Muniru A Asiru, Aug 16 2018
    
  • Haskell
    a001318 n = a001318_list !! n
    a001318_list = scanl1 (+) a026741_list -- Reinhard Zumkeller, Nov 15 2015
    
  • Magma
    [(6*n^2 + 6*n + 1 - (2*n + 1)*(-1)^n)/16 : n in [0..50]]; // Wesley Ivan Hurt, Nov 03 2014
    
  • Magma
    [(3*n^2 + 2*n + (n mod 2) * (2*n + 1)) div 8: n in [0..70]]; // Vincenzo Librandi, Nov 04 2014
    
  • Maple
    A001318 := -(1+z+z**2)/(z+1)**2/(z-1)**3; # Simon Plouffe in his 1992 dissertation; gives sequence without initial zero
    A001318 := proc(n) (6*n^2+6*n+1)/16-(2*n+1)*(-1)^n/16 ; end proc: # R. J. Mathar, Mar 27 2011
  • Mathematica
    Table[n*(n+1)/6, {n, Select[Range[0, 100], Mod[#, 3] != 1 &]}]
    Select[Accumulate[Range[0,200]]/3,IntegerQ] (* Harvey P. Dale, Oct 12 2014 *)
    CoefficientList[Series[x (1 + x + x^2) / ((1 + x)^2 (1 - x)^3), {x, 0, 70}], x] (* Vincenzo Librandi, Nov 04 2014 *)
    LinearRecurrence[{1,2,-2,-1,1},{0,1,2,5,7},70] (* Harvey P. Dale, Jun 05 2017 *)
    a[ n_] := With[{m = Quotient[n + 1, 2]}, m (3 m + (-1)^n) / 2]; (* Michael Somos, Jun 02 2018 *)
  • PARI
    {a(n) = (3*n^2 + 2*n + (n%2) * (2*n + 1)) / 8}; /* Michael Somos, Mar 24 2011 */
    
  • PARI
    {a(n) = if( n<0, n = -1-n); polcoeff( x * (1 - x^3) / ((1 - x) * (1-x^2))^2 + x * O(x^n), n)}; /* Michael Somos, Mar 24 2011 */
    
  • PARI
    {a(n) = my(m = (n+1) \ 2); m * (3*m + (-1)^n) / 2}; /* Michael Somos, Jun 02 2018 */
    
  • Python
    def a(n):
        p = n % 2
        return (n + p)*(3*n + 2 - p) >> 3
    print([a(n) for n in range(60)])  # Peter Luschny, Jul 15 2022
    
  • Python
    def A001318(n): return n*(n+1)-(m:=n>>1)*(m+1)>>1 # Chai Wah Wu, Nov 23 2024
  • Sage
    @CachedFunction
    def A001318(n):
        if n == 0 : return 0
        inc = n//2 if is_even(n) else n
        return inc + A001318(n-1)
    [A001318(n) for n in (0..59)] # Peter Luschny, Oct 13 2012
    

Formula

Euler: Product_{n>=1} (1 - x^n) = Sum_{n=-oo..oo} (-1)^n*x^(n*(3*n - 1)/2).
A080995(a(n)) = 1: complement of A090864; A000009(a(n)) = A051044(n). - Reinhard Zumkeller, Apr 22 2006
Euler transform of length-3 sequence [2, 2, -1]. - Michael Somos, Mar 24 2011
a(-1 - n) = a(n) for all n in Z. a(2*n) = A005449(n). a(2*n - 1) = A000326(n). - Michael Somos, Mar 24 2011. [The extension of the recurrence to negative indices satisfies the signature (1,2,-2,-1,1), but not the definition of the sequence m*(3*m -1)/2, because there is no m such that a(-1) = 0. - Klaus Purath, Jul 07 2021]
a(n) = 3 + 2*a(n-2) - a(n-4). - Ant King, Aug 23 2011
Product_{k>0} (1 - x^k) = Sum_{k>=0} (-1)^k * x^a(k). - Michael Somos, Mar 24 2011
G.f.: x*(1 + x + x^2)/((1 + x)^2*(1 - x)^3).
a(n) = n*(n + 1)/6 when n runs through numbers == 0 or 2 mod 3. - Barry E. Williams
a(n) = A008805(n-1) + A008805(n-2) + A008805(n-3), n > 2. - Ralf Stephan, Apr 26 2003
Sequence consists of the pentagonal numbers (A000326), followed by A000326(n) + n and then the next pentagonal number. - Jon Perry, Sep 11 2003
a(n) = (6*n^2 + 6*n + 1)/16 - (2*n + 1)*(-1)^n/16; a(n) = A034828(n+1) - A034828(n). - Paul Barry, May 13 2005
a(n) = Sum_{k=1..floor((n+1)/2)} (n - k + 1). - Paul Barry, Sep 07 2005
a(n) = A000217(n) - A000217(floor(n/2)). - Pierre CAMI, Dec 09 2007
If n even a(n) = a(n-1) + n/2 and if n odd a(n) = a(n-1) + n, n >= 2. - Pierre CAMI, Dec 09 2007
a(n)-a(n-1) = A026741(n) and it follows that the difference between consecutive terms is equal to n if n is odd and to n/2 if n is even. Hence this is a self-generating sequence that can be simply constructed from knowledge of the first term alone. - Ant King, Sep 26 2011
a(n) = (1/2)*ceiling(n/2)*ceiling((3*n + 1)/2). - Mircea Merca, Jul 13 2012
a(n) = (A008794(n+1) + A000217(n))/2 = A002378(n) - A085787(n). - Omar E. Pol, Jan 12 2013
a(n) = floor((n + 1)/2)*((n + 1) - (1/2)*floor((n + 1)/2) - 1/2). - Wesley Ivan Hurt, Jan 26 2013
From Oskar Wieland, Apr 10 2013: (Start)
a(n) = a(n+1) - A026741(n),
a(n) = a(n+2) - A001651(n),
a(n) = a(n+3) - A184418(n),
a(n) = a(n+4) - A007310(n),
a(n) = a(n+6) - A001651(n)*3 = a(n+6) - A016051(n),
a(n) = a(n+8) - A007310(n)*2 = a(n+8) - A091999(n),
a(n) = a(n+10)- A001651(n)*5 = a(n+10)- A072703(n),
a(n) = a(n+12)- A007310(n)*3,
a(n) = a(n+14)- A001651(n)*7. (End)
a(n) = (A007310(n+1)^2 - 1)/24. - Richard R. Forberg, May 27 2013; corrected by Zak Seidov, Mar 14 2015; further corrected by Jianing Song, Oct 24 2018
a(n) = Sum_{i = ceiling((n+1)/2)..n} i. - Wesley Ivan Hurt, Jun 08 2013
G.f.: x*G(0), where G(k) = 1 + x*(3*k + 4)/(3*k + 2 - x*(3*k + 2)*(3*k^2 + 11*k + 10)/(x*(3*k^2 + 11*k + 10) + (k + 1)*(3*k + 4)/G(k+1))); (continued fraction). - Sergei N. Gladkovskii, Jun 16 2013
Sum_{n>=1} 1/a(n) = 6 - 2*Pi/sqrt(3). - Vaclav Kotesovec, Oct 05 2016
a(n) = Sum_{i=1..n} numerator(i/2) = Sum_{i=1..n} denominator(2/i). - Wesley Ivan Hurt, Feb 26 2017
a(n) = A000292(A001651(n))/A001651(n), for n>0. - Ivan N. Ianakiev, May 08 2018
a(n) = ((-5 + (-1)^n - 6n)*(-1 + (-1)^n - 6n))/96. - José de Jesús Camacho Medina, Jun 12 2018
a(n) = Sum_{k=1..n} k/gcd(k,2). - Pedro Caceres, Apr 23 2019
Quadrisection. For r = 0,1,2,3: a(r + 4*k) = 6*k^2 + sqrt(24*a(r) + 1)*k + a(r), for k >= 1, with inputs (k = 0) {0,1,2,5}. These are the sequences A049453(k), A033570(k), A033568(k+1), A049452(k+1), for k >= 0, respectively. - Wolfdieter Lang, Feb 12 2021
a(n) = a(n-4) + sqrt(24*a(n-2) + 1), n >= 4. - Klaus Purath, Jul 07 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = 6*(log(3)-1). - Amiram Eldar, Feb 28 2022
a(n) = A002620(n) + A008805(n-1). Gary W. Adamson, Dec 10 2022
E.g.f.: (x*(7 + 3*x)*cosh(x) + (1 + 5*x + 3*x^2)*sinh(x))/8. - Stefano Spezia, Aug 01 2024

A049453 Second pentagonal numbers with even index: a(n) = n*(6*n+1).

Original entry on oeis.org

0, 7, 26, 57, 100, 155, 222, 301, 392, 495, 610, 737, 876, 1027, 1190, 1365, 1552, 1751, 1962, 2185, 2420, 2667, 2926, 3197, 3480, 3775, 4082, 4401, 4732, 5075, 5430, 5797, 6176, 6567, 6970, 7385, 7812, 8251, 8702, 9165, 9640, 10127, 10626, 11137, 11660, 12195
Offset: 0

Views

Author

Joe Keane (jgk(AT)jgk.org)

Keywords

Comments

Number of edges in the join of the complete tripartite graph of order 3n and the cycle graph of order n, K_n,n,n * C_n. - Roberto E. Martinez II, Jan 07 2002
Sequence found by reading the line (one of the diagonal axes) from 0, in the direction 0, 7, ..., in the square spiral whose vertices are the generalized pentagonal numbers A001318. - Omar E. Pol, Sep 08 2011
First bisection of A036498. - Bruno Berselli, Nov 25 2012

Crossrefs

Programs

Formula

G.f.: x*(7+5*x)/(1-x)^3.
a(n) = 12*n + a(n-1) - 5 with n > 0, a(0)=0. - Vincenzo Librandi, Aug 06 2010
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - G. C. Greubel, Jun 07 2017
From Amiram Eldar, Feb 18 2022: (Start)
Sum_{n>=1} 1/a(n) = 6 - sqrt(3)*Pi/2 - 2*log(2) - 3*log(3)/2.
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi + log(2) + sqrt(3)*log(2 + sqrt(3)) - 6. (End)
E.g.f.: exp(x)*x*(7 + 6*x). - Elmo R. Oliveira, Dec 12 2024

A331022 Numbers k such that the number of strict integer partitions of k is a power of 2.

Original entry on oeis.org

0, 1, 2, 3, 4, 6, 9, 16, 20, 29, 34, 45
Offset: 0

Views

Author

Gus Wiseman, Jan 10 2020

Keywords

Comments

An integer partition of n is a finite, nonincreasing sequence of positive integers (parts) summing to n. It is strict if the parts are all different. Integer partitions and strict integer partitions are counted by A000041 and A000009 respectively.
Conjecture: This sequence is finite.
Conjecture: The analogous sequence for non-strict partitions is: 0, 1, 2.
Next term > 5*10^4 if it exists. - Seiichi Manyama, Jan 12 2020

Examples

			The strict integer partitions of the initial terms:
  (1)  (2)  (3)    (4)    (6)      (9)
            (2,1)  (3,1)  (4,2)    (5,4)
                          (5,1)    (6,3)
                          (3,2,1)  (7,2)
                                   (8,1)
                                   (4,3,2)
                                   (5,3,1)
                                   (6,2,1)
		

Crossrefs

The version for primes instead of powers of 2 is A035359.
The version for factorizations instead of strict partitions is A330977.
Numbers whose number of partitions is prime are A046063.

Programs

  • Mathematica
    Select[Range[0,1000],IntegerQ[Log[2,PartitionsQ[#]]]&]

A036499 Numbers of the form k*(k+1)/6 for k = 2 or 3 modulo 6.

Original entry on oeis.org

1, 2, 12, 15, 35, 40, 70, 77, 117, 126, 176, 187, 247, 260, 330, 345, 425, 442, 532, 551, 651, 672, 782, 805, 925, 950, 1080, 1107, 1247, 1276, 1426, 1457, 1617, 1650, 1820, 1855, 2035, 2072, 2262, 2301, 2501, 2542, 2752, 2795, 3015, 3060, 3290, 3337, 3577, 3626
Offset: 1

Views

Author

Keywords

Comments

Numbers with an odd number of partitions with an extra odd partition; coefficient of z^p in Product_{n >= 1}(1-z^n) has coefficient (-1).
n such that the number of partitions of n into distinct parts with an odd number of parts exceed by 1 the number of partitions of n into distinct parts with an even number of parts. [Euler's 1754/55 pentagonal number theorem, see, e.g., the Freitag-Busam reference (in German). This reference is from Wolfdieter Lang, Jan 18 2016]
In formal power series, A010815=(product(1-x^k),k>0), ranks of coefficients -1. (A001318=ranks of nonzero (1 or -1) in A010815=ranks of odds terms in A000009).
Quasipolynomial of order 2. - Charles R Greathouse IV, Dec 08 2011
Union of A033568 and A033570. - Ray Chandler, Dec 09 2011

References

  • Eberhard Freitag and Rolf Busam, Funktionentheorie 1, Springer, Vierte Auflage, 2006, p. 410.

Crossrefs

Programs

  • Magma
    [(3*n*n-5*n+2)/2+(2*n-1)*(n mod 2): n in [1..50]]; // Vincenzo Librandi, Jan 19 2016
  • Maple
    seq(seq((6*k+i)*(6*k+i+1)/6,i=2..3),k=0..50); # Robert Israel, Jan 18 2016
  • Mathematica
    Table[ 1/8*(3 + (-1)^k - 6*k)*(1 + (-1)^k - 2*k), {k, 64} ]
    LinearRecurrence[{1,2,-2,-1,1},{1,2,12,15,35},50] (* or *)
    CoefficientList[Series[(1+x+8x^2+x^3+x^4)/((1-x)^3(1+x)^2),{x,0,100}],x] (* or *)
    Table[(2n+1)(3n+{1,2}),{n,0,24}]//Flatten (* Ray Chandler, Dec 09 2011 *)
  • PARI
    a(n)=n*(3*n-5)/2+1+n%2*(2*n-1) \\ Charles R Greathouse IV, Dec 08 2011
    

Formula

a(n) = (3*n*n-5*n+2)/2 + (2*n-1)*(n mod 2). - Frank Ellermann, Mar 16 2002
G.f.: (1+x+8*x^2+x^3+x^4)/((1-x)^3*(1+x)^2). - Ray Chandler, Dec 09 2011
Bisection: a(2*k+1) = A001318(1+4*k) = (2*k+1)*(3*k+1) = A033570(k), a(2*(k+1)) = A001318(2+4*k) = (2*k+1)*(3*k+2) = A033568(k+1), k >= 0. - Wolfdieter Lang, Jan 18 2016
a(n) = a(n-1)+2*a(n-2)-2*a(n-3)-a(n-4)+a(n-5) for n>5. - Wesley Ivan Hurt, Jan 18 2016
From Amiram Eldar, Feb 22 2022: (Start)
Sum_{n>=1} 1/a(n) = Pi/sqrt(3).
Sum_{n>=1} (-1)^(n+1)/a(n) = 3*log(3) - 4*log(2). (End)

Extensions

Better description from Claude Lenormand (claude.lenormand(AT)free.fr), Feb 12 2001
Edited by Ray Chandler, Dec 09 2011

A219191 Numbers of the form k*(7*k+1), where k = 0,-1,1,-2,2,-3,3,...

Original entry on oeis.org

0, 6, 8, 26, 30, 60, 66, 108, 116, 170, 180, 246, 258, 336, 350, 440, 456, 558, 576, 690, 710, 836, 858, 996, 1020, 1170, 1196, 1358, 1386, 1560, 1590, 1776, 1808, 2006, 2040, 2250, 2286, 2508, 2546, 2780, 2820, 3066, 3108, 3366, 3410, 3680, 3726, 4008
Offset: 1

Views

Author

Bruno Berselli, Nov 14 2012

Keywords

Comments

Equivalently, numbers m such that 28*m+1 is a square.
Also, integer values of h*(h+1)/7.
Let F(r) = Product_{n >= 1} 1 - q^(14*n-r). The sequence terms are the exponents in the expansion of F(0)*F(6)*F(8) = 1 - q^6 - q^8 + q^26 + q^30 - q^60 - q^66 + + - - ... (by the triple product identity).- Peter Bala, Dec 25 2024

Crossrefs

Cf. numbers of the form k*(i*k+1) with k in A001057: i=0, A001057; i=1, A110660; i=2, A000217; i=3, A152749; i=4, A074378; i=5, A219190; i=6, A036498; i=7, this sequence; i=8, A154260.
Cf. A113801 (square roots of 28*a(n)+1, see the comment).
Cf. similar sequences listed in A219257.
Subsequence of A011860.

Programs

  • Magma
    k:=7; f:=func; [0] cat [f(n*m): m in [-1,1], n in [1..25]];
    
  • Magma
    I:=[0,6,8,26,30]; [n le 5 select I[n] else Self(n-1)+2*Self(n-2)-2*Self(n-3)-Self(n-4)+Self(n-5): n in [1..50]]; // Vincenzo Librandi, Aug 18 2013
  • Maple
    A := proc (q) local n; for n from 0 to q do if type(sqrt(28*n+1), integer) then print(n) fi; od; end: A(4100); # Peter Bala, Dec 25 2024
  • Mathematica
    Rest[Flatten[{# (7 # - 1), # (7 # + 1)} & /@ Range[0, 25]]]
    CoefficientList[Series[2 x (3 + x + 3 x^2) / ((1 + x)^2 (1 - x)^3), {x, 0, 50}], x] (* Vincenzo Librandi, Aug 18 2013 *)
    LinearRecurrence[{1,2,-2,-1,1},{0,6,8,26,30},50] (* Harvey P. Dale, Sep 14 2022 *)

Formula

G.f.: 2*x^2*(3+x+3*x^2)/((1+x)^2*(1-x)^3).
a(n) = a(-n+1) = (14*n*(n-1)+5*(-1)^n*(2*n-1)+5)/8.
a(n) = 2*A057570(n) = (1/7)*A047335(n)*A047274(n+1).
Sum_{n>=2} 1/a(n) = 7 - cot(Pi/7)*Pi. - Amiram Eldar, Mar 17 2022

A219190 Numbers of the form k*(5*k+1), where k = 0,-1,1,-2,2,-3,3,...

Original entry on oeis.org

0, 4, 6, 18, 22, 42, 48, 76, 84, 120, 130, 174, 186, 238, 252, 312, 328, 396, 414, 490, 510, 594, 616, 708, 732, 832, 858, 966, 994, 1110, 1140, 1264, 1296, 1428, 1462, 1602, 1638, 1786, 1824, 1980, 2020, 2184, 2226, 2398, 2442, 2622, 2668, 2856, 2904, 3100
Offset: 1

Views

Author

Bruno Berselli, Nov 14 2012

Keywords

Comments

Equivalently, numbers m such that 20*m+1 is a square.
Also, integer values of h*(h+1)/5.
More generally, for the numbers of the form n*(k*n+1) with n in A001057, we have:
. generating function (offset 1): x^2*(k-1+2*x+(k-1)*x^2)/((1+x)^2*(1-x)^3);
. n-th term: b(n) = (2*k*n*(n-1)+(k-2)*(-1)^n*(2*n-1)+k-2)/8;
. first differences: (n-1)*((-1)^n*(k-2)+k)/2;
. b(2n+1)-b(2n) = 2*n (independent from k);
. (4*k)*b(n)+1 = (2*k*n+(k-2)*(-1)^n-k)^2/4.

Crossrefs

Subsequence of A011858.
Cf. A090771: square roots of 20*a(n)+1 (see the first comment).
Cf. numbers of the form n*(k*n+1) with n in A001057: k=0, A001057; k=1, A110660; k=2, A000217; k=3, A152749; k=4, A074378; k=5, this sequence; k=6, A036498; k=7, A219191; k=8, A154260.
Cf. similar sequences listed in A219257.

Programs

  • Magma
    k:=5; f:=func; [0] cat [f(n*m): m in [-1,1], n in [1..25]];
    
  • Magma
    I:=[0,4,6,18,22]; [n le 5 select I[n] else Self(n-1)+2*Self(n-2)-2*Self(n-3)-Self(n-4)+Self(n-5): n in [1..50]]; // Vincenzo Librandi, Aug 18 2013
  • Mathematica
    Rest[Flatten[{# (5 # - 1), # (5 # + 1)} & /@ Range[0, 25]]]
    CoefficientList[Series[2 x (2 + x + 2 x^2) / ((1 + x)^2 (1 - x)^3), {x, 0, 50}], x] (* Vincenzo Librandi, Aug 18 2013 *)
    LinearRecurrence[{1,2,-2,-1,1},{0,4,6,18,22},50] (* Harvey P. Dale, Jan 21 2015 *)

Formula

G.f.: 2*x^2*(2 + x + 2*x^2)/((1 + x)^2*(1 - x)^3).
a(n) = a(-n+1) = (10*n*(n-1) + 3*(-1)^n*(2*n - 1) + 3)/8.
a(n) = 2*A057569(n) = A008851(n+1)*A047208(n)/5.
a(n) = a(n-1) + 2*a(n-2) - 2*a(n-3) - a(n-4) + a(n-5). - Harvey P. Dale, Jan 21 2015
Sum_{n>=2} 1/a(n) = 5 - sqrt(1+2/sqrt(5))*Pi. - Amiram Eldar, Mar 15 2022
a(n) = A132356(n-1)/2, n >= 1. - Bernard Schott, Mar 15 2022

A129356 G.f.: A(x) = Product_{n>=1} [ (1-x)^3*(1 + 3x + 6x^2 +...+ n(n+1)/2*x^(n-1)) ].

Original entry on oeis.org

1, -3, -3, 15, -15, 66, -261, 618, -1155, 1040, 2361, -11616, 23733, -27027, 29394, -132318, 545790, -1383459, 2418896, -3383679, 4278462, -3127320, -8332866, 42021990, -99069516, 160683318, -200247795, 214883010, -345461022, 1184850729, -3966311448, 9899287254, -18787986009
Offset: 0

Views

Author

Paul D. Hanna, Apr 10 2007

Keywords

Comments

a(k) != 0 (mod 3) at k = 9*A001318(n) for n>=0, where A001318 are the generalized pentagonal numbers: m(3m-1)/2, m=0,+-1,+-2,...; a(k) == 1 (mod 3) at k = 9*A036498(n) (n>=0); a(k) == -1 (mod 3) at k = 9*A036499(n) (n>=0).

Examples

			A(x) = (1-3x+3x^2-x^3)(1-6x^2+8x^3-3x^4)(1-10x^3+15x^4-6x^5)*...
*( 1 - (n+1)(n+2)/2*x^n + n(n+2)*x^(n+1) - n(n+1)/2*x^(n+2) )*...
Terms are divisible by 3 except at positions given by:
a(n) == 1 (mod 3) at n = [0, 45, 63, 198, 234, 459,...,9*A036498(k),..];
a(n) == -1 (mod 3) at n = [9, 18, 108, 135, 315, 360,..,9*A036499(k),..].
		

Crossrefs

Programs

  • PARI
    {a(n)=if(n==0,1,polcoeff(prod(k=1,n,(1-x)^3*sum(j=1,k,j*(j+1)/2*x^(j-1)) +x*O(x^n)),n))}

Formula

G.f.: A(x) = Product_{n>=1} [ 1 - (n+1)(n+2)/2*x^n + n(n+2)*x^(n+1) - n(n+1)/2*x^(n+2) ].

A129358 G.f.: A(x) = Product_{n>=1} [ (1-x)^5*(1 + 5x + 15x^2 +...+ n(n+1)(n+2)(n+3)/4!*x^(n-1)) ].

Original entry on oeis.org

1, -5, -5, 70, -180, 770, -4760, 20840, -68085, 147890, -795, -1679855, 8378195, -25065005, 56439545, -145200415, 612604910, -2764023765, 10020060660, -28723695265, 67618167310, -128945409045, 137921330680, 375948665405, -3167538981120, 12823443150644, -38103903888575
Offset: 0

Views

Author

Paul D. Hanna, Apr 11 2007

Keywords

Comments

a(k) != 0 (mod 5) at k = 25*A001318(n) for n>=0, where A001318 are the generalized pentagonal numbers: m(3m-1)/2, m=0,+-1,+-2,...; a(k) == 1 (mod 5) at k = 25*A036498(n) (n>=0); a(k) == -1 (mod 5) at k = 25*A036499(n) (n>=0).

Examples

			A(x) = (1-5x+10x^2-10x^3+5x^4-x^5)*(1-15x^2+40x^3-45x^4+24x^5-5x^6)*(1-35x^3+105x^4-126x^5+70x^6-15x^7)*(1-70x^4+224x^5-280x^6+160x^7-35x^8)*...
Terms are divisible by 5 except at positions given by 25*A001318(n):
a(n) == 1 (mod 5) at n = [0, 125, 175, 550, 650,...,25*A036498(k),...];
a(n) == -1 (mod 5) at n = [25, 50, 300, 375, 875,...,25*A036499(k),...].
		

Crossrefs

Programs

  • PARI
    {a(n)=if(n==0,1,polcoeff(prod(k=1,n,(1-x)^5*sum(j=1,k,binomial(j+3,4)*x^(j-1)) +x*O(x^n)),n))}

Formula

G.f.: A(x) = Product_{n>=1} [ 1 - (n+1)(n+2)(n+3)(n+4)/4!*x^n + 4n(n+2)(n+3)(n+4)/4!*x^(n+1) - 6n(n+1)(n+3)(n+4)/4!*x^(n+2) + 4n(n+1)(n+2)(n+4)/4!*x^(n+3) - n(n+1)(n+2)(n+3)/4!*x^(n+4) ].

A186741 Expansion of f(x^5, x^7) in powers of x where f(, ) is Ramanujan's general theta function.

Original entry on oeis.org

1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0
Offset: 0

Views

Author

Michael Somos, Jan 21 2012

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

Examples

			G.f. = 1 + x^5 + x^7 + x^22 + x^26 + x^51 + x^57 + x^92 + x^100 + x^145 + ...
G.f. = q + q^121 + q^169 + q^529 + q^625 + q^1225 + q^1369 + q^2209 + q^2401 + ...
		

Crossrefs

Programs

  • Mathematica
    a[n_]:= SeriesCoefficient[ QPochhammer[-q^5,q^12]*QPochhammer[-q^7,q^12] *QPochhammer[q^12,q^12], {q, 0, n}]; (* G. C. Greubel, Dec 08 2017 *)
  • PARI
    {a(n) = my(m); if( !issquare( 24*n + 1, &m), 0, m%12 == 1 || m%12 == 11)};

Formula

Euler transform of period 24 sequence [ 0, 0, 0, 0, 1, 0, 1, 0, 0, -1, 0, -1, 0, -1, 0, 0, 1, 0, 1, 0, 0, 0, 0, -1, ...].
a(n) is the characteristic function of A036498. a(n) = max( 0, A010815(n)).
G.f.: Sum_{k in Z} x^(6*k^2 - k) = Product_{k>0} (1 + x^(12*k - 7)) * (1 + x^(12*k - 5)) * (1 - x^(12*k)).
Sum_{k=1..n} a(k) ~ sqrt(2*n/3). - Amiram Eldar, Jan 13 2024

A308400 Expansion of 1 / Sum_{k=-oo..oo} (-x)^(k*(6*k + 1)).

Original entry on oeis.org

1, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 0, 2, 0, 1, 1, 0, 3, 0, 3, 1, 1, 3, 0, 6, 1, 3, 3, 1, 8, 1, 8, 3, 3, 9, 2, 14, 3, 9, 9, 4, 19, 4, 19, 9, 10, 21, 6, 32, 10, 22, 22, 12, 42, 12, 43, 23, 25, 48, 18, 67, 25, 51, 51, 31, 88, 31, 90, 54, 59, 101, 44, 137, 60, 108, 109, 73, 177, 73
Offset: 0

Views

Author

Ilya Gutkovskiy, May 24 2019

Keywords

Comments

Number of partitions of n into parts congruent to {0, 5, 7} mod 12.
Convolution inverse of A247223.

Crossrefs

Programs

  • Mathematica
    nmax = 78; CoefficientList[Series[1/Sum[(-x)^(k (6 k + 1)), {k, -nmax, nmax}], {x, 0, nmax}], x]
    nmax = 78; CoefficientList[Series[Product[1/((1 - x^(12 k - 7)) * (1 - x^(12 k - 5)) * (1 - x^(12 k))), {k, 1, nmax}], {x, 0, nmax}], x]

Formula

G.f.: 1 / Sum_{k>=1} (-x)^A036498(k).
G.f.: Product_{k>=1} 1 / ((1 - x^(12*k - 7)) * (1 - x^(12*k - 5)) * (1 - x^(12*k))).
a(n) ~ (sqrt(3) - 1) * exp(sqrt(n/6)*Pi) / (2^(5/2)*n). - Vaclav Kotesovec, May 25 2019
Showing 1-10 of 10 results.