#include<fstream>
#include<cstring>
#define MOD 3210121
#define N 510
using namespace std;
ifstream f("iv.in");
ofstream g("iv.out");
int n,m,i,j,k,l,x,sol,d[2][N][N];
char a[N],b[N];
int main()
{
    f>>(a+1)>>(b+1);
    n=strlen(a+1);
    m=strlen(b+1);
    d[0][0][0]=1;
    for(i=0,l=1;i<=n;++i,l^=1)
    for(j=0;i+j<=n;++j)
    for(k=0;k<=m;++k)
    {///////////////
        x=i+k-j;///
        if(x<0||k+x>m)
        continue;
        if(x+i+j+k==n+m||x+i+j+k==n+m-1)
        {
            sol+=d[l^1][j][k];
            if(sol>MOD)
            sol-=MOD;
            continue;
        }
        if(a[i+1]==a[n-j])
        {
            d[l][j+1][k]+=d[l^1][j][k];
            if(d[l][j+1][k]>MOD)
            d[l][j+1][k]-=MOD;
        }
        if(m-x>=1&&a[i+1]==b[m-x])
        {
            d[l][j][k]+=d[l^1][j][k];
            if(d[l][j][k]>MOD)
            d[l][j][k]-=MOD;
        }
        if(b[k+1]==a[n-j])
        {
            d[l^1][j+1][k+1]+=d[l^1][j][k];
            if(d[l^1][j+1][k+1]>MOD)
            d[l^1][j+1][k+1]-=MOD;
        }
        if(m-x>=1&&b[k+1]==b[m-x])
        {
            d[l^1][j][k+1]+=d[l^1][j][k];
            if(d[l^1][j][k+1]>MOD)
            d[l^1][j][k+1]-=MOD;
        }
        d[l^1][j][k]=0;
    }
    g<<sol<<'\n';
    return 0;
}