a=int(input()) b=list(map(int,input())) mod=10**9+7 ans=0 n=len(b) for i in range(1,n): cnt=(pow(a,i,mod)-pow(a,i-1,mod))%mod ans+=(i-1)*cnt ans%=mod x=0 for y in b: x=(x*a+y)%mod ans+=(n-1)*(x-pow(a,n-1,mod)+1) ans%=mod print(ans)