Деление в двоичной системе счисления с остатком

#include <iostream>
#include <bitset> 
using namespace std;
unsigned short div(unsigned int a,unsigned short b)
{
    return (unsigned short)(a%b);
}
int main()
{
   unsigned int x;
   unsigned short y;
   cin>>x>>y; // ввод десятичного представления чисел
   cin.get();
   x &= 0x007FFFFF; // обнулить разряды старше 23-го
   y &= 0x0FFF; // обнулить разряды старше 12-го
   bitset<12> b( div(x,y) );
   cout<<b<<endl; // вывод двоичного представления b
   cin.get();
   return 0;
}

#include <iostream>
#include <bitset> 
#include <cstdlib> 
using namespace std;
bitset<12> div(bitset<23> a,bitset<12> b)
{
    return bitset<12>(a.to_ulong()%b.to_ulong());
}
int main()
{
   bitset<23> x;
   bitset<12> y;
   cin>>x; // ввод двоичного представления x
   cin>>y;  // ввод двоичного представления y
   cin.get();
   bitset<12> b = div(x,y);
   cout<<b<<endl;  // вывод двоичного представления b
   cin.get();
   return 0;
}

#include <iostream>
#include <cassert>
 
unsigned poly_div_23_by_12(unsigned a, unsigned b)
{
    if (b == 0)
        return 0;
    assert(b > 0);
    assert((a & 0x7FFFFF) == a); // a - 23-bit value
    assert((b & 0xFFF) == b); // b - 12-bit value
    
    int cnt = 0;
    unsigned mask = 1 << 22;
    while (not (b & mask))
    {
        b <<= 1;
        ++cnt;
    }
 
    while (cnt-- >= 0)
    {
        if (a & mask)
            a ^= b;
        b >>= 1;
        mask >>= 1;
    }
    
    return a;
}
 
 
int main()
{
    unsigned a = 4837261;
    unsigned b = 3189;
 
    std::cout << poly_div_23_by_12(a, b) << std::endl;
}

#include <bitset>
#include <iostream>
#include <cassert>
 
struct poly_div_result
{
    unsigned quotient;
    unsigned remainder;
};
 
poly_div_result poly_div_23_by_12(unsigned a, unsigned b)
{
    poly_div_result res;
    res.quotient = 0;
    res.remainder = 0;
    if (b == 0)
        return res;
    assert(b > 0);
    assert((a & 0x7FFFFF) == a); // a - 23-bit value
    assert((b & 0xFFF) == b); // b - 12-bit value
 
    int cnt = 0;
    unsigned mask = b;
    while ((mask & mask-1) > 0)
    {
        mask &= mask-1;
    }
 
    while (2*mask <= a)
    {
        mask <<= 1;
        b <<= 1;
        ++cnt;
    }
 
    while (cnt-- >= 0)
    {
        res.quotient <<= 1;
        if (a & mask)
        {
            a ^= b;
            res.quotient += 1;
        }
        b >>= 1;
        mask >>= 1;
    }
    res.remainder = a;
    return res;
}
 
 
int main()
{
    unsigned a = 0b10010011100111110001101; // 4837261
    unsigned b = 0b110001110101; // 3189
 
    poly_div_result r = poly_div_23_by_12(a, b);
    std::cout << "quotient:  0b" << std::bitset<23>(r.quotient) << std::endl;
    std::cout << "remainder: 0b" << std::bitset<12>(r.remainder) << std::endl;
}

unsigned poly_div(unsigned a, unsigned b)
{
    if (b == 0)
        return 0;
    assert(b > 0);
    assert((a & 0x7FFFFF) == a); // a - 23-bit value
    assert((b & 0xFFF) == b); // b - 12-bit value
    
    int cnt = 0;
    unsigned mask = 1 << 22;
    while (not (b & mask))
    {
        b <<= 1;
        ++cnt;
    }
 
    unsigned quotient = 0;
    while (cnt-- >= 0)
    {
        quotient <<= 1;
        if (a & mask)
        {
            a ^= b;
            quotient += 1;
        }
        b >>= 1;
        mask >>= 1;
    }
    
    return quotient;
}