#include<bits/stdc++.h>
#define Gc() getchar()
#define Me(x,y) memset(x,y,sizeof(x))
#define Mc(x,y) memcpy(x,y,sizeof(x))
#define d(x,y) ((m)*(x-1)+(y))
#define R(n) (rnd()%(n)+1)
#define Pc(x) putchar(x)
#define LB lower_bound
#define UB upper_bound
#define fi first
#define se second
#define eb emplace_back
using namespace std;using ll=long long;using db=double;using lb=long db;using ui=unsigned;using ull=unsigned long long;using pii=pair<int,int>;
const int N=6e6+5,M=N*10+5,K=1000+5,mod=998244353,Mod=+-1;const db eps=1e-9;const int INF=1e9+7;mt19937 rnd(263082);
#define Tp template<typename T>
#define Ts template<typename T,typename... Ar>
namespace Debug{
	Tp void _debug(char* f,T t){cerr<<f<<'='<<t<<endl;}
	Ts void _debug(char* f,T x,Ar... y){while(*f!=',') cerr<<*f++;cerr<<'='<<x<<",";_debug(f+1,y...);}
	#ifdef LOCAL
	#define gdb(...) _debug((char*)#__VA_ARGS__,__VA_ARGS__)
	#else 
	#define gdb(...) void()
	#endif
}using namespace Debug;
int n,siz[N];char s[N];vector<int> S[N];
namespace PAM{
	int son[N][10],link[N],cnt=1,len[N];
	int find(int x,int id){
		while(x^1&&s[id]^s[id-len[x]-1]) x=link[x];
		return x;
	}
	ll ans;
	void dfs(int x){
		int w=0;
		for(int i:S[x]) dfs(i),siz[x]+=siz[i];
		if(len[x]<=n&&len[x]>0){
			gdb(siz[x],x,len[x]);
			ans+=1ll*siz[x]*siz[x]%mod*len[x]%mod;
		}
	}
	void build(){
		link[0]=1;len[1]=-1;int x=0;
		for(int i=1;i<=2*n;i++){
			x=find(x,i);
			if(!son[x][s[i]-'0']){
				len[++cnt]=len[x]+2;
				link[cnt]=son[find(link[x],i)][s[i]-'0'];
				son[x][s[i]-'0']=cnt;
			}
			x=son[x][s[i]-'0'];
			if(i>n) siz[x]++;
		}
		for(int i=0;i<=cnt;i++) if(i^1) S[link[i]].emplace_back(i);
		dfs(1);
		printf("%lld\n",ans%mod);
	}
}
void Solve(){
	int i,j;scanf("%d",&n);
	scanf("%s",s+1);for(i=n+1;i<=2*n;i++) s[i]=s[i-n];
	PAM::build();
}
int main(){
	int t=1;
	// scanf("%d",&t);
	while(t--) Solve();
	cerr<<clock()*1.0/CLOCKS_PER_SEC<<'\n';
}