The dependence on the structure functions and Z, N numbers of the nuclear binding energy is investigated within the inverse problem(IP) approach. This approach allows us to infer the underlying model parameters from experimental observation, rather than to predict the observations from the model parameters. The IP was formulated for the numerical generalization of the semi-empirical mass formula of BW. It was solved in step by step way based on the AME2012 nuclear database. The established parametrization describes the measured nuclear masses of 2564 isotopes with a maximum deviation less than 2.6 MeV, starting from the number of protons and number of neutrons equal to 1. The set of parameters ${a_{i}}$, $i=1,dots, {mathcal{N}}_{rm{param}}$ of our fit represent the solution of an overdetermined system of nonlinear equations, which represent equalities between the binding energy $E_{B,j}^{rm{Expt}}(A,Z)$ and its model $E_{B,j}^{rm{Th}}(A,Z,{a_{i}})$, where $j$ is the index of the given isotope. The solution of the overdetermined system of nonlinear equations has been obtained with the help of the Aleksandrovs auto-regularization method of Gauss-Newton(GN) type for ill-posed problems. The efficiency of the above methods was checked by comparing relevant results with the results obtained independently. The explicit form of unknown functions was discovered in a step-by-step way using the modified least $chi^{2}$ procedure, that realized in the algorithms which were developed by Aleksandrov to solve nonlinear systems of equations via the GN method, lets us to choose between two functions with same $chi^{2}$ the better one. In the obtained generalized model the corrections to the binding energy depend on nine proton (2, 8, 14, 20, 28, 50, 82, 108, 124) and ten neutron (2, 8, 14, 20, 28, 50, 82, 124, 152, 202) magic numbers as well on the asymptotic boundaries of their influence.