#include <bits/stdc++.h>
using namespace std;

const int N=5e5+5,H=998244353;
int n,m,fac[N],ifac[N],ans=0;

int adc(int a,int b){return a+b>=H?a+b-H:a+b;}
int dec(int a,int b){return a<b?a-b+H:a-b;}
int mul(int a,int b){return 1ll*a*b%H;}
void add(int &a,int b){a=adc(a,b);}
void del(int &a,int b){a=dec(a,b);}

int qpow(int a,int b=H-2){
  int res=1;
  while(b){if(b&1) res=mul(res,a);a=mul(a,a),b>>=1;}
  return res;
}

namespace Poly{
  #define poly vector<int>
  int U[N<<2],P,w[N<<2];
  const int g=3,ig=qpow(g);
  
  void Init(int n){
  	for(int i=1;i<=n;i++) U[i]=(U[i>>1]>>1)|((i&1)?n>>1:0);
  }
  
  void Cro(poly &F,poly G){
  	for(int i=0;i<(int)F.size();i++) F[i]=mul(F[i],G[i]);
  }
  
  void NTT(poly &p,int len,int op){
  	for(int i=0;i<len;i++) if(i<U[i]) swap(p[i],p[U[i]]);
  	int u,tmp;w[0]=1;
  	for(int k=2,lst=1;k<=len;lst=k,k<<=1){
  	  u=qpow(op==1?g:ig,(H-1)/k);
	  for(int i=1;i<=lst;i++) w[i]=mul(w[i-1],u);
	  for(int i=0;i<len;i+=k)
	    for(int l=i;l<i+lst;l++)
	      tmp=mul(p[l+lst],w[l-i]),p[l+lst]=dec(p[l],tmp),add(p[l],tmp);
	}
	if(op==-1){
	  int I=qpow(len);
	  for(int i=0;i<len;i++) p[i]=mul(p[i],I);
	}
  }
  
  poly operator *(poly F,poly G){
  	int n=(int)F.size(),m=(int)G.size();
  	for(P=1;P<n+m;P<<=1);
  	Init(P),F.resize(P),G.resize(P);
  	
  	NTT(F,P,1),NTT(G,P,1),Cro(F,G),NTT(F,P,-1);
  	
  	F.resize(n+m-1);return F;
  }
}
using namespace Poly;

void init(int n){
  fac[0]=1;
  for(int i=1;i<=n;i++) fac[i]=mul(fac[i-1],i);
  ifac[n]=qpow(fac[n]);
  for(int i=n;i>=1;i--) ifac[i-1]=mul(ifac[i],i);
}

int calc(int x,int y){
  return mul(mul(fac[x+y-1],fac[x+y-2]),mul(mul(ifac[x],ifac[y]),mul(ifac[x-1],ifac[y-1])));
}

int main(){
  ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
  cin>>n>>m;
  init(n+m);
  
  poly F(n+m-1),G(n+m-1);
  
  for(int i=1;i<=n-1;i++) F[i]=mul(ifac[i],ifac[i-1]);
  for(int i=1;i<=m-1;i++) G[i]=mul(ifac[i],ifac[i-1]);
  
  F=F*G;
  
  for(int i=2;i<=n+m-2;i++) add(ans,mul(mul(fac[i-1],fac[i-2]),F[i]));
  ans=adc(mul(ans,2),2);
  for(int j=1;j<=m;j++) add(ans,mul(calc(n,j),m-j+1));
  for(int i=1;i<n;i++) add(ans,mul(calc(i,m),n-i+1));
  
  cout<<ans;
  return 0;
}