k=int(input())
P=1000000007
T=(((1,1),(0,1)),((1,0),(1,1)))
s=input()
def m(a,b):
 return(
  ((a[0][0]*b[0][0]+a[0][1]*b[1][0])%P,(a[0][0]*b[0][1]+a[0][1]*b[1][1])%P),
  ((a[1][0]*b[0][0]+a[1][1]*b[1][0])%P,(a[1][0]*b[0][1]+a[1][1]*b[1][1])%P)
 )
X=((1,0),(0,1))
for c in s:
 X=m((T[0]if c=='0'else T[1]),X)
Y=((1,0),(0,1))
while k:
 if k&1:Y=m(Y,X)
 X=m(X,X)
 k=k//2
a=Y[0][0]+Y[0][1]+Y[1][0]+Y[1][1]
print(a%P)