Standard Monte Carlo techniques were used to compute the potential of mean force between pairs and between triplets of freely jointed hard-sphere polymers in dilute solutions. Segment-segment interactions at poor solvent conditions were represented by square-well potentials. Well width equaled half a segment diameter and well depth was either zero or -0.30 k(B)T. Polymer chains contained 25 segments. For polymer triplets at a set of selected two-body distances, the pair-wise additivity of the potential of mean force provides a reasonable approximation for the three-body potential of mean force. At athermal conditions, the error introduced by assuming additivity is generally less than 10-15% of the total three-body interaction, while for well depth -0.30 kBT, the error rises, but is still generally less than 20-30%. Deviations from the calculated three-body potential of mean force are a function of solvent conditions and of relative positions of the interacting polymers. For polymer chains containing 15, 25 or 30 segments, simulation results do not depend significantly on polymer length.