|
Search: id:A058360
|
|
|
| A058360 |
|
Number of partitions of n whose reciprocal sum is an integer. |
|
+0 2
|
|
| 1, 1, 1, 2, 2, 2, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 17, 19, 23, 25, 31, 33, 38, 42, 51, 57, 66, 75, 86, 97, 109, 122, 138, 155, 177, 200, 230, 253, 287, 320, 363, 405, 456, 507, 572, 639, 707, 785, 877, 971, 1079, 1198, 1334, 1476, 1635, 1802, 2002, 2213, 2445, 2700
(list; graph; listen)
|
|
|
OFFSET
|
1,4
|
|
|
REFERENCES
|
From a question posted to the news group comp.soft-sys.math.mathematica by "Juan" (erfa11(AT)hotmail.com) at Steven M. Christensen and Associates, Inc and MathTensor, Inc. Jan 22, 2002 08:46:57 +0000 (UTC).
|
|
EXAMPLE
|
a(12) = 7 because the partitions of 12 whose reciprocal sum is an integer are: {{6, 3, 2, 1}, {4, 4, 2, 1, 1}, {3, 3, 3, 1, 1, 1}, {2, 2, 2, 2, 2, 2}, {2, 2, 2, 2, 1, 1, 1, 1}, {2, 2, 1, 1, 1, 1, 1, 1, 1, 1}, {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1}}. Individually their reciprocal sums are: 2, 3, 4, 3, 6, 9, and 12.
|
|
MATHEMATICA
|
<< DiscreteMath`Combinatorica`; f[n_] := (p = Partitions[n]; is = Compile[ {{x, _Integer, 1}}, Plus @@ (1/x)]; ans = p[[ Flatten[ Position[ FractionalPart[ is /@ p], x_ /; x < .000001 || x > 0.999999]]]]); Table[ Length[ f[n]], {n, 1, 65} ]
|
|
CROSSREFS
|
Cf. A066824.
Adjacent sequences: A058357 A058358 A058359 this_sequence A058361 A058362 A058363
Sequence in context: A029115 A029101 A029080 this_sequence A098527 A035635 A029100
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Robert G. Wilson v (rgwv(AT)rgwv.com), Jan 25 2002
|
|
|
Search completed in 0.002 seconds
|