-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathE.cpp
83 lines (79 loc) · 1.69 KB
/
E.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
#include<bits/stdc++.h>
using namespace std;
#define fo(i,n) for(i=0;i<n;i++)
#define fo2(i,start,end) for(i=start;i<=end;i++)
#define rfo(i,n) for(i=n-1;i>=0;i--)
#define ll long long int
#define deb(x) cout<<#x<<"="<<x<<endl
#define deb2(x,y) cout<<#x<<"="<<x<<","<<#y<<"="<<y<<endl
#define pb push_back
#define mp make_pair
#define F first
#define S second
#define all(x) x.begin(), x.end()
#define sortall(x) sort(all(x))
#define tr(it,a) for(auto it=a.begin(); it!=a.end();it++)
#define nl '\n'
#define PI 3.1415926535897932384626
#define mod 998244353
typedef pair<int, int> pii;
typedef pair<ll, ll> pl;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef vector<pii> vpii;
typedef vector<pl> vpl;
const int N=1e6+10;
ll fact[N];
ll power(ll a,ll n){
ll res=1;
while(n){
if(n%2){
res=(res*a)%mod;
}
a=(a*a)%mod;
n/=2;
}
return res;
}
ll ncr(ll n,ll r){
ll ans,num,denom;
num=fact[n];
denom=(fact[r]*fact[n-r])%mod;
denom=power(denom,mod-2);
ans=(num*denom)%mod;
return ans;
}
void solve(){
int i,j;
int n,k;
cin>>n>>k;
// if(k>n){
// cout<<"0\n";
// return;
// }
ll count,ans=0,temp;
for(i=1;i<=n;i++){
count=n/i-1;
if(count<k-1){
break;
}
temp=ncr(count,k-1);
ans=(ans+temp)%mod;
}
cout<<ans<<nl;
}
int main(){
// freopen("input.txt", "r", stdin);
// freopen("output.txt", "w", stdout);
ios_base::sync_with_stdio(0), cin.tie(0), cout.tie(0);
int t=1;
fact[0]=1;
fact[1]=1;
for(int i=2;i<N;i++){
fact[i]=(fact[i-1]*i)%mod;
}
// cin>>t;
while(t--){
solve();
}
}