// >>> TEMPLATES #include using namespace std; using ll = long long; using ld = long double; using i32 = int32_t; using i64 = int64_t; using u32 = uint32_t; using u64 = uint64_t; #define int ll #define rep(i, n) for (int i = 0; i < (int)(n); i++) #define rep1(i, n) for (int i = 1; i <= (int)(n); i++) #define repR(i, n) for (int i = (int)(n)-1; i >= 0; i--) #define rep1R(i, n) for (int i = (int)(n); i >= 1; i--) #define loop(i, a, B) for (int i = a; i B; i++) #define loopR(i, a, B) for (int i = a; i B; i--) #define all(x) begin(x), end(x) #define allR(x) rbegin(x), rend(x) #define rng(x, l, r) begin(x) + (l), begin(x) + (r) #define pb push_back #define eb emplace_back #define fst first #define snd second template auto mp(A &&a, B &&b) { return make_pair(forward(a), forward(b)); } template auto mt(T&&... x) { return make_tuple(forward(x)...); } template auto constexpr inf_ = numeric_limits::max()/2-1; auto constexpr INF32 = inf_; auto constexpr INF64 = inf_; auto constexpr INF = inf_; #ifdef LOCAL #include "debug.hpp" #else #define dump(...) (void)(0) #define say(x) (void)(0) #define debug if (0) #endif template struct pque : priority_queue, Comp> { vector &data() { return this->c; } void clear() { this->c.clear(); } }; template using pque_max = pque>; template using pque_min = pque>; template ::value, int> = 0> ostream& operator<<(ostream& os, T const& a) { bool f = true; for (auto const& x : a) os << (f ? "" : " ") << x, f = false; return os; } template ::value, int> = 0> ostream& operator<<(ostream& os, const T (&a)[N]) { bool f = true; for (auto const& x : a) os << (f ? "" : " ") << x, f = false; return os; } template ())), class = typename enable_if::value>::type> istream& operator>>(istream& is, T &a) { for (auto& x : a) is >> x; return is; } template ostream& operator<<(ostream& os, pair const& p) { return os << p.first << " " << p.second; } template istream& operator>>(istream& is, pair& p) { return is >> p.first >> p.second; } struct IOSetup { IOSetup() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(15); } } iosetup; template struct FixPoint : private F { constexpr FixPoint(F&& f) : F(forward(f)) {} template constexpr auto operator()(T&&... x) const { return F::operator()(*this, forward(x)...); } }; struct MakeFixPoint { template constexpr auto operator|(F&& f) const { return FixPoint(forward(f)); } }; #define MFP MakeFixPoint()| #define def(name, ...) auto name = MFP [&](auto &&name, __VA_ARGS__) template struct vec_impl { using type = vector::type>; template static type make_v(size_t n, U&&... x) { return type(n, vec_impl::make_v(forward(x)...)); } }; template struct vec_impl { using type = T; static type make_v(T const& x = {}) { return x; } }; template using vec = typename vec_impl::type; template auto make_v(Args&&... args) { return vec_impl::make_v(forward(args)...); } template void quit(T const& x) { cout << x << endl; exit(0); } template constexpr bool chmin(T& x, U const& y) { if (x > y) { x = y; return true; } return false; } template constexpr bool chmax(T& x, U const& y) { if (x < y) { x = y; return true; } return false; } template constexpr auto sumof(It b, It e) { return accumulate(b, e, typename iterator_traits::value_type{}); } template int sz(T const& x) { return x.size(); } template int lbd(C const& v, T const& x) { return lower_bound(begin(v), end(v), x)-begin(v); } template int ubd(C const& v, T const& x) { return upper_bound(begin(v), end(v), x)-begin(v); } int64_t mod(int64_t x, int64_t m) { assert(m != 0); return (x %= m) < 0 ? x+m : x; } const int dx[] = { 1, 0, -1, 0, 1, -1, -1, 1 }; const int dy[] = { 0, 1, 0, -1, 1, 1, -1, -1 }; constexpr int popcnt(ll x) { return __builtin_popcountll(x); } template pair operator+(pair const& p, pair const& q) { return { p.first + q.first, p.second + q.second }; } template pair operator-(pair const& p, pair const& q) { return { p.first - q.first, p.second - q.second }; } mt19937_64 seed_{random_device{}()}; template Int rand(Int a, Int b) { return uniform_int_distribution(a, b)(seed_); } i64 irand(i64 a, i64 b) { return rand(a, b); } // [a, b] u64 urand(u64 a, u64 b) { return rand(a, b); } // template void shuffle(It l, It r) { shuffle(l, r, seed_); } // <<< constexpr int mul(int x, int y) { if (x == 0 || y == 0) return 0; int sgn = 1; if (x < 0) sgn = -sgn, x = -x; if (y < 0) sgn = -sgn, y = -y; int p = (x > INF/y ? INF : x*y); if (sgn < 0) p = -p; return p; } constexpr int add(int x, int y) { x = clamp(x, -INF, +INF); y = clamp(y, -INF, +INF); return clamp(x + y, -INF, +INF); } int32_t main() { int n, k, m; cin >> n >> k >> m; map cnt; rep1 (a, int(1e6)) { int prod = mul(a, a+k); for (int b = 2; ; b++) { int x = add(a, mul(b, k)); prod = mul(prod, x); if (prod > n) break; cnt[prod]++; } } dump(cnt); auto g = [&](int x) { int l = 0, r = INF32; while (l+1 < r) { auto m = (l+r)/2; (mul(m, m+k) <= x ? l : r) = m; } return mp(l, mul(l, l+k)); }; if (m == 1) { int ans = 0; ans += g(n).fst; dump(ans); for (auto [x, num] : cnt) { if (g(x).snd == x) { ans--; num++; } if (num == m) ans++; } cout << ans << "\n"; } else { int ans = 0; for (auto [x, num] : cnt) { if (num == m or num == m-1) { if (g(x).snd == x) num++; if (num == m) ans++; } } cout << ans << "\n"; } }