Reduced density matrix functional theory for superconductors