#include<bits/stdc++.h>
#define N 3000000
using namespace std;
const int mod=998244353;
inline int fmod(const int &x) {return x<mod?x:x-mod;}
inline int qmod(const int &x) {return x<0?x+mod:x;}
inline void add(int &x,const int &y) {x=fmod(x+y);}
inline void sub(int &x,const int &y) {x=qmod(x-y);}
inline int quick_pow(int x,int b)
{
	int ret=1;
	for(;b;b>>=1,x=1ll*x*x%mod)
		if(b&1) ret=1ll*ret*x%mod;
	return ret;
}
inline int inv(int x) {return quick_pow(x,mod-2);}
struct node
{
	int trans[10],fail,num,len;
}t[N*2+5];
int n,ans,cnt,last=1;
char s[N*2+5];
inline void insert(int p)
{
	while(s[p]!=s[p-t[last].len-1]) last=t[last].fail;
	int ch=s[p]^'0';
	if(t[last].trans[ch]) last=t[last].trans[ch];
	else
	{
		t[last].trans[ch]=++cnt;
		t[cnt].len=t[last].len+2;
		last=t[last].fail;
		while(last && s[p]!=s[p-t[last].len-1]) last=t[last].fail;
		t[cnt].fail=t[last].trans[ch];
		last=cnt;
	}
	if(p>=n) ++t[last].num;
}
int main()
{
	scanf("%d",&n);
	scanf("%s",s+1);
	cnt=2,t[1].len=-1,t[2].fail=1;
	for(int i=0;i<10;++i) t[0].trans[i]=2;
	for(int i=1;i<n;++i) s[i+n]=s[i];
	for(int i=1;i<2*n;++i) insert(i);
	for(int i=cnt;i>1;--i)
	{
		t[t[i].fail].num+=t[i].num;
		if(t[i].len<=n) add(ans,1ll*t[i].num*t[i].num%mod*t[i].len%mod);
	}
	printf("%d\n",ans);
	return 0;
}