Ackermann function

#include <cstddef>

template< std::size_t m, std::size_t n >
constexpr std::size_t a = a< (m - 1), a< m, (n - 1) > >;

template< std::size_t m >
constexpr std::size_t a< m, 0 > = a< (m - 1), 1 >;

template< std::size_t n >
constexpr std::size_t a< 0, n > = (n + 1);

static_assert(a< 5, 3 > == 253);