This article concerns an optimal control problem of initial distributions for a class of hierarchical two age-structured population model, in which the state system is described by integro-partial differential equations, and the object functional presents the net profits from the exploitation of the final population distributions. With the proper construction of an adjoint system, the optimal controllers are exactly characterized via normal cones, and the optimal initial distributions must be structured by the truncated function of the adjoint vector. Furthermore, the existence and uniqueness of the optimal strategies are established by means of Ekeland variational principle and fixed point approach. In addition to theoretic analysis, some numerical experiments are conducted to examine the effects of the price functions on the optimal benefits.