#include<bits/stdc++.h>
using namespace std;
using ll=long long;
template<typename T>
ostream& operator << (ostream &out,const vector<T>&x){
	if(x.empty())return out<<"[]";
	out<<'['<<x[0];
	for(int len=x.size(),i=1;i<len;i++)out<<','<<x[i];
	return out<<']';
}
template<typename T>
vector<T> ary(const T *a,int l,int r){
	return vector<T>{a+l,a+1+r};
}
template<typename T>
void debug(T x){
	cerr<<x<<'\n';
}
template<typename T,typename ...S>
void debug(T x,S ...y){
	cerr<<x<<' ',debug(y...);
}
using LL=__int128;
const int N=1e6+10,Mod=998244353;
const ll mod=(ll)1e18+3,base=1e9+7;
int n,m,k;
char a[N],b[N];
ll f[N],g[N],pw[N];
ll F(int l,int r){
	return (f[r]+LL(mod-f[l-1])*pw[r-l+1])%mod;
}
ll G(int l,int r){
	return (g[r]+LL(mod-g[l-1])*pw[r-l+1])%mod;
}
int LCP(int i,int j){
	int l=0,r=min(n-i+1,m-j+1)+1,mid;
	for(;l+1<r;){
		mid=(l+r)>>1;
		if(F(i,i+mid-1)==G(j,j+mid-1))l=mid;
		else r=mid;
	}
	return l;
}
int c[3][N],dp[3][N];
int calc(int i){
	int len=LCP(i,1);
	if(len==m||len==n-i+1)return len;
	return len+1+LCP(i+len+1,1+len+1);
}
void add(int i,int l,int r,int x){
	(c[i][l]+=x)%=Mod,(c[i][r+1]+=Mod-x)%=Mod;
}
int main(){
	scanf("%s%s",a+1,b+1),n=strlen(a+1),m=strlen(b+1);
	for(int i=1;i<=n;i++)f[i]=((LL)f[i-1]*base+a[i])%mod;
	for(int i=1;i<=m;i++)g[i]=((LL)g[i-1]*base+b[i])%mod;
	for(int i=pw[0]=1;i<=n||i<=m;i++)pw[i]=(LL)pw[i-1]*base%mod;
	dp[0][0]=1;
//	debug(F(3,5),G(1,3));
	for(int i=1;i<=n;i++){
		int len=calc(i);
//		debug(i,len);
		add(0,i,i+len-1,dp[0][i-1]);
		add(1,i,i+len-1,(dp[0][i-1]+dp[1][i-1])%Mod);
		add(2,i,i+len-1,(dp[0][i-1]+2ll*dp[1][i-1]+dp[2][i-1])%Mod);
		for(int j=0;j<3;j++){
			dp[j][i]=(c[j][i]+=c[j][i-1])%=Mod;
		}
	}
	cout<<dp[2][n]<<endl;
	return 0;
}