Реализация Random(a,b) вызовами Random(0,1) (задача 5.1-2 из Кормен изд.2)

typedef unsigned __int8 uint_small;
typedef unsigned __int64 uint_large;
 
typedef std::map<uint_small, int> distribution_map_t;
typedef distribution_map_t::iterator iter_distribution_map_t;
 
bool Random01()
{
    return rand() > (RAND_MAX / 2);
}
 
uint_small Random(uint_small a, uint_small b)
{
    assert( b > a );
    unsigned __int64 result = 0;
 
    ++b;
 
    for ( int j = 0; j < sizeof(uint_large) * 8; ++j )
    {
        result <<= 1;
        if ( Random01() )
            result += 1;
    }
 
    //if ( result == 0 ) // Не наблюдал, поэтому и двигаю ++b, --b и  a + res - 1;
        //return a;
 
    uint_small res = (uint_small)(result % (uint_large)(b - a));
 
    if ( res == 0 )
        return --b;
 
    return a + res - 1;
}
 
int main(int argc, char* argv[])
{
    uint_small uint_small_max = 0;
    memset((void*)&uint_small_max, 0xFF, sizeof(uint_small));
 
    distribution_map_t mp;
 
    while ( true )
    {
        int rand1 = rand() % uint_small_max;
        int rand2 = rand() % uint_small_max;
 
        assert(rand1 >= 0 && rand1 < uint_small_max);
        assert(rand2 >= 0 && rand2 < uint_small_max);
 
        if ( rand1 == rand2 )
            continue;
 
        if ( rand1 > rand2 )
            std::swap(rand1, rand2);
 
        mp.clear();
 
        {
            for ( int i = rand1; i <= rand2; ++i )
                mp.insert(std::make_pair(i, 0));
        }
 
        printf("range: %d...%d\n", rand1, rand2);
 
        const int howmanycalls = (rand2 - rand1) * 100; // это неважно, чем больше, тем точнее
 
        for ( int i = 0; i < howmanycalls; ++i )
        {
            uint_small res = Random((uint_small)rand1, (uint_small)rand2);
            assert(res >= (uint_small)rand1);
            assert(res <= (uint_small)rand2);
 
            iter_distribution_map_t it = mp.find(res);
            if ( it != mp.end() )
                ++(it->second);
            else
                assert(false);
        }
 
        { // Вывод результатов распределения
            iter_distribution_map_t it = mp.begin();
            iter_distribution_map_t itEnd = mp.end();
 
            for ( ; it != itEnd; ++it )
                printf("num: %d, count: %d\n", it->first, it->second);  
        }
    } // while ( true )
 
    return 0;
}