In this paper, we consider the problem of minimizing the function f ( x ) = g 1 ( x ) + g 2 ( x ) − h ( x ) over R n , where g 1 is a proper and lower semicontinuous function, g 2 is continuously differentiable with a Hölder continuous gradient and h is a convex function that may be nondifferentiable. This problem has important practical applications but is challenging to solve due to the presence of nonconvexities and nonsmoothness. To address this issue, we propose an algorithm based on a proximal gradient method that uses a quadratic approximation of the function g 2 and a nonconvex regularization term. We show that the number of iterations required to reach our stopping criterion is O ( max { ϵ − β + 1 β , η 2 β ϵ − 2 ( β + 1 ) β } ) . Our approach offers a promising strategy for solving this challenging optimization problem and has potential applications in various fields. Numerical examples are provided to illustrate the theoretical results.