#include <iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iomanip>
#include<algorithm>
#include<deque>
#include<unordered_map>
#include<map>
#define P 31525197391593473ll
#define ll long long
#define l(p) t[p].l
#define sum(p) t[p].sum
#define cnt(p) t[p].cnt
#define r(p) t[p].r
#define mod 998244353
#define sc(A) scanf("%d",&A)
#define rep(A,B,C) for(int C=A;C<=B;++C)
#define fep(A,B,C) for(int C=A;C>=B;--C)
#define put(x) printf("%d\n",x)
#define db double
using namespace std;
const int MAXN=5010;
int n,m,q,x,maxx=62,top;

char a[MAXN],b[MAXN];
int nex[MAXN];
int w[MAXN];
int c1[MAXN][MAXN];
int c2[MAXN][MAXN];
struct wy
{
	struct jl
	{
		int ch[26];
		int len,fa;
		jl(){memset(ch,0,sizeof(ch));len=fa=0;}
	}t[MAXN<<2];
	int len;
	int cnt=1,last=1;
	int lin[MAXN<<2],ver[MAXN<<2],nex[MAXN<<2];
	ll sz[MAXN<<2],ans;
	inline void add(int x,int y)
	{
		ver[++len]=y;
		nex[len]=lin[x];
		lin[x]=len;
	}
	inline void insert(int x)
	{
		int p=last;
		int np=last=++cnt;sz[cnt]=1;
		t[np].len=t[p].len+1;
		while(p&&!t[p].ch[x])
		{
			t[p].ch[x]=np;
			p=t[p].fa;
		}
		if(!p)t[np].fa=1;
		else//此时存在一个子串出现过了扩充right集
		{
			int q=t[p].ch[x];
			if(t[q].len==t[p].len+1)t[np].fa=q;
			else//分裂出节点
			{
				int nq=++cnt;//代替q 成为一个right集最大的
				t[nq]=t[q];t[nq].len=t[p].len+1;
				t[q].fa=t[np].fa=nq;
				while(p&&t[p].ch[x]==q)
				{
					t[p].ch[x]=nq;
					p=t[p].fa;
				}
			}
		}
	}
	inline void dfs(int x)
	{
		for(int i=lin[x];i;i=nex[i])
		{
			int tn=ver[i];
			dfs(tn);
			sz[x]+=sz[tn];
		}
		if(sz[x]>1)ans=max(ans,sz[x]*t[x].len);
	}
}H,H1;


signed main() 
{
	//freopen("1.in","r",stdin);
	scanf("%s",a+1);
	scanf("%s",b+1);
	int n=strlen(a+1);
	int m=strlen(b+1);

	for(int i=1;i<=n;++i)H.insert(a[i]-'a');
	for(int i=2;i<=H.cnt;++i)H.add(H.t[i].fa,i);
	H.dfs(1);

	for(int i=1;i<=m;++i)H1.insert(b[i]-'a');
	for(int i=2;i<=H1.cnt;++i)H1.add(H1.t[i].fa,i);
	H1.dfs(1);


	ll ans=0;
	rep(1,m,i)
	{
		int now=1;
		rep(i,m,j)
		{
			if(!H.t[now].ch[b[j]-'a'])break;
			now=H.t[now].ch[b[j]-'a'];
			ans+=H.sz[now];
			c2[i][j]=H.sz[now];
		}
	}

	rep(1,n,i)
	{
		int now=1;
		rep(i,n,j)
		{
			if(!H1.t[now].ch[a[j]-'a'])break;
			now=H1.t[now].ch[a[j]-'a'];
			c1[i][j]=H1.sz[now];
		}
	}

	rep(2,n,j)
	{
		int now=j;
		w[j+1]=1;
		w[j]=0;
		nex[j+1]=j;
		rep(j+2,n,i)
		{
			while(now!=j&&a[i]!=a[now+1])now=nex[now];
			if(a[i]==a[now+1])++now;
			nex[i]=now;
			w[i]=w[now]+1;
		}
		now=j;
		rep(1,j-1,i)
		{
			while(now!=j&&a[i]!=a[now+1])now=nex[now];
			if(a[i]==a[now+1])++now;
			ans+=(ll)w[now]*c1[i+1][j];
		}
	}

	rep(2,m,j)
	{
		int now=j;
		w[j+1]=1;
		w[j]=0;
		nex[j+1]=j;
		rep(j+2,m,i)
		{
			while(now!=j&&b[i]!=b[now+1])now=nex[now];
			if(b[i]==b[now+1])++now;
			nex[i]=now;
			w[i]=w[now]+1;
		}
		now=j;
		rep(1,j-1,i)
		{
			while(now!=j&&b[i]!=b[now+1])now=nex[now];
			if(b[i]==b[now+1])++now;
			ans+=(ll)w[now]*c2[i+1][j];
		}
	}

	printf("%lld\n",ans);
	return 0;
}