#include<bits/stdc++.h>
using namespace std;
using ll=long long;
const int N=1<<19;
const int mod=998244353;
using arr=int[N+5];
int add(int x,int y){ return (x+=y)>=mod?x-mod:x; }
int sub(int x,int y){ return (x-=y)<0?x+mod:x; }
void red(int &x){ x+=x>>31&mod;}
int ksm(ll x,int tp,int s=1){
    for(;tp;x=x*x%mod,tp>>=1) if(tp&1) s=x*s%mod;
    return s;
}
int fac[N+5],ifac[N+5];
void prep(){
    fac[0]=1;
    for(int i=1;i<=N;i++) fac[i]=1ll*fac[i-1]*i%mod;
    ifac[N]=ksm(fac[N],mod-2);
    for(int i=N;i;i--) ifac[i-1]=1ll*ifac[i]*i%mod;
}
int C(int n,int m){ return n<0||m<0||n<m?0:1ll*fac[n]*ifac[m]%mod*ifac[n-m]%mod; }
namespace poly_base{
	int l,n,iv; arr g;
	void r_prep(int tl){
		int i;
		l=(int)log2(tl-1)+1; n=1<<l;
		g[0]=1; g[n]=ksm(31,1<<21-l);
		for(i=l;i;i--) g[1<<i-1]=1ll*g[1<<i]*g[1<<i]%mod;
		for(i=0;i<n;i++) g[i]=1ll*g[i&i-1]*g[i&-i]%mod;
	}
	void prep(int tl){ n=1<<(l=(int)log2(tl-1)+1); iv=mod-(mod-1)/n; }
	void DIT(arr T){
		int i,*t,*j,*k,v;
		for(i=n>>1;i;i>>=1)
			for(t=g,j=T;j!=T+n;j+=i<<1,++t)
				for(k=j;k!=j+i;++k)
					v=1ll**t*k[i]%mod,red(k[i]=*k-v),red(*k+=v-mod);
	}
	void DIF(arr T){
		int i,*t,*j,*k,v;
		for(i=1;i<n;i<<=1)
			for(t=g,j=T;j!=T+n;j+=i<<1,++t)
				for(k=j;k!=j+i;++k)
					red(v=*k+k[i]-mod),k[i]=1ll**t*(*k-k[i]+mod)%mod,*k=v;
		reverse(T+1,T+n);
		for(int i=0;i<n;i++) T[i]=1ll*T[i]*iv%mod;
	}
	void NTT(arr a){ DIT(a); }
	void INTT(arr a){ DIF(a); }
	void NTT(arr a,arr b){ memcpy(b,a,n*sizeof(int)); DIT(b); }
	void INTT(arr a,arr b){ memcpy(b,a,n*sizeof(int)); DIF(b); }
}
namespace poly{
	using namespace poly_base;
	arr pa,pb,pc,pd,pe,pf;
	#define szf sizeof(int)
	void cl(arr a){ memset(a,0,n*szf); }
	void Mult(int t,arr a,arr b,arr c){
		if(t==1){ *c=1ll**a**b%mod; return; }
		prep(t); NTT(a,c); NTT(b,pa);
		for(int i=0;i<n;i++) c[i]=1ll*c[i]*pa[i]%mod;
		INTT(c); cl(pa);
	}
}
int h,w;
int ans;
arr f,g,res;
int work(int a,int b){
    return ((1ll*C(a+b-2,a-1)*C(a+b-2,b-1)-1ll*C(a+b-2,a)*C(a+b-2,b))%mod+mod)%mod;
}
void calc(int h,int w){
    for(int i=1;i<h;i++)
        ans=(ans+1ll*(h-i-1)*work(w,i))%mod;
}
int main(){
    prep();
    poly_base::r_prep(N);
    scanf("%d %d",&h,&w);
    for(int i=0;i<h;i++)
        f[i]=1ll*ifac[i]*ifac[i]%mod;
    for(int i=0;i<w;i++)
        g[i]=1ll*ifac[i]*ifac[i]%mod;
    poly::Mult(h+w,f,g,res);
    for(int i=0;i<h+w;i++)
        ans=(ans+1ll*fac[i]%mod*fac[i]%mod*res[i])%mod;
    memset(f,0,sizeof f);
    memset(g,0,sizeof g);
    memset(res,0,sizeof res);
    for(int i=1;i<h;i++)
        f[i]=1ll*ifac[i+1]*ifac[i-1]%mod;
    for(int i=1;i<w;i++)
        g[i]=1ll*ifac[i+1]*ifac[i-1]%mod;
    poly::Mult(h+w,f,g,res);
    for(int i=0;i<h+w;i++)
        ans=(ans-1ll*fac[i]%mod*fac[i]%mod*res[i])%mod;
    red(ans);
    // for(int i=1;i<=h;i++)
    //     for(int j=1;j<=w;j++)
    //         ans=add(ans,work(i,j));
    ans=add(ans,ans);
    calc(h,w);
    calc(w,h);
    ans=sub(ans,work(h,w));
    ans=add(ans,2);
    printf("%d\n",ans);
    return 0;
}