#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define MP make_pair
#define mod MOD
mt19937 rnd(time(0));
const int N=1<<19,MOD=998244353;
inline void add(int& x,ll y){ x=(x+y)%mod; }
int inv[N],fac[N],ifac[N];
namespace P {
const int G=3;
int rev[N],w[N<<1];
int ksm(int a,int b=MOD-2) {
	int ret=1;
	for(;b;a=1ll*a*a%MOD,b=b>>1) if(b&1) ret=1ll*ret*a%MOD;
	return ret;
}
void poly_init() {
	inv[1]=1;
	for(int i=2;i<N;++i) inv[i]=1ll*(MOD-MOD/i)*inv[MOD%i]%MOD;
	fac[0]=ifac[0]=1;
	for(int i=1;i<N;++i) fac[i]=1ll*fac[i-1]*i%MOD,ifac[i]=1ll*ifac[i-1]*inv[i]%MOD;
	for(int k=1;k<=N;k<<=1) {
		int x=ksm(G,(MOD-1)/k); w[k]=1;
		for(int i=1;i<k;++i) w[i+k]=1ll*x*w[i+k-1]%MOD;
	}
}
int plen(int x) { int y=1; for(;y<x;y<<=1); return y;  }
void ntt(int *f,bool idft,int n) {
	for(int i=0;i<n;++i) {
		rev[i]=(rev[i>>1]>>1);
		if(i&1) rev[i]|=n>>1;
	}
	for(int i=0;i<n;++i) if(rev[i]<i) swap(f[i],f[rev[i]]);
	for(int k=2,x,y;k<=n;k<<=1) {
		for(int i=0;i<n;i+=k) {
			for(int j=i;j<i+k/2;++j) {
				x=f[j],y=1ll*f[j+k/2]*w[k+j-i]%MOD;
				f[j]=(x+y>=MOD)?x+y-MOD:x+y,f[j+k/2]=(x>=y)?x-y:x+MOD-y;
			}
		}
	}
	if(idft) {
		reverse(f+1,f+n);
		for(int i=0,x=ksm(n);i<n;++i) f[i]=1ll*f[i]*x%MOD;
	}
}
void poly_mul(const int *f,const int *g,int *h,int n,int m) {
	static int a[N],b[N];
	for(int i=0;i<n;++i) a[i]=f[i];
	for(int i=0;i<m;++i) b[i]=g[i];
	int len=plen(n+m-1);
	ntt(a,0,len),ntt(b,0,len);
	for(int i=0;i<len;++i) h[i]=1ll*a[i]*b[i]%MOD;
	ntt(h,1,len);
	memset(a,0,sizeof(int)*len);
	memset(b,0,sizeof(int)*len);
}
}
int n,m;
inline constexpr int sqr(int x){ return (ll)x*x%mod; }
inline int binom(int u,int d){
	if(u<0||d>u) return 0;
	return (ll)fac[u]*ifac[d]%mod*ifac[u-d]%mod;
}
inline int calc(int x,int y){
	return (sqr(binom(x+y-2,x-1))-(ll)binom(x+y-2,x)*binom(x+y-2,x-2)%mod+mod)%mod;
}
inline int calc2(int i,int j){
	return (ll)sqr(fac[i+j-2])%mod*(i+j-1)%mod*
	sqr(ifac[i-1])%mod*sqr(ifac[j-1])%mod*P::ksm(i*j)%mod;
}
int f[N],g[N];
signed main(){
	P::poly_init();
	ios::sync_with_stdio(false);
	// freopen("Otomachi_Una.in","r",stdin);
	// freopen("Otomachi_Una.out","w",stdout);
	int n,m;
	cin>>n>>m;
	int ans=2;
	for(int i=1;i<n||i<m;i++){
		f[i]=(ll)sqr(ifac[i-1])*inv[i]%mod;
	}
	P::poly_mul(f,f,g,n,m);
	for(int i=2;i<=n+m-2;i++){
		add(ans,(ll)sqr(fac[i-2])*(i-1)%mod*g[i]%mod*2%mod);
	}
	// for(int i=1;i<n;i++){
	// 	for(int j=1;j<m;j++){
	// 		assert(calc(i,j)==calc2(i,j));
	// 		add(ans,2*calc2(i,j));
	// 	}
	// }
	for(int i=1;i<n;i++){
		add(ans,(ll)calc(m,i)*(n-i+1)%mod);
	}
	for(int i=1;i<m;i++){
		add(ans,(ll)calc(n,i)*(m-i+1)%mod);
	}
	add(ans,(ll)calc(n,m));
	cout<<ans<<endl;
	return 0;
}