#include <stdio.h>
#define maxn 202
using namespace std;

long d[maxn][maxn][maxn];
long n, i, j, k, l1, l2, sol, l, p1, p2, q1, q2;
char a[maxn], b[maxn];
const long co=3210121;

int main()
{
    freopen("iv.in", "r", stdin);
    freopen("iv.out", "w", stdout);
    scanf("%c", &a[1]);
    l1=1;
    while(a[l1]>='A')
    {
        l1++;
        scanf("%c", &a[l1]);
    }
    l1--;
    scanf("%c", &b[1]);
    l2=1;
    while(b[l2]>='A')
    {
        l2++;
        scanf("%c", &b[l2]);
    }
    l2--;
    d[0][0][l1+1]=1;
    for(i=0; i<=(l1+l2)/2; i++)
    {
        for(p1=0; p1<=l1; p1++)
        {
            for(q1=p1+1; q1<=l1+1; q1++)
            {
                p2=i-p1;
                q2=l1-q1+1+l2-i+1;
                if(p2<q2)
                {
                  //  printf("*");
                    if(a[p1+1]==a[q1-1])
                    {
                        d[i+1][p1+1][q1-1]+=d[i][p1][q1];
                        if(d[i+1][p1+1][q1-1]>co) d[i+1][p1+1][q1-1]-=co;
                    }
                    if(a[p1+1]==b[q2-1])
                    {
                        d[i+1][p1+1][q1]+=d[i][p1][q1];
                        if(d[i+1][p1+1][q1]>co) d[i+1][p1+1][q1]-=co;
                    }
                    if(b[p2+1]==a[q1-1])
                    {
                        d[i+1][p1][q1-1]+=d[i][p1][q1];
                        if(d[i+1][p1][q1-1]>co) d[i+1][p1][q1-1]-=co;
                    }
                    if(b[p2+1]==b[q2-1])
                    {
                        d[i+1][p1][q1]+=d[i][p1][q1];
                        if(d[i+1][p1][q1]>co) d[i+1][p1][q1]-=co;
                    }
                }
           //     printf("%d ", d[i][p1][q1]);
            }
      //      printf("\n");
        }
     //   printf("\n");
    }
    i=(l1+l2)/2+(l1+l2)%2;
    for(p1=0; p1<=l1; p1++)
    {
        for(q1=p1; q1<=p1+1; q1++)
        {
            p2=i-p1;
            q2=l1-q1+1+l2-i+1;
            if(q2-p2<=1)
            {
                //if(q2-p2 == 1 || q1-p1 == 1)
                {
                    sol+=d[i][p1][q1];
                }
            }
        }
    }
    printf("%d\n", sol);
    return 0;
}