A hybrid multilevel-active set method for large box-constrained discrete ill-posed problems