本文版权归ljh2000和博客园共有，欢迎转载，但须保留此声明，并给出原文链接，谢谢合作。

//It is made by ljh2000
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <ctime>
#include <vector>
#include <queue>
#include <map>
#include <set>
#include <string>
#include <complex>
using namespace std;
typedef long long LL;
typedef long double LB;
typedef complex<double> C;
const double pi = acos(-1);
const int MAXN = 5011;
const int mod = 100000000;
char ch[MAXN],s[MAXN];
int n,m,f[2][MAXN],g[2][MAXN];inline int getint(){int w=0,q=0; char c=getchar(); while((c<'0'||c>'9') && c!='-') c=getchar();if(c=='-') q=1,c=getchar(); while (c>='0'&&c<='9') w=w*10+c-'0',c=getchar(); return q?-w:w;
}inline void work(){scanf("%s",ch+1); scanf("%s",s+1);n=strlen(ch+1); m=strlen(s+1); n--; m--;for(int i=0;i<=m;i++) g[0][i]=1;//边界!int tag=0; g[1][0]=1;for(int i=1;i<=n;i++) {tag^=1; //memset(f[tag],0,sizeof(f[tag]));//memset(g[tag],0,sizeof(g[tag]));for(int j=1;j<=m;j++) {if(ch[i]==s[j]) {f[tag][j]=f[tag^1][j-1]+1;g[tag][j]=g[tag^1][j-1];g[tag][j]+=(f[tag][j]==f[tag^1][j])*g[tag^1][j];g[tag][j]+=(f[tag][j]==f[tag][j-1])*g[tag][j-1];}else {f[tag][j]=max(f[tag][j-1],f[tag^1][j]);g[tag][j]=(f[tag][j]==f[tag^1][j])*g[tag^1][j];g[tag][j]+=(f[tag][j]==f[tag][j-1])*g[tag][j-1];g[tag][j]-=(f[tag][j]==f[tag^1][j-1])*g[tag^1][j-1];}g[tag][j]%=mod;}}printf("%d\n",f[tag][m]);g[tag][m]+=mod; g[tag][m]%=mod;printf("%d",g[tag][m]);
}int main()
{work();return 0;
}