The mathematical essence of structural topology optimization is large-scale nonlinear integer programming. To overcome its huge computational burden, a popular way is to relax the 0–1 variable constraints and transform the integer programming problem into a continuous variable programming problem which can be solved by using gradient based mathematical programming methods. To cope with the variable transformation, the well-known SIMP (Solid Isotropic Material with Penalty) method introduces interpolation schemes for the material properties versus design variables with penalty and achieves great success and popularity. However, there is no doubt that directly tackling the large-scale nonlinear integer programming is very important. This paper solves the structural topology optimization problems with single or multiple constraints by applying the Canonical Dual Theory (CDT) by Gao and Ruan (2010) together with Sequential Approximate Programming approach under the classic structural topology optimization formulations.Following the Sequential Approximate Programming (SAP) frame from the structural optimization, the present paper firstly utilizes sensitivity information to construct the explicit and separable approximate Sequential Quadratic Integer Programming (SQIP) or Sequential Linear Integer Programming (SLIP) subproblems for the topology optimization. And then, the subproblems are solved by applying the Canonical relaxation algorithm based on CDT theory. Their special mathematical structures are exploited to develop analytic solution of Kuhn–Tucker condition of the dual programming. Numerical experiments of two linear and quadratic integer programming problems with random coefficients assert that the Canonical relaxation algorithm can obtain approximate solutions with good properties very efficiently and the dual gap is negligible when the number of design variables increases.Because move limit strategies play a key role in many search algorithms of structural optimization, this paper combines one of two different move limit strategies within the new method. The new method first solves a set of classic topology optimization problems with only material usage constraint, including minimum structural compliance design under constant load, maximum heat transfer efficiency for the heat conduction problem. And then we apply the method to the topology optimization problems with multiple constraints, including minimum structural compliance design under an additional local displacement constraint and minimum structural compliance design under infill constraints. The results of these problems demonstrate that the new method can efficiently solve the discrete variable structural topology optimization problems with multiple nonlinear constraints or many local linear constraints in a unified and systematic way. Beyond that, the new method can achieve integer solutions when combined with the move limit strategy of controlling volume fraction parameter. It can deal with much more design variables than the general branch-and-bound method and makes no use of any sensitivity threshold or heuristic stabilization scheme during the iterative process in comparison with BESO method. Finally, the new method in this paper can be further developed as a general solver for these large-scale discrete variable structural topology optimization problems.