In this paper, the lifetime model based on series systems with a random number of components from the family of power series distributions has been considered. First, some basic theoretical results have been obtained, which have been used to optimize the number of components in series systems. The average lifetime of the system, the cost function, and the total time on test have been used as an objective function in optimization. The issue has been investigated in detail when the lifetimes of system components have Weibull distribution, and the number of components has geometric, logarithmic, or zero-truncated Poisson distributions. The results have been given analytically and numerically. Finally, a real data set has been used to illustrate the obtained results.