A295281 Number of complete strict tree-factorizations of n > 1.
1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 4, 1, 0, 1, 1, 1, 3, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 9, 1, 1, 1, 0, 1, 4, 1, 1, 1, 4, 1, 6, 1, 1, 1, 1, 1, 4, 1, 1, 0, 1, 1, 9, 1, 1, 1, 1, 1, 9, 1
Offset: 2
Keywords
Examples
The a(72) = 6 complete strict tree-factorizations are: 2*3*(2*(2*3)), 2*(2*3*(2*3)), 2*(2*(3*(2*3))), 2*(3*(2*(2*3))), 3*(2*(2*(2*3))), (2*3)*(2*(2*3)).
Links
- Andrew Howroyd, Table of n, a(n) for n = 2..10000
Programs
-
Mathematica
postfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[postfacs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]]; sftc[n_]:=Prepend[Join@@Function[fac,Tuples[sftc/@fac]]/@Select[postfacs[n],And[Length[#]>1,UnsameQ@@#]&],n]; Table[Length[Select[sftc[n],FreeQ[#,_Integer?(!PrimeQ[#]&)]&]],{n,2,100}]
-
PARI
seq(n)={my(v=vector(n), w=vector(n)); v[1]=1; for(k=2, n, w[k]=v[k]+isprime(k); forstep(j=n\k*k, k, -k, v[j]+=w[k]*v[j/k])); w[2..n]} \\ Andrew Howroyd, Nov 18 2018
Comments