#include <bits/stdc++.h>
// #include <boost/multiprecision/cpp_int.hpp>
using namespace std;

#define all(a) (a).begin(), (a).end()
#define rall(a) (a).rbegin(), (a).rend()
#define MIN(a) *std::min_element((a).begin(), (a).end())
#define MAX(a) *std::max_element((a).begin(), (a).end())
#define SUM(a) std::accumulate((a).begin(), (a).end(), 0LL)
#define REP(i, n) for (std::size_t (i) = 0; (i) < (n); (i)++)
#define RREP(i, n) for (std::size_t (i) = (n)-1; (i) >= 0; (i)--)
#define FOR(i, m, n) for (std::size_t (i) = (m); (i) < (n); i++)
#define FORR(i, m, n) for (std::size_t (i) = (n)-1; (i) >= (m); i--)
#define debug(x) std::cerr << #x << " = " << x << '\n'
// #define int long long

using ll = long long;
using ull = unsigned long long;
using ld = long double;
using vl = std::vector<long long>;
using vvl = std::vector<std::vector<long long>>;
// using bint = boost::multiprecision::cpp_int;
// using rat = boost::rational<bint>;

constexpr int inf = (1 << 30) - 1;
constexpr long long INF = 1LL << 60;

inline void ios_() {std::cin.tie(nullptr); std::ios::sync_with_stdio(false);}
template<typename T> inline T Div(const T& a, const T& b) {return (a + b - 1) / b;}
template<typename T> bool chmin(T& a, const T& b) {if(a > b){a = b; return true;} return false;}
template<typename T> bool chmax(T& a, const T& b) {if(a < b){a = b; return true;} return false;}

template<typename T>
inline void finish(const T& msg) {
    std::cout << msg << std::endl;
    std::exit(0);
}

template<typename T>
long long binarySearch(long long ng, long long ok, const T& fn) {
    while (std::abs(ok - ng) > 1) {
        T mid = (ok + ng) / 2;
        if (fn(mid)) ok = mid;
        else ng = mid;
    }
    return ok;
}

struct ModInt {
    long long value;
    static const long long mod; // mod must be a prime number

    ModInt(long long v = 0) : value(((v % mod) + mod) % mod) {}

    ModInt pow(long long n) const {
        long long x = value, res = 1;
        while(n > 0) {
            if(n & 1) {
                res = res * x % mod;
            }
            x = x * x % mod;
            n >>= 1;
        }
        return ModInt(res);
    }

    ModInt inverse() const {
        return pow(mod - 2);
    }

    ModInt operator-() const {
        return ModInt((mod - value) % mod);
    }

    ModInt& operator+=(const ModInt& a) {
        value += a.value;
        value %= mod;
        return *this;
    }

    ModInt& operator-=(const ModInt& a) {
        value += (-a).value;
        value %= mod;
        return *this;
    }

    ModInt& operator*=(const ModInt& a) {
        value *= a.value;
        value %= mod;
        return *this;
    }

    ModInt& operator/=(const ModInt& a) {
        return (*this) *= a.inverse();
    }

    ModInt operator+(const ModInt& a) const {
        return ModInt(*this) += a;
    }

    ModInt operator-(const ModInt& a) const {
        return ModInt(*this) -= a;
    }

    ModInt operator*(const ModInt& a) const {
        return ModInt(*this) *= a;
    }

    ModInt operator/(const ModInt& a) const {
        return ModInt(*this) /= a;
    }

    bool operator==(const ModInt& a) const {
        return value == a.value;
    }

    bool operator!=(const ModInt& a) const {
        return value != a.value;
    }

    friend std::istream& operator>>(std::istream& is, ModInt& a) {
        long long v;
        is >> v;
        a = v;
        return is;
    }

    friend std::ostream& operator<<(std::ostream& os, const ModInt& a) {
        os << a.value;
        return os;
    }

    static ModInt factorial(int n) {
        static std::vector<ModInt> fact(1 << 18);
        if(fact.size() <= n) {
            fact.resize(n + 1);
        }
        if(n < 0) {
            return 0;
        }
        if(n == 0) {
            return 1;
        }
        if(fact[n] != 0) {
            return fact[n];
        }
        return fact[n] = factorial(n - 1) * n;
    }

    static ModInt combin(int n, int k) {
        if(n < 0 || k < 0 || n < k) {
            return 0;
        }
        return factorial(n) / (factorial(k) * factorial(n - k));
    }
};
const long long ModInt::mod = 1e9+7;
// const long long ModInt::mod = 998244353;
using mint = ModInt;



signed main() {
    ll n, k;
    cin >> n >> k;

    ll m = 1;
    while((m+1)*(m+1) <= n) m++;

    vector<vector<vector<mint>>> dp(2, vector<vector<mint>>(k, vector<mint>(m+2, 0)));

    FOR(i, 1, m+1){
        dp[0][0][i] = 1;
        dp[1][0][i] = n/i - max(m, (ll)(n/(i+1)));
    }

    REP(i, k-1){
        FOR(j, 1, m+1){
            dp[0][i+1][1] += dp[0][i][j];
            dp[0][i+1][m+1] -= dp[0][i][j];
            dp[1][i+1][j] += dp[0][i][j];
            dp[1][i+1][m+1] -= dp[0][i][j];

            dp[0][i+1][1] += dp[1][i][j];
            dp[0][i+1][j+1] -= dp[1][i][j];
        }

        REP(j, m+1){
            dp[0][i+1][j+1] += dp[0][i+1][j];
            dp[1][i+1][j+1] += dp[1][i+1][j];
        }

        REP(j, m+2){
            dp[1][i+1][j] *= dp[1][0][j];
        }
    }

    mint ans = accumulate(all(dp[0][k-1]), mint(0)) + accumulate(all(dp[1][k-1]), mint(0));
    cout << ans << endl;

    return 0;
}


Measured memory_bytes: 74880