A002129 -id:A002129 - cp's OEIS frontend

A008441 Number of ways of writing n as the sum of 2 triangular numbers.

1, 2, 1, 2, 2, 0, 3, 2, 0, 2, 2, 2, 1, 2, 0, 2, 4, 0, 2, 0, 1, 4, 2, 0, 2, 2, 0, 2, 2, 2, 1, 4, 0, 0, 2, 0, 4, 2, 2, 2, 0, 0, 3, 2, 0, 2, 4, 0, 2, 2, 0, 4, 0, 0, 0, 4, 3, 2, 2, 0, 2, 2, 0, 0, 2, 2, 4, 2, 0, 2, 2, 0, 3, 2, 0, 0, 4, 0, 2, 2, 0, 6, 0, 2, 2, 0, 0, 2, 2, 0, 1, 4, 2, 2, 4, 0, 0, 2, 0, 2, 2, 2, 2, 0, 0

Offset: 0

N. J. A. Sloane

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700). The present sequence gives the expansion coefficients of psi(q)^2.

Also the number of positive odd solutions to equation x^2 + y^2 = 8*n + 2. - Seiichi Manyama, May 28 2017

			G.f. = 1 + 2*x + x^2 + 2*x^3 + 2*x^4 + 3*x^6 + 2*x^7 + 2*x^9 + 2*x^10 + 2*x^11 + ...
G.f. for B(q) = q * A(q^4) = q + 2*q^5 + q^9 + 2*q^13 + 2*q^17 + 3*q^25 + 2*q^29 + 2*q^37 + 2*q^41 + ...

B. C. Berndt, Ramanujan's Notebooks Part III, Springer-Verlag. See p. 139 Example (iv).
J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 102.
R. W. Gosper, Strip Mining in the Abandoned Orefields of Nineteenth Century Mathematics, in Computers in Mathematics (Ed. D. V. Chudnovsky and R. D. Jenks). New York: Dekker, 1990. See p. 279.
R. W. Gosper, Experiments and discoveries in q-trigonometry, in Symbolic Computation, Number Theory, Special Functions, Physics and Combinatorics. Editors: F. G. Garvan and M. E. H. Ismail. Kluwer, Dordrecht, Netherlands, 2001, pp. 79-105. [See Pi_q.]
P. A. MacMahon, Combinatory Analysis, Cambridge Univ. Press, London and New York, Vol. 1, 1915 and Vol. 2, 1916. See vol. 2, p 31, Article 272.
Ivan Niven, Herbert S. Zuckerman and Hugh L. Montgomery, An Introduction to the Theory Of Numbers, Fifth Edition, John Wiley and Sons, Inc., NY 1991, p. 165.

T. D. Noe, Table of n, a(n) for n = 0..10000
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, Tenth Printing, 1972, p. 575, 16.23.1 and 16.23.2.
R. P. Agarwal, Lambert series and Ramanujan, Prod. Indian Acad. Sci. (Math. Sci.), v. 103, n. 3, 1993, pp. 269-293 (see p. 285).
N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 72, Eq (31.2); p. 78, Eq. following (32.25).
J. W. L. Glaisher, On the function chi(n), Quarterly Journal of Pure and Applied Mathematics, 20 (1884), 97-167. See p. 108.
R. W. Gosper, Experiments and discoveries in q-trigonometry, Preprint.
R. W. Gosper, q-Trigonometry: Some Prefatory Afterthoughts.
Michael D. Hirschhorn, The number of representations of a number by various forms, Discrete Mathematics 298 (2005), 205-211.
Nadia Heninger, E. M. Rains and N. J. A. Sloane, On the Integrality of n-th Roots of Generating Functions, arXiv:math/0509316 [math.NT], 2005-2006; J. Combinatorial Theory, Series A, 113 (2006), 1732-1745.
Nicco, Conjectured identity of the product of two theta functions, Math Stackexchange, Sep 09 2015.
K. Ono, S. Robins and P. T. Wahl, On the representation of integers as sums of triangular numbers, Aequationes mathematicae, August 1995, Volume 50, Issue 1-2, pp. 73-94.
Hjalmar Rosengren, Sums of triangular numbers from the Frobenius determinant, arXiv:math/0504272 [math.NT], 2005.
James A. Sellers and Roberto Tauraso, Arithmetic properties of MacMahon-type sums of divisors: the odd case, arXiv:2505.12813 [math.NT], 2025.
Michael Somos, Introduction to Ramanujan theta functions.
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions.

A002654, A008442, A035154, A035181, A035184, A112301, A113406, A113414, A113446, A113652 all satisfy A(4n+1)=a(n).
Cf. A004020, A005883, A104794, A052343, A199015 (partial sums).
Number of ways of writing n as a sum of k triangular numbers, for k=1,...: A010054, A008441, A008443, A008438, A008439, A008440, A226252, A007331, A226253, A226254, A226255, A014787, A014809.
Cf. A002175, A053692, A113407, A121444.
Cf. A000217, A001938, A004018.
Cf. A274621 (reciprocal series).
Cf. A001826, A001842.
Cf. A019669, A079006.

a052343 = (flip div 2) . (+ 1) . a008441
-- Reinhard Zumkeller, Jul 25 2014

A := Basis( ModularForms( Gamma1(8), 1), 420); A[2]; /* Michael Somos, Jan 31 2015 */

sigmamr := proc(n,m,r) local a,d ; a := 0 ; for d in numtheory[divisors](n) do if modp(d,m) = r then a := a+1 ; end if; end do: a; end proc:
A002654 := proc(n) sigmamr(n,4,1)-sigmamr(n,4,3) ; end proc:
A008441 := proc(n) A002654(4*n+1) ; end proc:
seq(A008441(n),n=0..90) ; # R. J. Mathar, Mar 23 2011

