#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,N=40000007;
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
{
	int lin[N],nex[MAXN*MAXN],e[MAXN*MAXN];
	ll ver[MAXN*MAXN];
	int len=0;
	inline void add(ll y)
	{
		int x=y%N;
		ver[++len]=y;
		nex[len]=lin[x];
		lin[x]=len;
		e[len]=1;
	}
	int find(ll y,int w)
	{
		int x=y%N;
		for(int i=lin[x];i;i=nex[i])
		{
			ll tn=ver[i];
			if(tn==y){e[i]+=w;return e[i];}
		}
		return 0;
	}
	void clear(){
		fill(lin,lin+N,0);len=0;
		/*
		fill(nex,nex+MAXN*MAXN,0);
		fill(e,e+MAXN*MAXN,0);
		fill(ver,ver+MAXN*MAXN,0);*/ 
	} 
}H;

signed main() {
	//freopen("1.in", "r", stdin);
	//freopen("1.out", "w", stdout);
	scanf("%s", a + 1);
	scanf("%s", b + 1);
	int n = strlen(a + 1);
	int m = strlen(b + 1);
	rep(1, n, i) {
		int base1 = 29;
		int base2 = 114514;
		unsigned ll s2 = 0;
		ll s1 = 0;
		rep(i, n, j) {
			s1 = ((ll)s1 * base1 + a[j]) % P;
			//s2=s2*base2+a[j];
			if(H.find(s1,1)==0)H.add(s1);
			//++H[make_pair(s2,s1)];
		}
	}

	ll ans = 0;
	rep(1, m, i) {
		int base1 = 29;
		int base2 = 114514;
		unsigned ll s2 = 0;
		ll s1 = 0;
		rep(i, m, j) {
			s1 = ((ll)s1 * base1 + b[j]) % P;
			c2[i][j]=H.find(s1,0);
			ans +=c2[i][j];
		}
	}
	H.clear();
	rep(1, m, i) {
		int base1 = 29;
		int base2 = 114514;
		unsigned ll s2 = 0;
		ll s1 = 0;
		rep(i, m, j) {
			s1 = ((ll)s1 * base1 + b[j]) % P;
			if(H.find(s1,1)==0)H.add(s1);
		}
	}

	rep(1, n, i) {
		int base1 = 29;
		int base2 = 114514;
		unsigned ll s2 = 0;
		ll s1 = 0;
		rep(i, n, j) {
			s1 = ((ll)s1 * base1 + a[j]) % P;
			c1[i][j] = H.find(s1,0);
		}
	}



	rep(2, n-1, j) {
		int now = j;
		w[j + 1] = 1;
		w[j] = 0;
		nex[j + 1] = j;
		nex[j]=j;
		rep(j + 2, min(2 * j, 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];

			if (now == n)
				now = nex[now];
		}
	}
	memset(nex, 0, sizeof(nex));

	rep(2, m-1, j) {
		int now = j;
		w[j + 1] = 1;
		w[j] = 0;
		nex[j + 1] = j;
		nex[j]=j;
		rep(j + 2, min(2 * j, 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];

			if (now == m)
				now = nex[now];
		}
	}

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