#include<bits/stdc++.h> using namespace std; const int N = 510,M = 1e5+10,md = 998244353; int n; long long int m; long long int inv[N+M],pre[N+M]; long long int quickpow(long long int x,int y){ long long int tot=1; while(y){ if(y&1) tot=tot%md*x%md; x=x%md*x%md; y=y>>1; } return tot; } long long int C(int x,int y){ if(x<y) return 0; else return (pre[x]%md*inv[y]%md*inv[x-y]%md+md)%md; } int main(){ cin>>n>>m; pre[0]=1; for(int i=1;i<=n+m;i++){ pre[i]=pre[i-1]*i%md; } for(int i=n+m;i>=0;i--){ inv[i]=quickpow(pre[i],md-2); } long long int ANS=0; for(int i=0;i<=m;i++){ for(int j=0;j<=n;j++){ long long int z=(2*m-2*i)%md*(2*m-2*i)%md; long long int p=C(n+i-1,i)*C(n,j)%md; long long int q=C(m-i-1,j-1)*C(n-j+m-i-1,m-i)%md; ANS=(ANS%md+z%md*p%md*q%md)%md; } } cout<<ANS%md<<endl; return 0; }