Plus@@((-1)^(1/2 (Divisors[4#+1]-1)))& /@ Range[0, 104] (* Ant King, Dec 02 2010 *)
a[ n_] := SeriesCoefficient[ (1/2) EllipticTheta[ 2, 0, q] EllipticTheta[ 3, 0, q], {q, 0, n + 1/4}]; (* Michael Somos, Jun 19 2012 *)
a[ n_] := SeriesCoefficient[ (1/4) EllipticTheta[ 2, 0, q]^2, {q, 0, 2 n + 1/2}]; (* Michael Somos, Jun 19 2012 *)
a[ n_] := If[ n < 0, 0, DivisorSum[ 4 n + 1, (-1)^Quotient[#, 2] &]];  (* Michael Somos, Jun 08 2014 *)
QP = QPochhammer; s = QP[q^2]^4/QP[q]^2 + O[q]^100; CoefficientList[s, q] (* Jean-François Alcover, Nov 27 2015, adapted from PARI *)
TriangleQ[n_] := IntegerQ@Sqrt[8n +1]; Table[Count[FrobeniusSolve[{1, 1}, n], {?TriangleQ}], {n, 0, 104}] (* Robert G. Wilson v, Apr 15 2017 *)

{a(n) = if( n<1, n==0, polcoeff( sum(k=0, (sqrtint(8*n + 1) - 1)\2, x^(k * (k+1)/2), x * O(x^n))^2, n) )};

{a(n) = if( n<0, 0, n = 4*n + 1; sumdiv(n, d, (-1)^(d\2)))}; /* Michael Somos, Sep 02 2005 */

{a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^4 / eta(x + A)^2, n))};

{a(n) = if( n<0, 0, n = 4*n + 1; sumdiv( n, d, (d%4==1) - (d%4==3)))}; /* Michael Somos, Sep 14 2005 */

{ my(q='q+O('q^166)); Vec(eta(q^2)^4 / eta(q)^2) } \\ Joerg Arndt, Apr 16 2017

ModularForms( Gamma1(8), 1, prec=420).1; # Michael Somos, Jun 08 2014

This sequence is the quadrisection of many sequences. Here are two examples:
a(n) = A002654(4n+1), the difference between the number of divisors of 4*n+1 of form 4*k+1 and the number of form 4*k-1. - David Broadhurst, Oct 20 2002
a(n) = b(4*n + 1), where b(n) is multiplicative and b(2^e) = 0^e, b(p^e) = (1 + (-1)^e) / 2 if p == 3 (mod 4), b(p^e) = e+1 if p == 1 (mod 4). - Michael Somos, Sep 14 2005
G.f.: (Sum_{k>=0} x^((k^2 + k)/2))^2 = (Sum_{k>=0} x^(k^2 + k)) * (Sum_{k in Z} x^(k^2)).
Expansion of Jacobi theta (theta_2(0, sqrt(q)))^2 / (4 * q^(1/4)).
Sum[d|(4n+1), (-1)^((d-1)/2) ].
Given g.f. A(x), then B(q) = q * A(q^4) satisfies 0 = f(B(q), B(q^2), B(q^4)) where f(u, v, w) = v^3 + 4 * v * w^2 - u^2 * w. - Michael Somos, Sep 14 2005
Given g.f. A(x), then B(q) = q * A(q^4) satisfies 0 = f(B(q), B(q^2), B(q^3), B(q^6)) where f(u1, u2, u3, u6) = u1 * u3 - (u2 - u6) * (u2 + 3*u6). - Michael Somos, Sep 14 2005
Expansion of Jacobi k/(4*q^(1/2)) * (2/Pi)* K(k) in powers of q^2. - Michael Somos, Sep 14 2005. Convolution of A001938 and A004018. This appears in the denominator of the Jacobi sn and cn formula given in the Abramowitz-Stegun reference, p. 575, 16.23.1 and 16.23.2, where m=k^2. - Wolfdieter Lang, Jul 05 2016
G.f.: Sum_{k>=0} a(k) * x^(2*k) = Sum_{k>=0} x^k / (1 + x^(2*k + 1)).
G.f.: Sum_{k in Z} x^k / (1 - x^(4*k + 1)). - Michael Somos, Nov 03 2005
Expansion of psi(x)^2 = phi(x) * psi(x^2) in powers of x where phi(), psi() are Ramanujan theta functions.
Moebius transform is period 8 sequence [ 1, -1, -1, 0, 1, 1, -1, 0, ...]. - Michael Somos, Jan 25 2008
G.f. is a period 1 Fourier series which satisfies f(-1 / (8 t)) = 1/2 (t/i) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A104794.
Euler transform of period 2 sequence [ 2, -2, ...].
G.f.: q^(-1/4) * eta(q^2)^4 / eta(q)^2. See also the Fine reference.
a(n) = Sum_{k=0..n} A010054(k)*A010054(n-k). - Reinhard Zumkeller, Nov 03 2009
A004020(n) = 2 * a(n). A005883(n) = 4 * a(n).
Convolution square of A010054.
G.f.: Product_{k>0} (1 - x^(2*k))^2 / (1 - x^(2*k-1))^2.
a(2*n) = A113407(n). a(2*n + 1) = A053692(n). a(3*n) = A002175(n). a(3*n + 1) = 2 * A121444(n). a(9*n + 2) = a(n). a(9*n + 5) = a(9*n + 8) = 0. - Michael Somos, Jun 08 2014
G.f.: exp( Sum_{n>=1} 2*(x^n/n) / (1 + x^n) ). - Paul D. Hanna, Mar 01 2016
a(n) = A001826(2+8*n) - A001842(2+8*n), the difference between the number of divisors 1 (mod 4) and 3 (mod 4) of 2+8*n. See the Ono et al. link, Corollary 1, or directly the Niven et al. reference, p. 165, Corollary (3.23). - Wolfdieter Lang, Jan 11 2017
Expansion of continued fraction 1 / (1 - x^1 + x^1*(1 - x^1)^2 / (1 - x^3 + x^2*(1 - x^2)^2 / (1 - x^5 + x^3*(1 - x^3)^2 / ...))) in powers of x^2. - Michael Somos, Apr 20 2017
Given g.f. A(x), and B(x) is the g.f. for A079006, then B(x) = A(x^2) / A(x) and B(x) * B(x^2) * B(x^4) * ... = 1 / A(x). - Michael Somos, Apr 20 2017
a(0) = 1, a(n) = (2/n)*Sum_{k=1..n} A002129(k)*a(n-k) for n > 0. - Seiichi Manyama, May 06 2017
From Paul D. Hanna, Aug 10 2019: (Start)
G.f.: Sum_{n>=0} x^n * Sum_{k=0..n} binomial(n,k) * (x^(2*n+1) - x^(2*k))^(n-k) = Sum_{n>=0} a(n)*x^(2*n).
G.f.: Sum_{n>=0} x^n * Sum_{k=0..n} binomial(n,k) * (x^(2*n+1) + x^(2*k))^(n-k) * (-1)^k = Sum_{n>=0} a(n)*x^(2*n). (End)
From Peter Bala, Jan 05 2021: (Start)
G.f.: Sum_{n = -oo..oo} x^(4*n^2+2*n) * (1 + x^(4*n+1))/(1 - x^(4*n+1)). See Agarwal, p. 285, equation 6.20 with i = j = 1 and mu = 4.
For prime p of the form 4*k + 3, a(n*p^2 + (p^2 - 1)/4) = a(n).
If n > 0 and p are coprime then a(n*p + (p^2 - 1)/4) = 0. The proofs are similar to those given for the corresponding results for A115110. Cf. A000729.
For prime p of the form 4*k + 1 and for n not congruent to (p - 1)/4 (mod p) we have a(n*p^2 + (p^2 - 1)/4) = 3*a(n) (since b(n), where b(4*n+1) = a(n), is multiplicative). (End)
From Peter Bala, Mar 22 2021: (Start)
G.f. A(q) satisfies:
A(q^2) = Sum_{n = -oo..oo} q^n/(1 - q^(4*n+2)) (set z = q, alpha = q^2, mu = 4 in Agarwal, equation 6.15).
A(q^2) = Sum_{n = -oo..oo} q^(2*n)/(1 - q^(4*n+1)) (set z = q^2, alpha = q, mu = 4 in Agarwal, equation 6.15).
A(q^2) = Sum_{n = -oo..oo} q^n/(1 + q^(2*n+1))^2 = Sum_{n = -oo..oo} q^(3*n+1)/(1 + q^(2*n+1))^2. (End)
G.f.: Sum_{k>=0} a(k) * q^k = Sum_{k>=0} (-1)^k * q^(k*(k+1)) + 2 * Sum_{n>=1, k>=0} (-1)^k * q^(k*(k+2*n+1)+n). - Mamuka Jibladze, May 17 2021
G.f.: Sum_{k>=0} a(k) * q^k = Sum_{k>=0} (-1)^k * q^(k*(k+1)) * (1 + q^(2*k+1))/(1 - q^(2*k+1)). - Mamuka Jibladze, Jun 06 2021
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Pi/2 (A019669). - Amiram Eldar, Oct 15 2022

More terms and information from Michael Somos, Mar 23 2003

A321543 a(n) = Sum_{d|n} (-1)^(d-1)*d^2.

Original entry on oeis.org

1, -3, 10, -19, 26, -30, 50, -83, 91, -78, 122, -190, 170, -150, 260, -339, 290, -273, 362, -494, 500, -366, 530, -830, 651, -510, 820, -950, 842, -780, 962, -1363, 1220, -870, 1300, -1729, 1370, -1086, 1700, -2158, 1682, -1500, 1850, -2318, 2366, -1590, 2210, -3390, 2451, -1953, 2900, -3230, 2810, -2460, 3172

Offset: 1

N. J. A. Sloane, Nov 23 2018

Seiichi Manyama, Table of n, a(n) for n = 1..10000
J. W. L. Glaisher, On the representations of a number as the sum of two, four, six, eight, ten, and twelve squares, Quart. J. Math. 38 (1907), 1-62 (see p. 4 and p. 8).
Index entries for sequences mentioned by Glaisher.

Apart from signs, same as A064027.
Glaisher's zeta_i, i=0..12: A048272, A002129, A321543, A138503, A279395, A321544, A321545, A321546, A321547, A321548, A321549, A321550, A321551.
Cf. A321552 - A321565, A321807 - A321836 for similar sequences.

with(numtheory):
a := n -> add( (-1)^(d-1)*d^2, d in divisors(n) ): seq(a(n), n = 1..40);
#  Peter Bala, Jan 11 2021

f[p_, e_] := (p^(2*e + 2) - 1)/(p^2 - 1); f[2, e_] := 2 - (2^(2*e + 2) - 1)/3; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 55] (* Amiram Eldar, Nov 04 2022 *)

apply( a(n)=sumdiv(n, d, (-1)^(d-1)*d^2), [1..30]) \\ M. F. Hasler, Nov 26 2018

G.f.: Sum_{k>=1} (-1)^(k-1)*k^2*x^k/(1 - x^k). - Ilya Gutkovskiy, Dec 23 2018
G.f.: Sum_{n >= 1} x^n*(1 - x^n)/(1 + x^n)^3. - Peter Bala, Jan 11 2021
Multiplicative with a(2^e) = 2 - (2^(2*e + 2) - 1)/3, and a(p^e) = (p^(2*e + 2) - 1)/(p^2 - 1) for p > 2. - Amiram Eldar, Nov 04 2022

A008443 Number of ordered ways of writing n as the sum of 3 triangular numbers.

Original entry on oeis.org

1, 3, 3, 4, 6, 3, 6, 9, 3, 7, 9, 6, 9, 9, 6, 6, 15, 9, 7, 12, 3, 15, 15, 6, 12, 12, 9, 12, 15, 6, 13, 21, 12, 6, 15, 9, 12, 24, 9, 18, 12, 9, 18, 15, 12, 13, 24, 9, 15, 24, 6, 18, 27, 6, 12, 15, 18, 24, 21, 15, 12, 27, 9, 13, 18, 15, 27, 27, 9, 12, 27, 15, 24, 21, 12, 15, 30, 15, 12

Offset: 0

N. J. A. Sloane

Fermat asserted that every number is the sum of three triangular numbers. This was proved by Gauss, who recorded in his Tagebuch entry for Jul 10 1796 that: EYPHEKA! num = DELTA + DELTA + DELTA. See also Gauss, DA, art. 293.
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Andrews (2016), Theorem 2, shows that A008443(n) = A290735(n) + A290737(n) + A290739(n). = N. J. A. Sloane, Aug 10 2017

			5 can be written as 3+1+1, 1+3+1, 1+1+3, so a(5) = 3.
G.f. = 1 + 3*x + 3*x^2 + 4*x^3 + 6*x^4 + 3*x^5 + 6*x^6 + 9*x^7 + 3*x^8 + ...
G.f. = q^3 + 3*q^11 + 3*q^19 + 4*q^27 + 6*q^35 + 3*q^43 + 6*q^51 + 9*q^59 + 3*q^67 + ...

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 102.
C. F. Gauss, Disquisitiones Arithmeticae, Yale University Press, 1966, New Haven and London, p. 342, art. 293.
M. Nathanson, Additive Number Theory: The Classical Bases, Graduate Texts in Mathematics, Volume 165, Springer-Verlag, 1996. See Chapter 1.

N. J. A. Sloane, Table of n, a(n) for n = 0..20000 (first 5050 terms from T. D. Noe)
George E. Andrews, EYPHKA! num = Delta + Delta + Delta, J. Number Theory 23 (1986), 285-293. [The Y in the title is really the Greek letter Upsilon and Delta is really the Greek letter of that name.]
George E. Andrews, The Bhargava-Adiga Summation and Partitions, 2016.
M. Doring, J. Haidenbauer, U.-G. Meissner, and A. Rusetsky, Dynamical coupled-channel approaches on a momentum lattice, arXiv preprint arXiv:1108.0676 [hep-lat], 2011.
M. D. Hirschhorn and J. A. Sellers, Partitions into three triangular numbers, Australasian Journal of Combinatorics, Volume 30 (2004), Pages 307-318; Submission.
M. D. Hirschhorn and J. A. Sellers, On Representations Of A Number As A Sum Of Three Triangles, Acta Arithmetica 77 (1996), 289-301.
K. Ono, S. Robins, and P. T. Wahl, On the representation of integers as sums of triangular numbers, Aequationes mathematicae, August 1995, Volume 50, Issue 1-2, pp 73-94.
Michael Somos, Introduction to Ramanujan theta functions.
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions.

Number of ways of writing n as a sum of k triangular numbers, for k=1,...: A010054, A008441, A008443, A008438, A008439, A008440,A226252, A007331, A226253, A226254, A226255, A014787, A014809.
Cf. A053604, A002636.
Partial sums are in A038835.
Cf. A005869, A005875, A005878, A005886.
Cf. A290733-A290740.

Basis( ModularForms( Gamma0(16), 3/2), 630)[4]; /* Michael Somos, Aug 26 2015 */

s1 := sum(q^(n*(n+1)/2), n=0..30): s2 := series(s1^3, q, 250): for i from 0 to 200 do printf(`%d,`,coeff(s2, q, i)) od:

s1 = Sum[q^(n (n + 1)/2), {n, 0, 12}]; s2 = Series[s1^3, {q, 0, 80}]; CoefficientList[s2, q] (* Jean-François Alcover, Oct 04 2011, after Maple *)
a[ n_] := SeriesCoefficient[ (1/8) EllipticTheta[ 2, 0, q]^3, {q, 0, 2 n + 3/4}]; (* Michael Somos, May 29 2012 *)
QP = QPochhammer; CoefficientList[(QP[q^2]^2/QP[q])^3 + O[q]^80, q] (* Jean-François Alcover, Nov 24 2015 *)

{a(n) = if( n<0, 0, polcoeff( sum(k=0, (sqrtint(8*n + 1) - 1)\2, x^((k^2 + k)/2), x * O(x^n))^3, n))}; /* Michael Somos, Oct 25 2006 */

{a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x^2 + A)^2 / eta(x + A))^3, n))}; /* Michael Somos, Oct 25 2006 */

