In this paper, we consider the shape inverse problem of discontinuity identification for the heat conduction. Based on the continuous adjoint method, we derive the structure of the derivative for the cost functional by employing the differentiability of a minimax formulation involving a Lagrange functional with the function space parametrization technique. Moreover, we present a gradient-type algorithm to the shape inverse problem. Finally, numerical results demonstrate the feasibility and effectiveness of the proposed algorithm.