#include <iostream>
#include <fstream>
#include <cstring>
#include <algorithm>
#define mod 666013

using namespace std;
ifstream fin("subsir.in");
ofstream fout("subsir.out");
int n,m,dp[1200][1200],nr[1200][1200];
char a[1200],b[1200];
void citire()
{
    fin>>a;
    fin>>b;
    n=strlen(a);
    m=strlen(b);
}
void cmlsc()
{
    for(int i=0; i<n; i++)
        if(a[i]==b[0]||dp[i-1][0]) dp[i][0]=1;
    for(int i=0; i<m; i++)
        if(b[i]==a[0]||dp[0][i-1]) dp[0][i]=1;
    for(int i=1; i<n; i++)
        for(int j=1; j<m; j++)
        {
            if(a[i]==b[j]) dp[i][j]=dp[i-1][j-1]+1;
            else
                dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
        }
}
void numar()
{
    for(int i=0; i<n; i++)
        nr[i][0]=1;
        for(int i=0; i<m; i++)
            nr[0][i]=1;
        for(int i=1; i<n; i++)
            for(int j=1; j<m; j++)
            {
                if(a[i]==b[j]) nr[i][j]=nr[i-1][j-1]%mod;
                else
                {
                    if(dp[i-1][j]==dp[i][j]) nr[i][j]+=nr[i-1][j];
                    if(dp[i][j-1]==dp[i][j]) nr[i][j]+=nr[i][j-1];
                    if(dp[i-1][j-1]==dp[i][j]) nr[i][j]-=nr[i-1][j-1];
                }
            nr[i][j]+=mod;
            nr[i][j]=nr[i][j]%mod;
            }
}
int main()
{
    citire();
    cmlsc();
    numar();
    fout<<nr[n-1][m-1]<<"\n";
    return 0;
}