#include<bits/stdc++.h>
#define clr(a) memset(a,0,sizeof(a))
using namespace std;
const int N=2.1e6+10,mod=998244353,inv3=(mod+1)/3;
inline int Mod(int a){return a<mod?a:a-mod;}
int l;
typedef long long ll;
ll fastpow(ll a,int b)
{
	ll s=1;
	for(;b;b>>=1,a=a*a%mod)
		if(b&1) s=s*a%mod;
	return s;
}
void dif(int *a)
{
	ll x,y,wn;
	for(int i=l>>1;i;i>>=1)
	{
		wn=fastpow(3,(mod-1)/i/2);
		for(int j=0;j<l;j+=(i<<1))
			for(int k=0,w=1;k<i;k++,w=w*wn%mod)
			{
				x=a[j+k],y=a[i+j+k];
				a[j+k]=Mod(x+y);
				a[i+j+k]=(x-y+mod)*w%mod;
			}
	}
}
void dit(int *a)
{
	ll x,y,wn;
	for(int i=1;i<l;i<<=1)
	{
		wn=fastpow(inv3,(mod-1)/i/2);
		for(int j=0;j<l;j+=(i<<1))
			for(int k=0,w=1;k<i;k++,w=w*wn%mod)
			{
				x=a[j+k],y=1ll*a[i+j+k]*w%mod;
				a[j+k]=Mod(x+y),a[i+j+k]=Mod(x-y+mod);
			}
	}
	x=fastpow(l,mod-2);
	for(int i=0;i<l;i++) a[i]=a[i]*x%mod;
}
void mul(ll *a,ll *b,int n,int m,ll *c)
{
	static int A[N],B[N];
	for(l=1;l<=n+m;l<<=1);
	for(int i=0;i<l;i++) A[i]=B[i]=0;
	for(int i=0;i<=n;i++) A[i]=a[i];
	for(int i=0;i<=m;i++) B[i]=b[i];
	dif(A),dif(B);
	for(int i=0;i<l;i++) A[i]=1ll*A[i]*B[i]%mod;
	dit(A);
	for(int i=0;i<=n+m;i++) c[i]=A[i];
	for(int i=n+m+1;i<l;i++) c[i]=0;
}
ll a[N],b[N],c[N],fac[N],ifac[N];
int n,m;
void init(int n)
{
	fac[0]=ifac[0]=1;
	for(int i=1;i<=n;i++) fac[i]=fac[i-1]*i%mod;
	ifac[n]=fastpow(fac[n],mod-2);
	for(int i=n-1;i;i--) ifac[i]=ifac[i+1]*(i+1)%mod;
}
void pls(ll &a,ll b){a=(a+b)%mod;}
ll F(int a,int b){return a<0||b<0?0:ifac[a]*ifac[b]%mod*fac[a+b]%mod;}
ll sq(ll x){return x*x%mod;}
int main()
{
	ll ans=2;
	scanf("%d%d",&n,&m);
	init(n+m);
	for(int i=0;i<n;i++) a[i]=sq(ifac[i]);
	for(int i=0;i<m;i++) b[i]=sq(ifac[i]);
	mul(a,b,n-1,m-1,c);
	for(int i=0;i<n+m-1;i++)
		pls(ans,c[i]*sq(fac[i]));
	a[0]=b[0]=0;
	for(int i=1;i<n;i++) a[i]=ifac[i-1]*ifac[i+1]%mod;
	for(int i=1;i<m;i++) b[i]=ifac[i-1]*ifac[i+1]%mod;
	mul(a,b,n-1,m-1,c);
	for(int i=2;i<n+m-1;i++)
		pls(ans,(mod-c[i])*sq(fac[i]));
	clr(a),clr(b);
	for(int i=0;i<n-1;i++) a[i]=sq(ifac[i]);
	for(int i=0;i<m-1;i++) b[i]=sq(ifac[i]);
	mul(a,b,n-1,m-1,c);
	for(int i=0;i<n+m-1;i++)
		pls(ans,c[i]*sq(fac[i]));
	a[0]=b[0]=0;
	for(int i=1;i<n-1;i++) a[i]=ifac[i-1]*ifac[i+1]%mod;
	for(int i=1;i<m-1;i++) b[i]=ifac[i-1]*ifac[i+1]%mod;
	mul(a,b,n-1,m-1,c);
	for(int i=2;i<n+m-1;i++)
		pls(ans,(mod-c[i])*sq(fac[i]));
	clr(a),clr(b);
	for(int i=1;i<n;i++)
	{
		pls(ans,(n-i)*sq(F(i-1,m-1)));
		pls(ans,(mod-(n-i))*F(i-2,m)%mod*F(i,m-2));
	}
	for(int i=1;i<m;i++)
	{
		pls(ans,(m-i)*sq(F(i-1,n-1)));
		pls(ans,(mod-(m-i))*F(i-2,n)%mod*F(i,n-2));
	}
	printf("%lld\n",ans);
}