Researchers have proposed a polynomial-time approximation algorithm that achieves a factor of \(2+\epsilon\) for this problem, which includes various geometric shapes such as disks and rectangles [4].