[Solution] For any natural number n and A={a1,a2,⋯,an}. The sum of the products of the elements of all non-empty subsets of A is
a1=+a2+⋯+an+(a1a2+⋯+a1an+⋯+an−1an)+(a1a2a3+⋯+an−1ana1)+⋯+a1a2⋯an(1+a1)(1+a2)⋯(1+an)−1.
The number of all non-empty subsets of A is 2n−1. Therefore,
m(A)=2n−1(1+a1)⋯(1+an)−1,
which implies
(1+a1)⋯(1+an)−1=m(A)(2n−1).
From this equation, when n=r+1,
==(2r+1−1)m(S∪{ar+1})+1(1+a1)⋯(1+ar)(1+ar+1)[m(S)(2r−1)+1](1+ar+1).
Given m(S)=13,m(S∪{ar+1})=49,
49(2r+1−1)+1=[13(2r−1)+1](1+ar+1).
From this equation, we can derive
2r=13ar+1−8512(ar+1−3).
Notice that the right-hand side of the above equation is a function of ar+1, and it is a decreasing function in the intervals (−∞,1385) and (1385, +∞).
When ar+1=1, the right-hand side of (1) equals 7224=31. Therefore, when ar+1=1,2,3,4,5,6, the right-hand side of (1) is less than 1. Thus, we only need to consider the case where ar+1⩾7. When ar+1⩾8, the right-hand side of (1) is not an integer, so we must have ar+1=7.
From (1), we get
2r=91−8512×4=8
Thus, r=3.
Therefore,
=(1+a1)(1+a2)(1+a3)=(23−1)m(S)+17×13+1=92.
Hence,
(1+a1)(1+a2)(1+a3)=92=2×2×23.
From this, we obtain the unique positive integer solution.