//code by Emissary
#pragma GCC optimize(2)
#pragma GCC optimize(3)
#pragma GCC optimize("Ofast")
#include<bits/stdc++.h>

#define fi first
#define se second
#define vc vector
#define db double
#define ll long long
#define mk make_pair
#define pb push_back
#define PI pair<int,int>
#define ull unsigned long long
#define err cerr << "   -_-   " << endl
#define debug cerr << " ------------------- " << endl

#define input(x) freopen(#x".in","r",stdin)
#define output(x) freopen(#x".out","w",stdout)

#define NO puts("No")
#define YES puts("Yes")

//#define int long long

using namespace std;

namespace IO{
	inline int read(){
		int X=0, W=0; char ch=getchar();
		while(!isdigit(ch)) W|=ch=='-', ch=getchar();
		while(isdigit(ch)) X=(X<<1)+(X<<3)+(ch^48), ch=getchar();
		return W?-X:X;
	}
	inline void write(int x){
		if(x<0) x=-x, putchar('-');
		if(x>9) write(x/10);
		putchar(x%10+'0');
	}
	inline void sprint(int x){write(x), putchar(32);}
	inline void eprint(int x){write(x), putchar(10);}
}using namespace IO;

const int MAXN = 5005;
const int base = 13331;
const int mod = 998244353;

int lcp[MAXN][MAXN], add[MAXN];
int n, m, pres[MAXN], pret[MAXN];
int S[MAXN][MAXN], T[MAXN][MAXN], powc[MAXN];

ll ans;

char s[MAXN], t[MAXN];

int ss[MAXN], tt[MAXN], sss, ttt;

unordered_map<int,int> total;

inline int hs(int l, int r){return (pres[r]-1ll*pres[l-1]*powc[r-l+1]%mod+mod)%mod;}

inline int ht(int l, int r){return (pret[r]-1ll*pret[l-1]*powc[r-l+1]%mod+mod)%mod;}

inline ll solve1(){
	ll sum=0;
	for(int i=1;i<=n;++i)
		for(int j=1;j<=m;++j) lcp[i][j]=(s[i]==t[j]?1:0);
	for(int i=n;i>=2;--i)
		for(int j=m;j>=2;--j)
			if(s[i-1]==t[j-1]) lcp[i-1][j-1]=max(lcp[i-1][j-1],lcp[i][j]+1);
	for(int l=1;l<=n;++l){
		memset(add,0,sizeof add);
		for(int j=1;j<=m;++j) add[lcp[l][j]]++;
		for(int j=n;j>=1;--j) add[j]+=add[j+1];
		for(int len=1;l+len-1<=n;++len) S[l][l+len-1]=add[len], sum+=S[l][l+len-1];
	} 
	for(int l=1;l<=m;++l){
		memset(add,0,sizeof add);
		for(int j=1;j<=n;++j) add[lcp[j][l]]++;
		for(int j=m;j>=1;--j) add[j]+=add[j+1];
		for(int len=1;l+len-1<=m;++len) T[l][l+len-1]=add[len];
	} 
	return sum;
}

inline ll solve2(){
	ll sum=0;
	for(int i=1;i<=n;++i)
		for(int j=1;j<=n;++j)
			lcp[i][j]=(s[i]==s[j]?j:0);
	for(int i=n;i>=2;--i)
		for(int j=n;j>=2;--j)
			if(s[i-1]==s[j-1]) lcp[i-1][j-1]=max(lcp[i-1][j-1],lcp[i][j]);
	for(int r=n-1;r>=1;--r){
		memset(add,0,sizeof add); int now=0;
		for(int i=1;i<r;++i) add[lcp[r+1][i]]++;
		for(int i=n;i>=r;--i) now+=add[i];
		for(int l=r;l>=2;--l){
			now+=add[l-1];
			if(s[l]==s[r+1] && l!=r) now--;
			sum+=1ll*S[l][r]*now; 
			//cerr << "range:" << l << ' ' << r << ' ' << S[l][r] << ' ' << now << endl;
		}
	}
	return sum;
}

signed main(){
	scanf("%s%s",s+1,t+1);
	n=strlen(s+1), m=strlen(t+1);
	//n=m=5000;
	//for(int i=1;i<=n;++i) s[i]='a'+rand()%26;
	//for(int i=1;i<=n;++i) t[i]='a'+rand()%26;
	ans+=solve1();
	//cerr << "ans=" << ans << endl;
	ans+=solve2();
//	cerr << "ans=" << ans << endl;
	//return 0;
	swap(n,m); swap(s,t); swap(S,T);
	ans+=solve2();
	cout << ans << endl;
	return 0;
}