Efficient Gradient Approximation Method for Constrained Bilevel Optimization