#include<iostream>
#include<vector>
#include<cstring>
using namespace std;
typedef long long ll;
const int mod = 998244353;
int fac[1 << 19], inv[1 << 19], Inv[1 << 19];
int omg[1 << 20], rev[1 << 20];
int F[1 << 19], X[1 << 19], Y[1 << 19], Z[1 << 19];
int ans[1 << 19];
int qpow(int x, int y)
{
    int ans = 1;
    while(y)
    {
        if(y & 1)
            ans = (ll)ans * x % mod;
        x = (ll)x * x % mod;
        y >>= 1;
    }
    return ans;
}
void prep(int n)
{
    int i, j, w, m = (1 << n) - 1;
    fac[0] = 1;
    for(i = 1; i <= m; i++)
        fac[i] = (ll)fac[i - 1] * i % mod;
    inv[m] = qpow(fac[m], mod - 2);
    for(i = m; i >= 1; i--)
        inv[i - 1] = (ll)inv[i] * i % mod;
    for(i = 1; i <= m; i++)
        Inv[i] = (ll)inv[i] * fac[i - 1] % mod;
    for(i = 0; i <= n; i++)
    {
        for(j = 1; j < (1 << i); j++)
            rev[(1 << i) + j] = (rev[(1 << (i - 1)) + (j >> 1)] | ((j & 1) << (i - 1)));
        omg[1 << i] = 1;
        w = qpow(3, (mod - 1) >> i);
        for(j = 1; j < (1 << i); j++)
            omg[(1 << i) + j] = (ll)omg[(1 << i) + j - 1] * w % mod;
    }
}
void ntt(int *A, int n, bool t)
{
    int i, j, k, w;
    for(i = 0; i < (1 << n); i++)
        if(i < rev[(1 << n) + i])
            swap(A[i], A[rev[(1 << n) + i]]);
    for(i = 0; i < n; i++)
        for(j = 0; j < (1 << n); j += (2 << i))
            for(k = 0; k < (1 << i); k++)
            {
                w = (ll)omg[(2 << i) + k] * A[(1 << i) + j + k] % mod;
                A[(1 << i) + j + k] = (A[j + k] + mod - w) % mod;
                A[j + k] = (A[j + k] + w) % mod;
            }
    if(t)
    {
        for(i = 1; i < (1 << n); i++)
            if(i < (1 << n) - i)
                swap(A[i], A[(1 << n) - i]);
        w = qpow(1 << n, mod - 2);
        for(i = 0; i < (1 << n); i++)
            A[i] = (ll)A[i] * w % mod;
    }
}
struct poly
{
    vector<int> c;
    poly()
    {
        c.clear();
    }
    friend poly operator * (poly f, poly g)
    {
        if(f.c.empty() || g.c.empty())
            return poly();
        poly h;
        int i, l = 0;
        while((f.c.size() + g.c.size()) >> l)
            l++;
        memset(X, 0, 4 << l);
        memset(Y, 0, 4 << l);
        for(i = 0; i < f.c.size(); i++)
            X[i] = f.c[i];
        for(i = 0; i < g.c.size(); i++)
            Y[i] = g.c[i];
        ntt(X, l, 0);
        ntt(Y, l, 0);
        for(i = 0; i < (1 << l); i++)
            X[i] = (ll)X[i] * Y[i] % mod;
        ntt(X, l, 1);
        h.c.resize(f.c.size() + g.c.size() - 1);
        for(i = 0; i < h.c.size(); i++)
            h.c[i] = X[i];
        return h;
    }
    friend poly operator / (poly f, poly g)
    {
        if(f.c.empty() || g.c.empty())
            return poly();
        poly h;
        int i, l = 0;
        while((f.c.size() + g.c.size()) >> l)
            l++;
        memset(X, 0, 4 << l);
        memset(Y, 0, 4 << l);
        for(i = 0; i < f.c.size(); i++)
            X[i] = f.c[i];
        for(i = 0; i < g.c.size(); i++)
            Y[i] = g.c[g.c.size() - 1 - i];
        ntt(X, l, 0);
        ntt(Y, l, 0);
        for(i = 0; i < (1 << l); i++)
            X[i] = (ll)X[i] * Y[i] % mod;
        ntt(X, l, 1);
        h.c.resize(f.c.size());
        for(i = 0; i < h.c.size(); i++)
            h.c[i] = X[i + g.c.size() - 1];
        return h;
    }
};
int n,m,s[400010];
inline int C(int n,int m)
{
	return n<m?0:(long long)
	s[n]*inv[m]%mod*inv[n-m]%mod;
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	prep(19);
	cin>>n>>m;
	for(int i=s[0]=inv[0]=1;i<=n+m;i++)
	{
		s[i]=(long long)s[i-1]*i%mod;
		inv[i]=qpow(s[i],mod-2);
	}
	__int128_t ans=0;
	ans+=(long long)C(n+m-2,n-1)*C(n+m-2,n-1);
	ans-=(long long)C(n+m-2,n)*C(n+m-2,m);
	for(int i=1;i<n;i++)
	{
		ans+=2ll*C(m+i-2,m-1)*C(m+i-2,m-1);
		ans-=2ll*C(m+i-2,i)*C(m+i-2,m);
	}
	for(int i=1;i<m;i++)
	{
		ans+=2ll*C(n+i-2,n-1)*C(n+i-2,n-1);
		ans-=2ll*C(n+i-2,i)*C(n+i-2,n);
	}
	vector<int> v1(n-1),v2(m-1);
	for(int i=0;i<n-1;i++)
	{
		v1[i]=(long long)
		inv[i]*inv[i]%mod;
	}
	for(int i=0;i<m-1;i++)
	{
		v2[i]=(long long)
		inv[i]*inv[i]%mod;
	}
	poly a,b,c;
	a.c=v1,b.c=v2,c=a*b;
	for(int i=0;i<c.c.size();i++)
	{
		ans+=2ll*s[i]*s[i]%mod*c.c[i];
	}
	for(int i=0;i<n-1;i++)
	{
		v1[i]=i?(long long)
		inv[i+1]*inv[i-1]%mod:0;
	}
	for(int i=0;i<m-1;i++)
	{
		v2[i]=i?(long long)
		inv[i+1]*inv[i-1]%mod:0;
	}
	a.c=v1,b.c=v2,c=a*b;
	for(int i=0;i<c.c.size();i++)
	{
		ans-=2ll*s[i]*s[i]%mod*c.c[i];
	}
	for(int i=1;i<=n-2;i++)
	{
		ans+=(long long)(n-1-i)*C(i+m-2,m-1)%mod*C(i+m-2,m-1);
		ans-=(long long)(n-1-i)*C(i+m-2,i)%mod*C(i+m-2,m);
	}
	for(int i=1;i<=m-2;i++)
	{
		ans+=(long long)(m-1-i)*C(i+n-2,n-1)%mod*C(i+n-2,n-1);
		ans-=(long long)(m-1-i)*C(i+n-2,i)%mod*C(i+n-2,n);
	}
	cout<<int((ans+2)%mod+mod)%mod<<endl;
	return 0;
}