radix packing

radix.cpp

template<class U>
constexpr U static_pow(U base, U exponent)
{
	return exponent == U(0) ? U(1) : base * static_pow(base, exponent - 1);
}

template<size_t N, size_t I, class U>
U radix_put(U w, U x)
{
	constexpr U M = static_pow(N, I), MN = M * N;
	return w % M + x * M + (w / MN) * MN;
}

template<size_t N, size_t I, class U>
U radix_get(U w)
{
	return (w / static_pow(N, I)) % N;
}