#pragma GCC optimize("Ofast")
#include<bits/stdc++.h>
using namespace std;
#define pb push_back
const int maxn=1e5+5,M=5e3+5;
char a[maxn],b[maxn];
long long ans=0;
int sum[M][M],lcs[M][M];
struct SAM{
    int len[maxn],tag[maxn],ch[maxn][26],fa[maxn],las=1,tot=1,sz[maxn];
    vector<int>e[maxn];
    void clear(){
        //cout<<tot<<endl;
        for(int i=1;i<=tot;i++)memset(ch[i],0,sizeof ch[i]),fa[i]=0,tag[i]=0,len[i]=0,sz[i]=0,e[i].clear();
        tot=las=1;
    }
    void ins(int x){
        //cout<<x<<endl;
        int p=las,t=++tot;las=tot;len[t]=len[p]+1;sz[t]=1;
        for(;p&&!ch[p][x];p=fa[p])ch[p][x]=t;
        if(!p){fa[t]=1;return;}int q=ch[p][x];
        if(len[p]+1==len[q]){fa[t]=q;return ;}
        len[++tot]=len[p]+1;fa[tot]=fa[q];fa[q]=fa[t]=tot;
        memcpy(ch[tot],ch[q],sizeof ch[q]);
        for(;ch[p][x]==q;p=fa[p])ch[p][x]=tot;
    }
    void build(){
        for(int i=2;i<=tot;i++)e[fa[i]].pb(i);
    }
    void dfs(int x){
        for(auto i:e[x]){
            dfs(i);sz[x]+=sz[i];
        }
    }
}S;
void solve(char *a,char *b ,int op) {
    int n=strlen(a+1),m=strlen(b+1);
    //int n=5000,m=5000;
    S.clear();//cout<<"A";cout<<"A";
    //return;
    for(int i=1;i<=n;i++){
        S.ins(a[i]-'a');
        //if(i%100==0)cout<<"B";
    }//return;
    S.build();S.dfs(1);//return ;
    for(int pos=1;pos<=m;pos++){
        lcs[pos][pos]=m-pos+1;int id=pos-1,r=0;
        for(int i=pos+1;i<=m;i++){
            if(i<=r)lcs[pos][i]=min(lcs[pos][i-id+pos],r-i+1);
            while(i+lcs[pos][i]<=m&&b[i+lcs[pos][i]]==b[pos+lcs[pos][i]])lcs[pos][i]++;
            if(i+lcs[pos][i]-1>r)id=i,r=i+lcs[pos][i]-1;
            //cout<<r<<endl;
        }
    }//return;
    for(int j=1;j<=m;j++){
        for(int i=1;i<j;i++){
            int x=lcs[i][j];
            sum[i+1][j-1]++,sum[x+i+1][j-1]--;
        }
        for(int i=1;i<=j;i++){
            sum[i][j-1]+=sum[i-1][j-1];
        }
    }
    
   // cout<<S.tot<<endl;
    for(int i=1;i<=m;i++){
        int p=1;
        for(int j=i;j<=m;j++){
            int x=b[j]-'a';
            if(!S.ch[p][x])break;
            p=S.ch[p][x];
           // cout<<sum[i][j]<<" "<<S.sz[p]<<endl;
            ans=ans+(op+sum[i][j])*S.sz[p];
            
        }
    }
    for(int i=1;i<=m;i++)
        for(int j=1;j<=m;j++)sum[i][j]=lcs[i][j]=0;
  //  cout<<ans<<endl;
}

int main() {
    scanf("%s",a+1),scanf("%s",b+1);
    //for(int i=1;i<=5000;i++)if(i%2==1)a[i]='a',b[i]='b';else a[i]='b',b[i]='a';
    solve(a,b,1);solve(b,a,0);
    cout<<ans<<endl;
    return 0;
}