%I #97 Jan 22 2024 07:30:15

%S 1,7,13,35,31,91,57,155,130,217,133,455,183,399,403,651,307,910,381,

%T 1085,741,931,553,2015,806,1281,1210,1995,871,2821,993,2667,1729,2149,

%U 1767,4550,1407,2667,2379,4805,1723,5187,1893,4655,4030,3871,2257,8463,2850,5642,3991,6405,2863

%N Number of sublattices of index n in generic 3-dimensional lattice.

%C These sublattices are in 1-1 correspondence with matrices

%C [a b d]

%C [0 c e]

%C [0 0 f]

%C with acf = n, b = 0..c-1, d = 0..f-1, e = 0..f-1. The sublattice is primitive if gcd(a,b,c,d,e,f) = 1.

%C Total area of all distinct rectangles whose side lengths are divisors of n, and whose length is an integer multiple of the width. - _Wesley Ivan Hurt_, Aug 23 2020

%D Richard P. Stanley, Enumerative Combinatorics, Cambridge, Vol. 2, 1999; see Problem 5.13(d), pp. 76 and 113.

%H Amiram Eldar, <a href="/A001001/b001001.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from T. D. Noe)

%H Michael Baake, <a href="https://arxiv.org/abs/math/0605222">Solution of the coincidence problem in dimensions d <= 4</a>, in R. V. Moody, ed., Math. of Long-Range Aperiodic Order, Kluwer 1997, pp. 9-44.; arXiv:math/0605222 [math.MG], 2006.

%H M. Baake and N. Neumarker, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL12/Baake/baake7.html">A Note on the Relation Between Fixed Point and Orbit Count Sequences</a>, JIS 12 (2009) 09.4.4, Section 3.

%H J. Liouville, <a href="http://sites.mathdoc.fr/JMPA/PDF/JMPA_1857_2_2_A4_0.pdf">Théorème concernant les sommes de diviseurs des nombres</a>, Journal de mathématiques pures et appliquées 2e série, tome 2 (1857), p. 56.

%H V. A. Liskovets and A. Mednykh, <a href="https://doi.org/10.1080/00927870008826924">Enumeration of subgroups in the fundamental groups of orientable circle bundles over surfaces</a>, Commun. in Algebra, 28, No. 4 (2000), 1717-1738.

%H Y. Puri and T. Ward, <a href="http://www.cs.uwaterloo.ca/journals/JIS/VOL4/WARD/short.html">Arithmetic and growth of periodic orbits</a>, J. Integer Seqs., Vol. 4 (2001), #01.2.1.

%H J. S. Rutherford, <a href="https://dx.doi.org/10.1107/S0108767392000898">The enumeration and symmetry-significant properties of derivative lattices</a>, Act. Cryst. (1992) A48, 500-508

%H J. S. Rutherford, <a href="https://doi.org/10.1107/S0108767392007657">The enumeration and symmetry-significant properties of derivative lattices II</a>, Acta Cryst. A49 (1993), 293-300. [_N. J. A. Sloane_, Mar 14 2009]

%H Tad White, <a href="http://arxiv.org/abs/1304.2830">Counting Free Abelian Actions</a>, arXiv:1304.2830 [math.CO], 2013.

%H <a href="/index/Su#sublatts">Index entries for sequences related to sublattices</a>.

%F If n = Product p^m, a(n) = Product (p^(m + 1) - 1)(p^(m + 2) - 1)/(p - 1)(p^2 - 1). Or, a(n) = Sum_{d|n} sigma(n/d)*d^2, Dirichlet convolution of A000290 and A000203.

%F a(n) = Sum_{d|n} d*sigma(d). - _Vladeta Jovovic_, Apr 06 2001

%F Multiplicative with a(p^e) = ((p^(e+1)-1)(p^(e+2)-1))/((p-1)(p^2-1)). - _David W. Wilson_, Sep 01 2001

%F Dirichlet g.f.: zeta(s)*zeta(s-1)*zeta(s-2).

%F L.g.f.: -log(Product_{k>=1} (1 - x^k)^sigma(k)) = Sum_{n>=1} a(n)*x^n/n. - _Ilya Gutkovskiy_, May 23 2018

%F a(n) = Sum_{d1|n, d2|n, d1|d2} d1*d2. - _Wesley Ivan Hurt_, Aug 23 2020

%F Sum_{k=1..n} a(k) ~ c * n^3, where c = Pi^2*zeta(3)/18 = 0.659101... . - _Amiram Eldar_, Oct 19 2022

%p nmax := 100:

%p L12 := [seq(1,i=1..nmax) ];

%p L27 := [seq(i,i=1..nmax) ];

%p L290 := [seq(i^2,i=1..nmax) ];

%p DIRICHLET(L12,L27) ;

%p DIRICHLET(%,L290) ; # _R. J. Mathar_, Sep 25 2017

%t a[n_] := Sum[ d*DivisorSigma[1, d], {d, Divisors[n]}]; Table[ a[n], {n, 1, 42}] (* _Jean-François Alcover_, Jan 20 2012, after _Vladeta Jovovic_ *)

%t f[p_, e_] := Product[(p^(e + k) - 1)/(p^k - 1), {k, 1, 2}]; a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* _Amiram Eldar_, Aug 29 2019 *)

%o (PARI)

%o N=17; default(seriesprecision,N); x=z+O(z^(N+1))

%o c=sum(j=1,N,j*x^j);

%o t=1/prod(j=1,N, eta(x^(j))^j)

%o t=log(t)

%o t=serconvol(t,c)

%o Vec(t)

%o /* _Joerg Arndt_, May 03 2008 */

%o (PARI) a(n)=sumdiv(n,d, d * sumdiv(d,t, t ) ); /* _Joerg Arndt_, Oct 07 2012 */

%o (PARI) a(n)=sumdivmult(n,d, sigma(d)*d) \\ _Charles R Greathouse IV_, Sep 09 2014

%Y Column 3 of A160870.

%Y Cf. A060983, A064987 (Mobius transform).

%Y Cf. A061256, A226313, A301777.

%Y Primes in this sequence are in A053183.

%Y Cf. A038991, A038992, A038993, A038994, A038995, A038996, A038997, A038998, A038999.

%K nonn,easy,nice,mult

%O 1,2

%A _N. J. A. Sloane_