#include<bits/stdc++.h>
#define MOD 666013
using namespace std;
char a[505],b[505];
int d[505][505],t[505][505],n,m;
int main()
{
    freopen("subsir.in","r",stdin);
    freopen("subsir.out","w",stdout);
    gets(a+1);
    gets(b+1);
    n=strlen(a+1);
    m=strlen(b+1);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            if(a[i]==b[j])
            {
                d[i][j]=d[i-1][j-1]+1;
            }
                else
            {
                d[i][j]=max(d[i][j-1],d[i-1][j]);
            }
        }
    }
    for(int i=1;i<=n;i++)
        t[i][1]=1;
    for(int i=1;i<=m;i++)
        t[1][i]=1;
    for(int i=2;i<=n;i++)
    {
        for(int j=2;j<=m;j++)
        {
            if (a[i]==b[j])
            {
                t[i][j]=t[i-1][j-1];
                continue;
            }
            if (d[i][j]==d[i-1][j])
            {
                t[i][j]=(t[i][j]+t[i-1][j])%MOD;
            }
            if (d[i][j]==d[i][j-1])
            {
                t[i][j]=(t[i][j]+t[i][j-1])%MOD;
            }
            if (d[i][j]==d[i-1][j-1])
            {
                t[i][j]=t[i][j]-t[i-1][j-1];
                if (t[i][j]<0) t[i][j]+=MOD;
            }
        }
    }
    printf("%d\n",t[n][m]);
    return 0;
}