It is considered jump boundary value problem for multiply connected domains with unknown internal domains (inclusions). Distribution of the inclusions is determined by the optimality of the functional of effective conductivity and the constancy of their complete area. The solution is given in form of sums depending on centers of inclusions. For two and three inclusions optimal configurations are described, and the closed form solution is found in the important for application particular case.