#include <fstream>
#include<unordered_map>
using namespace std;
ifstream in("subsir.in");
ofstream out("subsir.out");
#define int long long
int dp[1001][1001];
int previ[1001][1001];
string a,b;
int memo[1001][1001];
int cnt=0;
string s;
unordered_map<string,bool> ma;
void solve(int n,int m)
{
    memo[n][m]=1;
    if(n==0 || m==0)
    {
        return;
    }
    int i=n,j=m;
    if(previ[i][j]==1)
    {
        solve(i-1,j);
    }
    else if(previ[i][j]==2)
    {
        solve(i,j-1);
    }
    else
    {
        cnt++;
        solve(i-1,j-1);
        s+=a[i-1];
    }
}
signed main()
{
    in>>a>>b;
    int n=a.size();
    int m=b.size();
    for(int i=1 ; i<=n; i++)
    {
        for(int j=1; j<=m; j++)
        {
            dp[i][j]=max(dp[i-1][j],dp[i][j-1]);
            if(a[i-1]==b[j-1])
            {
                dp[i][j]=max(dp[i][j],dp[i-1][j-1]+1);
            }
            else
            {
                dp[i][j]=max(dp[i][j],dp[i-1][j-1]);
            }
            if(dp[i][j]==dp[i-1][j])
            {
                previ[i][j]=1;
            }
            else if(dp[i][j]==dp[i][j-1])
            {
                previ[i][j]=2;
            }
            else
            {
                previ[i][j]=3;
            }
        }
    }
    int ans=dp[n][m],ras=0;
    for(int i=n; i>=1; i--)
    {
        for(int j=m; j>=1; j--)
        {
            if(memo[i][j]==0)
            {
                cnt=0;
                s="";
                solve(i,j);
                if(cnt==ans)
                {
                    if(ma[s]==0)
                    {
                        ma[s]=1;
                        ras++;
                    }
                }
            }
        }
    }
    out<<ras%666013;
    return 0;
}