In this paper, we first prove that the existence of a solution of SDEs under the assumptions that the drift coefficient is of linear growth and path--dependent, and diffusion coefficient is bounded, uniformly elliptic and Hölder continuous. We apply Gaussian upper bound for a probability density function of a solution of SDE without drift coefficient and local Novikov condition, in order to use Maruyama--Girsanov transformation. The aim of this paper is to prove the existence with explicit representations (under linear/super--linear growth condition), Gaussian two--sided bound and Hölder continuity (under sub--linear growth condition) of a probability density function of a solution of SDEs with path--dependent drift coefficient. As an application of explicit representation, we provide the rate of convergence for an Euler--Maruyama (type) approximation, and an unbiased simulation scheme.