Expansion of Jacobi theta constant theta_2^3 /8. G.f. is cube of g.f. for A010054.
Expansion of psi(q)^3 in powers of q where psi() is a Ramanujan theta function (A010054). - Michael Somos, Oct 25 2006
Expansion of q^(-3/8) * (eta(q^2)^2 / eta(q))^3 in powers of q. - Michael Somos, May 29 2012
Euler transform of period 2 sequence [ 3, -3, ...]. - Michael Somos, Oct 25 2006
G.f. is a period 1 Fourier series which satisfies f(-1 / (16 t)) = 2^(-3/2) (t/i)^(3/2) g(t) where q = exp(2 Pi i t) and g() is the g.f. for A213384. - Michael Somos, Jun 23 2012
a(3*n) = A213627(n). a(3*n + 1) = 3 * A213617(n). a(3*n + 2) = A181648(n). - Michael Somos, Jun 23 2012
G.f.: (Sum_{k>0} x^((k^2 - k)/2))^3 = (Product_{k>0} (1 + x^k) * (1 - x^(2*k)))^3. - Michael Somos, May 29 2012
a(n) = A005869(n)/2 = A005886(n)/4 = A005878(n)/8.
a(n) = A005875(8*n+3)/8. See, e.g., the Ono et al. link: The case k=3. - Wolfdieter Lang, Jan 12 2017
a(0) = 1, a(n) = (3/n)*Sum_{k=1..n} A002129(k)*a(n-k) for n > 0. - Seiichi Manyama, May 06 2017

