#include<bits/stdc++.h>
using namespace std;
const long long mod=998244353;
long long T,a,b,c,d[1000001],v[1000001],o,h[1000001],fa[1000001],q,w,e,an,cn,fac[1000001],inv[1000001],st[1000001],u[1000001],t[1000001][3],f[101][101];
char s[1000001];
struct p{long long q,w;}l[1000001],g[1000001];
long long tmp[1000001],tmp1[1000001],tmp2[1000001],id[1000001],lim=(1<<19),G=3,tmp3[1000001],lg;
long long pow_(long long qq,long long ww){long long ee=1;while(ww){if(ww&1) ee*=qq,ee%=mod;qq*=qq,qq%=mod,ww>>=1;}return ee%mod;}
inline int read(){int x=0,f=1;char c=getchar();while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+c-'0';c=getchar();}return x*f;}
void add(long long qq,long long ww){l[++o].q=ww,l[o].w=h[qq],h[qq]=o;}
long long gcd(long long qq,long long ww){return !ww?qq:gcd(ww,qq%ww);}
long long find(long long qq){return qq==fa[qq]?qq:fa[qq]=find(fa[qq]);}
void merge(long long qq,long long ww){long long f1=find(qq),f2=find(ww);if(f1==f2) return;fa[f1]=f2;}
long long C(long long qq,long long ww){if(qq<0||ww<0||qq<ww) return 0;return fac[qq]*inv[ww]%mod*inv[qq-ww]%mod;}
void ntt(long long *A,long long *B,long long fl)
{
	for(int i=0;i<lim;i++) B[i]=A[i];
	for(int i=0;i<lim;i++) if(i<id[i]) swap(B[i],B[id[i]]);
	for(int i=1;i<lim;i*=2)
	{
		long long un=pow_(G,(mod-1)/(i*2));
		if(fl==-1) un=pow_(un,mod-2);
		for(int j=0;j<lim;j+=i*2)
		{
			long long om=1;
			for(int k=0;k<i;k++)
			{
				long long x=B[j+k],y=B[i+j+k]*om%998244353;
				B[j+k]=(x+y)%998244353,B[i+j+k]=(x-y+998244353)%998244353;om=om*un%998244353;
			}
		}
	}
}
int main()
{
//	freopen("1.in","r",stdin);
	srand((unsigned)(time(0)^(*new int)));
	fac[0]=1;for(int i=1;i<=1000000;i++) fac[i]=fac[i-1]*i%mod;
	inv[1000000]=pow_(fac[1000000],mod-2);for(int i=999999;i>=0;i--) inv[i]=inv[i+1]*(i+1)%mod;
	scanf("%lld%lld",&a,&b);
	for(int i=a-1;i>=1;i--)
	{
		l[++cn]=p{0,i};
	}
	for(int i=0;i<b;i++) l[++cn]=p{i,0};
	cn=0;
	for(int i=1;i<=b;i++) g[++cn]=p{i,a};
	for(int i=a-1;i>=1;i--) g[++cn]=p{b,i};
	for(int i=1;i<=100;i++)
	{
		for(int j=1;j<=100;j++)
		{
//			0,1~i-1,j
//			1,0~i,j-1
			f[i][j]=(C(i+j-2,i-1)*C(i+j-2,i-1)-C(i+j-2,i-2)*C(i+j-2,i))%mod;
		}
	}
	lg=log2(lim);for(int i=0;i<lim;i++) id[i]=(id[i>>1]>>1)+((i&1)<<(lg-1));
	for(int i=0;i<=a-2;i++) tmp[i]=inv[i]*inv[i]%mod;
	for(int i=0;i<=b-1;i++) tmp1[i]=inv[i]*inv[i]%mod;
	ntt(tmp,tmp2,1);ntt(tmp1,tmp3,1);
	for(int i=0;i<lim;i++) tmp2[i]=tmp2[i]*tmp3[i]%mod;
	ntt(tmp2,tmp3,-1);long long tt=pow_(lim,mod-2);
	for(int i=0;i<lim;i++) tmp3[i]=tmp3[i]*tt%mod;
	for(int i=0;i<lim;i++) an=(an+tmp3[i]*fac[i]%mod*fac[i]%mod)%mod;
	memset(tmp,0,sizeof(tmp));
	memset(tmp1,0,sizeof(tmp1));
	for(int i=1;i<=a-2;i++) tmp[i]=inv[i-1]*inv[i+1]%mod;
	for(int i=1;i<=b-1;i++) tmp1[i]=inv[i-1]*inv[i+1]%mod;
	ntt(tmp,tmp2,1);ntt(tmp1,tmp3,1);
	for(int i=0;i<lim;i++) tmp2[i]=tmp2[i]*tmp3[i]%mod;
	ntt(tmp2,tmp3,-1);tt=pow_(lim,mod-2);
	for(int i=0;i<lim;i++) tmp3[i]=tmp3[i]*tt%mod;
	for(int i=0;i<lim;i++) an=(an-tmp3[i]*fac[i]%mod*fac[i]%mod)%mod;
//	for(int i=0;i<=a-2;i++)
//	{
//		for(int j=0;j<=b-1;j++)
//		{
//			an=(an+C(i+j,j)*C(i+j,i)-C(i+j,j-1)*C(i+j,i-1))%mod;
//		}
//	}
	an=an*2%mod;
	for(int i=1;i<=b;i++)
	{
		an=(an+(C(a+i-2,i-1)*C(a+i-2,i-1)-C(a+i-2,i-2)*C(a+i-2,i))%mod*(b-i+1))%mod;
	}
	for(int i=1;i<=a-2;i++)
	{
		an=(an+(C(b+i-2,i-1)*C(b+i-2,i-1)-C(b+i-2,i-2)*C(b+i-2,i))%mod*(a-1-i))%mod;
	}
	an=(an+2)%mod;
	printf("%lld",(an%mod+mod)%mod);
	return 0;
}