Branching diffusion representation of semilinear PDEs and Monte Carlo approximation