An analytical iterative algorithm for fast computation of scattering from multiple conductive cylinders is developed. It takes account of the independent 1-st order, coupling 2-nd and higher order scattering of multiple objects. The 1-st order scattering, as independent scattering of a single object self, is derived by using the expansion of cylindrical waves in local coordinate of the object, which are solved by the boundary conditions. Exciting by the 1-st order scattering from one object, the 2-nd order scattering of multi-objects are obtained iteratively using the addition theorem of Hankel function and the boundary conditions. The same approach is applicable to calculations of the 3-rd order and higher order scattering. It is found that the phase differences caused by the cylinder centers indicate the scattering wave interferences of multi-objects. Comparison with numerical MoM results well vailidates the analytical iterative algorithm, which significantly accelerates the scattering computation for wide-band frequency f and overall azimuth angles theta over 360 degrees. In order to depress the ringing effect caused by the abrupt change of scattering field in computed region to zero field in uncomputed region, the Hamming window is used to make gradual change of the scattering field to zero. Using FFT and 2D 3-order spline interpolation, discrete scattering data over f - theta plane yield uniform discrete data over X - Y plane. The image reconstruction is finally performed for multiple conductive cylinders. The image can well identify the positions and sizes of multi-objects.