#include<bits/stdc++.h>
using namespace std;
#define mp make_pair
#define db double
#define endl '\n'
#define lowbit(x) x&-x
#define intz(x,a) memset(x,a,sizeof(x));
const int N=5e3+10;
inline int min(int x,int y){return x<y?x:y;}
int lcp[N][N],lens,lent,dp[N][N],cnt[N][N],ct[N][N],sum[N][N];char s[N],t[N];
signed main(){
//	freopen("tst.out","r",stdin);
	scanf("%s%s",s+1,t+1);lens=strlen(s+1),lent=strlen(t+1);
	for(int j=lent;j>0;j--)
		for(int i=lens;i>0;i--)dp[i][j]=(s[i]^t[j]?0:dp[i+1][j+1]+1),cnt[i][dp[i][j]]++;
	for(int j=1;j<=lens;j++){
		for(int i=lens-j+1;i>0;i--)cnt[j][i]+=cnt[j][i+1];
		for(int i=j;i<=lens;i++)ct[j][i]=cnt[j][i-j+1],cnt[j][i-j+1]=0;;
	}
	for(int j=lens;j>0;j--)
		for(int i=j;i>0;i--)lcp[i][j]=(s[i]^s[j]?0:lcp[i+1][j+1]+1);
	for(int i=lens;i>0;i--)
		for(int j=i;j>0;j--)sum[i][j]=sum[i][j+1]+ct[j][i];
	long long ans=0;int len;
	for(int i=1;i<=lens;i++)
		for(int j=i+1;j<=lens;j++){len=lcp[i][j];lcp[i][j]=0;ans+=sum[j-1][i+1]-sum[j-1][min(i+len,j-1)+1];}
	for(int i=1;i<=min(lens,lent);i++)swap(s[i],t[i]);
	for(int i=min(lens,lent)+1;i<=max(lens,lent);i++)if(lens>lent)t[i]=s[i],s[i]='.';else s[i]=t[i],t[i]='.';swap(lens,lent);intz(dp,0);
	for(int j=lent;j>0;j--)
		for(int i=lens;i>0;i--)dp[i][j]=(s[i]^t[j]?0:dp[i+1][j+1]+1),cnt[i][dp[i][j]]++;
	for(int j=1;j<=lens;j++){
		for(int i=lens-j+1;i>0;i--)cnt[j][i]+=cnt[j][i+1];
		for(int i=j;i<=lens;i++)ct[j][i]=cnt[j][i-j+1];
	}
	for(int j=lens;j>0;j--)
		for(int i=j;i>0;i--)lcp[i][j]=(s[i]^s[j]?0:lcp[i+1][j+1]+1);
	for(int i=lens;i>0;i--)
		for(int j=i;j>0;j--)sum[i][j]=sum[i][j+1]+ct[j][i];
	for(int i=1;i<=lens;i++)
		for(int j=i+1;j<=lens;j++){len=lcp[i][j];ans+=sum[j-1][i+1]-sum[j-1][min(i+len,j-1)+1];}
	for(int i=1;i<=lens;i++)
		for(int j=i+1;j<=lens+1;j++)ans+=ct[i][j-1];
	cout<<ans;
	return 0;
}