#include <cstdio>
#include <vector>
using namespace std;
const int maxn=1e5+10;
const int Mod=998244353;
int dp[maxn][31];
int jc[maxn];
int ny[maxn];
int jcny[maxn];
vector<int> ve[maxn];
int C(int x,int y){
	if(x<y||y<0) return 0;
	return 1ll*jc[x]*jcny[y]%Mod*jcny[x-y]%Mod;
}
int main(){
	int m,n;
	scanf("%d%d",&m,&n);
	jc[0]=jc[1]=jcny[0]=jcny[1]=ny[1]=1;
	for(int i=2;i<=m;i++){
		jc[i]=1ll*jc[i-1]*i%Mod;
		ny[i]=1ll*(Mod-Mod/i)*ny[Mod%i]%Mod;
		jcny[i]=1ll*jcny[i-1]*ny[i]%Mod;
	}
	for(int i=1;i<=m;i++)
		for(int j=i;j<=m;j+=i)
			ve[j].push_back(i);
	int ans=1ll*C(m,n)*n%Mod;
	for(int i=1;i<=m;i++){
		dp[i][1]=1;
		for(int j:ve[i]){
			if(i!=j){
				for(int k=2;k<=30;k++)
					dp[i][k]=(dp[i][k]+dp[j][k-1])%Mod;
			}
		}
		for(int k=1;k<=30;k++){
			if(n<k) continue;
			ans=(ans-1ll*C(m/i-1,n-k)*dp[i][k]+Mod)%Mod;
		}
	}
	printf("%d",ans);
	return 0;
}