More terms from James Sellers, Feb 07 2001

A007331 Fourier coefficients of E_{infinity,4}.

Original entry on oeis.org

0, 1, 8, 28, 64, 126, 224, 344, 512, 757, 1008, 1332, 1792, 2198, 2752, 3528, 4096, 4914, 6056, 6860, 8064, 9632, 10656, 12168, 14336, 15751, 17584, 20440, 22016, 24390, 28224, 29792, 32768, 37296, 39312, 43344, 48448, 50654, 54880, 61544, 64512

Offset: 0

N. J. A. Sloane, Mira Bernstein

E_{infinity,4} is the unique normalized weight-4 modular form for Gamma_0(2) with simple zeros at i*infinity. Since this has level 2, it is not a cusp form, in contrast to A002408.
a(n+1) is the number of representations of n as a sum of 8 triangular numbers (from A000217). See the Ono et al. link, Theorem 5.
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
a(n) gives the sum of cubes of divisors d of n such that n/d is odd. This is called sigma^#3(n) in the Ono et al. link. See a formula below. - _Wolfdieter Lang, Jan 12 2017

			G.f. = q + 8*q^2 + 28*q^3 + 64*q^4 + 126*q^5 + 224*q^6 + 344*q^7 + 512*q^8 + ...

B. C. Berndt, Ramanujan's Notebooks Part III, Springer-Verlag, see p. 139, Ex (ii).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1001 from T. D. Noe)
B. Brent, Quadratic Minima and Modular Forms, Experimental Mathematics, v.7 no.3, 257-274.
H. H. Chan and C. Krattenthaler, Recent progress in the study of representations of integers as sums of squares, arXiv:math/0407061 [math.NT], 2004.
J. H. Conway and N. J. A. Sloane, Sphere Packings, Lattices and Groups, Springer-Verlag, p. 187.
J. W. L. Glaisher, On the representations of a number as the sum of two, four, six, eight, ten, and twelve squares, Quart. J. Math. 38 (1907), 1-62 (see p. 4 and p. 8).
Masao Koike, Modular forms on non-compact arithmetic triangle groups, Unpublished manuscript [Extensively annotated with OEIS A-numbers by N. J. A. Sloane, Feb 14 2021. I wrote 2005 on the first page but the internal evidence suggests 1997.]
K. Ono, S. Robins and P. T. Wahl, On the representation of integers as sums of triangular numbers, Aequationes mathematicae, August 1995, Volume 50, Issue 1-2, pp 73-94. Theorem 5.
H. Rosengren, Sums of triangular numbers from the Frobenius determinant, arXiv:math/0504272 [math.NT], 2005.
Index entries for sequences mentioned by Glaisher.

