Fast QR iterations for unitary plus low rank matrices