We derive the joint density of a Skew Brownian motion, its last visit to the origin, its local and occupation times. The result enables us to obtain explicit analytical formulas for pricing European options under both a two‐valued local volatility model and a displaced diffusion model with constrained volatility.