#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=5005,P1=998244353,P2=1000000007;
int n,m;
char S[N],T[N];
struct Hash{
    ll v1,v2;
    Hash(int v=0):v1(v),v2(v){}
    Hash(int v1,int v2):v1(v1),v2(v2){}
    Hash& operator+=(const Hash &o){(v1+=o.v1)%=P1;return *this;}
    Hash& operator-=(const Hash &o){(v1+=P1-o.v1)%=P1;return *this;}
    Hash& operator*=(const Hash &o){(v1*=o.v1)%=P1;return *this;}
    friend Hash operator+(Hash a,Hash b){return a+=b;}
    friend Hash operator-(Hash a,Hash b){return a-=b;}
    friend Hash operator*(Hash a,Hash b){return a*=b;}
    friend bool operator==(Hash a,Hash b){return a.v1==b.v1 and a.v2==b.v2;}
    friend bool operator<(Hash a,Hash b){return a.v1<b.v1 or (a.v1==b.v1 and a.v2<b.v2);}
}base(131);
// map<Hash,int>Ma[N],Mb[N];
// Hash ps[N],pt[N],pw[N];
int sa[N][N],sb[N][N];
struct Sam{
    struct node{
        int len,fa;
        int ch[30];
    }data[N<<2];
    int last,tot;
    int fa[N<<2][20],dep[N<<2];
    int va[N<<2],vb[N<<2];
    void init(){
        tot=1;
        data[1].fa=data[1].len=0;
        last=1;
    }
    int add(int c){
        int now=++tot;
        int f=last;last=now;
        data[now].len=data[f].len+1;
        for(;f and !data[f].ch[c];f=data[f].fa){
            data[f].ch[c]=now;
        }
        if(!f){
            data[now].fa=1;
            return now;
        }   
        int p=data[f].ch[c];
        if(data[p].len==data[f].len+1){
            data[now].fa=p;
            return now;
        }
        int q=++tot;
        data[q]=data[p];
        data[q].len=data[f].len+1;
        data[p].fa=data[now].fa=q;
        for(;f and data[f].ch[c]==p;f=data[f].fa){
            data[f].ch[c]=q;
        }
        return now;
    }
    vector<int>G[N<<2];
    bool cmp(int a,int b){
        return dep[a]<dep[b];
    }
    void ddfs(int x){
        // printf("%d %d\n",x,dep[x]);
        fa[x][0]=data[x].fa;
        dep[x]=dep[fa[x][0]]+1;
        for(int p=1;p<=15;p++){
            fa[x][p]=fa[fa[x][p-1]][p-1];
            // printf("%d ",fa[x][p-1]);
        }
        // puts("");
        for(auto y:G[x]) ddfs(y);
        va[fa[x][0]]+=va[x];
        vb[fa[x][0]]+=vb[x];
    }
    void dfs(){
        for(int i=2;i<=tot;i++){
            G[data[i].fa].push_back(i);
        }
        ddfs(1);
    }
    int lca(int x,int y){
        if(dep[x]<dep[y]){
            swap(x,y);
        }
        // printf("%d %d\n",x,y);
        int d=dep[x]-dep[y];
        for(int p=15;p>=0;p--){
            if((d>>p)&1) x=fa[x][p];
        }
        // printf("%d %d\n",x,y);
        if(x==y) return x;
        for(int p=15;p>=0;p--){
            if(fa[x][p]!=fa[y][p]){
                x=fa[x][p];
                y=fa[y][p];
            }
        }
        return fa[x][0];
    }
    int query(int x,int y){
        return data[lca(x,y)].len;
    }
}SAM;
int ra[N],rb[N];
int main(){
    // freopen("ex.in","r",stdin);
    // freopen("ex.our","w",stdout);
    scanf("%s%s",S+1,T+1);
    n=strlen(S+1);
    m=strlen(T+1);
    SAM.init();
    for(int i=1;i<=n;i++){
        ra[i]=SAM.add(S[i]-'a');
        SAM.va[ra[i]]=1;
        // printf("%d %d\n",ra[i],SAM.data[ra[i]].fa);
    }
    SAM.add(26);
    for(int i=1;i<=m;i++){
        rb[i]=SAM.add(T[i]-'a');
        SAM.vb[rb[i]]=1;
    }

    SAM.dfs();

    // pw[0]=Hash(1);
    // for(int i=1;i<N;i++){
    //     pw[i]=pw[i-1]*base;
    // }
    // ps[0]=Hash(0);
    // for(int i=1;i<=n;i++){
    //     ps[i]=ps[i-1]*base+Hash(S[i]);
    // }
    // pt[0]=Hash(0);
    // for(int i=1;i<=m;i++){
    //     pt[i]=pt[i-1]*base+Hash(T[i]);
    // }
    for(int i=1;i<=n;i++){
        for(int j=i+2;j<=n;j++){
            if(S[i]!=S[j]) continue;
            int L=min(j-i,SAM.query(ra[i],ra[j]));
            // printf("%d %d %d\n",ra[i],ra[j],L);
            sa[i+1][j-L]++;
            sa[i+1][j]--;
        }
    }
    for(int i=1;i<=m;i++){
        for(int j=i+2;j<=m;j++){
            if(T[i]!=T[j]) continue;
            int L=min(j-i,SAM.query(rb[i],rb[j]));
            sb[i+1][j-L]++;
            sb[i+1][j]--;
        }
    }

    // for(int len = 1;len<=n;len++){
    //     for(int i=1;i+len-1<=n;i++){
    //         int j=i+len-1;
    //         Ma[len][ps[j]-ps[i-1]*pw[len]]++;
    //     }
    // }
    // for(int len = 1;len<=m;len++){
    //     for(int i=1;i+len-1<=m;i++){
    //         int j=i+len-1;
    //         Mb[len][pt[j]-pt[i-1]*pw[len]]++;
    //     }
    // }
    ll ans=0;
    for(int i=1;i<=n;i++){
        int p = 1;
        for(int j=i;j<=n;j++){
            sa[i][j]+=sa[i][j-1];
            p = SAM.data[p].ch[S[j]-'a'];
            ans += 1ll * (sa[i][j]+1) * SAM.vb[p];
        }
    }
    for(int i=1;i<=m;i++){
        int p=1;
        for(int j=i;j<=m;j++){
            sb[i][j]+=sb[i][j-1];
            p = SAM.data[p].ch[T[j]-'a'];
            ans += 1ll * sb[i][j] * SAM.va[p];
        }
    }
    printf("%lld",ans);
    return 0;
}