#include <bits/stdc++.h>

using namespace std;

#include <bits/stdc++.h>

using namespace std;

const int N = 6e6 + 10;
const int mod = 998244353;
void AddMod(int &p, int k) { p = (p + k) % mod; }
int n; char s[N]; 

struct Node {
	int len;
	int fail, ch[10];
}nd[N];
int last, tot = 1;
vector<int> g[N];
int f[N];
int ans;

void add(int x, int y) {
	g[x].push_back(y);
}

int getfail(int x, int i) {
	while (i - nd[x].len - 1 < 1 || s[i] != s[i - nd[x].len - 1]) x = nd[x].fail; 
	return x;
}

void insert(int i) {
	int p = getfail(last, i);
	if (!nd[p].ch[s[i] - '0']) {
		nd[++tot].len = nd[p].len + 2;
		nd[tot].fail = nd[getfail(nd[p].fail, i)].ch[s[i] - '0'];
		nd[p].ch[s[i] - '0'] = tot;
		// 先连fail边, 再连扩展状态的边. 
	}
	last = nd[p].ch[s[i] - '0'];
	if (i > n) {
		++f[last];
	}
} 

void dfs(int x) {
	for (auto y : g[x]) {
		dfs(y);
		AddMod(f[x], f[y]);
	}
	if (nd[x].len > 0 && nd[x].len <= n) {
		int v = 1ll * f[x] * f[x] % mod;
		AddMod(ans, 1ll * nd[x].len * v % mod);
	}
}

int main() {
	scanf("%d%s", &n, s + 1);
//	n = strlen(s + 1);
	nd[0].fail = 1; nd[0].len = 0; add(1, 0); // 偶根 
	nd[1].len = -1; // 奇根
	last = 0; tot = 1; 
	for (int i = 1; i <= n; ++i) s[i + n] = s[i];
	for (int i = 1; i <= (n << 1); ++i) insert(i);
	for (int i = 2; i <= tot; ++i) add(nd[i].fail, i);
	dfs(1); printf("%d\n", ans);
}