from sys import stdin
input=lambda :stdin.readline()[:-1]

mod=998244353

n=int(input())
N=1<<n
S=input()
m=n//2
node=[[0]*N for i in range(m+1)]
node=[0]*N*(m+1)

pow11=[1]*(N+1)
pow2=[1]*(N+1)
for i in range(1,N+1):
  pow11[i]=pow11[i-1]*11%mod
  pow2[i]=pow2[i-1]*2%mod

for i in range(N):
  node[i]=int(S[i])

def calc(fL,Lsize,fR,Rsize):
  num=fL*pow11[Rsize]+fR*pow2[Lsize]
  return num%mod

def merge(i,j):
  L=j
  R=L^(1<<(i-1))
  node[(i<<n)|L]=calc(node[((i-1)<<n)|L],1<<(i-1),node[((i-1)<<n)|R],1<<(i-1))

for i in range(1,m+1):
  for j in range(N):
    merge(i,j)

q=int(input())
for _ in range(q):
  query=list(map(int,input().split()))
  
  if query[0]==1:
    x,y=query[1:]
    node[x]=y
    
    for i in range(1,m+1):
      l=(x>>i)<<i
      r=l+(1<<i)
      for j in range(l,r):
        merge(i,j)
  
  else:
    L,R,x=query[1:]
    R+=1
    res=0
    size=0
    
    s=(1<<m)
    while L<R and L%s!=0:
      res=calc(res,size,node[L^x],1)
      L+=1
      size+=1
          
    while L+s<=R:
      res=calc(res,size,node[N*m|L^x],s)
      L+=s
      size+=s

    while L<R:
      res=calc(res,size,node[L^x],1)
      L+=1
      size+=1
    
    print(res)