OFFSET
1,2
COMMENTS
The rule of building products is (a,b)*(x,y) = (a*x,b*y).
a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24=2^3*3 and 375=3*5^3 both have prime signature (3,1).
LINKS
FORMULA
EXAMPLE
From Alois P. Heinz and Antti Karttunen, Nov 24 2017: (Start)
a(4) = 9 because for pair (4,4) there are nine factorizations:
(4,4)
(1,4)*(4,1)
(1,2)*(4,2)
(2,1)*(2,4)
(2,2)*(2,2)
(1,2)*(2,1)*(2,2)
(1,4)*(2,1)*(2,1)
(4,1)*(1,2)*(1,2)
(1,2)*(1,2)*(2,1)*(2,1)
(End)
a(pq) = 15 for primes p<>q: (pq,pq); (p,1)(q,pq); (p,1)(q,1)(1,pq); (p,1)(q,1)(1,p)(1,q); (p,1)(q,q)(1,p); (p,1)(q,p)(1,q); (p,q)(q,p); (p,q)(q,1)(1,p); (p,p)(q,q) ; (p,p)(q,1)(1,q); (p,pq)(q,1); (pq,1)(1,pq); (pq,1)(1,p)(1,q); (pq,q)(1,p); (pq,p)(1,q). - R. J. Mathar, Nov 30 2017
PROG
(PARI) a(n) = if(n==1, return(1)); my(b, c, r, x, y, v=List([]), w=List([[n]])); while(#w>r, c++; for(k=r+1, r=#w, y=w[k]; if(!isprime(x=y[c]), fordiv(x, d, if(d!=1&&d!=x, listput(w, concat([y[1..c-1], d, x/d]))))))); for(i=1, #w, x=w[i]; r=#x; for(j=1, #w, y=w[j]; for(k=0, 2^r-1, b=concat(b=binary(k), vector(r-#b)); if(#y>=t=vecsum(b), c=0; listput(v, vecsort(vector(r+#y-t, m, if(m>r, [1, y[m-r+t]], if(b[m], [x[m], y[c++]], [x[m], 1]))))))))); #Set(v); \\ Jinyuan Wang, Jan 17 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Christian G. Bower, Jun 03 2005
STATUS
approved