#include <bits/stdc++.h>

using namespace std ;
const int N = 5100 ;
char S[N] , T[N] ;
int f[N][N] , g[N][N] , s[3][N] , n , m ;

int main()
{
 scanf("%s%s", S + 1 , T + 1) ;
 long long ans = 0 ;
 n = strlen(S + 1) ; m = strlen(T + 1) ;
 for(int i = 1 ; i <= n ; i ++)
 {
    for(int j = 1 ; j <= m ; j ++)
    {
        if(S[i] == T[j]) g[i][j] = g[i - 1][j - 1] + 1 ;
        else g[i][j] = 0 ;
    }
 }
 for(int l2 = n ; l2 >= 1 ; l2 --)
 {
    for(int i = 0 ; i <= n ; i ++) s[0][i] = s[1][i] = s[2][i] = 0 ;
    for(int l1 = l2 - 1 ; l1 >= 1 ; l1 --)
    {
        if(S[l1] == S[l2]) f[l1][l2] = f[l1 + 1][l2 + 1] + 1 ;
        else f[l1][l2] = 0 ;

        if(f[l1][l2])
        {
            int len = min(f[l1][l2] , l2 - l1 - 1) ;
            s[0][l1] ++ ;
            s[0][l1 + len] -- ; 
        }
    }
    for(int i = 1 ; i <= n ; i ++) s[1][i] = s[1][i - 1] + s[0][i] ;
    for(int i = 1 ; i <= n ; i ++) s[2][i] = s[2][i - 1] + s[1][i] ;
    for(int i = 1 ; i <= m ; i ++)
    {
        if(g[l2 - 1][i] != 0)
        {
            int l = l2 - 1 - g[l2 - 1][i] ;
            int r = l2 - 2 ;
            ans += s[2][r] - s[2][l - 1] ;
        }
        ans += g[l2][i] ;
    }
 }
 swap(n , m) ; swap(S , T) ;
 memset(g , 0 , sizeof g) ; memset(f , 0 , sizeof f) ;
 for(int i = 1;  i <= n ; i ++)
 {
    for(int j = 1 ; j <= m ; j ++)
    {
        if(S[i] == T[j]) g[i][j] = g[i - 1][j - 1] + 1 ;
        else g[i][j] = 0 ;
    }
 }
 for(int l2 = n ; l2 >= 1 ; l2 --)
 {
    for(int i = 0 ; i <= n ; i ++) s[0][i] = s[1][i] = s[2][i] = 0 ;
    for(int l1 = l2 - 1 ; l1 >= 1 ; l1 --)
    {
        if(S[l1] == S[l2]) f[l1][l2] = f[l1 + 1][l2 + 1] + 1 ;
        else f[l1][l2] = 0 ;
        if(f[l1][l2])
        {
            int len = min(f[l1][l2] , l2 - l1 - 1) ;
            s[0][l1] ++ ;
            s[0][l1 + len] -- ;
        }
    }
    for(int i = 1 ; i <= n ; i ++) s[1][i] = s[1][i - 1] + s[0][i] ;
    for(int i = 1 ; i <= n ; i ++) s[2][i] = s[2][i - 1] + s[1][i] ;

    for(int i = 1 ; i <= m ; i ++)
    {
        if(g[l2 - 1][i] != 0)
        {
            int l = l2 - 1 - g[l2 - 1][i] ;
            int r = l2 - 2 ;
            ans += s[2][r] - s[2][l - 1] ;
        }
    }
 }
 cout << ans << "\n" ;
 return 0 ;
}