OFFSET
0,3
COMMENTS
See a comment on A254811 about multiset partitions and the Knuth reference. - Wolfdieter Lang, Mar 26 2015
LINKS
Brian Chen, Table of n, a(n) for n = 0..64
EXAMPLE
For n = 2 the a(2) = 5 ways to put 2 red balls and 2 blue balls into 2 indistinguishable boxes are (RRBB)(), (RRB)(B), (RBB)(R), (RR)(BB), (RB)(RB).
MAPLE
with(numtheory):
b:= proc(n, k, i) option remember;
`if`(n>k, 0, 1) +`if`(isprime(n) or i<2, 0, add(
`if`(d>k, 0, b(n/d, d, i-1)), d=divisors(n) minus {1, n}))
end:
a:= n-> b(6^n$2, n):
seq(a(n), n=0..20); # Alois P. Heinz, Mar 26 2015
MATHEMATICA
b[n_, k_, i_] := b[n, k, i] = If[n > k, 0, 1] + If[PrimeQ[n] || i < 2, 0, Sum[If[d > k, 0, b[n/d, d, i - 1]], {d, Divisors[n] [[2 ;; -2]]}]]; a[n_] := b[6^n, 6^n, n]; Table[a[n], {n, 0, 20}] (* Jean-François Alcover, Jan 08 2016, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Brian Chen, Feb 08 2015
STATUS
approved