To present a prediction of the multi-particle production and multi-fragment emission in electron-nucleus ($eA$) collisions at the forthcoming Electron-Ion Collider (EIC), a simple hybrid model which is based on the multi-source thermal model and the ideal gas model is proposed in this article. According to the hybrid model, some statistical laws such as the two-component Erlang distribution and others are presented, which means a two-source production. These statistical laws are hopeful to describe the bulk properties of multiple particles produced in the scattering of electron-nucleon ($eN$) and multiple fragments emitted in the fragmentation of excited residual nucleus. Although both the scattering and fragmentation can occur in $eA$ collisions at the EIC, their two-sources are different. In $eN$ scattering, the multiple particles come from the soft excitation and hard scattering processes respectively, which are classified as two different types of events. In nuclear fragmentation, the multiple fragments come from the cold and hot sources which exist in the same excited nucleus.

Utilizing a partitioning method based on equal (or unequal) probabilities—without incorporating the alpha-cluster (α-cluster) model—allows for the derivation of diverse topological configurations of nuclear fragments resulting from fragmentation. Subsequently, we predict the multiplicity distribution of nuclear fragments for specific excited nuclei, such as ${}^9\mathrm{Be}^{*}$, ${}^{12}\mathrm{C}^{*}$, and ${}^{16}\mathrm{O}^{*}$, which can be formed as nuclear remnants in electron–nucleus ($eA$) collisions at high energy. Based on the $\alpha$-cluster model, an $\alpha$-cluster structure may result in deviations in the multiplicity distributions of nuclear fragments with charge $Z=2$, compared to those predicted by the partitioning methods. Furthermore, in the framework of Tsallis statistics, the nonextensive generalized temperature, entropy index, and q-entropy are obtained from the multiplicity distribution of nuclear fragments with a given charge number. Our work shows that fragmentation of nuclear remnants in electron–nucleus collisions at high energy is a nonextensive process.