#include <algorithm>
#include <bitset>
#include <cassert>
#include <cctype>
#include <cmath>
#include <cstdint>
#include <cstdio>
#include <cstring>
#include <deque>
#include <fstream>
#include <functional>
#include <iostream>
#include <limits>
#include <map>
#include <numeric>
#include <queue>
#include <set>
#include <vector>
using namespace std;
using ll = long long;
#define fst first
#define snd second

/* clang-format off */
template <class T, size_t D> struct _vec { using type = vector<typename _vec<T, D - 1>::type>; };
template <class T> struct _vec<T, 0> { using type = T; };
template <class T, size_t D> using vec = typename _vec<T, D>::type;
template <class T> vector<T> make_v(size_t size, const T& init) { return vector<T>(size, init); }
template <class... Ts> auto make_v(size_t size, Ts... rest) { return vector<decltype(make_v(rest...))>(size, make_v(rest...)); }
template <class T> inline void chmin(T &a, const T& b) { if (b < a) a = b; }
template <class T> inline void chmax(T &a, const T& b) { if (b > a) a = b; }
/* clang-format on */

template <std::uint_fast64_t Modulus>
class modint {
  using u64 = std::uint_fast64_t;

 public:
  u64 a;

  constexpr modint(const u64 x = 0) noexcept
      : a(x % Modulus) {
  }
  constexpr u64 &value() noexcept {
    return a;
  }
  constexpr const u64 &value() const noexcept {
    return a;
  }
  constexpr modint operator+(const modint rhs) const noexcept {
    return modint(*this) += rhs;
  }
  constexpr modint operator-(const modint rhs) const noexcept {
    return modint(*this) -= rhs;
  }
  constexpr modint operator*(const modint rhs) const noexcept {
    return modint(*this) *= rhs;
  }
  constexpr modint operator/(const modint rhs) const noexcept {
    return modint(*this) /= rhs;
  }
  constexpr modint &operator+=(const modint rhs) noexcept {
    a += rhs.a;
    if (a >= Modulus) {
      a -= Modulus;
    }
    return *this;
  }
  constexpr modint &operator-=(const modint rhs) noexcept {
    if (a < rhs.a) {
      a += Modulus;
    }
    a -= rhs.a;
    return *this;
  }
  constexpr modint &operator*=(const modint rhs) noexcept {
    a = a * rhs.a % Modulus;
    return *this;
  }
  constexpr modint &operator/=(const modint rhs) noexcept {
    return *this *= ~rhs;
  }
  constexpr modint operator^(u64 exp) const noexcept {
    modint v = 1, x = *this;
    while (exp) {
      if (exp & 1) {
        v *= x;
      }
      x *= x;
      exp >>= 1;
    }
    return v;
  }
  constexpr modint operator~() const noexcept {
    return *this ^ (Modulus - 2);
  }
};

using mint = modint<1000000007>;

struct DivisorRange {
  ll divisor;
  ll begin, end;  // [begin, end)
};

// [begin, end)
vector<DivisorRange> splitIntoDivisorRanges(ll N, ll begin, ll end) {
  vector<DivisorRange> res;
  for (ll x = begin; x < end;) {
    ll divisor = N / x;
    ll e;
    if (divisor == 0) {
      e = end;
    } else {
      e = min(end, N / divisor + 1);
    }
    res.push_back({divisor, x, e});
    x = e;
  }
  return res;
}

int main() {
#ifdef DEBUG
  ifstream ifs("in.txt");
  cin.rdbuf(ifs.rdbuf());
#endif
  int N, K;
  while (cin >> N >> K) {
    auto ranges = splitIntoDivisorRanges(N, 1, N + 1);
    int M = ranges.size();
    vector<mint> dp(M, mint(0));
    dp[0] = 1;
    for (int t = 1; t <= K; t++) {
      vector<mint> sum(M + 1, mint(0));
      for (int i = 0; i < M; i++) sum[i + 1] = sum[i] + dp[i] * (ranges[i].end - ranges[i].begin);
      vector<mint> ndp(M, mint(0));
      int j = 0;
      for (int i = M - 1; i >= 0; i--) {
        while (j < M && ranges[j].end <= ranges[i].divisor) {
          ++j;
        }
        if (j < M) {
          ndp[i] = sum[j] + dp[j] * (ranges[i].divisor - ranges[j].begin + 1);
        } else {
          ndp[i] = sum[M];
        }
      }
      dp = ndp;
    }
    mint res = 0;
    for (int i = 0; i < M; i++) res += dp[i] * (ranges[i].end - ranges[i].begin);
    cout << res.a << endl;
  }
  return 0;
}


Measured memory_bytes: 3188