Problem: Kavi has 2n points lying on the OX axis, i-th of which is located at x = i.

Kavi considers all ways to split these 2n points into n pairs. Among those, he is interested in good pairings, which are defined as follows:

Consider n segments with ends at the points in correspondent pairs. The pairing is called good, if for every 2 different segments A and B among those, at least one of the following holds:

  * One of the segments A and B lies completely inside the other. 
  * A and B have the same length. 



Consider the following example:

<image>

A is a good pairing since the red segment lies completely inside the blue segment.

B is a good pairing since the red and the blue segment have the same length.

C is not a good pairing since none of the red or blue segments lies inside the other, neither do they have the same size.

Kavi is interested in the number of good pairings, so he wants you to find it for him. As the result can be large, find this number modulo 998244353.

Two pairings are called different, if some two points are in one pair in some pairing and in different pairs in another.

Input

The single line of the input contains a single integer n (1≤ n ≤ 10^6).

Output

Print the number of good pairings modulo 998244353.

Examples

Input


1


Output


1

Input


2


Output


3

Input


3


Output


6

Input


100


Output


688750769

Note

The good pairings for the second example are: 

<image>

In the third example, the good pairings are: 

<image>

Solution (in CPP):

Solution:
// g++-10 a.cpp -o a && ./a  use endl in interactive
#include <bits/stdc++.h>
    using namespace std;
    // #include <ext/pb_ds/assoc_container.hpp>
    // #include <ext/pb_ds/tree_policy.hpp>
    // using namespace __gnu_pbds;
    // template <class T>
    // using oset = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;
    /*
        to use random numbers use first line random idk the use second line main stuff
        mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
        int rnd = uniform_int_distribution<int>(l,r)(rng);
    */
    #define int long long 
    #define float long double
    #define pb push_back
    #define mp make_pair
    #define F first
    #define S second
    #define B begin()
    #define E end()
    #define all(x) x.begin(), x.end()
    #define sortall(x) sort(all(x))
    #define rall(x) x.rbegin(), x.rend()
    #define revsortall(x) sort(rall(x))
    #define clr(x) memset(x, 0, sizeof(x))
    #define clr_dp(x) memset(x, -1, sizeof(x))
    #define loop(n) for(int i=0;i<n;i++)
    #define revloop(a,n) for(int i=n-1;i>=a;i--)
    #define Loop(k, n) for(int i=k;k<n?i<n:i>=n;k<n?i+=1:i-=1)
    #define len length()
    #define si size()
    #define tr(x) for(auto it = x.begin(); it != x.end(); it++)
    #define trrev(x) for(auto it = x.rbegin(); it != x.rend(); it++)
    typedef pair<int, int> pii;
    typedef vector<int> vi;
    typedef vector<pii> vpii;
    typedef vector<vi> vvi;
    const int mod = 1000000007;
    const int inf=LLONG_MAX;
    void IO(){
        ios_base::sync_with_stdio(0), cin.tie(0), cout.tie(0);
        cout.setf(ios::fixed);
        srand(chrono::high_resolution_clock::now().time_since_epoch().count());
            #ifndef ONLINE_JUDGE
        freopen("input.txt", "r", stdin); 
        freopen("output.txt", "w", stdout);
            #endif
    }
    int pw(int x, int y, int p=inf) {
        int res = 1;
        x = x % p;
        if(x == 0) 
            return 0;

        while(y > 0) {
            if(y & 1)
                res = (res * x) % p;

            y = y >> 1;
            x = (x * x) % p;
        }
        return res;
    }
    int invmod(int a, int m) {
        return pw(a, m - 2, m);
    }

    int cl(int a, int x) {
        return a % x == 0 ? a / x : a / x + 1;
    }
    int gcd(int a, int b)
    {
        if (a == 0)
            return b;
        return gcd(b % a, a);
    }

int m=998244353;

void run_time_terror()
{
    int dp[1000000+10];
    clr(dp);
    for (int i = 1; i <= 1000000+10; ++i) {
        for (int j = 2*i; j <= 1000000+10; j += i)
            dp[j]++;
    }
    int n;
    cin>>n;
    int s=1;
    dp[0]=s;
    for (int i = 1; i <=n; ++i)
    {
        dp[i]+=s%m;
        s+=dp[i];
        s=s%m;
    }
    cout<<dp[n];
}


int32_t main(){
    cout << setprecision(0);
    IO();
    int tt=1;
    // cin>>tt;
    for(int case_no=1;case_no<=tt;case_no++){
        // cout<<"Case #"<<case_no<<": ";
        run_time_terror();
    }
    cerr<<(float)clock()/CLOCKS_PER_SEC*1000<<" ms"<<endl;
    return 0;
}