OFFSET
1,24
COMMENTS
LINKS
Felix Huber, Table of n, a(n) for n = 1..10000 (terms 1..2048 from Antti Karttunen, terms 2049..8192 from Seiichi Manyama)
Wikipedia, Burnside's lemma
FORMULA
EXAMPLE
a(24) = 2 because 24 = 2*2*6 = 2*3*4, two products of three proper factors of 24.
MAPLE
A122179 := proc(n)
local i, l, x;
l := [seq(i[2], i in ifactors(n)[2])];
(mul(binomial(x + 2, 2), x in l) + 3*mul(floor(x/2) + 1, x in l) + 2*`if`(andmap(x -> x mod 3 = 0, l), 1, 0))/6 - ceil(mul(x + 1, x in l)/2);
end proc:
seq(A122179(n), n = 1 .. 88); # Felix Huber, May 27 2026
PROG
(PARI) for(n=1, 105, t=0; for(x=2, n-1, for(y=x, n-1, for(z=y, n-1, if(x*y*z==n, t++)))); print1(t, ", "))
(PARI) A122179(n) = { my(s=0); fordiv(n, x, if((x>1)&&(x<n), for(y=x, n-1, for(z=y, n-1, if(x*y*z==n, s++))))); (s); }; \\ Antti Karttunen, Aug 24 2017
(Python)
from math import prod, comb
from sympy import factorint
def A122179(n):
f = factorint(n).values()
return (prod(comb(e+2, 2) for e in f)+3*prod((e>>1)+1 for e in f)+((not any(e%3 for e in f))<<1))//6-(prod(e+1 for e in f)+1>>1) # Chai Wah Wu, May 28 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Rick L. Shepherd, Aug 23 2006
STATUS
approved
