OFFSET
1,2
COMMENTS
a(n) is also the number of lattices in R^4 with determinant n.
LINKS
Wikipedia, Hermite normal form
FORMULA
a(n) = Sum_{ d_1 * d_2 * d_3 * d_4 = n ; d_1 <= d_2 <= d_3 <= d_4 } d_2 * d_3^2 * d_4^3.
For p prime: a(p) = p^3.
EXAMPLE
For n=2, there are 8 4 X 4 matrices M in Hermite normal form (rows convention) with det(M) = 2:
{{1, 0, 0, 0}, {0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 2}},
{{1, 0, 0, 0}, {0, 1, 0, 0}, {0, 0, 1, 1}, {0, 0, 0, 2}},
{{1, 0, 0, 0}, {0, 1, 0, 1}, {0, 0, 1, 0}, {0, 0, 0, 2}},
{{1, 0, 0, 0}, {0, 1, 0, 1}, {0, 0, 1, 1}, {0, 0, 0, 2}},
{{1, 0, 0, 1}, {0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 2}},
{{1, 0, 0, 1}, {0, 1, 0, 0}, {0, 0, 1, 1}, {0, 0, 0, 2}},
{{1, 0, 0, 1}, {0, 1, 0, 1}, {0, 0, 1, 0}, {0, 0, 0, 2}},
{{1, 0, 0, 1}, {0, 1, 0, 1}, {0, 0, 1, 1}, {0, 0, 0, 2}}
MATHEMATICA
a[n_] :=Module[{factors = Union[Sort/@Select[Tuples[Divisors[n], 4], Times @@ # == n &]]}, Sum[Product[ factors[[i, j]]^(j - 1) , {j, 4 }] , {i, Length[factors]}]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Andrew Ortegaray, Sep 23 2025
STATUS
approved
