The BCS equations are the centerpiece of the microscopic description of superconductivity. Their space discretization yields a system of coupled ordinary differential equations. In this work, we come up with fast time evolution schemes based on a splitting approach. One of the schemes only requires basic operations. For the physically important case of the BCS equations for a contact interaction potential, the computational cost of the schemes increases only linearly with the dimension of the space discretization. Their accuracy is demonstrated in extensive numerical experiments. These experiments also show that the physical energy of the system is preserve up to very small errors.