A033685 -id:A033685 - cp's OEIS frontend

A004016 Theta series of planar hexagonal lattice A_2.

1, 6, 0, 6, 6, 0, 0, 12, 0, 6, 0, 0, 6, 12, 0, 0, 6, 0, 0, 12, 0, 12, 0, 0, 0, 6, 0, 6, 12, 0, 0, 12, 0, 0, 0, 0, 6, 12, 0, 12, 0, 0, 0, 12, 0, 0, 0, 0, 6, 18, 0, 0, 12, 0, 0, 0, 0, 12, 0, 0, 0, 12, 0, 12, 6, 0, 0, 12, 0, 0, 0, 0, 0, 12, 0, 6, 12, 0, 0, 12, 0

Offset: 0

N. J. A. Sloane

The hexagonal lattice is the familiar 2-dimensional lattice in which each point has 6 neighbors. This is sometimes called the triangular lattice.
a(n) is the number of integer solutions to x^2 + x*y + y^2 = n (or equivalently x^2 - x*y + y^2 = n). - Michael Somos, Sep 20 2004
a(n) is the number of integer solutions to x^2 + y^2 + z^2 = 2*n where x + y + z = 0. - Michael Somos, Mar 12 2012
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
Cubic AGM theta functions: a(q) (the present sequence), b(q) (A005928), c(q) (A005882).
a(n) = 6*A002324(n) if n>0, and A002324 is multiplicative, thus a(1)*a(m*n) = a(n)*a(m) if n>0, m>0 are relatively prime. - Michael Somos, Mar 17 2019
The first occurrence of a(n)= 6, 12, 18, 24, ... (multiples of 6) is at n= 1, 7, 49, 91, 2401, 637, 117649, ... (see A002324). - R. J. Mathar, Sep 21 2024

			G.f. = 1 + 6*x + 6*x^3 + 6*x^4 + 12*x^7 + 6*x^9 + 6*x^12 + 12*x^13 + 6*x^16 + ...
Theta series of A_2 on the standard scale in which the minimal norm is 2:
1 + 6*q^2 + 6*q^6 + 6*q^8 + 12*q^14 + 6*q^18 + 6*q^24 + 12*q^26 + 6*q^32 + 12*q^38 + 12*q^42 + 6*q^50 + 6*q^54 + 12*q^56 + 12*q^62 + 6*q^72 + 12*q^74 + 12*q^78 + 12*q^86 + 6*q^96 + 18*q^98 + 12*q^104 + 12*q^114 + 12*q^122 + 12*q^126 + 6*q^128 + 12*q^134 + 12*q^146 + 6*q^150 + 12*q^152 + 12*q^158 + ...

B. C. Berndt, Ramanujan's Notebooks Part IV, Springer-Verlag, see p. 171, Entry 28.
Harvey Cohn, Advanced Number Theory, Dover Publications, Inc., 1980, p. 89. Ex. 18.
J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 111.
M. N. Huxley, Area, Lattice Points and Exponential Sums, Oxford, 1996; p. 236.
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..1000 from T. D. Noe)
S. Ahlgren, The sixth, eighth, ninth and tenth powers of Ramanujan's theta function, Proc. Amer. Math. Soc., 128 (1999), 1333-1338; F_3(q).
J. M. Borwein and P. B. Borwein, A cubic counterpart of Jacobi's identity and the AGM, Trans. Amer. Math. Soc., 323 (1991), no. 2, 691-701. MR1010408 (91e:33012) see page 695.
G. L. Hall, Comment on the paper "Theta series and magic numbers for diamond and certain ionic crystal structures" [J. Math. Phys. 28, 1653 (1987)]. Journal of Mathematical Physics; Sep. 1988, Vol. 29 Issue 9, pp. 2090-2092. - From _N. J. A. Sloane_, Dec 18 2012
M. D. Hirschhorn, Three classical results on representations of a number, Séminaire Lotharingien de Combinatoire, B42f (1999), 8 pp.
Christian Kassel and Christophe Reutenauer, The zeta function of the Hilbert scheme of n points on a two-dimensional torus, arXiv:1505.07229v3 [math.AG], 2015. [Note that a later version of this paper has a different title and different contents, and the number-theoretical part of the paper was moved to the publication which is next in this list.]
Christian Kassel and Christophe Reutenauer, Complete determination of the zeta function of the Hilbert scheme of n points on a two-dimensional torus, arXiv:1610.07793 [math.NT], 2016.
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.]
G. Nebe and N. J. A. Sloane, Home page for hexagonal (or triangular) lattice A2.
Helmut Ruhland, A family of lattices with an unbounded number of unit vectors, arXiv:2410.16172 [math.MG], 2024. See p. 2.
N. J. A. Sloane et al., Binary Quadratic Forms and OEIS (Index to related sequences, programs, references).
N. J. A. Sloane, Tables of Sphere Packings and Spherical Codes, IEEE Trans. Information Theory, vol. IT-27, 1981 pp. 327-338.
N. J. A. Sloane, Theta series and magic numbers for diamond and certain ionic crystal structures, J. Math. Phys. 28 (1987), 1653-1657.
Michael Somos, Introduction to Ramanujan theta functions.
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions.
Index entries for sequences related to A2 = hexagonal = triangular lattice.

Cf. A002324, A003051, A003215, A005881, A005882, A005928, A008458, A033685, A033687, A038587-A038591, etc.
See also A035019.
Cf. A000007, A000122, A004015, A008444, A008445, A008446, A008447, A008448, A008449 (Theta series of lattices A_0, A_1, A_3, A_4, ...), A186706.

Basis( ModularForms( Gamma1(3), 1), 81) [1]; /* Michael Somos, May 27 2014 */

L := Lattice("A",2); A := ThetaSeries(L, 161); A; /* Michael Somos, Nov 13 2014 */

A004016 := proc(n)
    local a,j ;
    a := A033716(n) ;
    for j from 0 to n/3 do
        a := a+A089800(n-1-3*j)*A089800(j) ;
    end do:
    a;
end proc:
seq(A004016(n),n=0..49) ; # R. J. Mathar, Feb 22 2021

a[ n_] := If[ n < 1, Boole[ n == 0 ], 6 DivisorSum[ n, KroneckerSymbol[ #, 3] &]]; (* Michael Somos, Nov 08 2011 *)
a[ n_] := SeriesCoefficient[ (QPochhammer[ q]^3 + 9 q QPochhammer[ q^9]^3) / QPochhammer[ q^3], {q, 0, n}]; (* Michael Somos, Nov 13 2014 *)
a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q] EllipticTheta[ 3, 0, q^3] + EllipticTheta[ 2, 0, q] EllipticTheta[ 2, 0, q^3], {q, 0, n}]; (* Michael Somos, Nov 13 2014 *)
a[ n_] := Length @ FindInstance[ x^2 + x y + y^2 == n, {x, y}, Integers, 10^9]; (* Michael Somos, Sep 14 2015 *)
terms = 81; f[q_] = LatticeData["A2", "ThetaSeriesFunction"][-I Log[q]/Pi]; s = Series[f[q], {q, 0, 2 terms}]; CoefficientList[s, q^2][[1 ;; terms]] (* Jean-François Alcover, Jul 04 2017 *)

{a(n) = my(A, p, e); if( n<1, n==0, A = factor(n); 6 * prod( k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==3, 1, p%3==1, e+1, 1-e%2)))}; /* Michael Somos, May 20 2005 */ /* Editor's note: this is the most efficient program */

{a(n) = if( n<0, 0, polcoeff( 1 + 6 * sum( k=1,n, x^k / (1 + x^k + x^(2*k)), x * O(x^n)), n))}; /* Michael Somos, Oct 06 2003 */

{a(n) = if( n<1, n==0, 6 * sumdiv( n,d, kronecker( d, 3)))}; /* Michael Somos, Mar 16 2005 */

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

{a(n) = my(A); if( n<0, 0, n*=3; A = x * O(x^n); polcoeff( (eta(x + A)^3  + 3 * x * eta(x^9 + A)^3) / eta(x^3 + A), n))}; /* Michael Somos, May 20 2005 */

{a(n) = if( n<1, n==0, qfrep([ 2, 1; 1, 2], n, 1)[n] * 2)}; /* Michael Somos, Jul 16 2005 */

{a(n) = if( n<0, 0, polcoeff( 1 + 6 * sum( k=1, n, x^(3*k - 2) / (1 - x^(3*k - 2)) - x^(3*k - 1) / (1 - x^(3*k - 1)), x * O(x^n)), n))} /* Paul D. Hanna, Jul 03 2011 */

from math import prod
from sympy import factorint
def A004016(n): return 6*prod(e+1 if p%3==1 else int(not e&1) for p, e in factorint(n).items() if p != 3) if n else 1 # Chai Wah Wu, Nov 17 2022

ModularForms( Gamma1(3), 1, prec=81).0 ; # Michael Somos, Jun 04 2013

Expansion of a(q) in powers of q where a(q) is the first cubic AGM theta function.
Expansion of theta_3(q) * theta_3(q^3) + theta_2(q) * theta_2(q^3) in powers of q.
Expansion of phi(x) * phi(x^3) + 4 * x * psi(x^2) * psi(x^6) in powers of x where phi(), psi() are Ramanujan theta functions.
Expansion of (1 / Pi) integral_{0 .. Pi/2} theta_3(z, q)^3 + theta_4(z, q)^3 dz in powers of q^2. - Michael Somos, Jan 01 2012
Expansion of coefficient of x^0 in f(x * q, q / x)^3 in powers of q^2 where f(,) is Ramanujan's general theta function. - Michael Somos, Jan 01 2012
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = u^2 - 3*v^2 - 2*u*w + 4*w^2. - Michael Somos, Jun 11 2004
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = (u1-u3) * (u3-u6) - (u2-u6)^2. - Michael Somos, May 20 2005
G.f. is a period 1 Fourier series which satisfies f(-1 / (3 t)) = 3^(1/2) (t/i) f(t) where q = exp(2 Pi i t). - Michael Somos, Sep 11 2007
G.f. A(x) satisfies A(x) + A(-x) = 2 * A(x^4), from Ramanujan.
G.f.: 1 + 6 * Sum_{k>0} x^k / (1 + x^k + x^(2*k)). - Michael Somos, Oct 06 2003
G.f.: Sum_( q^(n^2+n*m+m^2) ) where the sum (for n and m) extends over the integers. - Joerg Arndt, Jul 20 2011
G.f.: theta_3(q) * theta_3(q^3) + theta_2(q) * theta_2(q^3) = (eta(q^(1/3))^3 + 3 * eta(q^3)^3) / eta(q).
G.f.: 1 + 6*Sum_{n>=1} x^(3*n-2)/(1-x^(3*n-2)) - x^(3*n-1)/(1-x^(3*n-1)). - Paul D. Hanna, Jul 03 2011
a(3*n + 2) = 0, a(3*n) = a(n), a(3*n + 1) = 6 * A033687(n). - Michael Somos, Jul 16 2005
a(2*n + 1) = 6 * A033762(n), a(4*n + 2) = 0, a(4*n) = a(n), a(4*n + 1) = 6 * A112604(n), a(4*n + 3) = 6 * A112595(n). - Michael Somos, May 17 2013
a(n) = 6 * A002324(n) if n>0. a(n) = A005928(3*n).
Euler transform of A192733. - Michael Somos, Mar 12 2012
a(n) = (-1)^n * A180318(n). - Michael Somos, Sep 14 2015
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2*Pi/sqrt(3) = 3.627598... (A186706). - Amiram Eldar, Oct 15 2022

A005882 Theta series of planar hexagonal lattice (A2) with respect to deep hole.

Original entry on oeis.org

3, 3, 6, 0, 6, 3, 6, 0, 3, 6, 6, 0, 6, 0, 6, 0, 9, 6, 0, 0, 6, 3, 6, 0, 6, 6, 6, 0, 0, 0, 12, 0, 6, 3, 6, 0, 6, 6, 0, 0, 3, 6, 6, 0, 12, 0, 6, 0, 0, 6, 6, 0, 6, 0, 6, 0, 9, 6, 6, 0, 6, 0, 0, 0, 6, 9, 6, 0, 0, 6, 6, 0, 12, 0, 6, 0, 6, 0, 0, 0, 6, 6, 12, 0, 0, 3, 12, 0, 0, 6, 6, 0, 6, 0, 6, 0, 3, 6, 0, 0, 12

Offset: 0

N. J. A. Sloane

nonn

The hexagonal lattice is the familiar 2-dimensional lattice in which each point has 6 neighbors. This is sometimes called the triangular lattice.
Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
On page 111 of Conway and Sloane is "If the origin is moved to a deep hole the theta series is Theta_{hex+[1]}(z) = theta_2(z) psi_6(3z) + theta_3(z) psi_3(3z) = 3 q^{1/3} + 3 q^{4/3} + 6 q^{7/3} + 6 q^{13/3} + ... (63)" where the psi_k() for integer k is defined on page 103 equation (11) as psi_k(z) = e^{Pi i/z^2} theta_3(Pi z/k | z) = Sum_{m in Z} q^{(m + 1/k)^2}. - Michael Somos, Sep 10 2018

			G.f. = 3 + 3*x + 6*x^2 + 6*x^4 + 3*x^5 + 6*x^6 + 3*x^8 + 6*x^9 + 6*x^10 + ...
G.f. = 3*q + 3*q^4 + 6*q^7 + 6*q^13 + 3*q^16 + 6*q^19 + 3*q^25 + 6*q^28 + ...

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 111.
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..1000 from T. D. Noe)
J. M. Borwein and P. B. Borwein, A cubic counterpart of Jacobi's identity and the AGM, Trans. Amer. Math. Soc., 323 (1991), no. 2, 691-701. MR1010408 (91e:33012) see page 695.
Christian Kassel and Christophe Reutenauer, The zeta function of the Hilbert scheme of n points on a two-dimensional torus, arXiv:1505.07229v3 [math.AG], 2015. [Note that a later version of this paper has a different title and different contents, and the number-theoretical part of the paper was moved to the publication which is next in this list.]
Christian Kassel and Christophe Reutenauer, Complete determination of the zeta function of the Hilbert scheme of n points on a two-dimensional torus, arXiv:1610.07793 [math.NT], 2016.
G. Nebe and N. J. A. Sloane, Home page for hexagonal (or triangular) lattice A2
N. J. A. Sloane and B. K. Teo, Theta series and magic numbers for close-packed spherical clusters, J. Chem. Phys. 83 (1985) 6520-6534.
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions

Essentially same as A033685 and A033687.

Cf. A113062, A273845.

Basis( ModularForms( Gamma1(9), 1), 302)[2] * 3; /* Michael Somos, Jul 19 2014 */

a[ n_] := SeriesCoefficient[ 3 QPochhammer[ q^3]^3 / QPochhammer[ q], {q, 0, n}]; (* Michael Somos, Jul 19 2014 *)

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

Expansion of q^(-1/3) * 3 * eta(q^3)^3 / eta(q) in powers of q.
Expansion of q^(-1/3) * c(q) in powers of q where c(q) is the third cubic AGM theta function.
Given g.f. A(x), then B(x) = x*A(x^3) satisfies 0 = f(B(x), B(x^2), B(x^4)) where f(u, v, w) = v^3 + 2*u*w^2 - u^2*w. - Michael Somos, Aug 15 2006
G.f.: 3 Product_{k>0} (1-q^(3k))^3/(1-q^k).
G.f.: Sum_{u,v in Z} x^(u*u + u*v + v*v + u + v). - Michael Somos, Jul 19 2014
a(n) = 3 * A033687(n). a(n) = A113062(3*n + 1) = A033685(3*n + 1).
Expansion of 2 * psi(x^2) * f(x^2, x^4) + phi(x) * f(x^1, x^5) in powers of x where phi(), psi() are Ramanujan theta functions and f(, ) is Ramanujan's general theta function. - Michael Somos, Sep 07 2018
Sum_{k=1..n} a(k) ~ 2*Pi*n/sqrt(3). - Vaclav Kotesovec, Dec 17 2022

A033687 Theta series of hexagonal lattice A_2 with respect to deep hole divided by 3.

Original entry on oeis.org

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

Offset: 0

N. J. A. Sloane

nonn

The hexagonal lattice is the familiar 2-dimensional lattice in which each point has 6 neighbors. This is sometimes called the triangular lattice.
a(n)=0 if and only if A000731(n)=0 (see the Han-Ono paper). - Emeric Deutsch, May 16 2008
Number of 3-core partitions of n (denoted c_3(n) in Granville and Ono, p. 340). - Brian Hopkins, May 13 2008
Denoted by g_1(q) in Cynk and Hulek in Remark 3.4 on page 12 (but not explicitly listed).
This is a member of an infinite family of integer weight modular forms. g_1 = A033687, g_2 = A030206, g_3 = A130539, g_4 = A000731. - Michael Somos, Aug 24 2012
Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).
Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).

			G.f. = 1 + x + 2*x^2 + 2*x^4 + x^5 + 2*x^6 + x^8 + 2*x^9 + 2*x^10 + 2*x^12 + 2*x^14 + ...
G.f. = q + q^4 + 2*q^7 + 2*q^13 + q^16 + 2*q^19 + q^25 + 2*q^28 + 2*q^31 + 2*q^37 + ...

J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 111.
N. J. Fine, Basic Hypergeometric Series and Applications, Amer. Math. Soc., 1988; p. 79, Eq. (32.35) and (32.351).

Seiichi Manyama, Table of n, a(n) for n = 0..10000 (terms 0..1000 from T. D. Noe)
Johnathan M. Borwein and Peter B. Borwein, A cubic counterpart of Jacobi's identity and the AGM, Trans. Amer. Math. Soc., 323 (1991), no. 2, 691-701. MR1010408 (91e:33012) see page 697.
Slawomir Cynk and Klaus Hulek, Construction and examples of higher-dimensional modular Calabi-Yau manifolds, arXiv:math/0509424 [math.AG], 2005-2006.
Andrew Granville and Ken Ono, Defect Zero p-blocks for Finite Simple Groups, Transactions of the American Mathematical Society, Vol. 348 (1996), pp. 331-347.
Guo-Niu Han and Ken Ono, Hook Lengths and 3-Cores, Ann. Comb. (2011) Vol. 15, 305-312. See also arXiv:0805.2461, [math.NT], 2008.
Jeremy Lovejoy and Olivier Mallet, n-color overpartitions, twisted divisor functions, and Rogers-Ramanujan identities, South East Asian J. Math. Math. Sci., 6 (2008), 23-36. [From _Jeremy Lovejoy_, Jun 12 2009]
Gabriele Nebe and N. J. A. Sloane, Home page for hexagonal (or triangular) lattice A2
Michael Somos, Introduction to Ramanujan theta functions
Eric Weisstein's World of Mathematics, Ramanujan Theta Functions

Cf. A000731, A002324, A005882, A006938, A033685.

Cf. A045831, A053723, A081622.

Basis( ModularForms( Gamma1(9), 1), 316) [2]; /* Michael Somos, May 06 2015 */

a[ n_] := If[ n < 0, 0, DivisorSum[ 3 n + 1, KroneckerSymbol[ -3, #] &]]; (* Michael Somos, Sep 23 2013 *)
a[ n_] := SeriesCoefficient[ QPochhammer[ x^3]^3 / QPochhammer[ x], {x, 0, n}]; (* Michael Somos, Sep 01 2015 *)

{a(n) = if( n<0, 0, sumdiv( 3*n + 1, d, kronecker( -3, d)))}; /* Michael Somos, Nov 03 2005 */

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

{a(n) = my(A, p, e); if( n<0, 0, A = factor( 3*n + 1); prod( k=1, matsize(A)[1], [p, e] = A[k, ]; if( p==3, 0, p%6==1, e+1, 1-e%2)))}; /* Michael Somos, May 06 2015 */

Euler transform of period 3 sequence [1, 1, -2, ...].
Expansion of q^(-1/3) * eta(q^3)^3 / eta(q) in powers of q.
a(4*n + 1) = a(n). - Michael Somos, Dec 06 2004
a(n) = b(3*n + 1) where b(n) is multiplicative and b(p^e) = 0^e if p = 3, b(p^e) = e+1 if p == 1 (mod 6), b(p^e) = (1+(-1)^e)/2 if p == 2, 5 (mod 6). - Michael Somos, May 20 2005
Given g.f. A(x), B(q) = q * A(q^3) satisfies 0 = f(B(q), B(q^2), B(q^4)) where f(u, v, w) = u^2*w - 2*u*w^2 - v^3. - Michael Somos, Dec 06 2004
Given g.f. A(x), B(q)= q * A(q^3) satisfies 0 = f(B(q), B(q^2), B(q^3), B(q^6)) where f(u1, u2, u3, u6) = u1*u3^2 + u1*u6^2 - u1*u3*u6 - u2^2*u3. - Michael Somos, May 20 2005
Given g.f. A(x), B(q) = q * A(q^3) satisfies 0 = f(B(q), B(q^2), B(q^3), B(q^6)) where f(u1, u2, u3, u6) = u2*u3^2 + 2*u2*u3*u6 + 4*u2*u6^2 - u1^2*u6. - Michael Somos, May 20 2005
G.f.: Product_{k>0} (1 - x^(3*k))^3 / (1 - x^k).
G.f.: Sum_{k in Z} x^k / (1 - x^(3*k + 1)) = Sum_{k in Z} x^k / (1 - x^(6*k + 2)). - Michael Somos, Nov 03 2005
Expansion of q^(-1) * c(q^3) / 3 = q^(-1) * (a(q) - b(q)) / 9 in powers of q^3 where a(), b(), c() are cubic AGM theta functions. - Michael Somos, Dec 25 2007
G.f. is a period 1 Fourier series which satisfies f(-1 / (27 t)) = 3^(1/2) (t/i) g(t) where q = exp(2 Pi i t) and g(t) is the g.f. for A005928.
a(n) = Sum_{d|3n+1} LegendreSymbol{d,3} - Brian Hopkins, May 13 2008
q-series for a(n): Sum_{n >= 0} q^(n^2+n)(1-q)(1-q^2)...(1-q^n)/((1-q^(n+1))(1-q^(n+2))...(1-q^(2n+1))). [From Jeremy Lovejoy, Jun 12 2009]
a(n) = A002324(3*n + 1). 3*a(n) = A005882(n) = A033685(3*n + 1). - Michael Somos, Apr 04 2003
G.f.: (2 * psi(x^2) * f(x^2, x^4) + phi(x) * f(x^1, x^5)) / 3 where phi(), psi() are Ramanujan theta functions and f(, ) is Ramanujan's general theta function. - Michael Somos, Sep 07 2018
Sum_{k=1..n} a(k) ~ 2*Pi*n/3^(3/2). - Vaclav Kotesovec, Dec 17 2022

A202822 Numbers of the form 3*(x^2 + xy + y^2 + x + y) + 1 where x and y are integers.

Original entry on oeis.org

1, 4, 7, 13, 16, 19, 25, 28, 31, 37, 43, 49, 52, 61, 64, 67, 73, 76, 79, 91, 97, 100, 103, 109, 112, 121, 124, 127, 133, 139, 148, 151, 157, 163, 169, 172, 175, 181, 193, 196, 199, 208, 211, 217, 223, 229, 241, 244, 247, 256, 259, 268, 271, 277, 283, 289, 292

Offset: 1

Michael Somos, Dec 25 2011

nonn

Closed under multiplication.

Löschian numbers of the form 3*k+1. - Altug Alkan, Nov 18 2015

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Joerg Arndt, Plane-filling curves on all uniform grids, arXiv preprint arXiv:1607.02433 [math.CO], 2016.

Cf. A033685, A045833.

Subsequence of A003136, A260682 (subsequence).

a202822 n = a202822_list !! (n-1)
a202822_list = filter ((== 1) . flip mod 3) a003136_list
-- Reinhard Zumkeller, Nov 16 2015

nf[{i_,j_}]:=3(i^2+i*j+j^2+i+j)+1; Union[nf/@Tuples[Range[-10,10],2]] (* Harvey P. Dale, Dec 31 2011 *)

isA(n) = if(n%3 == 0, 0, 0 != sumdiv( n, d, kronecker( -3, d)))

x='x+O('x^500); p=eta(x)^3/eta(x^3); for(n=0, 499, if(polcoeff(p, n) != 0 && n%3==1, print1(n, ", "))) \\ Altug Alkan, Nov 18 2015

is(n)=(n%3==1) && #bnfisintnorm(bnfinit(z^2+z+1), n); \\ Joerg Arndt, Jan 04 2016

list(lim)=my(v=List(), y, t); for(x=0, sqrtint(lim\3), my(y=x, t); while((t=x^2+x*y+y^2)<=lim, if((x-y)%3, listput(v, t)); y++)); Set(v) \\ Charles R Greathouse IV, Jul 05 2017

A033685(n) != 0 if and only if n is in the set.

A113062 Expansion of theta series of hexagonal net with respect to a node.

Original entry on oeis.org

1, 3, 0, 6, 3, 0, 0, 6, 0, 6, 0, 0, 6, 6, 0, 0, 3, 0, 0, 6, 0, 12, 0, 0, 0, 3, 0, 6, 6, 0, 0, 6, 0, 0, 0, 0, 6, 6, 0, 12, 0, 0, 0, 6, 0, 0, 0, 0, 6, 9, 0, 0, 6, 0, 0, 0, 0, 12, 0, 0, 0, 6, 0, 12, 3, 0, 0, 6, 0, 0, 0, 0, 0, 6, 0, 6, 6, 0, 0, 6, 0, 6, 0, 0, 12, 0, 0, 0, 0, 0, 0, 12, 0, 12, 0, 0, 0, 6, 0, 0

Offset: 0

Michael Somos, Oct 13 2005

The hexagonal net is the familiar 2-dimensional honeycomb (not a lattice) in which each node has 3 neighbors.

Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).

			G.f. = 1 + 3*q + 6*q^3 + 3*q^4 + 6*q^7 + 6*q^9 + 6*q^12 + 6*q^13 + 3*q^16 + ...

A. F. Wells, Structural Inorganic Chemistry, Oxford, 5th ed., 1984; see Fig. 3.9(a.1).

T. D. Noe, Table of n, a(n) for n = 0..1000
George L. Hall, Comment on the paper "Theta series and magic numbers for diamond and certain ionic crystal structures" [J. Math. Phys. 28, 1653 (1987)]. Journal of Mathematical Physics; Sep. 1988, Vol. 29 Issue 9, pp. 2090-2092. - From _N. J. A. Sloane_, Dec 18 2012
N. J. A. Sloane, Theta series and magic numbers for diamond and certain ionic crystal structures, J. Math. Phys. 28 (1987), 1653-1657.

Cf. A004016, A005882, A005928, A033685, A113063, A275486.

Cf. A000122, A000700, A004016, A005882, A005928, A010054, A121373.

a[0] = 1; a[n_] := 3*DivisorSum[n, {0, 1, -1, 1, 1, -1, -1, 1, -1}[[Mod[#, 9]+1]]&]; Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Nov 04 2015, after 1st PARI script *)

{a(n) = if( n<1, n==0, 3 * sumdiv(n, d, [ 0, 1, -1, 1, 1, -1, -1, 1, -1][d%9+1]))};

{a(n) = local(A, p, e); if( n<1, n==0, A = factor(n); 3 * prod(k=1, matsize(A)[1], if(p=A[k,1], e=A[k,2]; if(p==3, 2, if(p%6==1, e+1, !(e%2))))))};

{a(n) = local(A); if( n<0, 0, A = x * O(x^n); abs( polcoeff( eta(x + A)^3 / eta(x^3 + A), n)))};

{a(n) = local(A); if( n<0, 0, A = x * O(x^n); polcoeff( (eta(x + A)^3 + 6 * x * eta(x^9 + A)^3) / eta(x^3 + A), n))}; /* Michael Somos, Aug 15 2006 */

Moebius transform is period 9 sequence [ 3, -3, 3, 3, -3, -3, 3, -3, 0, ...].
Expansion of a(q^3) + c(q^3) in powers of q where a(), c() are cubic AGM theta functions. - Michael Somos, Aug 15 2006
For n>0, a(n) = 3*b(n) where b(n)=A113063(n) is multiplicative and b(p^e) = 2 if p = 3 and e>0, b(p^e) = e+1 if p == 1 (mod 6), b(p^e) = (1+(-1)^e)/2 if p == 2, 5 (mod 6).
a(3*n + 2) = 0. a(3*n + 1) = A005882(n) = A033685(3*n + 1) = -A005928(3*n + 1). a(3*n) = A004016(n) = A005928(3*n).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 4*Pi/(3*sqrt(3)) = 2.418399... (A275486). - Amiram Eldar, Dec 28 2023

Definition corrected Michael Somos, Oct 17 2005

A045833 Expansion of eta(q^9)^3 / eta(q^3) in powers of q.

Original entry on oeis.org

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

Offset: 0

N. J. A. Sloane

			G.f. = q + q^4 + 2*q^7 + 2*q^13 + q^16 + 2*q^19 + q^25 + 2*q^28 + 2*q^31 + ...

Amiram Eldar, Table of n, a(n) for n = 0..10000

Cf. A033685, A033687, A097195, A217219.

a[ n_] := SeriesCoefficient[ q QPochhammer[ q^9]^3 / QPochhammer[ q^3], {q, 0, n}]; (* Michael Somos, Feb 22 2015 *)
f[p_, e_] := If[Mod[p, 3] == 1, e + 1, (1 + (-1)^e)/2]; f[3, e_] := 0; a[0] = 0; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100, 0] (* Amiram Eldar, Oct 13 2022 *)

{a(n) = local(A, p, e); if( n<0, 0, A=factor(n); prod(k=1, matsize(A)[1], if( p=A[k,1], e=A[k,2]; if( p!=3, if( p%3==1, e+1, !(e%2))))))}; \\ Michael Somos, May 25 2005

{a(n) = local(A); if( (n<1) || (n%3!=1), 0, n = (n-1)/3; A = x * O(x^n); polcoeff( eta(x^3 + A)^3 / eta(x + A), n))}; \\ Michael Somos, May 25 2005

From Michael Somos, May 25 2005: (Start)
Euler transform of period 9 sequence [ 0, 0, 1, 0, 0, 1, 0, 0, -2, ...].
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = u^2*w - 2*u*w^2 - v^3.
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = u1*u3^2 + u1*u6^2 - u1*u3*u6 - u2^2*u3.
a(3*n) = a(3*n + 2) = 0. a(3*n + 1) = A033687(n). a(6*n + 1) = A097195(n). 3*a(n) = A033685(n).
Multiplicative with a(3^e) = 0^e, a(p^e) = e+1 if p == 1 (mod 3), a(p^e) = (1+(-1)^e)/2 if p == 2 (mod 3).
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = u2*u3^2 + 2*u2*u3*u6 + 4*u2*u6^2 - u1^2*u6. (End)
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 2*Pi/(9*sqrt(3)) = 0.403066... . - Amiram Eldar, Oct 13 2022
Dirichlet g.f.: L(chi_1,s)*L(chi_{-1},s), where chi_1 = A011655 and chi_{-1} = A102283 are respectively the principal and the non-principal Dirichlet character modulo 3. For the formula of the sequence whose Dirichlet g.f. is Product_{chi} L(chi,s), where chi runs through all Dirichlet characters modulo k, see A378006. This sequence is the case k = 3. - Jianing Song, Nov 13 2024

cp's OEIS Frontend

A004016 Theta series of planar hexagonal lattice A_2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

PARI

PARI

Python

Sage

Formula

A005882 Theta series of planar hexagonal lattice (A2) with respect to deep hole.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A033687 Theta series of hexagonal lattice A_2 with respect to deep hole divided by 3.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

PARI

Formula

A202822 Numbers of the form 3*(x^2 + xy + y^2 + x + y) + 1 where x and y are integers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

PARI

PARI

Formula

A113062 Expansion of theta series of hexagonal net with respect to a node.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs