OFFSET
1,4
COMMENTS
Given a finite multiset S of positive integers greater than one, let G(S) be the simple labeled graph with vertex set S and edges between any two vertices with a common divisor greater than 1. For example, G({6,14,15,35}) is a 4-cycle. This sequence counts factorizations S such that G(S) is a connected graph.
a(n) depends only on prime signature of n (cf. A025487). - Antti Karttunen, Nov 07 2018
LINKS
EXAMPLE
The a(72) = 10 factorizations:
(72),
(2*2*18), (2*3*12), (2*6*6), (3*4*6),
(2*36), (3*24), (4*18), (6*12),
(2*2*3*6).
MATHEMATICA
zsm[s_]:=With[{c=Select[Tuples[Range[Length[s]], 2], And[Less@@#, GCD@@s[[#]]]>1&]}, If[c=={}, s, zsm[Union[Append[Delete[s, List/@c[[1]]], LCM@@s[[c[[1]]]]]]]]];
facs[n_]:=If[n<=1, {{}}, Join@@Table[Map[Prepend[#, d]&, Select[facs[n/d], Min@@#>=d&]], {d, Rest[Divisors[n]]}]];
Table[Length[Select[facs[n], Length[zsm[#]]==1&]], {n, 100}]
PROG
(PARI)
is_connected(facs) = { my(siz=length(facs)); if(1==siz, 1, my(m=matrix(siz, siz, i, j, (gcd(facs[i], facs[j])!=1))^siz); for(n=1, siz, if(0==vecmin(m[n, ]), return(0))); (1)); };
A305193aux(n, m, facs) = if(1==n, is_connected(Set(facs)), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs, d); s += A305193aux(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Nov 07 2018
A305193(n) = if(1==n, 0, A305193aux(n, n, List([]))); \\ Antti Karttunen, Nov 07 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, May 27 2018
EXTENSIONS
More terms from Antti Karttunen, Nov 07 2018
STATUS
approved