Разложение натурального числа на сумму 5 слагаемых, где каждое слагаемое является степенью двойки

//////////////////////////////////////////////////////////////////////////////////////
#include <iostream>
#include <set>
//////////////////////////////////////////////////////////////////////////////////////
typedef std::multiset<unsigned>  T_multiset;
//////////////////////////////////////////////////////////////////////////////////////
int main()
{
    for(;;)
    {
        std::cout << std::endl
                  << std::endl
                  << "n = ";
        unsigned    n = 0;
        std::cin >> n;
        T_multiset  multiset;
        
        unsigned    nn  = n;
        unsigned    k   = 1;
        while(nn)
        {
            if(nn & 1)
            {
                multiset.insert(k);
            }
            k   <<= 1;
            nn  >>= 1;
        }
 
        T_multiset::iterator  lower_bound_of_2;
        while(multiset.size() < 5
              && (lower_bound_of_2 = multiset.lower_bound(2)) != multiset.end())
        {
            unsigned  first_not_unit_val = *lower_bound_of_2;            
            multiset.insert(first_not_unit_val / 2);
            multiset.insert(first_not_unit_val / 2);
            multiset.erase(multiset.find(first_not_unit_val));
        }        
 
        if(multiset.size() != 5)
        {
            std::cout << "No"
                      << std::endl;
        }
        else
        {
            std::cout << n
                      <<" -> ";
                std::copy(multiset.begin(), multiset.end(), 
                          std::ostream_iterator<unsigned>(std::cout, "  "));
            std::cout << std::endl;
        }    
    }
}