OFFSET
1,2
EXAMPLE
For p = 5, p^2 - 1 = 24 = 2^3 * 3^1. The number of abelian groups of order 24 = (the number of partitions of 3)*(the number of partitions of 1) = 3*1 = 3.
MAPLE
a:= n-> mul(combinat[numbpart](i[2]), i=ifactors(ithprime(n)^2-1)[2]):
seq(a(n), n=1..79); # Alois P. Heinz, Dec 31 2022
PROG
(Python)
from sympy import factorint, npartitions
from math import prod
def A359421(n): return prod(npartitions(d) for d in factorint(prime(n)**2-1).values()) # Chai Wah Wu, Jan 12 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Ali Ramsey, Dec 31 2022
STATUS
approved