0,2

For the unlabeled case, the number of one-element transitions from the partitions of n to the partitions of n+1 is given by A000070. Example: There are A000070(4) = 12 transitions from n=4 to n=5: [1111] -> [11111], [1111] -> [1112], [112] -> [1112], [112] -> 113], [112] -> [122], [13] -> [113], [13] -> [14], [13] -> [23], [22] -> [23], [22] -> [122], [4] -> [14], [4] -> [5].

a(n) is also the total number of digits in all partitions of the integer n which contain at least one digit 1.

T. Wieder, Home Page.

T. Wieder, (old) Home Page.

a(n) = Sum_k=1^p(n) (nops(p(k, n)) + 1), where p(n) is the number of partitions of n and nops(p(k, n)) is the number of parts in the k-th partition p(n, k) of n.

a(n) = Sum_k=1^p(n) nops(p(k, n)[subject to: at least one p(l, k, n) = 1]; p(n) = number of partitions of n, p(k, n) = k-th partition, p(l, k, n) = l-th part in the k-th partition p(k, n) of integer n.

In the labeled case, we have 9 one-element transitions from all partitions of n=3 to the partitions of n+1=4: [1,1,1] -> [1,1,1,1]; [1,1,1] -> [1,1,2]; [1,1,1] -> [1,1,2]; [1,1,1] -> [1,1,2]; [1,2] -> [1,1,2]; [1,2] -> [1,3]; [1,2] -> [2,2]; [3] -> [1,3]; [3] -> [4].

For n=4 we have the following partitions which contain at least one digit 1: [1111], [112], [13] and these partitions contain 4 + 3 + 2 = 9 = a(4) digits.

main := proc(n::integer) local a, ndxp, ListOfPartitions; with(combinat): with(ListTools): ListOfPartitions:=partition(n-1); a:=0; for ndxp from 1 to nops(ListOfPartitions) do if Occurrences(1, ListOfPartitions[ndxp]) > 0 then a:=a+nops(Flatten(ListOfPartitions[ndxp])); print("ndxp, Flatten(ListOfPartitions[ndxp]):", ndxp, Flatten(ListOfPartitions[ndxp])); print("ndxp, ListOfPartitions[ndxp], a:", ndxp, ListOfPartitions[ndxp], a); # End of if-clause *** Occurrences(1, ListOfPartitions[ndxp]) *** fi; end do; print("n, a(n):", n, a); end proc;

(* first do *) Needs["DiscreteMath`Combinatorica`"] (* then *) f[n_] := Block[{l = Sort[ Flatten[ Partitions[n]]]}, Length[l] - Count[l, 1]]; g[n_] := (f[n] + Sum[PartitionsP[k], {k, 0, n}]); Table[ g[n], {n, 0, 40}] (from Robert G. Wilson v Jul 13 2004)

Cf. A000070, A093695, A089378.

Sequence in context: A062492 A165271 A139672 this_sequence A068006 A000097 A081996

Adjacent sequences: A093691 A093692 A093693 this_sequence A093695 A093696 A093697

nonn

Thomas Wieder (wieder.thomas(AT)t-online.de), Apr 10 2004

More terms from Robert G. Wilson v (rgwv(AT)rgwv.com), Jul 13 2004