#include<bits/stdc++.h>
using namespace std;
#define int long long
using pii=pair<int,int>;

const int N=100005,M=100000;
const int mod=998244353;

int n,m;
int f[22][N],cnt[N];

int fac[N],infac[N];

int quick_pow(int a, int k, int p)
{
    int res = 1;
    while (k)
    {
        if (k & 1) res = res * a % p;
        a = a * a % p;
        k >>= 1;
    }
    return res;
}

int cc(int a,int b){
    if(b>a) return 0;
    return fac[a] * infac[b] % mod * infac[a - b] % mod;
}

void pre(){
     fac[0]=infac[0]=1;
     for(int i=1;i<=M;i++)
     {
         fac[i]=fac[i-1]*i%mod;
         infac[i]=infac[i - 1] * quick_pow(i,mod-2,mod)%mod;
     }

     //求cnt[i]

     for(int i=1; i<=M; i++) {
         int w=i;
         for(int j=2; j<=w/j; j++) {
             while(w%j==0) {
                 w/=j;
                 cnt[i]++;
             }
         }
         if(w>1) cnt[i]++;

         cnt[i]++;

     }
     cnt[1]=1;

     // 求f[i][j]
      for(int j=1; j<=M; j++) {
          // cout<<j<<":";
          for(int i=1; i<=cnt[j]; i++) {
              f[i][j]=cc(cnt[j]-1,i-1);
              // cout<<f[i][j]<<" \n"[i==cnt[j]+1];
          }
      }
}

//f[i][j] 表示长度为i，当前数为j的方案数

void solve(){
    cin>>m>>n;

    int ans=0;

    for(int j=1; j<=m; j++) {
        for(int i=1; i<=cnt[j]; i++) {
            ans+=f[i][j]*cc(m/j-1,n-i);
            ans%=mod;
        }
    }

    // cout<<cc(m,n)*n<<" "<<ans<<endl;
    cout<<(cc(m,n)*n-ans)%mod;
}

signed main(){
    ios::sync_with_stdio(false),cin.tie(nullptr);

    pre();

    solve();

}