#include <bits/stdc++.h>
using namespace std;


#define ls(x) ((x) << 1)
#define rs(x) ((x) << 1 | 1)
#define lowbit(x) ((x) & (-x))

const int N = 10010007;
const int M = 120;
const int INF = (1ull << 60) - 1;
const int mod = 998244353;
using pii = pair<int, int>;
using ll = long long;

struct node{
    int l;int r;
    int num;
}b[N];
int a[N];
int p[N];
int ans[N];
void solve(){   
    int n;
    cin >> n;
    cin >> a[0];
    const int P = 1e9 + 7;
    for (int i = 1;i <= n;i++)a[i] = (1LL * a[i - 1] * mod + P) % (P + 2);
    for (int i = 1;i <= n;i++)p[i] = i;
    for (int i = 1;i <= n;i++)swap(p[i], p[a[i] % i + 1]);
    for (int i = 1;i <= n;i++){
        b[i].l = i - 1;
        b[i].r = i + 1;
        b[i].num = i;
    }

    int mid = (n + 1)/2;
    b[0].r = 1;
    b[n + 1].l = n;
    for (int i = n;i >= 1;i--){
        ans[i] = b[mid].num;
        //cout << ans[i] << '\n';
        int L = b[p[i]].l;
        int R = b[p[i]].r;
        b[L].r = R;
        b[R].l = L;
        if (i & 1){
            if (p[i] == b[mid].num)mid = b[mid].l;
            else if (p[i] > b[mid].num)mid = b[mid].l;
        }
        else{
            
            if (p[i] == b[mid].num)mid = b[mid].r;
            else if (p[i] < b[mid].num)mid = b[mid].r;
        }
    }
    ll sum = 0;
    ll res = 19;
    for (int i = 1;i <= n;i++){
        
        sum = (1LL * ans[i] * res % mod + 1LL * sum ) % mod;
        //cout << ans[i] << " " << sum << '\n';
        res *= 19;
        res %= mod;
    }
    cout << sum << '\n';
}
void init(int n,int k) {
    
}
signed main()
{
    ios::sync_with_stdio(false);
    cin.tie(0); cout.tie(0);
    int T = 1;
    //cin >> T;

    
    while (T--) {
        solve();
    }
    return 0;
}