#include <bits/stdc++.h>
#define Mod 1000000007
#define getchar() ((p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<20,stdin),p1==p2))?EOF:*p1++)
char buf[1<<20],*p1=buf,*p2=buf;
using namespace std;
int get_str(char *S)
{
	int len=1;
	do S[1]=getchar();while(S[1]<'a'||S[1]>'z');
	do S[++len]=getchar();while(S[len]>='a'&&S[len]<='z');
	return len-1;
}
inline void add(int &a,int b){a+=b;if(a>=Mod) a-=Mod;}
inline void del(int &a,int b){a-=b;if(a<0) a+=Mod;}
inline int chk(int a){return a<0?a+Mod:a;}
char S[500010],T[500010];
int lenS,lenT;
int sum2[500010],ssum2[500010],ans;
namespace SAM{
	struct Node{
		int link,len,nxt[27];
		int cnt[2],sumy;
	}s[2000010];
	vector<int> son[2000010];
	int poi_cnt;
	void reset()
	{
		for(int i=0;i<=poi_cnt;i++) vector<int>().swap(son[i]);
		memset(s,0,sizeof(s));poi_cnt=0;
		s[0].link=-1;
	}
	int ins_nd(int las,int nw)
	{
		int cur=++poi_cnt;
		s[cur].len=s[las].len+1;
		int p=las;
		while(p!=-1&&!s[p].nxt[nw])
		{
			s[p].nxt[nw]=cur;
			p=s[p].link;
		}
		if(p==-1) s[cur].link=0;
		else
		{
			int q=s[p].nxt[nw];
			if(s[p].len+1==s[q].len) s[cur].link=q;
			else
			{
				int clo=++poi_cnt;
				s[clo].len=s[p].len+1;
				s[clo].link=s[q].link;
				memcpy(s[clo].nxt,s[q].nxt,27<<2);
				while(p!=-1&&s[p].nxt[nw]==q)
				{
					s[p].nxt[nw]=clo;
					p=s[p].link;
				}
				s[cur].link=s[q].link=clo;
			}
		}
		return cur;
	}
	void dfs(int a,bool op)
	{
		for(int v:son[a])
		{
			dfs(v,op);
			s[a].cnt[0]+=s[v].cnt[0];
			s[a].cnt[1]+=s[v].cnt[1];
			add(s[a].sumy,s[v].sumy);
		}
		if(a)
		{
			int tmp=0;
			if(op)
			{
				add(tmp,1ll*(s[a].len-s[s[a].link].len)*s[a].cnt[0]%Mod*sum2[lenS]%Mod);
				add(tmp,1ll*chk(ssum2[s[a].len]-ssum2[s[s[a].link].len])*s[a].cnt[0]%Mod);
			}
			del(tmp,1ll*(s[a].len-s[s[a].link].len)*s[a].sumy%Mod);
			add(ans,1ll*tmp*s[a].cnt[1]%Mod);
		}
	}
	void calc(bool op)
	{
		for(int i=1;i<=poi_cnt;i++) son[s[i].link].push_back(i);
		dfs(0,op);
	}
}
int cur;
int main()
{
	lenS=get_str(S);
	lenT=get_str(T);
	
	for(int i=1;i<=500000;i++) sum2[i]=1ll*i*i%Mod;
	for(int i=1;i<=500000;i++) add(sum2[i],sum2[i-1]);
	for(int i=1;i<=500000;i++) add(sum2[i],sum2[i-1]);
	
	for(int i=2;i<=500000;i++) ssum2[i]=sum2[i-2];
	for(int i=1;i<=500000;i++) add(ssum2[i],ssum2[i-1]);

	SAM::reset();cur=0;
	for(int i=1;i<=lenS;i++)
	{
		cur=SAM::ins_nd(cur,S[i]-'a');
		SAM::s[cur].cnt[0]=1,
		SAM::s[cur].sumy=sum2[i-1];
	}
	cur=SAM::ins_nd(cur,26);
	for(int i=1;i<=lenT;i++)
	{
		cur=SAM::ins_nd(cur,T[i]-'a');
		SAM::s[cur].cnt[1]=1;
	}
	SAM::calc(true);

	SAM::reset();cur=0;
	for(int i=lenS;i;i--)
	{
		cur=SAM::ins_nd(cur,S[i]-'a');
		SAM::s[cur].cnt[0]=1;
		SAM::s[cur].sumy=sum2[lenS-i];
	}
	cur=SAM::ins_nd(cur,26);
	for(int i=lenT;i;i--)
	{
		cur=SAM::ins_nd(cur,T[i]-'a');
		SAM::s[cur].cnt[1]=1;
	}
	SAM::calc(false);

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