We provide a tight-binding model parametrization for black phosphorus (BP) with an arbitrary number of layers. The model is derived from partially self-consistent GW0 approach, where the screened Coulomb interaction W0 is calculated within the random phase approximation on the basis of density functional theory. We thoroughly validate the model by performing a series of benchmark calculations, and determine the limits of its applicability. The application of the model to the calculations of electronic and optical properties of multilayer BP demonstrates good quantitative agreement with ab initio results in a wide energy range. We also show that the proposed model can be easily extended for the case of external fields, yielding the results consistent with those obtained from first principles. The model is expected to be suitable for a variety of realistic problems related to the electronic properties of multilayer BP including different kinds of disorder, external fields, and many-body effects.