#include <bits/stdc++.h>
#define int long long
using namespace std;
#define YES "YES"
#define NO "NO"
#define all(a) a.begin(), a.end()
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
typedef pair<int, int> pii;
const double eps = 1e-9;
const int N = 3e6 + 10;
const int M = 2e6 + 10;
const int INF = 1e16 + 100;
const int mod = 998244353;
const int base = 23333;
mt19937 rng(std::chrono::steady_clock::now().time_since_epoch().count());
int n, cnt1, cnt2, ans;
string s;
int d[N + N / 2];
int nxt[M][25];
int lef[N + N / 2];
struct node
{
    int son[10], num = 0, add = 0, h = 0;
} trie[M];
void dfs1(int u)
{
    for (int i = 0; i < 10; i++)
    {
        int v = trie[u].son[i];
        if (v)
        {
            dfs1(v);
            trie[u].add = (trie[u].add + trie[v].add) % mod;
        }
        trie[u].son[i] = 0;
    }
    trie[u].num = (trie[u].num + trie[u].add) % mod;
    if (trie[u].h)
        ans = (ans + trie[u].num * trie[u].num % mod * (2 * trie[u].h - 1) % mod) % mod;
}
void dfs2(int u)
{
    for (int i = 0; i < 10; i++)
    {
        int v = trie[u].son[i];
        if (v)
        {
            dfs2(v);
            trie[u].add = (trie[u].add + trie[v].add) % mod;
        }
    }
    trie[u].num = (trie[u].num + trie[u].add) % mod;

    ans = (ans + trie[u].num * trie[u].num % mod * (2 * trie[u].h) % mod) % mod;
}
int get(int u, int k)
{
    int i = 0;
    while (k)
    {
        if (k & 1)
            u = nxt[u][i];
        k = k >> 1;
        i++;
    }
    return u;
}
void solve()
{
    cin >> n;
    cin >> s;
    s = s + s;
    // 奇数
    for (int i = 0, l = 0, r = -1; i < n + n / 2; i++)
    {
        int k = 0, pos = 0;
        if (i <= r)
        {
            k = min(d[l + r - i], r - i + 1);
            pos = lef[l + r - i];
            if (d[l + r - i] > r - i + 1)
                pos = get(pos, trie[pos].h - (r - i + 1));
        }
        if (pos)
            trie[pos].add++;
        while (0 <= i - k && i + k < 2 * n && k < (n + 1) / 2 && s[i - k] == s[i + k])
        {
            int v = s[i + k] - '0';
            if (!trie[pos].son[v])
                trie[pos].son[v] = ++cnt1;
            int x = trie[pos].son[v];
            trie[x].num++;
            trie[x].h = trie[pos].h + 1;

            // update

            nxt[x][0] = pos;
            for (int j = 1; j < 25; j++)
            {
                if ((1 << j) > trie[x].h)
                    break;
                nxt[x][j] = nxt[nxt[x][j - 1]][j - 1];
            }

            pos = trie[pos].son[v];
            k++;
        }
        d[i] = k--;
        lef[i] = pos;
        if (i + k > r)
        {
            l = i - k;
            r = i + k;
        }
    }
    for (int i = 0; i < n / 2; i++)
    {
        trie[lef[i]].add = (trie[lef[i]].add + mod - 1) % mod;
    }
    dfs1(0);

    for (int i = 0; i <= cnt1; i++)
    {
        trie[i].add = trie[i].h = trie[i].num = 0;
    }
    // 偶数
    for (int i = 0, l = 0, r = -1; i < n + n / 2; i++)
    {
        int k = 0, pos = 0;
        if (i <= r)
        {
            k = min(d[l + r - i + 1], r - i + 1);
            pos = lef[l + r - i + 1];
            if (d[l + r - i + 1] > r - i + 1)
                pos = get(pos, trie[pos].h - (r - i + 1));
        }
        if (pos)
            trie[pos].add++;
        while (0 <= i - k - 1 && i + k < 2 * n && k < n / 2 && s[i - k - 1] == s[i + k])
        {
            int v = s[i + k] - '0';
            if (!trie[pos].son[v])
                trie[pos].son[v] = ++cnt2;
            int x = trie[pos].son[v];
            trie[x].num++;
            trie[x].h = trie[pos].h + 1;

            // update

            nxt[x][0] = pos;
            for (int j = 1; j < 25; j++)
            {
                if ((1 << j) > trie[x].h)
                    break;
                nxt[x][j] = nxt[nxt[x][j - 1]][j - 1];
            }

            pos = trie[pos].son[v];
            k++;
        }
        d[i] = k--;
        lef[i] = pos;
        if (i + k > r)
        {
            l = i - k - 1;
            r = i + k;
        }
    }
    for (int i = 0; i < n / 2; i++)
    {
        trie[lef[i]].add = (trie[lef[i]].add + mod - 1) % mod;
    }
    dfs2(0);
    cout << ans << '\n';
}
signed main()
{
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int t = 1;
    // cin >> t;
    while (t--)
        solve();
    return 0;
}