Efficient algorithms for solving three dimensional parabolic interface problem with variable coefficients