#include <bits/stdc++.h>
using namespace std;
#define int long long
const int N = 1e6 + 10;
#define pii pair<int, int>
#define lowbit(x) (x & (-x))
const int mod = 998244353;
int fact[N];
int infact[N];
int qmi(int a, int b)
{
    a %= mod;
    int res = 1;
    while (b)
    {
        if (b & 1)
            res = res * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return res;
}
void init()
{
    fact[0] = 1;
    for (int i = 1; i < N; i++)
        fact[i] = fact[i - 1] * i % mod;
    infact[N - 1] = qmi(fact[N - 1], mod - 2);
    for (int i = N - 2; i >= 0; i--)
        infact[i] = infact[i + 1] * (i + 1) % mod;
}
int C(int n, int m)
{
    if (m < 0 || m > n)
        return 0;
    return fact[n] * infact[n - m] % mod * infact[m] % mod;
}
void solve()
{
    int n, m;
    cin >> m >> n;
    unordered_map<int, unordered_map<int, int>> f;
    int res = C(m, n) % mod * n % mod;
    for (int i = 1; i <= m; i++)
        f[1][i] = 1;
    for (int i = 1; i <= n; i++)
    {
        for (auto t : f[i])
        {
            for (int j = 2; j * t.first <= m; j++)
                (f[i + 1][j * t.first] += f[i][t.first]) %= mod;
        }
    }
    for (int i = 1; i <= n; i++)
    {
        for (auto [x, y] : f[i])
        {
            res -= y % mod * C(m / x - 1, n - i) % mod;
            res %= mod;
        }
    }
    cout << (res + mod) % mod << '\n';
}
signed main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int t = 1;
    init();
    // cin >> t;
    while (t--)
        solve();
}