#include<bits/stdc++.h>
using namespace std;
#define endl '\n'
typedef long long LL;

const int N=100010,mod=998244353;

int primes[N],minp[N],idx;
bool st[N];

void is_primes(int n)
{
	for(int i=2;i<=n;i++)
	{
		if(!st[i])primes[idx++]=i,minp[i]=i;
		for(int j=0;i*primes[j]<=n;j++)
		{
			st[i*primes[j]]=true;
			minp[i*primes[j]]=primes[j];
			if(i%primes[j]==0)break;
		}
	}
}

int fact[N],infact[N];//fact表示阶乘，infact表示阶乘的逆元

int qmi(int a,int b)
{
    int res=1;
    while(b)
    {
        if(b&1)res=(LL)res*a%mod;
        a=(LL)a*a%mod;
        b>>=1;
    }
    return res;
}

void init()
{
	fact[0]=infact[0]=1;
    for(int i=1;i<N;i++)
    {
        fact[i]=(LL)fact[i-1]*i%mod;
        infact[i]=(LL)infact[i-1]*qmi(i,mod-2)%mod;
    }
}

LL query(int a,int b)
{
	LL res=(LL)fact[a]*infact[b]%mod*infact[a-b]%mod;
	return res;
}

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

	LL res=0;
	for(int i=1;i<=n;i++)
	{
		map<int,int>mp;
		for(int j=i;j<=n;j+=i)
		{
			int val=j/i;
			while(val!=minp[val]&&val!=1)
			{
				int t=minp[val];
				mp[t]++;
				while(val%t==0)val/=t;
			}
			if(val!=1)mp[val]++;
		} 

		LL t=0;
		if(n/i-1>=m)t=query(n/i-1,m);
		for(auto [u,v]:mp)
		{
			if(v>=m)t-=query(v,m);
		}
		res=(res+t*m%mod)%mod;

		int cnt=0;
		for(int j=1;j*j<=i;j++)
		{
			if(i%j==0)
			{
				if(j!=i/j)cnt+=2;
				else cnt++;
			}
		}

		for(int j=1;j<=cnt;j++)
		{
			LL val=0;
			if(n/i-1>=m-j)val=query(n/i-1,m-j);
			for(auto [u,v]:mp)
			{
				if(v>=m-j)val-=query(v,m-j);
			}
			res=(res+query(cnt,j)*val%mod*(m-j)%mod)%mod;
		}
	}
	cout<<res<<endl;
}

int main()
{
	ios::sync_with_stdio(0);
	cin.tie(0),cout.tie(0);

	int T;
	T=1;
	init();
	is_primes(1e5);
	while(T--)
	{
		solve();
	}
	return 0;
}