#include #pragma GCC diagnostic ignored "-Wsign-compare" #pragma GCC diagnostic ignored "-Wsign-conversion" #define NDEBUG #define SHOW(...) static_cast(0) //!===========================================================!// //! dP dP dP !// //! 88 88 88 !// //! 88aaaaa88a .d8888b. .d8888b. .d888b88 .d8888b. 88d888b. !// //! 88 88 88ooood8 88' '88 88' '88 88ooood8 88' '88 !// //! 88 88 88. ... 88. .88 88. .88 88. ... 88 !// //! dP dP '88888P' '88888P8 '88888P8 '88888P' dP !// //!===========================================================!// template T read() { T v; return std::cin >> v, v; } template std::vector readVec(const std::size_t l) { std::vector v(l); for (auto& e : v) { std::cin >> e; } return v; } using ld = long double; using uint = unsigned int; using ll = long long; using ull = unsigned long long; constexpr unsigned int MOD = 1000000007; template constexpr T INF = std::numeric_limits::max() / 4; template constexpr F PI = static_cast(3.1415926535897932385); std::mt19937 mt{std::random_device{}()}; template bool chmin(T& a, const T& b) { return (a > b ? a = b, true : false); } template bool chmax(T& a, const T& b) { return (a < b ? a = b, true : false); } template std::vector Vec(const std::size_t n, T v) { return std::vector(n, v); } template auto Vec(const std::size_t n, Args... args) { return std::vector(n, Vec(args...)); } template constexpr T popCount(const T u) { #ifdef __has_builtin return u == 0 ? T(0) : (T)__builtin_popcountll(u); #else unsigned long long v = static_cast(u); return v = (v & 0x5555555555555555ULL) + (v >> 1 & 0x5555555555555555ULL), v = (v & 0x3333333333333333ULL) + (v >> 2 & 0x3333333333333333ULL), v = (v + (v >> 4)) & 0x0F0F0F0F0F0F0F0FULL, static_cast(v * 0x0101010101010101ULL >> 56 & 0x7f); #endif } template constexpr T log2p1(const T u) { #ifdef __has_builtin return u == 0 ? T(0) : T(64 - __builtin_clzll(u)); #else unsigned long long v = static_cast(u); return v = static_cast(v), v |= (v >> 1), v |= (v >> 2), v |= (v >> 4), v |= (v >> 8), v |= (v >> 16), v |= (v >> 32), popCount(v); #endif } template constexpr T clog(const T v) { return v == 0 ? T(0) : log2p1(v - 1); } template constexpr T msbp1(const T v) { return log2p1(v); } template constexpr T lsbp1(const T v) { #ifdef __has_builtin return __builtin_ffsll(v); #else return v == 0 ? T(0) : popCount((v & (-v)) - T(1)) + T(1); #endif } template constexpr bool ispow2(const T v) { return popCount(v) == 1; } template constexpr T ceil2(const T v) { return v == 0 ? T(1) : T(1) << log2p1(v - 1); } template constexpr T floor2(const T v) { return v == 0 ? T(0) : T(1) << (log2p1(v) - 1); } //!====================================================================================================!// //! 888888ba oo 88888888b dP !// //! 88 '8b 88 88 !// //! a88aaaa8P' 88d888b. dP 88d8b.d8b. .d8888b. a88aaaa d8888b. .d8888b. d8888P .d8888b. 88d888b. !// //! 88 88' '88 88 88''88''88 88ooood8 88 88' '88 88' '"" 88 88' '88 88' '88 !// //! 88 88 88 88 88 88 88. ... 88 88. .88 88. ... 88 88. .88 88 !// //! dP dP dP dP dP dP '88888P' dP '88888P8 '88888P' dP '88888P' dP !// //!====================================================================================================!// template std::vector> PrimeFactors(T n) { std::vector> ans; for (T i = 2; i * i <= n; i++) { std::size_t cnt = 0; for (; n % i == 0; n /= i, cnt++) {} if (cnt > 0) { ans.push_back({i, cnt}); } } if (n > 1) { ans.push_back({n, 1}); } return ans; } //!=====================================!// //! 8888ba.88ba oo !// //! 88 '8b '8b !// //! 88 88 88 .d8888b. dP 88d888b. !// //! 88 88 88 88' '88 88 88' '88 !// //! 88 88 88 88. .88 88 88 88 !// //! dP dP dP '88888P8 dP dP dP !// //!=====================================!// int main() { const int N = read(), K = read(), M = read(); auto fac = PrimeFactors(N); const int L = fac.size(); std::vector p(L); std::vector ind(L); std::size_t i = 0; for (auto& pr : fac) { pr.second *= K; ll q = 1; for (; q <= M; ind[i]++, q *= pr.first) {} p[i] = pr.first, ind[i]--, chmin(ind[i], pr.second), i++; } SHOW(p, ind); auto dfs = [&](auto&& self, const std::size_t i, const ll P) -> ll { if (i == L) { return 1LL; } ll prod = 1, ans = 0; for (std::size_t j = 0; j <= ind[i] and prod * P <= M; j++, prod *= p[i]) { ans += self(self, i + 1, prod * P); } return ans; }; std::cout << dfs(dfs, 0, 1) << std::endl; return 0; }