In this paper we study the asymptotic behavior of the (skew) Macdonald and Jack symmetric polynomials as the number of variables grows to infinity. We characterize their limits in terms of certain variational problems. As an intermediate step, we establish a large deviation principle for the $θ$ analogue of non-intersecting Bernoulli random walks. When $θ=1$, these walks are equivalent to random Lozenges tilings of strip domains, where the variational principle (with general domains and boundary conditions) has been proven in the seminal work by Cohn, Kenyon, and Propp. Our result gives a new argument of this variational principle, and also extends it to non-intersecting $θ$-Bernoulli random walks for any $θ\in (0,\infty)$. Remarkably, the rate functions remain identical, differing only by a factor of $1/θ$.