#include <bits/stdc++.h>
using namespace std;
#define int long long
const int mod = 998244353;
int qmi(int a, int b)
{
    int res = 1;
    while (b)
    {
        if (b & 1)
            res = res * a % mod;
        b >>= 1;
        a = a * a % mod;
    }
    return res;
}
int prime[100010],cnt,phi[100010];
void euler(int n){
    for(int i=2;i<=n;i++){
        if(!phi[i]){
            prime[++cnt]=i;
            phi[i]=i;
        }
        for(int j=1;j<=cnt&&prime[j]<=n/i;j++){
            phi[i*prime[j]]=prime[j];
            if(i%prime[j]==0)break;
        }
    }
}
signed main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    euler(1e5);
    vector<int> f(1e5 + 10);
    vector<int> inv(1e5 + 10);
    f[0] = inv[0] = 1;
    for (int i = 1; i <= 1e5; i++)
        f[i] = f[i - 1] * i % mod;
    inv[(int)1e5] = qmi(f[(int)1e5], mod - 2);
    for(int i=1e5-1;i>=1;i--){
        inv[i]=inv[i+1]*(i+1)%mod;
    }
    int n, m;
    cin >> n >> m;
    auto c = [&](int a, int b) -> int
    {
        if (a < 0 || b < 0 || a < b)
            return 0;
        return f[a] * inv[b] % mod * inv[a - b] % mod;
    };
    vector<int> d(n + 1);
    vector<vector<int>>dp(n+1,vector<int>(60)); 
    for (int i = 1; i <= n; i++)
    {
        dp[i][1]=1;
        for (int j = 2*i; j <= n; j += i)
        {
            for(int k=2;k<=60;k++){
                dp[j][k]=(dp[j][k]+dp[i][k-1])%mod;
            }
        }
    }
    int ans = 0;
    for (int i = 1; i <= n; i++)
    {
        int tmp=0;
        for(int j=1;j<=60;j++){
            tmp=(tmp+dp[i][j]*c(n/i-1,m-j)%mod)%mod;
            ans=(ans+dp[i][j]*c(n/i-1,m-j)%mod)%mod;
        }
    }
    ans=(m*c(n,m)%mod-ans)%mod;
    ans = (ans % mod + mod) % mod;
    cout << ans << '\n';
    // vector<int>cnt(12);
    // for(int i=0;i<(1<<n);i++){
    //     vector<int>tmp;
    //     for(int j=0;j<n;j++){
    //         if(i>>j&1){
    //             tmp.push_back(j+1);
    //         }
    //     }
    //     if(tmp.size()!=m)continue;
    //     vector<int>gcd(tmp.size());
    //     gcd[tmp.size()-1]=tmp.back();
    //     for(int i=tmp.size()-2;i>=0;i--){
    //         gcd[i]=__gcd(gcd[i+1],tmp[i]);
    //     }
    //     for(int i=0;i<(int)tmp.size();i++){
    //         if(tmp[i]==gcd[i]){
    //             cnt[tmp[i]]++;
    //         }
    //         else break;
    //     }
    // }
}