#include <bits/stdc++.h>
#define dbg() cout<<"-------"
typedef  long long ll;
using namespace std;
typedef pair<int,int> PII;
typedef unsigned long long ull;
const int N=2e6+10;
const int mod=998244353;
ll n,m;
ll dp[N][20];
ll fact[N];
ll qmi(ll x,ll y,ll mod){
    ll res=1,t=x%mod;
    while(y){
        if(y&1) res=res*t%mod;
        t=t*t%mod;
        y>>=1;
    }
    return res;
}
ll c(int n,int m){
    if(m>n) return 0;
    return fact[n]*qmi(fact[m],mod-2,mod)%mod*qmi(fact[n-m],mod-2,mod)%mod;
    //return fact[n]*qmi(fact[m]*fact[n-m]%mod,mod-2,mod)%mod;
}
void solve(){
    cin>>m>>n;
    fact[0]=1;for(int i=1;i<=m;i++) fact[i]=fact[i-1]*i%mod;
    //cout<<c(3,0)<<'\n';
   //  cout<<"*****"<<'\n';
    for(int i=1;i<=m;i++) dp[i][1]=1;

    for(int len=1;len<=17;len++){
        for(int i=1;i<=m;i++){
            for(int j=2*i;j<=m;j+=i){
                dp[j][len+1]+=dp[i][len];
                dp[j][len+1]%=mod;
            }
        }
    }
   // cout<<dp[1][1]<<' '<<dp[1][2]<<'\n';
 //   cout<<"..."<<'\n';
    
    ll ans=0;
    for(int i=1;i<=m;i++){
        int cnt=m/i-1;
        ans+=(c(m-i,n-1)-c(cnt,n-1))*n%mod;
        //cout<<i<<"..."<<ans<<'\n';
        ans%=mod;
    }
    //cout<<ans<<'\n';
  //  exit(0);
  //ans=0;/////////////////
    for(int g=1;g<=m;g++){
        for(int len=1;len<=min(17ll,n-1);len++){
            if(!dp[g][len]) break;
            for(int t=2*g;t<=m;t+=g){
                int cnt=m/g-t/g;
                int cnt2=m/t-1;
                ans+=(c(cnt,n-len-1)-c(cnt2,n-len-1))%mod*dp[g][len]%mod*(n-len)%mod;
                // if((c(cnt,n-len-1)-c(cnt2,n-len-1))*dp[g][len]%mod*(n-len)%mod){
                //     cout<<g<<' '<<len<<" "<<t<<'\n';
                // }
                ans%=mod;
            }
        }
      //  cout<<g<<"..."<<ans<<'\n';
    }
    cout<<((ans%mod)+mod)%mod<<'\n';



}
int main(){
    
    ios::sync_with_stdio(false);
    cin.tie(0);
    int t;
    //cin>>t;
    t=1;
    while(t--){
        solve();
    }
}