#include<stdio.h>
#include<string.h>
#define M 666013
char a[502],b[502];
int n,m,i,j,A[502][502],Nr[502][502];
int main(){
freopen("subsir.in","r",stdin);
freopen("subsir.out","w",stdout);
scanf("%s",a+1);
n=strlen(a+1);
scanf("%s",b+1);
m=strlen(b+1);
for(i=0;i<=n;++i){Nr[i][0]=1;}
for(j=0;j<=m;++j){Nr[0][j]=1;}
for(i=1;i<=n;++i)
  for(j=1;j<=m;++j)
    {if(a[i]==b[j])
      {A[i][j]=A[i-1][j-1]+1;
      Nr[i][j]=Nr[i-1][j-1];}
     else
      if(A[i-1][j]>A[i][j-1])
       {A[i][j]=A[i-1][j];
    Nr[i][j]=Nr[i-1][j];}
      else
       if(A[i][j-1]>A[i-1][j])
    {A[i][j]=A[i][j-1];
    Nr[i][j]=Nr[i][j-1];}
       else
    {A[i][j]=A[i-1][j];
    Nr[i][j]=(Nr[i][j-1]+Nr[i-1][j])%M;
    if(A[i-1][j-1]==A[i-1][j])
     Nr[i][j]=(Nr[i][j]-Nr[i-1][j-1]+M)%M;}}
printf("%d",Nr[n][m]);
return 0;}