#include<bits/stdc++.h>
#define x first
#define y second
#define bug(x) cerr<<#x<<'='<<x<<' '
#define debug(x) cerr<<#x<<'='<<x<<'\n'
#define FOR(a,b,c) for(int a=(b),a##_end=(c);a<=a##_end;++a)
#define ROF(a,b,c) for(int a=(b),a##_end=(c);a>=a##_end;--a)
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int INF=0x3f3f3f3f,N=1010,P=998244353;
template<class T>inline bool chkmin(T &x,const T &y){return y<x?x=y,1:0;}
template<class T>inline bool chkmax(T &x,const T &y){return x<y?x=y,1:0;}
template<class T>inline void Add(T &x,const T &y){x=(x+y+P)%P;}
int n,m;
char a[N],b[N];
ll dp[N][N][2],ans,res1,res2;
int main(){
	scanf("%s%s",a+1,b+1);
	n=strlen(a+1);
	m=strlen(b+1);
	dp[n+1][m+1][0]=dp[n+1][m+1][1]=1;
	ROF(i,n+1,1)ROF(j,m+1,1)FOR(k,0,1){
		if((k==0&&i==1)||(k==1&&j==1))continue;
		int las=k==0?a[i-1]:b[j-1];
		dp[i][j][k]=1;
		if(i<=n&&las!=a[i])
			Add(dp[i][j][k],dp[i+1][j][0]);
		if(j<=m&&las!=b[j])
			Add(dp[i][j][k],dp[i][j+1][1]);
	}
	dp[1][1][0]=(dp[2][1][0]+dp[1][2][1])%P;
	FOR(i,1,n)FOR(j,1,m){
		Add(ans,dp[i+1][j][0]);
		Add(ans,dp[i][j+1][1]);
	}
	FOR(i,1,n){
		int j=i+1;
		for(;j<=n;++j)
			if(a[j-1]==a[j])
				break;
		res1+=j-i;
	}
	FOR(i,1,m){
		int j=i+1;
		for(;j<=m;++j)
			if(b[j-1]==b[j])
				break;
		res2+=j-i;
	}
	ans=ans-n*res2-m*res1;
	printf("%lld\n",(ans%P+P)%P);
	return 0;
}
