To answer your first question, if this is what you mean, try this.
Clear[a, b, c]; p =
Solve[{a^2 + b^2 == c^2, a > 0, b > 0, c > 0, a < b, c <= 50}, {a, b,
c}, Integers]
and to then store the results, I use this method though there are other ways to achieve this.
q = {a, b, c} /. p
Hope this helps.