#include<bits/stdc++.h>
typedef long long  LL;
using namespace std;
const int N = 5005;
int lcp[N][N];
int n,m;
vector<int> e[N*2];
int tot = 1, np = 1;
int fa[N * 2] = { 0 }, ch[N * 2][27] = { 0 }, len[N * 2] = { 0 }, cnt[N * 2] = { 0 };
int F[N * 2] = { 0 }, s[N * 2] = { 0 }, sp = 0;
class SAM {
public:
    void extend(int c) {
        int p = np; np = ++tot;
        len[np] = len[p] + 1; 
        cnt[np] = 1;
        for (; p && !ch[p][c]; p = fa[p])ch[p][c] = np;
        if (p == 0)fa[np] = 1;
        else {
            int q = ch[p][c];
            if (len[q] == len[p] + 1)fa[np] = q;
            else {
                int nq = ++tot;
                len[nq] = len[p] + 1;
                fa[nq] = fa[q]; fa[q] = nq; fa[np] = nq;
                for (; p && ch[p][c] == q; p = fa[p])ch[p][c] = nq;
                memcpy(ch[nq], ch[q], sizeof(ch[q]));
            }
        }
    }
    void dfs(int u) {
        for (auto& v : e[u]) {
            dfs(v);
            cnt[u] += cnt[v];
        }
    }
    void create(string& a) {
        for (int i = 0; i < (int)a.size(); i++)extend(a[i] - 'a');
        for (int i = 2; i <= tot; i++)e[fa[i]].push_back(i);
        for (int i = tot; i >= 1; i--)F[i] = i;
        dfs(1);
    }
    int find(int u) {
        if (u == F[u])return u;
        return F[u] = find(F[u]);
    }
    void tarjan(int u) {
        for (auto v : e[u]) {
            tarjan(v);
            F[v] = u;
        }
        if (e[u].size() == 0) {
            for (int i = 1; i <= sp; i++)lcp[n - len[u] + 1][n - len[s[i]] + 1] = len[find(s[i])];
            s[++sp] = u;
        }
    }
    void init() {
        for (int i = 1; i <= tot; i++) {
            fa[i] = len[i] = cnt[i] = 0;
            e[i].clear();
            memset(ch[i], 0, sizeof(ch[i]));
        }
        tot = 1, np = 1;
        sp = 0;
    }
}SAM;
LL dp[N][N],cmo=0;
LL ans = 0;
void solve(string& a,string& b) {
    SAM.init();
    n = a.size(); m = b.size();
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n; j++) {
            lcp[i][j] = 0;
            dp[i][j] = 0;
        }
    }
    SAM.create(b);
    for (int i = 0; i < n; i++) {
        int p = 1;
        for (int j = i; j < n; j++) {
            if (ch[p][a[j] - 'a'])p = ch[p][a[j] - 'a'];
            else break;
            dp[i][j] = cnt[p];
        }
    }
    SAM.init();
    string t = a+'{';
    reverse(t.begin(),t.end());
    SAM.create(t);
    SAM.tarjan(1);
    for (int i = 0; i < n; i++)for (int j = i + 2; j < n; j++)lcp[i][j] = min(j - i - 1, max(lcp[i][j], lcp[j][i]));
    if(!cmo)
    for (int i = 0; i < n; i++)for (int j = i; j < n; j++)ans += dp[i][j];
    cmo = 1;
    for (int j = 0; j < n; j++)for (int i = 1; i <= j; i++)dp[i][j] += dp[i - 1][j];
    for (int i = 0; i < n; i++) {
        for (int j = i + 2; j < n; j++) {//add i+1...j-1,i+2...j-1,..,i+lcp[i][j]-1...j-1
            ans += dp[i + lcp[i][j]][j - 1] - dp[i][j - 1];
        }
    }
}
signed main() {
    ios::sync_with_stdio(0);
    cin.tie(0);
    string a, b;
    cin >> a>>b;
    solve(a, b);
    solve(b, a);
    cout << ans << '\n';
    return 0;
}