OFFSET
1,2
COMMENTS
a(n) is multiplicative: if m and n are relatively prime then a(m*n) = a(n) * a(m). For n >= 2, a(n)>=2 with equality iff n is prime.
LINKS
Eric M. Schmidt, Table of n, a(n) for n = 1..10000
Max Alekseyev, PARI scripts for various problems
G. A. Miller, On the Subgroups of an Abelian Group, The Annals of Mathematics, 2nd Ser., Vol. 6, No. 1. (1904), pp. 1-6. doi:10.2307/2007151 [See paragraph 4 entitled "Total number of subgroups in a group of order p^m". - M. F. Hasler, Dec 03 2007]
G. A. Miller, Determination of the number of subgroups of an abelian group, Bull. Amer. Math. Soc. 33 (1927), 192-194.
FORMULA
(C_2)^m has A006116(m) subgroups, so this is a lower bound if n is a power of 2 (e.g. a(16) >= 67). - N. J. A. Sloane, Dec 01 2007
EXAMPLE
a(16) = 67: C16 has 5 subgroups, C2 X C8 has 11 subgroups, (C2)^2 X C4 has 27 subgroups, (C2)^4 has 67 subgroups, (C4)^2 has 15 subgroups.
PROG
(PARI) { A061034(n) = my(f=factorint(n)); prod(i=1, #f~, vecmax( apply( x->numsubgrp(f[i, 1], x), partitions(f[i, 2]) ) ) ); } \\ See Alekseyev link for numsubgrp(), Max Alekseyev, 2008
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Ola Veshta (olaveshta(AT)my-deja.com), May 26 2001
EXTENSIONS
More terms from Victoria A Sapko (vsapko(AT)canes.gsw.edu), Jun 13 2003
STATUS
approved