Background: The inner crust of neutron stars consists of a Coulomb lattice of neutron-rich nuclei, immersed in a sea of superfluid neutrons with background relativistic electron gas. A proper quantum mechanical treatment for such a system under a periodic potential is the band theory of solids. The effect of band structure on the effective mass of dripped neutrons, the so-called entrainment effect, is currently in a debatable situation, and it has been highly desired to develop a nuclear band theory taking into account neutron superfluidity in a fully self-consistent manner.Purpose: The main purpose of the present work is twofold: 1) to develop a formalism of the time-dependent self-consistent band theory, taking full account of nuclear superfluidity, based on time-dependent density functional theory (TDDFT) extended for superfluid systems, and 2) to quantify the effects of band structure and superfluidity on crustal properties, applying the formalism to the slab phase of nuclear matter in the β equilibrium.Methods: The fully self-consistent time-dependent band theory, proposed in a previous work [K. Sekizawa, S. Kobayashi, and M. Matsuo, Phys. Rev. C 105, 045807 (2022)], is extended for superfluid systems. To this end, superfluid TDDFT with a local treatment of pairing, known as time-dependent superfluid local density approximation (TDSLDA), is formulated in the coordinate space with a Skyrme-type energy density functional and the Bloch boundary condition. A real-time method is employed to extract the collective mass of a slab and that of protons, which in turn quantifies conduction neutron number density and the neutron effective mass, i.e., the entrainment effect.Results: Static calculations have been performed for a range of baryon (nucleon) number density (nb = 0.04-0.07 fm −3 ) under the β-equilibrium condition with and without superfluidity, for various inter-slab spacings. From the results, we find that the system gains energy through the formation of Cooper pairs for all densities examined, which validates the existence of superfluidity in the inner crust of neutron stars. From a dynamic response to an external potential, we extract the collective mass of a slab and that of protons immersed in neutron superfluid. From the results, we find that the collective mass of a slab is substantially reduced by 57.5-82.5% for nb = 0.04-0.07 fm −3 , which corresponds to an enhancement of conduction neutron number density and, thus, to a reduction of the neutron effective mass, which we call the anti-entrainment effect. A comparison of results with and without superfluidity reveals that superfluid effects slightly enhance the anti-entrainment effects for the slab phase of neutron-star matter. We discuss novel phenomena associated with superfluidity, quasiparticle resonances in the inner crust, which are absent in normal systems.Conclusions: Our fully self-consistent, microscopic, superfluid band calculations based on (TD)DFT showed that the effective mass of dripped neutrons is reduced by about 20-40% for nb = 0.04-0...