#include <bits/stdc++.h>
 
//    #include <boost/multiprecision/cpp_int.hpp>
 #define int long long
 #define inf  1000000007
 #define pa pair<int,int>
 #define ll long long
 #define pal pair<double,double>
 #define ppap pair<pa,int>
  #define PI 3.14159265358979323846
  #define paa pair<int,char>
  #define  mp make_pair
  #define  pb push_back
  #define EPS (1e-8)
                                          
    int dx[8]={0,1,0,-1,1,1,-1,-1};
    int dy[8]={1,0,-1,0,-1,1,1,-1};
                                            using namespace std;
                                   			class pa3{
                                            	public:
                                            	int x;
                                   				int y,z;
                                            	pa3(int x=0,int y=0,int z=0):x(x),y(y),z(z) {}
                                            	bool operator < (const pa3 &p) const{
                                            		if(x!=p.x) return x<p.x;
                                            		if(y!=p.y) return y<p.y;
                                            		 return z<p.z;
                                            		//return x != p.x ? x<p.x: y<p.y;
                                            	}
                                   				bool operator > (const pa3 &p) const{
                                            		if(x!=p.x) return x>p.x;
                                            		if(y!=p.y) return y>p.y;
                                            		 return z>p.z;
                                            		//return x != p.x ? x<p.x: y<p.y;
                                            	}
                                            	bool operator == (const pa3 &p) const{
                                            		return x==p.x && y==p.y && z==p.z;
                                            	}
                                            		bool operator != (const pa3 &p) const{
                                            			return !( x==p.x && y==p.y && z==p.z);
                                            	}
                                            
                                            };
                                            
                                            class pa4{
                                            	public:
                                            	int x;
                                            	int y,z,w;
                                            	pa4(int x=0,int y=0,int z=0,int w=0):x(x),y(y),z(z),w(w) {}
                                            	bool operator < (const pa4 &p) const{
                                            		if(x!=p.x) return x<p.x;
                                            		if(y!=p.y) return y<p.y;
                                            		if(z!=p.z)return z<p.z;
                                            		return w<p.w;
                                            		//return x != p.x ? x<p.x: y<p.y;
                                            	}
                                            	bool operator > (const pa4 &p) const{
                                            		if(x!=p.x) return x>p.x;
                                            		if(y!=p.y) return y>p.y;
                                            		if(z!=p.z)return z>p.z;
                                            		return w>p.w;
                                            		//return x != p.x ? x<p.x: y<p.y;
                                            	}
                                            	bool operator == (const pa4 &p) const{
                                            		return x==p.x && y==p.y && z==p.z &&w==p.w;
                                            	}
                                            		
                                            
                                            };
                                            class pa2{
                                            	public:
                                            	int x,y;
                                            	pa2(int x=0,int y=0):x(x),y(y) {}
                                            	pa2 operator + (pa2 p) {return pa2(x+p.x,y+p.y);}
                                            	pa2 operator - (pa2 p) {return pa2(x-p.x,y-p.y);}
                                            	bool operator < (const pa2 &p) const{
                                            		return y != p.y ? y<p.y: x<p.x;
                                            	}
                                            	bool operator > (const pa2 &p) const{
                                            		return x != p.x ? x<p.x: y<p.y;
                                            	}
                                            	bool operator == (const pa2 &p) const{
                                            		return abs(x-p.x)==0 && abs(y-p.y)==0;
                                            	}
                                            	bool operator != (const pa2 &p) const{
                                            		return !(abs(x-p.x)==0 && abs(y-p.y)==0);
                                            	}
                                            		
                                            
                                            };
                                            
 
                      
                                string itos( int i ) {
                                ostringstream s ;
                                s << i ;
                                return s.str() ;
                                }
                                 
                                int gcd(int v,int b){
                                	if(v>b) return gcd(b,v);
                                	if(v==b) return b;
                                	if(b%v==0) return v;
                                	return gcd(v,b%v);
                                }
                 
                            
                                int mod;
