(defun power-set (s)
  (mappend* (op comb s) (range 0 (length s))))