uint64_t get_byte(size_t num, int ord_num, int bytes_to_take){
    return ((unsigned long long)num << ((sizeof (num) - bytes_to_take - ord_num * bytes_to_take) * 8)) >> (sizeof (num) * 8 - 8 * bytes_to_take);
}

void Radix_sort(uint64_t Temp[], int length, int bytes_to_take) {
     int cnt_bytes = sizeof (Temp[0]) / bytes_to_take;
     //std::cout << " bytes " << cnt_bytes << "\n";
     int cnt = pow(2, bytes_to_take * 8);
     uint64_t* Temp2 = new uint64_t[length];
     for (int byte = 0; byte < cnt_bytes; ++byte) {
         std::vector<int> byte_var(cnt,0);
         //std::cout << "byte_num: " << byte << "\n";
         for (int iter = 0; iter < length; ++iter) {
             byte_var[get_byte(Temp[iter], byte, bytes_to_take)]++;
             //std::cout << get_byte(Temp[iter], byte, bytes_to_take) << " " << Temp[iter] << "\n";
         }

         int count = 0;
         int temp;
         for (int iter = 0; iter < cnt; ++iter) {
             temp = byte_var[ iter ];
             byte_var[ iter ] = count;
             count += temp;
         }
         //std::cout<<"swaped\n";
         for (int iter = 0; iter < length; ++iter) {
             temp = get_byte(Temp[ iter ], byte, bytes_to_take);
             Temp2[byte_var[temp]] = Temp[ iter ];
             byte_var[ temp ]++;
         }
         std::swap(Temp, Temp2);
         //std::cout<<"swaped\n";
     }
     return;
}