int extgcd(int a, int b, int &x, int &y) {
    if (b == 0) {
        x = 1;
        y = 0;
        return a;
    }
    int d = extgcd(b, a%b, y, x);
    y -= a/b * x;
    return d;
}
pa operator+(const pa & l,const pa & r) {   
    return {l.first+r.first,l.second+r.second};                                    
}    
pa operator-(const pa & l,const pa & r) {   
    return {l.first-r.first,l.second-r.second};                                    
}  
                                
                int pr[1000010];
                int inv[1000010];
                
                int beki(int wa,int rr,int warukazu){
                	if(rr==0) return 1%warukazu;
                	if(rr==1) return wa%warukazu;
                	wa%=warukazu;
                	if(rr%2==1) return ((ll)beki(wa,rr-1,warukazu)*(ll)wa)%warukazu;
                	ll zx=beki(wa,rr/2,warukazu);
                	return (zx*zx)%warukazu;
                }
 
                
    			int comb(int nn,int rr){
    				if(rr<0 || rr>nn || nn<0) return 0;
    				int r=pr[nn]*inv[rr];
    				r%=mod;
    				r*=inv[nn-rr];
    				r%=mod;
    				return r;
    			}
               
                void gya(int ert){
                	pr[0]=1;
                	for(int i=1;i<=ert;i++){
                		pr[i]=(pr[i-1]*i)%mod;
                	}
                		inv[ert]=beki(pr[ert],mod-2,mod);
                	for(int i=ert-1;i>=0;i--){
                		inv[i]=inv[i+1]*(i+1)%mod;
                	}
                }
                
              //   cin.tie(0);
    		//	ios::sync_with_stdio(false);
    			//priority_queue<pa3,vector<pa3>,greater<pa3>> pq;            
                 //sort(ve.begin(),ve.end(),greater<int>());
 //   mt19937(clock_per_sec);
  //  mt19937_64 rng(chrono::steady_clock::now().time_since_epoch().count()) ;
      
       
       


                               //----------------kokomade tenpure------------
int n;
 	string s;
 	
 	vector<pa> ve;
 	
bool check(int e,int f){
	//cout<<e<<" "<<f<<endl;
	int t=ve.size();
	for(int i=0;i<t;i++){
		int p;
		if(i>0){
			p=f;
			
		}
		else p=e;
		
			if(ve[i].second>p){
				if(ve[i].first==1)return 1;
				if(ve[i].first==0){
					return 0;
				}
			}
			if(ve[i].second<p){
				if(ve[i].first==1){
					if(i==t-1) return 1;
					return 0;
				}
				if(ve[i].first==0){
					return 1;
				}
			}
	}
	return 1;
}
 	
int solve(int k){
	int e=n/k;
	
	int ok=0;
	
	if(ve[0].first==0){
	//	int si=min(e+1,ve[0].second+1);
		int si=e-ve[0].second;
		si=max(0ll,si);
		ok=si;
	//	cout<<"d"<<endl;
		if(ve[0].second<=e)if(check(ve[0].second,e)) ok++;
	}
	
	else{
		ok=e;
		ok+= min(e,ve[0].second-1);
		if(ve[0].second<=e)if(check(ve[0].second,e))ok++;
	}
	
	
	return ok;
	
}


int solve2(int k){
	int e=n/k;
	
	int ok=0;
	
	int ko=0;
	for(int i=0;i<e;i++){
		ko*=2;
		ko%=mod;
		if(s[i]=='1')ko++;
	}
	ko%=mod;
	ok+=ko;
	
	string t=s.substr(0,e);
	t+=t;
//	cout<<t<<endl;
	for(int i=e;i<2*e;i++){
		if(t[i]=='1')t[i]='0';
		else t[i]='1';
	}
	while(t.length()<s.length())t+=t.substr(0,2*e);
	
	while(t.length()>s.length())t.pop_back();
	//cout<<t<<endl;
	if(s>=t) ok++;
	return ok;
	
}



 signed main(){

   cin.tie(0);
	ios::sync_with_stdio(false);

 	cin>>n>>s;
 	mod= 998244353;
 	for(auto c:s){
 		if(ve.size()==0 || ve.back().first!=c-'0')ve.pb(mp(c-'0',1));
 		else ve.back().second++;
 	}
 	
 	int zen=0;
 	for(int i=0;i<n;i++){
 		zen*=2;
 		zen%=mod;
 		if(s[i]=='1')zen++;
 	
 	}
 	zen++;
 	zen%=mod;
 	//cout<<zen<<endl;
 	int ans=0;
 	vector<pa> F;
 	
 	for(int i=3;i<=n;i+=2)if(n%i==0){
 		int y=solve2(i);
 	//	cout<<n/i<<" "<<y<<endl;
 		F.pb(mp(n/i,y));
 	}
 	reverse(F.begin(),F.end());
 	int r=F.size();
 	
 	for(int i=0;i<r;i++)for(int j=i+1;j<r;j++)if(F[j].first%F[i].first==0){
 	//	cout<<i<<" "<<j<<endl;
 		F[j].second+=mod-F[i].second;
 		F[j].second%=mod;
 	}
 	
 	for(int i=0;i<r;i++){
 	//	cout<<F[i].first<<" "<<F[i].second<<endl;
 		zen+=mod-F[i].second;
 		ans+=2*F[i].first*F[i].second%mod;
 		ans%=mod;
 		zen%=mod;
 	}
 	
 	
 //	cout<<zen<<endl;
 	ans+=zen*2*n%mod;
 	
 	cout<<ans%mod<<endl;
 	
 	
 	
return 0;
 
 }
 

Measured memory_bytes: 6720