Cf. A002408, A004017, A035016, A045825, A076577, A096960, A222072.

Number of ways of writing n as a sum of k triangular numbers, for k=1,...: A010054, A008441, A008443, A008438, A008439, A008440, A226252, A007331, A226253, A226254, A226255, A014787, A014809, A076577.

Basis( ModularForms( Gamma0(2), 4), 10) [2]; /* Michael Somos, May 27 2014 */

nmax:=40: seq(coeff(series(x*(product((1-x^k)^8*(1+x^k)^16, k=1..nmax)), x, n+1), x, n), n=0..nmax); # Vaclav Kotesovec, Oct 14 2015

Prepend[Table[Plus @@ (Select[Divisors[k + 1], OddQ[(k + 1)/#] &]^3), {k, 0, 39}], 0] (* Ant King, Dec 04 2010 *)
a[ n_] := SeriesCoefficient[ EllipticTheta[ 2, 0, q^(1/2)]^8 / 256, {q, 0, n}]; (* Michael Somos, Jun 04 2013 *)
a[ n_] := If[ n < 1, 0, Sum[ d^3 Boole[ OddQ[ n/d]], {d, Divisors[ n]}]]; (* Michael Somos, Jun 04 2013 *)
f[n_] := Total[(2n/Select[ Divisors[ 2n], Mod[#, 4] == 2 &])^3]; Flatten[{0, Array[f, 40] }] (* Robert G. Wilson v, Mar 26 2015 *)
nmax=60; CoefficientList[Series[x*Product[(1-x^k)^8 * (1+x^k)^16, {k,1,nmax}],{x,0,nmax}], x] (* Vaclav Kotesovec, Oct 14 2015 *)
QP = QPochhammer; s = q * (QP[-1, q]/2)^16 * QP[q]^8 + O[q]^50; CoefficientList[s, q] (* Jean-François Alcover, Dec 01 2015, adapted from PARI *)

{a(n) = if( n<1, 0, sumdiv( n, d, (n/d%2) * d^3))}; /* Michael Somos, May 31 2005 */

{a(n) = local(A); if( n<1, 0, n--; A = x * O(x^n); polcoeff( (eta(x^2 + A)^2 / eta(x + A))^8, n))}; /* Michael Somos, May 31 2005 */

a(n)=my(e=valuation(n,2)); 8^e * sigma(n/2^e, 3) \\ Charles R Greathouse IV, Sep 09 2014

from sympy import divisors
def a(n):
    return 0 if n == 0 else sum(((n//d)%2)*d**3 for d in divisors(n))
print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 24 2017

ModularForms( Gamma0(2), 4, prec=33).1; # Michael Somos, Jun 04 2013

G.f.: q * Product_{k>=1} (1-q^k)^8 * (1+q^k)^16. - corrected by Vaclav Kotesovec, Oct 14 2015

a(n) = Sum_{0

G.f.: Sum_{n>0} n^3*x^n/(1-x^(2*n)). - Vladeta Jovovic, Oct 24 2002

Expansion of Jacobi theta constant theta_2(q)^8 / 256 in powers of q.

Expansion of eta(q^2)^16 / eta(q)^8 in powers of q. - Michael Somos, May 31 2005

Expansion of x * psi(x)^8 in powers of x where psi() is a Ramanujan theta function. - Michael Somos, Jan 15 2012

Expansion of (Q(x) - Q(x^2)) / 240 in powers of x where Q() is a Ramanujan Lambert series. - Michael Somos, Jan 15 2012

Expansion of E_{gamma,2}^2 * E_{0,4} in powers of q.

Euler transform of period 2 sequence [8, -8, ...]. - Michael Somos, May 31 2005

G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = v^3 - u^2*w + 16*u*v*w - 32*v^2*w + 256*v*w^2. - Michael Somos, May 31 2005

G.f. is a period 1 Fourier series which satisfies f(-1 / (2 t)) = 16^(-1) (t / i)^4 g(t) where q = exp(2 Pi i t) and g() is the g.f. for A035016. - Michael Somos, Jan 11 2009

Multiplicative with a(2^e) = 2^(3e), a(p^e) = (p^(3(e+1))-1)/(p^3-1). - Mitch Harris, Jun 13 2005

Dirichlet convolution of A154955 by A001158. Dirichlet g.f. zeta(s)*zeta(s-3)*(1-1/2^s). - R. J. Mathar, Mar 31 2011

A002408(n) = -(-1)^n * a(n).

Convolution square of A008438. - Michael Somos, Jun 15 2014

a(1) = 1, a(n) = (8/(n-1))*Sum_{k=1..n-1} A002129(k)*a(n-k) for n > 0. - Seiichi Manyama, May 06 2017

Sum_{k=1..n} a(k) ~ c * n^4, where c = Pi^4/384 = 0.253669... (A222072). - Amiram Eldar, Oct 19 2022

Additional comments from Barry Brent (barryb(AT)primenet.com)
Wrong Maple program replaced by Vaclav Kotesovec, Oct 14 2015
a(0)=0 prepended by Vaclav Kotesovec, Oct 14 2015

A181983 a(n) = (-1)^(n+1) * n.

Original entry on oeis.org

0, 1, -2, 3, -4, 5, -6, 7, -8, 9, -10, 11, -12, 13, -14, 15, -16, 17, -18, 19, -20, 21, -22, 23, -24, 25, -26, 27, -28, 29, -30, 31, -32, 33, -34, 35, -36, 37, -38, 39, -40, 41, -42, 43, -44, 45, -46, 47, -48, 49, -50, 51, -52, 53, -54, 55, -56, 57, -58, 59

Offset: 0

Michael Somos, Apr 04 2012

This is the Lucas U(-2,1) sequence. - R. J. Mathar, Jan 08 2013
Apparently the Mobius transform of A002129. - R. J. Mathar, Jan 08 2013
For n>0, a(n) is also the determinant of the symmetric n X n matrix M defined by M(i,j) = max(i,j) for 1 <= i,j <= n. - Enrique Pérez Herrero, Jan 14 2013
The sums of the terms of this sequence is the divergent series 1 - 2 + 3 - 4 + ... . Euler summed it to 1/4 which was one of the first examples of summing divergent series. - Michael Somos, Jun 05 2013

			G.f. = x - 2*x^2 + 3*x^3 - 4*x^4 + 5*x^5 - 6*x^6 + 7*x^7 - 8*x^8 + 9*x^9 + ...

Paulo Ribenboim, My Numbers, My Friends: Popular Lectures on Number Theory, Springer-Verlag, NY, 2000, p. 6.

Enrique Pérez Herrero, Max Determinant, 2013.
Wikipedia, 1 - 2 + 3 - 4 + ....
Wikipedia, Lucas sequence.
Index entries for linear recurrences with constant coefficients, signature (-2,-1).
Index entries for Lucas sequences.

Cf. A000108, A000169, A001057, A001477, A001787, A002129, A038608, A049347, A154955, A274922.

a181983 = negate . a038608
a181983_list = [0, 1] ++ map negate
   (zipWith (+) a181983_list (map (* 2) $ tail a181983_list))
-- Reinhard Zumkeller, Mar 20 2013

[(-1)^(n+1)*n: n in [0..30]]; // G. C. Greubel, Aug 11 2018

A181983:=n->-(-1)^n * n; seq(A181983(n), n=0..100); # Wesley Ivan Hurt, Feb 26 2014

a[ n_] := -(-1)^n n;
a[ n_] := Sign[n] SeriesCoefficient[ x / (1 + x)^2, {x, 0, Abs @n}];
a[ n_] := Sign[n] (Abs @n)! SeriesCoefficient[ x / Exp[ x], {x, 0, Abs @n}];
CoefficientList[Series[x/(1+x)^2,{x,0,60}],x] (* or *) LinearRecurrence[{-2,-1},{0,1},60] (* or *) Table[If[OddQ[n],n,-n],{n,0,60}]  (* Harvey P. Dale, Apr 22 2022 *)

PARI
```
{a(n) = -(-1)^n * n};
```

G.f.: x / (1 + x)^2.
E.g.f.: x / exp(x).
a(n) = -a(-n) = -(-1)^n * A001477(n) for all n in Z.
a(n+1) = p(n+1) where p(x) is the unique degree-n polynomial such that p(k) = Bernoulli(k) for k = 0, 1, ..., n.
A001787(n) = p(0) where p(x) is the unique degree-n polynomial such that p(k) = a(k) for k = 1, ..., n+1. - Michael Somos, Jun 05 2013
Euler transform of length 2 sequence [-2, 2].
Series reversion of g.f. is A000108(n) (Catalan numbers) with a(0)=0.
Series reversion of e.g.f. is A000169. INVERT transform omitting a(0)=0 is A049347. PSUM transform is A001057. BINOMIAL transform is A154955. - Michael Somos, Jun 05 2013
n * a(n) = A162395(n). - Michael Somos, Jun 05 2013
a(n) = - A038608(n). - Reinhard Zumkeller, Mar 20 2013
a(n+2) = a(n) - 2*(-1)^n. - G. C. Greubel, Aug 11 2018
a(n) = - A274922(n) if n>0. - Michael Somos, Sep 24 2019
From Amiram Eldar, Oct 24 2023: (Start)
Multiplicative with a(2^e) = -2^e, and a(p^e) = p^e for an odd prime p.
Dirichlet g.f.: zeta(s-1) * (1-2^(2-s)). (End)

A059358 Coefficients in expansion of Sum_{n >= 1} x^n/(1-x^n)^4.

Original entry on oeis.org

0, 1, 5, 11, 25, 36, 71, 85, 145, 176, 260, 287, 455, 456, 649, 726, 961, 970, 1376, 1331, 1820, 1866, 2315, 2301, 3175, 2961, 3736, 3830, 4729, 4496, 5966, 5457, 6945, 6842, 8114, 7890, 10196, 9140, 11215, 11126, 13420, 12342, 15730, 14191, 17515, 17106, 19601

Offset: 0

N. J. A. Sloane, Jan 27 2001

nonn

G. C. Greubel, Table of n, a(n) for n = 0..1000

Cf. A000203, A000292, A001157, A001158, A002129, A007437, A013662, A048272, A073570, A101289.

a:= proc(n) option remember;
      add(d*(d+1)*(d+2)/6, d=numtheory[divisors](n))
    end:
seq(a(n), n=0..60);  # Alois P. Heinz, Jun 12 2023

With[{nn=50},CoefficientList[Series[Sum[x^n/(1-x^n)^4,{n,nn}],{x,0,nn}],x]] (* Harvey P. Dale, May 14 2013 *)

a(n) = if(n==0, 0, sumdiv(n, d, binomial(d+2, 3))); \\ Seiichi Manyama, Apr 19 2021

a(n) = if(n==0, 0, my(f = factor(n)); (sigma(f, 3) + 3*sigma(f, 2) + 2 * sigma(f)) / 6); \\ Amiram Eldar, Dec 29 2024

a(n) = (1/6)*(sigma_3(n) + 3*sigma_2(n) + 2*sigma_1(n)), i.e., this sequence is the inverse Möbius transform of tetrahedral (or pyramidal) numbers: n*(n+1)(n+2)/6 with g.f. 1/(1-x)^4 (cf. A000292). - Vladeta Jovovic, Aug 31 2002
L.g.f.: -log(Product_{k>=1} (1 - x^k)^((k+1)*(k+2)/6)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, May 21 2018
From Amiram Eldar, Dec 29 2024: (Start)
Dirichlet g.f.: zeta(s) * (zeta(s-3) + 3*zeta(s-2) + 2*zeta(s-1)) / 6.
Sum_{k=1..n} a(k) ~ (zeta(4)/24) * n^4. (End)

A008439 Expansion of Jacobi theta constant theta_2^5 /32.

Original entry on oeis.org

1, 5, 10, 15, 25, 31, 35, 55, 60, 60, 90, 90, 95, 135, 125, 126, 170, 180, 175, 215, 220, 195, 285, 280, 245, 340, 300, 320, 405, 350, 351, 450, 465, 415, 515, 480, 425, 620, 590, 505, 655, 625, 590, 755, 660, 650, 805, 770, 755, 855, 841, 730, 1045, 960, 770, 1100

Offset: 0

N. J. A. Sloane

nonn

Also number of ways of writing n as a sum of five triangular numbers. - N. J. A. Sloane, Jun 01 2013

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 102.

Seiichi Manyama, Table of n, a(n) for n = 0..10000
K. Ono, S. Robins and P. T. Wahl, On the representation of integers as sums of triangular numbers, Aequationes mathematicae, August 1995, Volume 50, Issue 1-2, pp 73-94.

Number of ways of writing n as a sum of k triangular numbers, for k=1,...: A010054, A008441, A008443, A008438, A008439, A008440, A226252, A007331, A226253, A226254, A226255, A014787, A014809.

a002129[n_]:=-Sum[(-1)^d*d, {d, Divisors[n]}]; a[n_]:=a[n]=If[n==0, 1, 5 Sum[a002129[k] a[n - k], {k, n}]/n]; Table[a[n], {n, 0, 100}] (* Indranil Ghosh, Aug 02 2017 *)

a(0) = 1, a(n) = (5/n)*Sum_{k=1..n} A002129(k)*a(n-k) for n > 0. - Seiichi Manyama, May 06 2017

G.f.: exp(Sum_{k>=1} 5*(x^k/k)/(1 + x^k)). - Ilya Gutkovskiy, Jul 31 2017

More terms from Seiichi Manyama, May 05 2017

A008440 Expansion of Jacobi theta constant theta_2^6 /(64q^(3/2)).

Original entry on oeis.org

1, 6, 15, 26, 45, 66, 82, 120, 156, 170, 231, 276, 290, 390, 435, 438, 561, 630, 651, 780, 861, 842, 1020, 1170, 1095, 1326, 1431, 1370, 1716, 1740, 1682, 2016, 2145, 2132, 2415, 2550, 2353, 2850, 3120, 2810, 3321, 3486, 3285, 3906, 4005, 3722, 4350

Offset: 0

N. J. A. Sloane

Number of representations of n as sum of 6 triangular numbers. - Michel Marcus, Oct 24 2012. See the Ono et al. link.

			G.f. = 1 + 6*x + 15*x^2 + 26*x^3 + 45*x^4 + 66*x^5 + 82*x^6 + ... - _Michael Somos_, Jun 25 2019
G.f. = q^3 + 6*q^7 + 15*q^11 + 26*q^15 + 45*q^19 + 66*q^23 + 82*q^27 + ...

B. C. Berndt, Fragments by Ramanujan on Lambert series, in Number Theory and Its Applications, K. Gyory and S. Kanemitsu, eds., Kluwer, Dordrecht, 1999, pp. 35-49, see Entry 6.
J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 102.

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from Vincenzo Librandi)
B. C. Berndt, Fragments by Ramanujan on Lambert series.
K. Ono, S. Robins and P. T. Wahl, On the representation of integers as sums of triangular numbers, Aequationes mathematicae, August 1995, Volume 50, Issue 1-2, pp 73-94. Theorem 4.

CoefficientList[(QPochhammer[q^2]^2 / QPochhammer[q])^6 + O[q]^50, q] (* Jean-François Alcover, Nov 05 2015 *)
a[ n_] := If[ n < 0, 0, -DivisorSum[ 4 n + 3, Re[I^(# - 1)] #^2 &] / 8]; (* Michael Somos, Jun 25 2019 *)

{a(n) = if( n<0, 0, polcoeff( sum(k=0, (sqrtint(8*n+1)-1)\2, x^((k^2+k)/2), x * O(x^n))^6, n))}; /* Michael Somos, May 23 2006 */

{a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x^2 + A)^2 / eta(x + A))^6, n))}; /* Michael Somos, May 23 2006 */

{a(n)= -sumdiv(4*n + 3, d, real(I^(d-1))*d^2)/8}; /* Michael Somos, Oct 24 2012 */

Expansion of Ramanujan phi^6(q) in powers of q.
Expansion of q^(-3/4)(eta(q^2)^2/eta(q))^6 in powers of q.
Euler transform of period 2 sequence [6, -6, ...]. - Michael Somos, May 23 2006
G.f.: (Sum_{n>=0} x^((n^2+n)/2))^6.
a(n) = (-1/8)*Sum_{d divides (4n+3)} Chi_2(4;d)*d^2. - Michel Marcus, Oct 24 2012. See the Ono et al. link. Theorem 4.
a(n) =(-1/8)*A002173(4*n+3). This is the preceding formula. - Wolfdieter Lang, Jan 12 2017
a(0) = 1, a(n) = (6/n)*Sum_{k=1..n} A002129(k)*a(n-k) for n > 0. - Seiichi Manyama, May 06 2017
G.f.: exp(Sum_{k>=1} 6*(x^k/k)/(1 + x^k)). - Ilya Gutkovskiy, Jul 31 2017

A113184 Absolute difference between sum of odd divisors of n and sum of even divisors of n.

Original entry on oeis.org

1, 1, 4, 5, 6, 4, 8, 13, 13, 6, 12, 20, 14, 8, 24, 29, 18, 13, 20, 30, 32, 12, 24, 52, 31, 14, 40, 40, 30, 24, 32, 61, 48, 18, 48, 65, 38, 20, 56, 78, 42, 32, 44, 60, 78, 24, 48, 116, 57, 31, 72, 70, 54, 40, 72, 104, 80, 30, 60, 120, 62, 32, 104, 125, 84, 48, 68, 90, 96, 48, 72

Offset: 1

Michael Somos, Oct 17 2005

The generating function equals 1/8 at q = Lambda = 0.1076539192... (A072558) the "One-Ninth" constant. - Michael Somos, Jul 21 2006
Absolute value of A002129. - John W. Layman, Sep 27 2012
The Möbius transform is 1, 0, 3, 4, 5, 0, 7, 8, 9, 0, 11, 12, 13, 0, 15, 16, 17, 0, 19, 20, 21, 0, 23, 24, 25, 0, 27, ... - R. J. Mathar, Jan 08 2013

			From _Peter Bala_, Dec 11 2020: (Start)
n = 15: n is a triangular number, so e(n) = (-1)^(n+1)*n = 15 and a(15) = 15 + a(14) + a(12) - a(9) - a(5) =  15 + 8 + 20 - 13 - 6 = 24;
n = 16: n is a not triangular number, so e(n) = 0 and a(16) = a(15) + a(13) - a(10) - a(6) + a(1) =  24 + 14 - 6 - 4 + 1 = 29. (End)

G. H. Hardy, Ramanujan: twelve lectures on subjects suggested by his life and work, AMS Chelsea Publishing, Providence, Rhode Island, 2002, p. 142.

G. C. Greubel, Table of n, a(n) for n = 1..1000
Eric Weisstein's World of Mathematics, One-Ninth Constant.

Cf. A002129(n) = -(-1)^n a(n).

Cf. A072558, A284362.

f[n_]:=Module[{dn=Divisors[n],odn,edn},odn=Select[dn,OddQ];edn=Select[dn,EvenQ];Abs[Total[odn]-Total[edn]]]
f/@Range[80]  (* Harvey P. Dale, Feb 25 2011 *)
max = 80; s = (1/x)*Sum[k*x^k/(1 - (-x)^k), {k, 1, max}] + O[x]^max; CoefficientList[s, x] (* Jean-François Alcover, Dec 04 2015 *)
f[p_, e_] := If[p == 2, 2^(e + 1) - 3, (p^(e + 1) - 1)/(p - 1)]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Jun 12 2022 *)

a(n)=if(n<1, 0, (-1)^n*sumdiv(n,d,(-1)^d*d))

{a(n)=local(A,p,e); if(n<1, 0, A=factor(n); prod(k=1,matsize(A)[1], if(p=A[k,1], e=A[k,2]; if(p==2, 2^(e+1)-3, (p^(e+1)-1)/(p-1)))))}

Multiplicative with a(2^e) = 2^(e+1)-3 if e>0, a(p^e) = (p^(e+1)-1)/(p-1) if p>2.
G.f.: Sum_{k>0} -(-x)^k/(1+(-x)^k)^2 = Sum_{k>0} k*x^k/(1-(-x)^k).
Expansion of (1-(2/Pi)^2(2E(k)-K(k))K(k))/8 in powers of nome q where E(k) and K(k) are complete elliptic integrals and q=exp(-Pi*K(k')/K(k)). - Michael Somos, Jul 21 2006
Bisection: a(2*k-1) = A000203(2*k-1), a(2*k) = A146076(2*k) - A000593(2*k), k >= 1. See the Hardy reference where a(n) = sigma^*1(n). - _Wolfdieter Lang, Jan 07 2017
From Peter Bala, Dec 11 2020: (Start)
a(n) = Sum_{d | n, d != 2 (mod 4)} d.
O.g.f.: Sum_{k >= 1, k != 2 (mod 4)} k*x^k/(1 - x^k). Cf. A284362.
Define a(n) = 0 for n < 1. Then a(n) = e(n) + a(n-1) + a(n-3) - a(n-6) - a(n-10) + + - -, where [1, 3, 6, 10, ...] is the sequence of triangular numbers A000217, and e(n) = (-1)^(n+1)*n if n is a triangular number; otherwise e(n) = 0. Examples of this recurrence are given below. (End)
Dirichlet g.f.: Sum_{n>0} a(n)/n^s = zeta(s) * zeta(s-1) * (1+2^(3-3*s)) / (1+2^(1-s)). - Werner Schulte, Jan 23 2021
Sum_{k=1..n} a(k) ~ Pi^2 * n^2 / 16. - Vaclav Kotesovec, Aug 20 2021

Name corrected by Wolfdieter Lang, Jan 07 2017

A226252 Number of ways of writing n as the sum of 7 triangular numbers.

Original entry on oeis.org

1, 7, 21, 42, 77, 126, 175, 253, 357, 434, 567, 735, 833, 1057, 1302, 1400, 1708, 2037, 2191, 2597, 3003, 3151, 3619, 4242, 4389, 4935, 5691, 5740, 6594, 7434, 7371, 8400, 9303, 9506, 10626, 11592, 11585, 12761, 14427, 14203, 15519, 17241, 16808, 18788, 20559, 19950, 21882, 23898, 23786

Offset: 0

N. J. A. Sloane, Jun 01 2013

nonn

Seiichi Manyama, Table of n, a(n) for n = 0..10000
K. Ono, S. Robins and P. T. Wahl, On the representation of integers as sums of triangular numbers, Aequationes mathematicae, August 1995, Volume 50, Issue 1-2, pp 73-94.

Number of ways of writing n as a sum of k triangular numbers, for k=1,...: A010054, A008441, A008443, A008438, A008439, A008440, A226252, A007331, A226253, A226254, A226255, A014787, A014809.

G.f. is 7th power of g.f. for A010054.
a(0) = 1, a(n) = (7/n)*Sum_{k=1..n} A002129(k)*a(n-k) for n > 0. - Seiichi Manyama, May 06 2017
G.f.: exp(Sum_{k>=1} 7*(x^k/k)/(1 + x^k)). - Ilya Gutkovskiy, Jul 31 2017

Previous Showing 11-20 of 84 results. Next

cp's OEIS Frontend

A008441 Number of ways of writing n as the sum of 2 triangular numbers.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

Sage

Formula

Extensions

A321543 a(n) = Sum_{d|n} (-1)^(d-1)*d^2.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A008443 Number of ordered ways of writing n as the sum of 3 triangular numbers.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

Formula

Extensions

A007331 Fourier coefficients of E_{infinity,4}.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

Sage

Formula

Extensions

A181983 a(n) = (-1)^(n+1) * n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs