A003050 -id:A003050 - cp's OEIS frontend

A001615 Dedekind psi function: n * Product_{p|n, p prime} (1 + 1/p).

1, 3, 4, 6, 6, 12, 8, 12, 12, 18, 12, 24, 14, 24, 24, 24, 18, 36, 20, 36, 32, 36, 24, 48, 30, 42, 36, 48, 30, 72, 32, 48, 48, 54, 48, 72, 38, 60, 56, 72, 42, 96, 44, 72, 72, 72, 48, 96, 56, 90, 72, 84, 54, 108, 72, 96, 80, 90, 60, 144, 62, 96, 96, 96, 84, 144, 68, 108, 96

Offset: 1

N. J. A. Sloane

Number of primitive sublattices of index n in generic 2-dimensional lattice; also index of Gamma_0(n) in SL_2(Z).
A generic 2-dimensional lattice L =  consists of all vectors of the form mV + nW, (m,n integers). A sublattice S =  has index |ad-bc| and is primitive if gcd(a,b,c,d) = 1. The generic lattice L has precisely a(2) = 3 sublattices of index 2, namely <2V,W>,  and  (which = ) and so on for other indices.
The sublattices of index n are in 1-to-1 correspondence with matrices [a b; 0 d] with a>0, ad=n, b in [0..d-1]. The number of these is Sum_{d|n} = sigma(n), which is A000203. A sublattice is primitive if gcd(a,b,d) = 1; the number of these is n * product_{p|n} (1+1/p), which is the present sequence.
SL_2(Z) = Gamma is the group of all 2 X 2 matrices [a b; c d] where a,b,c,d are integers with ad-bc = 1 and Gamma_0(N) is usually defined as the subgroup of this for which N|c. But conceptually Gamma is best thought of as the group of (positive) automorphisms of a lattice , its typical element taking V -> aV + bW, W -> cV + dW and then Gamma_0(N) can be defined as the subgroup consisting of the automorphisms that fix the sublattice  of index N. - J. H. Conway, May 05 2001
Dedekind proved that if n = k_i*j_i for i in I represents all the ways to write n as a product, and e_i=gcd(k_i,j_i), then a(n)= sum(k_i / (e_i * phi(e_i)), i in I ) [cf. Dickson, History of the Theory of Numbers, Vol. 1, p. 123].
Also a(n)= number of cyclic subgroups of order n in an Abelian group of order n^2 and type (1,1) (Fricke). - Len Smiley, Dec 04 2001
The polynomial degree of the classical modular equation of degree n relating j(z) and j(nz) is psi(n) (Fricke). - Michael Somos, Nov 10 2006; clarified by Katherine E. Stange, Mar 11 2022
The Mobius transform of this sequence is A063659. - Gary W. Adamson, May 23 2008
The inverse Mobius transform of this sequence is A060648. - Vladeta Jovovic, Apr 05 2009
The Dirichlet inverse of this sequence is A008836(n) * A048250(n). - Álvar Ibeas, Mar 18 2015
The Riemann Hypothesis is true if and only if a(n)/n - e^gamma*log(log(n)) < 0 for any n > 30. - Enrique Pérez Herrero, Jul 12 2011
The Riemann Hypothesis is also equivalent to another inequality, see the Sole and Planat link. - Thomas Ordowski, May 28 2017
An infinitary analog of this sequence is the sum of the infinitary divisors of n (see A049417). - Vladimir Shevelev, Apr 01 2014
Problem: are there composite numbers n such that n+1 divides psi(n)? - Thomas Ordowski, May 21 2017
The sum of divisors d of n such that n/d is squarefree. - Amiram Eldar, Jan 11 2019
Psi(n)/n is a new maximum for each primorial (A002110) [proof in link: Patrick Sole and Michel Planat, Proposition 1 page 2]. - Bernard Schott, May 21 2020
From Jianing Song, Nov 05 2022: (Start)
a(n) is the number of subgroups of C_n X C_n that are isomorphic to C_n, where C_n is the cyclic group of order n. Proof: the number of elements of order n in C_n X C_n is A007434(n) (they are the elements of the form (a,b) in C_n X C_n where gcd(a,b,n) = 1), and each subgroup isomorphic to C_n contains phi(n) generators, so the number of such subgroups is A007434(n)/phi(n) = a(n).
The total number of order-n subgroups of C_n X C_n is A000203(n). (End)

			Let L = <V,W> be a 2-dimensional lattice. The 6 primitive sublattices of index 4 are generated by <4V,W>, <V,4W>, <4V,W+-V>, <2V+W,2W>, <2V,2W+V>. Compare A000203.
G.f. = x + 3*x^2 + 4*x^3 + 6*x^4 + 6*x^5 + 12*x^6 + 8*x^7 + 12*x^8 + 12*x^9 + ...

Tom Apostol, Intro. to Analyt. Number Theory, page 71, Problem 11, where this is called phi_1(n).
David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989, p. 228.
R. Fricke, Die elliptischen Funktionen und ihre Anwendungen, Teubner, 1922, Vol. 2, see p. 220.
Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004. See Section B41, p. 147.
B. Schoeneberg, Elliptic Modular Functions, Springer-Verlag, NY, 1974, p. 79.
G. Shimura, Introduction to the Arithmetic Theory of Automorphic Functions, Princeton, 1971, see p. 25, Eq. (1).
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).

T. D. Noe and N. J. A. Sloane, Table of n, a(n) for n = 1..10000
O. Bordelles and B. Cloitre, An Alternating Sum Involving the Reciprocal of Certain Multiplicative Functions, J. Int. Seq. 16 (2013) #13.6.3.
Harriet Fell, Morris Newman, and Edward Ordman, Tables of genera of groups of linear fractional transformations, J. Res. Nat. Bur. Standards Sect. B 67B 1963 61-68.
M. Hampejs, N. Holighaus, L. Toth and C. Wiesmeyr, Representing and counting the subgroups of the group Z_m X Z_n, arXiv:1211.1797 [math.GR], 2012.
W. Hürlimann, Dedekind's arithmetic function and primitive four squares counting functions, Journal of Algebra, Number Theory: Advances and Applications 14:2 (2015), 73-88.
F. A. Lewis et al., Problem 4002, Amer. Math. Monthly, Vol. 49, No. 9, Nov. 1942, pp. 618-619.
E. Pérez Herrero, Recycling Hardy & Wright, Average Order of Dedekind Psi Function, Psychedelic Geometry Blogspot.
Michel Planat, Riemann hypothesis from the Dedekind psi function, arXiv:1010.3239 [math.GM], 2010.
Patrick Sole and Michel Planat, Extreme values of the Dedekind Psi function, to appear in Journal of Combinatorics and Number Theory, arXiv:1011.1825 [math.NT], 2010-2011.
Eric Weisstein's World of Mathematics, Dedekind Function
Wikipedia, Dedekind psi function
Index entries for "core" sequences
Index entries for sequences related to sublattices

Other sequences that count lattices/sublattices: A000203 (with primitive condition removed), A003050 (hexagonal lattice instead), A003051, A054345, A160889, A160891.
Cf. A002110, A027748, A124010, A019269.
Cf. A301594.
Cf. A063659 (Möbius transform), A082020 (average order), A156303 (Euler transform), A173290 (partial sums), A175836 (partial products), A203444 (range).
Cf. A210523 (record values).
Algebraic combinations with other core sequences: A000082, A033196, A175732, A291784, A344695.
Sequences of the form n^k * Product_ {p|n, p prime} (1 + 1/p^k) for k=0..10: A034444 (k=0), this sequence (k=1), A065958 (k=2), A065959 (k=3), A065960 (k=4), A351300 (k=5), A351301 (k=6), A351302 (k=7), A351303 (k=8), A351304 (k=9), A351305 (k=10).
Cf. A082695 (Dgf at s=3), A339925 (Dgf at s=4).

import Data.Ratio (numerator)
a001615 n = numerator (fromIntegral n * (product $
            map ((+ 1) . recip . fromIntegral) $ a027748_row n))
-- Reinhard Zumkeller, Jun 03 2013, Apr 12 2012

m:=75; R:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&+[MoebiusMu(k)^2*x^k/(1-x^k)^2: k in [1..2*m]]) )); // G. C. Greubel, Nov 23 2018

A001615 := proc(n) n*mul((1+1/i[1]),i=ifactors(n)[2]) end; # Mark van Hoeij, Apr 18 2012

Join[{1}, Table[n Times @@ (1 + 1/Transpose[FactorInteger[n]][[1]]), {n, 2, 100}]] (* T. D. Noe, Jun 11 2006 *)
Table[DirichletConvolve[j, MoebiusMu[j]^2, j, n], {n, 100}] (* Jan Mangaldan, Aug 22 2013 *)
a[n_] := n Sum[MoebiusMu[d]^2/d, {d, Divisors[n]}]; (* Michael Somos, Jan 10 2015 *)
Table[n Product[1 + 1/p, {p, Select[Divisors[n], PrimeQ]}], {n, 1, 100}] (* Vaclav Kotesovec, May 08 2021 *)
Table[n DivisorSum[n, MoebiusMu[#]^2/# &], {n, 20}] (* Eric W. Weisstein, Mar 09 2025 *)

{a(n) = if( n<1, 0, direuler( p=2, n, (1 + X) / (1 - p*X)) [n])};

{a(n) = if( n<1, 0, n * sumdiv( n, d, moebius(d)^2 / d))}; /* Michael Somos, Nov 10 2006 */

a(n)=my(f=factor(n)); prod(i=1,#f~, f[i,1]^f[i,2] + f[i,1]^(f[i,2]-1)) \\ Charles R Greathouse IV, Aug 22 2013

a(n) = n * sumdivmult(n, d, issquarefree(d)/d) \\ Charles R Greathouse IV, Sep 09 2014

from math import prod
from sympy import primefactors
def A001615(n):
    plist = primefactors(n)
    return n*prod(p+1 for p in plist)//prod(plist) # Chai Wah Wu, Jun 03 2021

def A001615(n) : return n*mul(1+1/p for p in prime_divisors(n))
[A001615(n) for n in (1..69)] # Peter Luschny, Jun 10 2012

Dirichlet g.f.: zeta(s) * zeta(s-1) / zeta(2*s). - Michael Somos, May 19 2000
Multiplicative with a(p^e) = (p+1)*p^(e-1). - David W. Wilson, Aug 01 2001
a(n) = A003557(n)*A048250(n) = n*A000203(A007947(n))/A007947(n). - Labos Elemer, Dec 04 2001
a(n) = n*Sum_{d|n} mu(d)^2/d, Dirichlet convolution of A008966 and A000027. - Benoit Cloitre, Apr 07 2002
a(n) = Sum_{d|n} mu(n/d)^2 * d. - Joerg Arndt, Jul 06 2011
From Enrique Pérez Herrero, Aug 22 2010: (Start)
a(n) = J_2(n)/J_1(n) = J_2(n)/phi(n) = A007434(n)/A000010(n), where J_k is the k-th Jordan Totient Function.
a(n) = (1/phi(n))*Sum_{d|n} mu(n/d)*d^(b-1), for b=3. (End)
a(n) = n / Sum_{d|n} mu(d)/a(d). - Enrique Pérez Herrero, Jun 06 2012
a(n^k)= n^(k-1) * a(n). - Enrique Pérez Herrero, Jan 05 2013
If n is squarefree, then a(n) = A049417(n) = A000203(n). - Vladimir Shevelev, Apr 01 2014
a(n) = Sum_{d^2 | n} mu(d) * A000203(n/d^2). - Álvar Ibeas, Dec 20 2014
The average order of a(n) is 15*n/Pi^2. - Enrique Pérez Herrero, Jan 14 2012. See Apostol. - N. J. A. Sloane, Sep 04 2017
G.f.: Sum_{k>=1} mu(k)^2*x^k/(1 - x^k)^2. - Ilya Gutkovskiy, Oct 25 2018
a(n) = Sum_{d|n} 2^omega(d) * phi(n/d), Dirichlet convolution of A034444 and A000010. - Daniel Suteu, Mar 09 2019
From Richard L. Ollerton, May 07 2021: (Start)
a(n) = Sum_{k=1..n} 2^omega(gcd(n,k)).
a(n) = Sum_{k=1..n} 2^omega(n/gcd(n,k))*phi(gcd(n,k))/phi(n/gcd(n,k)). (End)
a(n) = abs(A158523(n)) = A158523(n) * A008836(n). - Enrique Pérez Herrero, Nov 07 2022
a(n) = (1/n) * Sum_{d|n} mu(n/d)*sigma(d^2). - Ridouane Oudra, Mar 26 2025

More terms from Olivier Gérard, Aug 15 1997

A007997 a(n) = ceiling((n-3)(n-4)/6).

Original entry on oeis.org

0, 0, 1, 1, 2, 4, 5, 7, 10, 12, 15, 19, 22, 26, 31, 35, 40, 46, 51, 57, 64, 70, 77, 85, 92, 100, 109, 117, 126, 136, 145, 155, 166, 176, 187, 199, 210, 222, 235, 247, 260, 274, 287, 301, 316, 330, 345, 361, 376, 392, 409, 425, 442, 460, 477, 495, 514, 532, 551, 571, 590, 610

Offset: 3

N. J. A. Sloane

Number of solutions to x+y+z=0 (mod m) with 0<=x<=y<=z

Nonorientable genus of complete graph on n nodes.

Also (with different offset) Molien series for alternating group A_3.

(1+x^3 ) / ((1-x)*(1-x^2)*(1-x^3)) is the Poincaré series [or Poincare series] (or Molien series) for H^*(S_6, F_2).

a(n+5) is the number of necklaces with 3 black beads and n white beads.

The g.f./x^5 is Z(C_3,x), the 3-variate cycle index polynomial for the cyclic group C_3, with substitution x[i]->1/(1-x^i), i=1,2,3. Therefore by Polya enumeration a(n+5) is the number of cyclically inequivalent 3-necklaces whose 3 beads are labeled with nonnegative integers such that the sum of labels is n, for n=0,1,2,... . See A102190 for Z(C_3,x). - Wolfdieter Lang, Feb 15 2005

a(n+1) is the number of pairs (x,y) with x and y in {0,...,n}, x = (y mod 3), and x+y < n. - Clark Kimberling, Jul 02 2012

From Gus Wiseman, Oct 17 2020: (Start)

Also the number of 3-part integer compositions of n - 2 that are either weakly increasing or strictly decreasing. For example, the a(5) = 1 through a(13) = 15 compositions are:

(111) (112) (113) (114) (115) (116) (117) (118) (119)

(122) (123) (124) (125) (126) (127) (128)

(222) (133) (134) (135) (136) (137)

(321) (223) (224) (144) (145) (146)

(421) (233) (225) (226) (155)

(431) (234) (235) (227)

(521) (333) (244) (236)

(432) (334) (245)

(531) (532) (335)

(621) (541) (344)

(631) (542)

(721) (632)

(641)

(731)

(821)

(End)

			For m=7 (n=12), the 12 solutions are xyz = 000 610 520 511 430 421 331 322 662 653 644 554.

A. Adem and R. J. Milgram, Cohomology of Finite Groups, Springer-Verlag, 2nd. ed., 2004, p. 204.
D. J. Benson, Polynomial Invariants of Finite Groups, Cambridge, 1993, p. 105.
J. L. Gross and T. W. Tucker, Topological Graph Theory, Wiley, 1987; see \bar{I}(n) p. 221.
J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 740.
E. V. McLaughlin, Numbers of factorizations in non-unique factorial domains, Senior Thesis, Allegeny College, Meadville, PA, 2004.

T. D. Noe, Table of n, a(n) for n = 3..1000
Chwas Ahmed, Paul Martin, and Volodymyr Mazorchuk, On the number of principal ideals in d-tonal partition monoids, arXiv preprint arXiv:1503.06718 [math.CO], 2015-2019.
Alexander Taveira Blomenhofer, On Uniqueness of Power Sum Decomposition, SIAM J. Appl. Alg. Geom. (2025) Vol. 9, Iss. 1.
David Broadhurst and Xavier Roulleau, Number of partitions of modular integers, arXiv:2502.19523 [math.NT], 2025. See pp. 3, 19.
Magnus Rahbek Hansen, Invariant Theory and Graphs, Master's Thesis, Univ. Copenhagen (Denmark, 2023). See p. 38.
Mónica A. Reyes, Cristina Dalfó, Miguel Àngel Fiol, and Arnau Messegué, A general method to find the spectrum and eigenspaces of the k-token of a cycle, and 2-token through continuous fractions, arXiv:2403.20148 [math.CO], 2024. See p. 6.
Index entries for linear recurrences with constant coefficients, signature (2,-1,1,-2,1).
Index entries for Molien series.
Index entries for sequences related to necklaces.

Cf. A000031, A007998, A003050, A047996, A048259, A053618.
Apart from initial term, same as A058212.
Cf. A000212, A000217, A001840, A128422, A156040.
A001399(n-6)*2 = A069905(n-3)*2 = A211540(n-1)*2 counts the strict case.
A014311 intersected with A225620 U A333256 ranks these compositions.
A218004 counts these compositions of any length.
A000009 counts strictly decreasing compositions.
A000041 counts weakly increasing compositions.
A001523 counts unimodal compositions, with complement counted by A115981.
A007318 and A097805 count compositions by length.
A032020 counts strict compositions, ranked by A233564.
A333149 counts neither increasing nor decreasing strict compositions.
Cf. A014311, A332834, A333147, A337481, A337482-A337484.

a007997 n = ceiling $ (fromIntegral $ (n - 3) * (n - 4)) / 6
a007997_list = 0 : 0 : 1 : zipWith (+) a007997_list [1..]
-- Reinhard Zumkeller, Dec 18 2013

x^5*(1+x^3)/((1-x)*(1-x^2)*(1-x^3));
seq(ceil(binomial(n,2)/3), n=0..63); # Zerinvary Lajos, Jan 12 2009
a := n -> (n*(n-7)-2*([1,1,-1][n mod 3 +1]-7))/6;
seq(a(n), n=3..64); # Peter Luschny, Jan 13 2015

k = 3; Table[Apply[Plus, Map[EulerPhi[ # ]Binomial[n/#, k/# ] &, Divisors[GCD[n, k]]]]/n, {n, k, 30}] (* Robert A. Russell, Sep 27 2004 *)
Table[Ceiling[((n-3)(n-4))/6],{n,3,100}] (* or *) LinearRecurrence[ {2,-1,1,-2,1},{0,0,1,1,2},100] (* Harvey P. Dale, Jan 21 2014 *)

a(n)=(n^2-7*n+16)\6 \\ Charles R Greathouse IV, Sep 24 2015

a(n) = a(n-3) + n - 2, a(0)=0, a(1)=0, a(2)=1 [Offset 0]. - Paul Barry, Jul 14 2004
G.f.: x^5*(1+x^3)/((1-x)*(1-x^2)*(1-x^3)) = x^5*(1-x+x^2)/((1-x)^2*(1-x^3)).
a(n+5) = Sum_{k=0..floor(n/2)} C(n-k,L(k/3)), where L(j/p) is the Legendre symbol of j and p. - Paul Barry, Mar 16 2006
a(3)=0, a(4)=0, a(5)=1, a(6)=1, a(7)=2, a(n) = 2*a(n-1) - a(n-2) + a(n-3) - 2*a(n-4) + a(n-5). - Harvey P. Dale, Jan 21 2014
a(n) = (n^2 - 7*n + 14 - 2*(-1)^(2^(n + 1 - 3*floor((n+1)/3))))/6. - Luce ETIENNE, Dec 27 2014
a(n) = A001399(n-3) + A001399(n-6). Compare to A140106(n) = A001399(n-3) - A001399(n-6). - Gus Wiseman, Oct 17 2020
a(n) = (40 + 3*(n - 7)*n - 4*cos(2*n*Pi/3) - 4*sqrt(3)*sin(2*n*Pi/3))/18. - Stefano Spezia, Dec 14 2021
Sum_{n>=5} 1/a(n) = 6 - 2*Pi/sqrt(3) + 2*Pi*tanh(sqrt(5/3)*Pi/2)/sqrt(15). - Amiram Eldar, Oct 01 2022

A003051 Number of inequivalent sublattices of index n in hexagonal lattice, where two sublattices are equivalent if they are related by a rotation or reflection preserving the hexagonal lattice.

Original entry on oeis.org

1, 1, 2, 3, 2, 3, 3, 5, 4, 4, 3, 8, 4, 5, 6, 9, 4, 8, 5, 10, 8, 7, 5, 15, 7, 8, 9, 13, 6, 14, 7, 15, 10, 10, 10, 20, 8, 11, 12, 20, 8, 18, 9, 17, 16, 13, 9, 28, 12, 17, 14, 20, 10, 22, 14, 25, 16, 16, 11, 34, 12, 17, 21, 27, 16, 26, 13, 24, 18, 26, 13, 40, 14

Offset: 1

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.
From Andrey Zabolotskiy, Mar 10 2018: (Start)
If only primitive sublattices are considered, we get A003050.
Here only rotations and reflections preserving the parent hexagonal lattice are allowed. If reflections are not allowed, we get A145394. If any rotations and reflections are allowed, we get A300651.
In other words, the parent lattice of the sublattices under consideration has Patterson symmetry group p6mm, and two sublattices are considered equivalent if they are related via a symmetry from that group [Rutherford]. For other 2D Patterson groups, the analogous sequences are A000203 (p2), A069734 (p2mm), A145391 (c2mm), A145392 (p4), A145393 (p4mm), A145394 (p6).
Rutherford says at p. 161 that his sequence for p6mm differs from this sequence, but it seems that with the current definition and terms of this sequence, this actually is his p6mm sequence, and the sequence he thought to be this one is actually A300651. Also, he says that a(n) != A300651(n) only when A002324(n) > 2 (first time happens at n = 49), but actually these two sequences differ at other terms, too, for example, at n = 42 (see illustration). (End)

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
A. Altshuler, Construction and enumeration of regular maps on the torus, Discrete Math. 4 (1973), 201-217.
A. Altshuler, Construction and enumeration of regular maps on the torus, Discrete Math. 4 (1973), 201-217. [Annotated and corrected scanned copy]
M. Bernstein, N. J. A. Sloane and P. E. Wright, On Sublattices of the Hexagonal Lattice, Discrete Math. 170 (1997) 29-39 (Abstract, pdf, ps).
Amihay Hanany, Domenico Orlando, and Susanne Reffert, Sublattice counting and orbifolds, High Energ. Phys., 2010 (2010), 51, arXiv.org:1002.2981 [hep-th] (see Table 3).
Daejun Kim, Seok Hyeong Lee, and Seungjai Lee, Zeta functions enumerating subforms of quadratic forms, arXiv:2409.05625 [math.NT], 2024.
W. Kurth, Enumeration of Platonic maps on the torus, Discrete Math. 61 (1986), 71-83.
G. Nebe and N. J. A. Sloane, Home page for hexagonal (or triangular) lattice A2
John S. Rutherford, Sublattice enumeration. IV. Equivalence classes of plane sublattices by parent Patterson symmetry and colour lattice group type, Acta Cryst. (2009). A65, 156-163. [See Table 2].
Andrey Zabolotskiy, Sublattices of the hexagonal lattice (illustrations for n = 1..7, 14)
Andrey Zabolotskiy, Coweight lattice A^*_n and lattice simplices, arXiv:2003.10251 [math.CO], 2020.
Index entries for sequences related to sublattices
Index entries for sequences related to A2 = hexagonal = triangular lattice

Cf. A003050, A054384, A001615, A006984, A054345, A054346, A000203, A069734, A145391, A145392, A145393, A145394, A112689, A159842.

max = 73; A145390 = Drop[ CoefficientList[ Series[ Sum[(1 + Cos[n*Pi/2])*x^n/(1 - x^n), {n, 1, max}], {x, 0, max}], x], 1]; A002324[n_] := (dn = Divisors[n]; Count[dn, ?(Mod[#, 3] == 1 & )] - Count[dn, ?(Mod[#, 3] == 2 & )]); a[n_] := (DivisorSigma[1, n] + 2 A002324[n] + 3*A145390[[n]])/6; Table[a[n], {n, 1, max}] (* Jean-François Alcover, Oct 11 2011, after given formula *)

a(n) = Sum_{ m^2 | n } A003050(n/m^2).
a(n) = (A000203(n) + 2*A002324(n) + 3*A145390(n))/6. [Rutherford] - N. J. A. Sloane, Mar 13 2009
a(n) = Sum_{ d|n } A112689(d+1). - Andrey Zabolotskiy, Aug 29 2019
a(n) = Sum_{ d|n } floor(d/6) + 1 - 1*[d == 2 or 6 (mod 12)] + 1*[d == 4 (mod 12)]. [Kurth] - Brahadeesh Sankarnarayanan, Feb 24 2023

A145377 a(n) = A002324(n) mod 2.

Original entry on oeis.org

1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 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, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0

Offset: 1

N. J. A. Sloane, Mar 12 2009

Antti Karttunen, Table of n, a(n) for n = 1..65537
J. S. Rutherford, Generating functions for the cage isomers of the C_{20n} icosahedral fullerenes, J. Mathematical Chem., 14 (1993), 385-390. See Eq. (3).
J. S. Rutherford, The enumeration and symmetry-significant properties of derivative lattices, Act. Cryst. A48 (1992), 500-508. [See Sect. (7), p. 505.]
Index entries for characteristic functions

Essentially same as A195198.

Cf. A002324, A154272, A003050.

Array[Boole@ OddQ@ If[# < 1, 0, DivisorSum[#, KroneckerSymbol[-3, #] &]] &, 105] (* Michael De Vlieger, Nov 05 2017, after Michael Somos at A002324 *)

A002324(n) = if( n<1, 0, sumdiv(n, d, (d%3==1) - (d%3==2)));
A145377(n) = (A002324(n)%2); \\ Antti Karttunen, Nov 06 2017

a(n) = A195198(n) for n >= 1.

a(n) = Sum_{ m: m^2|n } A154272(n/m^2). - Andrey Zabolotskiy, May 07 2018

More terms from Antti Karttunen, Nov 05 2017

A300651 Number of inequivalent sublattices of index n in hexagonal lattice, where two sublattices are considered equivalent if they are related by any rotation or reflection.

Original entry on oeis.org

1, 1, 2, 3, 2, 3, 3, 5, 4, 4, 3, 8, 4, 5, 6, 9, 4, 8, 5, 10, 8, 7, 5, 15, 7, 8, 9, 13, 6, 14, 7, 15, 10, 10, 10, 20, 8, 11, 12, 20, 8, 17, 9, 17, 16, 13, 9, 28, 11, 17, 14, 20, 10, 22, 14, 25, 16, 16, 11, 34, 12, 17, 20, 27, 16, 26, 13, 24, 18, 24, 13, 40

Offset: 1

Andrey Zabolotskiy, Mar 10 2018

nonn

If we count sublattices as equivalent only if they are related by a rotation, we get A054384 instead of this sequence. If we only allow rotations and reflections that preserve the parent (hexagonal) lattice, we get A003051; the first discrepancy is at n = 42 (see illustration), the second is at n = 49. If both restrictions are applied, i.e., only rotations preserving the parent lattice are allowed, we get A145394. The analog for square lattice is A054346.

Although A003051 has its counterpart A003050 which counts primitive sublattices only, this sequence has no such counterpart sequence because a primitive sublattice can turn to a non-primitive one via a non-parent-lattice-preserving rotation, so the straightforward definition of primitiveness does not work in this case.

Andrey Zabolotskiy, Table of n, a(n) for n = 1..1000
Daejun Kim, Seok Hyeong Lee, and Seungjai Lee, Zeta functions enumerating subforms of quadratic forms, arXiv:2409.05625 [math.NT], 2024. See section 6.1 for the Dirichlet g.f. zeta^GL_{x^2+xy+y^2}(s).
Andrey Zabolotskiy, Sublattices of the hexagonal lattice (illustrations for n = 1..7, 14)
Index entries for sequences related to sublattices
Index entries for sequences related to A2 = hexagonal = triangular lattice

Cf. A003050, A003051, A054384, A145394, A054346.

# See A159842 and A054384 for the definitions of functions used here
def a_GL(n):
    return (a_SL(n) + dc(fin(1, -1, 0, 2), u, u, g2)(n)) / 2
print([a_GL(n) for n in range(1, 100)]) # Andrey Zabolotskiy, Sep 22 2024

A048259 Number of distinct solutions to x + y + z = 0 (mod n), where two solutions are equivalent if one can be obtained from the other by multiplying by units in Z/nZ and permuting x,y,z.

Original entry on oeis.org

1, 1, 2, 3, 4, 3, 7, 4, 8, 6, 8, 4, 15, 5, 10, 11, 14, 5, 17, 6, 18, 14, 12, 6, 31, 9, 14, 13, 22, 7, 33, 8, 24, 16, 16, 16, 39, 9, 18, 19, 38, 9, 41, 10, 28, 28, 20, 10, 57, 15, 30, 21, 32, 11, 43, 20, 46, 24, 24, 12, 77, 13, 26, 35, 42, 23, 53, 14, 38, 26, 52, 14, 83

Offset: 0

N. J. A. Sloane

nonn

			For n=6 the 7 solutions are (x,y,z) = (0,0,0), (5,1,0), (4,2,0), (4,1,1), (3,3,0), (3,2,1), (2,2,2).

Sean A. Irvine, Java program (github)

Cf. A007997, A007998, A003050, A003051.

iscanon(n,v)={for(d=1, n-1, if(gcd(n,d)==1 && lex(v,vecsort(v*d%n))>0, return(0))); 1}
a(n)={if(n==0, 1, sum(x=0, n-1, sum(y=x, n-1, my(z=(-x-y)%n); y<=z && iscanon(n,[x,y,z]) )))} \\ Andrew Howroyd, Jun 11 2021

a(42) onward corrected by Sean A. Irvine, Jun 10 2021

A006984 Greatest minimal norm of sublattice of index n in hexagonal lattice.

Original entry on oeis.org

1, 1, 3, 4, 3, 4, 7, 7, 9, 7, 7, 12, 13, 12, 13, 16, 13, 13, 19, 16, 21, 19, 19, 21, 25, 21, 27, 28, 21, 27, 31, 28, 27, 28, 31, 36, 37, 31, 39, 37, 37, 36, 43, 39, 39, 39, 39, 48, 49, 43, 43

Offset: 1

N. J. A. Sloane, Mira Bernstein

The hexagonal lattice is the familiar 2-dimensional lattice in which each point has 6 neighbors. This is sometimes called the triangular lattice.

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Andrey Zabolotskiy, Table of n, a(n) for n = 1..500
M. Bernstein, N. J. A. Sloane and P. E. Wright, On Sublattices of the Hexagonal Lattice, Discrete Math. 170 (1997) 29-39 (Abstract, pdf, ps).
G. Nebe and N. J. A. Sloane, Home page for hexagonal (or triangular) lattice A2.
N. J. A. Sloane, Computer printout with notes, Mar. 1994.

Cf. A003051, A003050, A001615.

A007998 Number of primitive solutions to x + y + z = 0 (mod n).

Original entry on oeis.org

1, 2, 3, 3, 4, 5, 5, 7, 8, 9, 9, 15, 11, 17, 17, 18, 18, 31, 21, 33, 30, 33, 30, 53, 34, 48, 46, 53, 43, 87, 47, 73, 69, 79, 66, 110, 65, 100, 94, 111, 80, 159, 86, 135, 130, 135, 108, 202, 116, 181, 153, 175, 133, 255, 150, 212

Offset: 1

N. J. A. Sloane

I can't reconstruct this. Possibly an erroneous version of A003050 or A048259. - N. J. A. Sloane

Cf. A007997.

A157235 Number of primitive inequivalent oblique sublattices of hexagonal (triangular) lattice of index n (equivalence and symmetry of sublattices are determined using only parent lattice symmetries).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 1, 2, 1, 2, 1, 3, 2, 2, 2, 5, 2, 4, 3, 5, 3, 4, 4, 6, 5, 6, 4, 10, 4, 6, 6, 8, 6, 10, 5, 9, 7, 8, 6, 14, 6, 10, 10, 11, 7, 12, 8, 14, 10, 12, 8, 17, 10, 12, 11, 14, 9, 20, 9, 15, 14, 14, 12, 22, 10, 16, 14, 22, 11, 20, 11, 18, 18, 18

Offset: 1

N. J. A. Sloane, Feb 25 2009

nonn

John S. Rutherford, Sublattice enumeration. IV. Equivalence classes of plane sublattices by parent Patterson symmetry and colour lattice group type, Acta Cryst. (2009). A65, 156-163. [See Table 5.]

Cf. A003051 (all sublattices), A003050 (all primitive sublattices), A154272 (primitive sublattices fully inheriting the parent lattice symmetry, inlcuding the orientation of the mirrors), A000086 (primitive rotation-symmetric sublattices, counting mirror images as distinct), A060594 (primitive mirror-symmetric sublattices), A145377 (all sublattices inheriting the parent lattice symmetry), A304182.

a(n) = A003050(n) - (A000086(n)-A154272(n))/2 - A060594(n). - Andrey Zabolotskiy, Mar 19 2021

New name and a(1)=0 prepended by Andrey Zabolotskiy, May 09 2018

Terms a(31) and beyond from Andrey Zabolotskiy, Mar 19 2021

Showing 1-9 of 9 results.

cp's OEIS Frontend

A001615 Dedekind psi function: n * Product_{p|n, p prime} (1 + 1/p).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

Python

Sage

Formula

Extensions

A007997 a(n) = ceiling((n-3)(n-4)/6).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A003051 Number of inequivalent sublattices of index n in hexagonal lattice, where two sublattices are equivalent if they are related by a rotation or reflection preserving the hexagonal lattice.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

Formula

A145377 a(n) = A002324(n) mod 2.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A300651 Number of inequivalent sublattices of index n in hexagonal lattice, where two sublattices are considered equivalent if they are related by any rotation or reflection.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

SageMath

A048259 Number of distinct solutions to x + y + z = 0 (mod n), where two solutions are equivalent if one can be obtained from the other by multiplying by units in Z/nZ and permuting x,y,z.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

PARI

Extensions

A006984 Greatest minimal norm of sublattice of index n in hexagonal lattice.